看到这种题,要么猥琐数学题,要么DP,还能搞搞什么矩乘什么的。
然后这题就是硬刚DP。很容易看出每行每列的棋子数都不超过2吧
f[i][j][k]表示枚举到第i行,有j列有1个棋子,有k列有2个棋子。然后m-j-k就可以的到没放的多少列吧。
枚举每一行,分六种情况。
1、这一行不放
2、放一个,放在没有棋子的列里
3、放一个,放在有棋子的列里
4、放两个,都放在没有棋子的列里
5、放两个,一个放在没有棋子的列里,一个放在有棋子的列里
6、放两个,都放在有棋子的列里
然后转移。方程不写了(其实代码你也看得到),推的要死要活的,丑。
无脑long long谢谢。
#include<cstdio> #include<iostream> #include<cstring> #include<cstdlib> #include<algorithm> #include<cmath> using namespace std; typedef long long LL; const LL mod=9999973; LL f[110][110][110]; int main() { LL n,m; scanf("%lld%lld",&n,&m); if(n==1){printf("%lld\n",1+m+m*(m-1)/2);return 0;} memset(f,0,sizeof(f)); f[1][0][0]=1; f[1][1][0]=m; f[1][2][0]=m*(m-1)/2; //init LL ans=0; for(int i=2;i<=n;i++) for(int j=0;j<=m;j++) for(int k=0;k+j<=m;k++) { f[i][j][k]=(f[i][j][k]+f[i-1][j][k])%mod;//不放 //一个 if(j>0) f[i][j][k]=(f[i][j][k]+f[i-1][j-1][k]*(m-k-(j-1)))%mod;//没有 if(k>0&&j<m) f[i][j][k]=(f[i][j][k]+f[i-1][j+1][k-1]*(j+1))%mod;//有 //两个 if(j>1) f[i][j][k]=(f[i][j][k]+f[i-1][j-2][k]*(m-k-(j-2))*(m-k-(j-2)-1)/2)%mod;//没有 if(j>0&&k>0) f[i][j][k]=(f[i][j][k]+f[i-1][j-1+1][k-1]*(m-(k-1)-j)*j)%mod;//一有一无 if(k>1&&j<m-1) f[i][j][k]=(f[i][j][k]+f[i-1][j+2][k-2]*(j+2)*(j+1)/2)%mod;//都有 if(i==n)ans=(ans+f[i][j][k])%mod; } printf("%lld\n",ans); return 0; }