题意概述:n行m列棋盘放若干个棋子每行每列最多两个求方案总数,答案对9999973取模。
可以比较容易看出这是个dp,设f[i][j][k]表示前i行j列放1个棋子k列放2个棋子的方案总数。转移时分类讨论+计数就好了。分类讨论比较麻烦需要仔细考虑一下。
吐槽:我把循环里的m写成n还有50,查了半天(真的半天)取模QAQ
#include<cstring> #include<iostream> #include<cctype> #include<cstdio> #include<algorithm> #define now f[i][j][k] #define ll unsigned long long using namespace std; inline int read() { register int X=0; register char ch=0; bool flag=0; for(; !isdigit(ch); ch=getchar()) if(ch=='-') flag=1; for(; isdigit(ch); ch=getchar()) X=(X<<3)+(X<<1)+ch-'0'; return (flag ? -X : X); } int N=101; const ll mod=9999973; int n,m; ll f[201][201][201],ans; int C(int x) { return x*(x-1)/2; } int main() { n=read(),m=read(),f[0][0][0]=1; for(int i=0;i<n;i++) for(int j=0;j<=m;j++) for(int k=0;k+j<=m;k++) { if(now) { (f[i+1][j][k]+=now)%=mod; if(m-j-k > 0) (f[i+1][j+1][k]+=now*(m-j-k))%=mod; if(j > 0) (f[i+1][j-1][k+1]+=now*j)%=mod; if(m-j-k > 1) (f[i+1][j+2][k]+=now*(C(m-j-k)%mod))%=mod; if(j > 0 && m-j-k > 0) (f[i+1][j][k+1]+=now*j%mod*(m-j-k))%=mod; if(j > 1) (f[i+1][j-2][k+2]+=now*(C(j)%mod))%=mod; } } for(int i=0; i<=m; i++) for(int j=0; j+i<=m; j++) (ans+=f[n][i][j])%=mod; printf("%llu\n",ans); }