题目链接
题解:很明显这是个dp,很明显每一行只能最多放两个炮,这样就很容易状态表示了,我们用dp[i][j][k]来表示前i行放了一个炮的列数有j个,放了两个的炮列数有k个,我们很容易可以得到放了0个炮的列数有多少个。具体怎么转移看代码注释。
#include<bits/stdc++.h>
#define ll long long
#define pr pair<ll,ll>
#define ios ios::sync_with_stdio(false)
#define CRL(a) memset(a,0,sizeof a)
#define endl "\n"
using namespace std;
const int maxn = 2e5 + 5;
const int mod = 9999973;
ll dp[105][105][105];
ll p(ll a)//C(a,2)
{
return ((a * (a - 1)) / 2) % mod;
}
int main()
{
int n, m;
cin >> n >> m;
dp[0][0][0] = 1;
for (int i = 1; i <= n;i++)
{
for (int j = 0; j <= m;j++)
{
for (int k = 0; k + j <= m;k++)
{
dp[i][j][k] += dp[i - 1][j][k];//这一行什么也不放
if(k>=1)
(dp[i][j][k] += dp[i - 1][j + 1][k - 1] * (j + 1)) %= mod;//从第i减一行里的j+1个放一个炮的里挑一个
if(j>=1)
(dp[i][j][k] += dp[i - 1][j - 1][k] * (m - k - j + 1)) %= mod;//从第i减一行里的m-k-j+1个放零炮的里挑一个
if(k>=2)
(dp[i][j][k] += dp[i - 1][j + 2][k - 2] * p(j + 2)) %= mod;//从第i减一行里的j+2个放一个炮的里挑两个
if(j>=2)
(dp[i][j][k] += dp[i - 1][j - 2][k] * p(m - k - j + 2)) %= mod;//同上 只不过是从放零个里挑
if(k>=1)
(dp[i][j][k] += dp[i - 1][j][k - 1] * j * (m - j - k + 1)) %= mod;//从放零个的行里面选一个,放一个的行里选一个,这样j是不变的。
}
}
}
ll ans = 0;
for (int i = 0; i <= m;i++)
{
for (int j = 0; i + j <= m;j++)
{
(ans += dp[n][i][j]) %= mod;
}
}
cout << ans;
}