(吐槽请无视:自从写STL卡到吐之后写题就一直狂WA,TLE,RE不止,woc我真是没救了)。
不难发现方案符合条件的充分必要条件是:对于任意的i,都有编号<=i的人的个数>i=。那么用f[i][j]表示有j个人,他们的编号都在1..i之间时的方案数。显然f[i][j]>0仅当j>=i,然后枚举编号为i的人的个数,显然可以得到下述转移方程:
f[i][j]=Σ(k=已经确定得编号为i的人的个数,j)f[i-1][j-k]*C(现在还可以选的人的总数,当前位置还没有确定的人的个数)
AC代码如下:
#include<iostream>
#include<cstdio>
#include<cstring>
#define N 305
#define ll long long
using namespace std;
int n,m,p,a[N],sum[N],f[N][N],c[N][N];
void pfs(){
int i,j; memset(c,0,sizeof(c)); c[0][0]=1;
memset(a,0,sizeof(a)); memset(f,0,sizeof(f));
for (i=1; i<=n; i++){
c[i][0]=1;
for (j=1; j<=i; j++) c[i][j]=(c[i-1][j]+c[i-1][j-1])%p;
}
}
int main(){
int cas; scanf("%d",&cas);
while (cas--){
scanf("%d%d%d",&n,&m,&p); int i,j,k; pfs();
for (i=1; i<=m; i++){
int x,y; scanf("%d%d",&x,&y);
a[y]++;
}
f[0][0]=1; sum[0]=n-m;
for (i=1; i<=n; i++){
sum[i]=sum[i-1]+a[i];
if (sum[i]<i){ puts("NO"); break; }
}
if (i<=n) continue;
for (i=1; i<=n; i++)
for (j=i; j<=sum[i]; j++)
for (k=a[i]; k<=j-i+1; k++)
f[i][j]=(f[i][j]+(ll)f[i-1][j-k]*c[sum[i-1]+k-j][k-a[i]]%p)%p;
printf("YES %d\n",f[n][n]);
}
return 0;
}
by lych
2016.2.26