题意
略。
题解
想到如果相对位置确定,求出具体方案数就是再乘一个简便的组合数即可。
所以将人数和分数分开来求。
观察题目的限制条件:
1.
1.
1. 有
k
k
k 名同学分数小于等于 B。
2.
2.
2. 剩下的
n
−
k
−
1
n-k-1
n−k−1 名同学至少有一门科目大于 B。
所以到这里就有两种可行的方向:
1.
1.
1. 通过
d
p
dp
dp 求出具体有几名同学被 B 碾压。
2.
2.
2. 通过容斥求出至少有
i
i
i 名同学被 B 碾压。
本题解采用第一种方法:
首先想到如果有同学前面已经有一科大于 B 了,那么接下来的学科他无论什么样的分数都不会影响 B 碾压的同学数。
而前面被 B 碾压的同学必须得在接下来的学科继续被 B 碾压才能算入 B 碾压的同学数。
综上,考虑 f [ i ] [ j ] f[i][j] f[i][j] 表示前 i i i 个学科,有 j j j 个同学被 B 碾压的方案数(不考虑具体分数)。
稍微思考可得到转移式: f [ i ] [ j ] = ∑ k = j n − 1 C k − j k × C r [ i ] − 1 − ( k − j ) n − k − 1 × p [ i ] f[i][j]=\sum_{k=j}^{n-1} C_{k-j}^{k} \times C^{n-k-1}_{r[i]-1-(k-j)} \times p[i] f[i][j]=∑k=jn−1Ck−jk×Cr[i]−1−(k−j)n−k−1×p[i], p [ i ] p[i] p[i] 是为第 i i i 门学科分配分数的方案数。
相当于前一门学科有 k − j k-j k−j 个同学被 B 碾压,然后在第 i i i 门学科分数大于 B,所以去枚举这 k − j k-j k−j 个同学即 C k − j k C^{k}_{k-j} Ck−jk,然后还有剩下 n − k − 1 n-k-1 n−k−1 名同学是放哪里都无所谓的,选择其中满足有 r [ i ] − 1 r[i]-1 r[i]−1 个同学大于 B 的剩下 r [ i ] − 1 − ( k − j ) r[i]-1-(k-j) r[i]−1−(k−j) 名,即 C r [ i ] − 1 − ( k − j ) n − k − 1 C^{n-k-1}_{r[i]-1-(k-j)} Cr[i]−1−(k−j)n−k−1。
然后关于 p [ i ] p[i] p[i] 的求法,先不考虑 u [ i ] u[i] u[i] 非常大的朴素做法:直接枚举 B 的位置,于是 p [ i ] = ∑ j = 1 u [ i ] j n − r [ i ] × ( u [ i ] − j ) r [ i ] − 1 p[i]= \sum_{j=1}^{u[i]} j^{n-r[i]} \times (u[i]-j)^{r[i]-1} p[i]=∑j=1u[i]jn−r[i]×(u[i]−j)r[i]−1。
可以发现,这是个关于 u [ i ] u[i] u[i] 的多项式,于是用拉格朗日把 u [ i ] 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;
}