hdu6377 度度熊看球赛

链接

http://acm.hdu.edu.cn/showproblem.php?pid=6377

题解

$f_{ij}$表示$i$对情侣进行排列，有$j$对坐在一起的方案数
新加入的一对，可能转移到$f_{i,j+1},f_{i,j},f_{i,j-1},f_{i,j-2}$

代码

#include<bits/stdc++.h>
#define maxn 1010
#define linf (1ll<<60)
#define iinf 0x3f3f3f3f
#define eps 1e-8
#define cl(x) memset(x,0,sizeof(x))
#define mod 998244353ll
using namespace std;
typedef long long ll;
ll N, D, dp[maxn][maxn];
ll read(ll x=0)
{
    ll c, f(1);
    for(c=getchar();!isdigit(c);c=getchar())if(c=='-')f=-f;
    for(;isdigit(c);c=getchar())x=x*10+c-48;
    return f*x;
}
ll fastpow(ll a, ll b)
{
    ll t(a), ans(1);
    for(;b;b>>=1,t=t*t%mod)if(b&1)ans=ans*t%mod;
    return ans;
}
int main()
{
    ll i, j, k, ans(0);
    dp[1][1]=2;
    for(i=1;i<=1000;i++)for(j=0;j<=i;j++)
    {
        auto f=dp[i], g=dp[i+1];
        (g[j+1]+=f[j]*(2*i+1-j)*2)%=mod;
        (g[j]+=f[j]*j*2+f[j]*(2*i+1-j)*(2*i-j))%=mod;
        (g[j-1]+=f[j]*j*(2*i+1-j)*2)%=mod;
        (g[j-2]+=f[j]*(j)*(j-1))%=mod;
    }
    while(~scanf("%lld%lld",&N,&D))
    {
        ans=0;
        for(i=0;i<=N;i++)(ans+=fastpow(D,i)*dp[N][i])%=mod;
        printf("%lld\n",ans);
    }
	return 0;
}