bzoj4806 炮

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=4806

这种题应该想状压的。

于是发现压不下，结合每一行每一列最多放两个炮想到记一下放炮的列就行。

**又考虑到当前行能怎么放只与之前放0/1/2个炮的列数有关，与具体在哪一列无关。

　　于是记录一下前 i-1 行有多少列放了1/2个炮，就能顺利转移了。

别忘了考虑全所有情况。

（转移条件只考虑数组不越界就行。像最后那个转移，j=0体现在*j上）

#include<iostream>
#include<cstdio>
#include<cstring>
#define ll long long
using namespace std;
const int N=;const ll mod=;
ll n,m,dp[N][N][N],ans;
ll c(ll a){return (a*(a-)/)%mod;}
int main()
{
    scanf("%lld%lld",&n,&m);
    dp[][][]=;
    for(int i=;i<=n;i++)
        for(int j=;j<=m;j++)
            for(int k=;k<=m-j;k++)
            {
                dp[i][j][k]+=(dp[i-][j][k])%=mod;        //
                if(j)(dp[i][j][k]+=dp[i-][j-][k]*(m-j+-k))%=mod;    //
                if(k&&j+<=m)(dp[i][j][k]+=dp[i-][j+][k-]*(j+))%=mod;    //
                if(j>=)(dp[i][j][k]+=dp[i-][j-][k]*c(m-j+-k))%=mod;    //
                if(k>=)(dp[i][j][k]+=dp[i-][j+][k-]*c(j+))%=mod;    //
                if(k)(dp[i][j][k]+=dp[i-][j][k-]*j*(m-j-k+))%=mod;    ////////
            }
    for(int j=;j<=m;j++) for(int k=;k<=m-j;k++) (ans+=dp[n][j][k])%=mod;
    printf("%lld",ans);
    return ;
}