【LOJ】#2026. 「JLOI / SHOI2016」成绩比较

题解

用\(f[i][j]\)表示考虑了前i个排名有j个人被碾压
\(f[i][j] = f[i - 1][k] \* C[k][j] \* C[N - k - 1][N - r[i] - j] \* P[i]\)
P[i]是成绩排列的方式，意义是在前面k个人里选了j个来碾压，并将人数空缺用上一次没有碾压的来填补

\(P[i]\)怎么求，对于一个i，考虑枚举B君的成绩，也就是
\(\sum_{j = 1}^{u_{i}} j^{N - r[i]}(u_{i} - j)^{r[i] - 1}\)
观察一下把\(u_{i}\)当成一个变量，这是一个n次多项式，插值就能做了，暴力插值即可
复杂度\(O(n^3)\)

代码

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <ctime>
#define pii pair<int,int>
#define fi first
#define se second
#define mo 974711
//#define ivorysi
using namespace std;
typedef long long int64;
const int MOD = 1000000007;
int N,M,K,r[105],u[105];
int64 f[105],P[105],inv[105];
int64 fpow(int64 x,int c) {
    int64 res = 1,t = x;
    while(c) {
        if(c & 1) res = res * t % MOD;
        t = t * t % MOD;
        c >>= 1;
    }
    return res;
}

void Solve() {
    scanf("%d%d%d",&N,&M,&K);
    for(int i = 1 ; i <= M ; ++i) scanf("%d",&u[i]);
    for(int i = 1 ; i <= M ; ++i) scanf("%d",&r[i]);
    for(int T = 1 ; T <= M ; ++i) {
        for(int i = 1 ; i <= r[T] ; ++i) {
            f[i] = 0;
            for(int j = 1 ; j < i ; ++j) {
                f[i] += fpow(j,n - r[T] + 1) * fpow(i - j,r[T] - 1) % MOD;
                f[i] %= MOD;
            }
        }
        for(int i = 1 ; i <= r[T] ; ++i) {
            for(int j = 1 ; j <= r[T] ; ++j) {
                if(i == j) continue;

            }
        }
    }   
}
int main() {
#ifdef ivorysi
    freopen("f1.in","r",stdin);
#endif
    Solve();
}