srm 711 DerangementsDiv2(dp)

johsnows

于 2017-07-02 09:31:14 发布

阅读量315

点赞数

CC 4.0 BY-SA版权

分类专栏： Topcoder acm studyproblem 动态规划文章标签： dp srm

本文链接：https://blog.csdn.net/johsnows/article/details/74090429

acm 同时被 3 个专栏收录

188 篇文章

订阅专栏

studyproblem

96 篇文章

订阅专栏

动态规划

42 篇文章

订阅专栏

本文介绍了一种使用动态规划解决错排问题的方法，具体为求1-n+m的数排列中，前m个数不等于其位置的排列数量。通过定义dp状态并给出转移方程，最终实现了对错排问题的有效求解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题意：

1-n+m的数排列，要求前m个数p[i]!=i，问有多少种排列方式，对1e9+7取模

解题思路：

不会做，最后问的p神才会了。

dp状态是dp[i][j]，表示i个数排列要求前j个数错排的方案数。

转移dp[i][j]=(i-j)*d[i-1][j]+j*dp[i-1][j-1];

分别表示最后一个数放在后n个位置还是前m个位置的两种方法。

i=j的时候就是错排公式，要单独写。

代码：

#include <bits/stdc++.h>
#define LL long long
using namespace std;
const int maxn=55;
const int mod=1e9+7;
LL fa[202];
int i, j;
LL dp[105][104];;
void init()
{
    fa[0]=1;
    for(int i=1; i<=200; i++)
    {
        fa[i]=(fa[i-1]*(LL)i)%mod;
    }
    return;
}

struct  DerangementsDiv2
{
    
    LL i, j;
    int count(int n, int m) 
    {
        init();
        memset(dp, 0, sizeof dp);
        dp[0][0]=1, dp[1][1]=0, dp[2][2]=1;

        for(i=3; i<=n+m; i++)
        {
            dp[i][i]=((i-1)*((dp[i-1][i-1]+dp[i-2][i-2]))%mod)%mod;
        }
        for(i=0; i<=n+m; i++)
        {
            dp[i][0]=fa[i];
            if(i)dp[i][1]=((i-1)*fa[i-1])%mod;
        }

        for(i=1; i<=n+m; i++)
        {
            for(j=2; j<=m; j++)
            {
                if(i==j)continue;
                dp[i][j]=(((i-j)*(dp[i-1][j]))%mod+ (j*dp[i-1][j-1])%mod)%mod;
            }
        }
        return dp[n+m][m];
        
    }
};
/*

int main()
{
    DerangementsDiv2 a;
    int x, y;
    while(~scanf("%d %d", &x, &y))
    a.count(x, y);
    
    
}
*/