luogu 2051 [AHOI2009]中国象棋

luogu 2051 [AHOI2009]中国象棋

---

真是一道令人愉♂悦丧心并框的好题。。。
首先“没有一个炮可以攻击到另一个炮”有个充分条件就是没有三个炮在同一行或同一列。证明：显然。
所以每行/每列最多2个炮。
还有一个点就是任意交换行和列不会改变这个方案的可行性。
那么可以脑补摆完前i行之后就把所有列排个序，会有摆了1个的列，摆了2个的列和没摆的列。
如果是桶排，就记录一下每个有多少就行了。
所以需要3个桶。可以省掉一个，因为知道了总共m列。
那么设f[i][j][k]为摆完了前i行，且j列摆了1个炮，k列摆了2个炮，你摆放的总方案数，然后滚动一下第一维
然后很多种情况，有注释。
PS.要多取膜！！！被坑到了

// It is made by XZZ
#include<cstdio>
#include<cstring>
#include<algorithm>
#define Fname "cchess"
using namespace std;
#define rep(a,b,c) for(rg int a=b;a<=c;a++)
#define drep(a,b,c) for(rg int a=b;a>=c;a--)
#define erep(a,b) for(rg int a=fir[b];a;a=nxt[a])
#define il inline
#define rg register
#define vd void
typedef long long ll;
#define mod 9999973
il int gi(){
    rg int x=0;rg bool flg=0;rg char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')flg=1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
    return flg?-x:x;
}
ll f[2][101][101];
il vd add(ll&a,ll b){a+=b;if(a>=mod)a-=mod;}
int main(){
#ifdef xzz
	freopen(Fname".in","r",stdin);
	freopen(Fname".out","w",stdout);
#endif
	int n=gi(),m=gi(),now=1;
	f[0][0][0]=1;
	rg ll t;
	rep(i,0,n-1){
		now^=1;
		memset(f[now^1],0,sizeof f[now^1]);
		drep(j,m,0)drep(k,m-j,0)if(f[now][j][k]){
			t=(f[now][j][k]%=mod);
			//第i行
			add(f[now^1][j][k],t);//不放
			add(f[now^1][j+1][k],t*(m-j-k)%mod);//放一个 新占一列
			if(j)add(f[now^1][j-1][k+1],t*j%mod);//放一个 与一列放了一个的放一起
			if(m-j-k>1)add(f[now^1][j+2][k],t*(m-j-k)*(m-j-k-1)/2%mod);//放两个 新占两列
			if(j&&m-j-k>0)add(f[now^1][j][k+1],t*(m-j-k)*j%mod);//放两个 一个新占一列 另一个与一列放了一个的放一起
			if(j>1)add(f[now^1][j-2][k+2],t*j*(j-1)/2%mod);//放两个 都与一列放了一个的放一起
		}
	}
	ll ans=0;
	now^=1;
	rep(i,0,m)rep(j,0,m-i)add(ans,f[now][i][j]);
	printf("%lld
",ans);
	return 0;
}