DTOJ #2672. 成绩比较 （mark）

题意

略。

题解

想到如果相对位置确定，求出具体方案数就是再乘一个简便的组合数即可。
所以将人数和分数分开来求。

观察题目的限制条件：
$1.$ 有 $k$ 名同学分数小于等于 B。
$2.$ 剩下的 $n-k-1$ 名同学至少有一门科目大于 B。

所以到这里就有两种可行的方向：
$1.$ 通过 $dp$ 求出具体有几名同学被 B 碾压。
$2.$ 通过容斥求出至少有 $i$ 名同学被 B 碾压。

本题解采用第一种方法：
首先想到如果有同学前面已经有一科大于 B 了，那么接下来的学科他无论什么样的分数都不会影响 B 碾压的同学数。
而前面被 B 碾压的同学必须得在接下来的学科继续被 B 碾压才能算入 B 碾压的同学数。

综上，考虑 $f[i][j]$ 表示前 $i$ 个学科，有 $j$ 个同学被 B 碾压的方案数（不考虑具体分数）。

稍微思考可得到转移式：$f[i][j]={\sum }_{k=j}^{n-1}{C}_{k-j}^k\times C_{r[i]-1-(k-j)}^{n-k-1}\times p[i]$，$p[i]$ 是为第 $i$ 门学科分配分数的方案数。

相当于前一门学科有 $k-j$ 个同学被 B 碾压，然后在第 $i$ 门学科分数大于 B，所以去枚举这 $k-j$ 个同学即 $C_{k-j}^k$，然后还有剩下 $n-k-1$ 名同学是放哪里都无所谓的，选择其中满足有 $r[i]-1$ 个同学大于 B 的剩下 $r[i]-1-(k-j)$ 名，即 $C_{r[i]-1-(k-j)}^{n-k-1}$。

然后关于 $p[i]$ 的求法，先不考虑 $u[i]$ 非常大的朴素做法：直接枚举 B 的位置，于是 $p[i]={\sum }_{j=1}^{u[i]}{j}^{n-r[i]}\times (u[i]-j{)}^{r[i]-1}$。

可以发现，这是个关于 $u[i]$ 的多项式，于是用拉格朗日把 $u[i]$ 处的值插出来即可。

代码

#include<bits/stdc++.h>
using namespace std;
const int N=505;
const int P=1e9+7;
const int Maxn=500;
int n,m,k,u[N],r[N],d[N],y[N],jc[N],f[N][N];
int ksm(int a,int b){
	if(!a&&!b)return 1;
	if(!a)return 0;
	int s=1;while(b){
		if(b&1)s=1ll*s*a%P;
		a=1ll*a*a%P;b>>=1;
	}
	return s;
}
inline int C(int x,int y){
	if(y>x)return 0;
	if(x<0||y<0)return 0;
	return 1ll*jc[x]*ksm(1ll*jc[x-y]*jc[y]%P,P-2)%P;
}
int main(){
	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]);
	jc[0]=1;for(int i=1;i<=Maxn;++i)jc[i]=1ll*jc[i-1]*i%P;
	for(int i=1;i<=m;++i){
		for(int j=1;j<=Maxn;++j)y[j]=0;
		for(int j=1;j<=Maxn;++j){
			for(int z=1;z<=j;++z)y[j]=(y[j]+1ll*ksm(z,n-r[i])*ksm(j-z,r[i]-1)%P)%P;
		}
		if(u[i]<=Maxn){
			d[i]=y[u[i]];
			continue;
		}
		for(int j=1;j<=Maxn;++j){
			int s1=y[j],s2=1;
			for(int z=1;z<=Maxn;++z){
				if(z==j)continue;
				s1=1ll*s1*(u[i]-z)%P;
				s2=1ll*s2*((j-z+P)%P)%P;
			}
			d[i]=(d[i]+1ll*s1*ksm(s2,P-2)%P+P)%P;
		}
	}
	f[0][n-1]=1;
	for(int i=1;i<=m;++i){
		for(int j=k;j<n;++j){
			for(int z=j;z<n;++z){
				f[i][j]=(f[i][j]+1ll*f[i-1][z]*C(z,z-j)%P*C(n-z-1,r[i]-1-z+j)%P*d[i]%P)%P;
			}
		}
	}
	printf("%d\n",f[m][k]);
	return 0;
}