【DP】BZOJ2302 [HAOI2011]Problem c

分析：

比较水的题。。

注意别去用欧拉定理算逆元+Lucas，尼玛N最多才300，直接杨辉三角算组合数。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define SF scanf
#define PF printf
#define MAXN 310
using namespace std;
int n,mod,t,m,val,num;
int a[MAXN];
int dp[MAXN][MAXN],C[MAXN][MAXN],dp1[MAXN][MAXN];
void prepare(){
	C[0][0]=1;
	for(int i=1;i<=n;i++){
		C[i][0]=1;
		for(int j=1;j<=i;j++)
			C[i][j]=(C[i-1][j-1]+C[i-1][j])%mod;
	}
}
int miv[MAXN];
int main(){
	SF("%d",&t);
	while(t--){
		SF("%d%d%d",&n,&m,&mod);
		prepare();
		memset(miv,0,sizeof miv);
		for(int i=1;i<=m;i++){
			SF("%d%d",&val,&num);
			miv[n-num+1]++;
		}
		memset(dp,0,sizeof dp);
		memset(dp1,0,sizeof dp1);
		dp[0][0]=1;
		dp1[0][0]=1;
		int tot=n-m;
		for(int i=0;i<n;i++){
			for(int j=0;j<=i;j++){
//				PF("[%d %d %d]\n",i,j,dp[i][j]);
				for(int k=miv[i+1];k+j<=i+1;k++){
					int add=k-miv[i+1];
					dp1[i+1][j+k]|=dp1[i][j];
					dp[i+1][j+k]=(dp[i+1][j+k]+1ll*dp[i][j]*C[tot-j][add]%mod)%mod;
				}
			}
			tot+=miv[i+1];
		}
		if(dp1[n][n]==0)
			PF("NO\n");
		else
			PF("YES %d\n",dp[n][n]);
	}
}