当前位置 博文首页 > 文章内容

    luogu P2051 [AHOI2009]中国象棋 dp 状态压缩+容斥

    作者: 栏目:未分类 时间:2020-07-03 14:01:23

    本站于2023年9月4日。收到“大连君*****咨询有限公司”通知
    说我们IIS7站长博客,有一篇博文用了他们的图片。
    要求我们给他们一张图片6000元。要不然法院告我们

    为避免不必要的麻烦,IIS7站长博客,全站内容图片下架、并积极应诉
    博文内容全部不再显示,请需要相关资讯的站长朋友到必应搜索。谢谢!

    另祝:版权碰瓷诈骗团伙,早日弃暗投明。

    相关新闻:借版权之名、行诈骗之实,周某因犯诈骗罪被判处有期徒刑十一年六个月

    叹!百花齐放的时代,渐行渐远!



    参考:https://www.luogu.com.cn/blog/RPdreamer/p2051

    #include<map>
    #include<queue>
    #include<time.h>
    #include<limits.h>
    #include<cmath>
    #include<ostream>
    #include<iterator>
    #include<set>
    #include<stack>
    #include<cstdio>
    #include<cstring>
    #include<iostream>
    #include<algorithm>
    using namespace std;
    #define rep_1(i,m,n) for(int i=m;i<=n;i++)
    #define mem(st) memset(st,0,sizeof st)
    inline int read()
    {
        int num=0, w=0;
        char ch=0;
        while (!isdigit(ch))
        {
            w|=ch=='-';
            ch = getchar();
        }
        while (isdigit(ch))
        {
            num = (num<<3) + (num<<1) + (ch^48);
            ch = getchar();
        }
        return w? -num: num;
    }
    typedef long long ll;
    typedef pair<int,int> pii;
    typedef unsigned long long ull;
    typedef pair<double,double> pdd;
    const int inf = 0x3f3f3f3f;
    const int mod=9999973;
    const int N=110;
    #define int ll
    int f[N][N][N];
    //f[i][j][k]代表放了前i行,有j列是有一个棋子,有k列是有2个棋子的合法方案数.
    inline int C(int x)
    {
        return ((x*(x-1))/2)%mod;
    }
    ll n,m,ans;
    void solve()
    {
        cin>>n>>m;
        f[0][0][0]=1;
        for(int i=1; i<=n; i++)
            for(int j=0; j<=m; j++)
                for(int k=0; k<=m-j; k++)
                {
                    //这一行不放棋子
                    f[i][j][k]=f[i-1][j][k];
                    //这一行放一个棋子,分两种:
                    //在一个有棋子的列放
                    if(k>=1)
                        f[i][j][k]+=f[i-1][j+1][k-1]*(j+1);
                    //在没有棋子的列放
                    if(j>=1)
                        f[i][j][k]+=f[i-1][j-1][k]*(m-j-k+1);
                    //放两个棋子:分三种
                    //都放在有一个棋子的列
                    if(k>=2)
                        f[i][j][k]+=f[i-1][j+2][k-2]*(((j+2)*(j+1))/2);
                    //一个放在有棋子的列,一个放在没有棋子的列
                    if(k>=1)
                        f[i][j][k]+=f[i-1][j][k-1]*j*(m-j-k+1);
                    //都放在没有棋子的列
                    if(j>=2)
                        f[i][j][k]+=f[i-1][j-2][k]*C(m-j-k+2);
                    f[i][j][k]%=mod;
                }
        for(int i=0; i<=m; i++)
            for(int j=0; j<=m; j++)
                (ans+=f[n][i][j])%=mod;
        cout<<(ans+mod)%mod<<endl;
    }
    signed main()
    {
        int t=1;
        while(t--)
            solve();
        return 0;
    }