题目
题解
不敢相信这道题居然是我想出来的qwq。
首先发现象是走斜线,所以先考虑把棋盘顺时针旋转45度,板正。然后稍微变化一下格子的形状,就成了这个样子(0表示空格)
0 0 0 0 0 0 0 0 0
然后你发现,,其实可以每隔一行等价拆成两个棋盘(也就是把黑格和白格拆开)。
0 0 0 0 0 0 0 0 0
然后就可以分别dp了,dp[i][j]表示前i行放了j个,可以枚举这一行放不放,要放的话有多少种可能的位置,可以试着写一写转移方程。
但是大小不一,从宽往窄转移的时候会有问题(如果你试了就会知道),但是我们发现,列的顺序与dp结果并没有关系,于是上面左边的图可以等价为.
0 0 0 0 0
这样再做就简单多了,先做完这个dp,然后处理每一组数据的时候,枚举黑格白格分别放多少个就行了
代码
//QWsin
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=30+10;
int n,k;
namespace Init{
ll C[maxn][maxn];
inline void init_C(){
C[0][0]=1;
for(int i=1;i<=30;++i)
for(int j=0;j<=i;++j)
if(!j||j==i) C[i][j]=1;
else C[i][j]=C[i-1][j]+C[i-1][j-1];
}
ll dp1[maxn][maxn*maxn],dp2[maxn][maxn*maxn];
//dp[i][j]表示前i行共放了j列的方案数
//分奇偶 1奇2偶 奇数第i层1+(i-1)/2*2 偶数第i层(i+1)/2*2
inline int count1(int x){return 1+(x-1)/2*2;}
inline int count2(int x){return (x+1)/2*2;}
inline void init_dp()
{
dp1[0][0]=dp2[0][0]=1;
for(int i=1;i<=30;++i)
{
//分别处理
for(int j=0,sz=count1(i);j<=sz;++j)
{
dp1[i][j]+=dp1[i-1][j];
if(j) dp1[i][j]+=dp1[i-1][j-1]*(sz-j+1);
}
for(int j=0,sz=count2(i);j<=sz;++j)
{
dp2[i][j]+=dp2[i-1][j];
if(j) dp2[i][j]+=dp2[i-1][j-1]*(sz-j+1);
}
}
}
inline void init(){init_C();init_dp();}
}
inline void solve()
{
ll ans=0;
using namespace Init;
for(int x=0;x<=k;++x)
ans+=dp1[n][x]*dp2[n-1][k-x];
printf("%lld\n",ans);
}
int main()
{
Init::init();
while(scanf("%d%d",&n,&k)==2&&(n||k)) solve();
return 0;
}