这是一道很简单的dp题,但是思路比较难想
主要方法就是枚举每个点走i次走回自己的方案数
然后略微容斥一下,就可以A掉了
#include<bits/stdc++.h>
using namespace std;
int mod,n,m,c[2005][2005],w[2005],dp[2005];
int moc(int x)
{
if(x>=mod) return x-mod;
if(x<0) return x+mod;
return x;
}
int fast(int x,int y)
{
int ans=1;
while(y)
{
if(y&1) ans=1ll*ans*x%mod;
y>>=1;
x=1ll*x*x%mod;
}
return ans;
}
int main()
{
cin>>n>>m>>mod;
for(int i=0;i<=n;i++)
{
c[i][0]=1;
for(int j=1;j<=i;j++)
{
c[i][j]=moc(c[i-1][j-1]+c[i-1][j]);
}
}
w[0]=1;
for(int i=1;i<=m;i++)
{
for(int j=n/2*2;j>=0;j-=2)
{
for(int k=2;k<=j;k+=2)
{
w[j]=moc(w[j]+1ll*w[j-k]*c[j][k]%mod*c[k][k>>1]%mod);
}
}
}
dp[0]=1;
for(int i=1;i<=n;i++)
{
dp[i]=1ll*fast((m<<1)%mod,i)*(i+1)%mod;
for(int j=0;j<i;j++)
{
dp[i]=moc(dp[i]-1ll*dp[j]*w[i-j]%mod);
}
}
cout<<dp[n];
return 0;
}
/*
4 2 100000
*/