hdu 6377 : 度度熊看球赛

题目链接

题解：

• 将原问题转换为 对于全部 (2n)! 种情况，每种情况对ans的贡献为 D^k，其中k表示该情况下有k对情侣座位相邻。

• 预处理好共有 i (1<=i<=N)对情侣时，出现 j (0<=j<=i) 对情侣坐在一起时情况数，用dp[i][j]记录
• 初始条件为dp[1][1]=2
• 当总情侣对数由 i 向 i+1 转移时，j 有四种变化的可能

1. j-1 → j : 新来的情侣“一起”插入到不打断先前情侣的“相邻”座位处 dp[i+1][j] += 2 * dp[i][j-1] * (2*i+1-(j-1))
2. j → j : 新来的情侣“一起”插入到打断先前某对情侣的“相邻”座位处 dp[i+1][j] += 2 * dp[i][j] * j ，或者“分开”插入到不打断先前情侣的“不相邻”座位处  dp[i+1][j] += 2 * dp[i][j] * C(2*i+1-j,2)
3. j+1 → j : 新来的情侣，一个插入到打断先前某对情侣的座位处，另一个插入到不打断先前情侣的座位处 dp[i+1][j] += 2 * dp[i][j+1] * (j+1)*(2*i+1-(j+1)
4. j+2 → j : 新来的情侣“分开”插入到打断先前某两对情侣的“相邻”座位处 dp[i+1][j] += 2 * dp[i][j+2] * C(j+2,2)

• 以上为预处理过程

• 最后 ans += d^i * dp[n][i] ，i=0→n

代码如下

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;

const LL mod=;

LL qpow(LL x,LL n)        //求x^n%mod
{
    LL ret=;
    for(; n; n>>=)
    {
        if(n&) ret=ret*x%mod;
        x=x*x%mod;
    }
    return ret;
}

int C(LL a,LL b=)
{
    return a*(a-)%mod*qpow(b,mod-)%mod;
}

const int N=1e3+;
LL dp[N][N];

void init()
{
    dp[][]=;
    for(int i=;i<N;i++)
    {
        dp[i+][]=*dp[i][]%mod*C(*i+,)%mod+
                   *dp[i][]%mod**i%mod+
                   *dp[i][]%mod;
        dp[i+][]%=mod;
        for(int j=;j<=i+;j++)
        {
            dp[i+][j]=*dp[i][j-]%mod*(*i+-(j-))%mod+
                       *dp[i][j]%mod*(j+C(*i+-j,))%mod+
                       *dp[i][j+]%mod*((j+)*(*i+-(j+))%mod)%mod+
                       *dp[i][j+]%mod*C(j+,)%mod;
            dp[i+][j]%=mod;
        }
    }
}

int main()
{
    init();
    LL n,d;
    while(~scanf("%lld%lld",&n,&d))
    {
        LL ret=;
        if(n==)
        {
            printf("%lld\n",*d%mod);
            continue;
        }
        for(int i=;i<=n;i++) ret=(ret+dp[n][i]*qpow(d,i)%mod)%mod;
        printf("%lld\n",ret);
    }
}