代码短的惨绝人寰,思维复杂度高的惨绝人寰。
首先很容易抽象出模型,求长度为n的抖动子序列个数,然后我们开始蛋疼。
状态 f[i][j]表示以j开头的长度为i的,并且第一位是下降的抖动子序列个数。显然的是,实际上的答案是个数*2,因为这样并没有算上第一位是上升的。
状态转移方程:
(1)首先我们考虑开头为j-1的抖动子序列个数,显然这样的个数为f[i][j-1];
(2)然后我们考虑开头为j的抖动子序列个数,那么我们就要求第二位x一定要比j小,即x=j-1,但显然这样的序列第一位上升,我们需要把它翻转过来,j-1就变成了(i-1)-(j-1)+1=i-j+1,这样的答案就是f[i-1][i-j+1]。
f[i][j]=f[i][j-1]+f[i-1][i-j+1]
初始化 f[2][2]=1;
蛋疼结束了么?远远没有,答案其实是f[n+1][n+1]*2,因为只有这样才能使第一位放小于n+1的[1,n]。同时必须使用滚动数组,否则MLE;code:
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n,p,f[2][5001];
int main()
{
int ans,i,j,x;
scanf("%d%d",&n,&p);
f[0][2]=1;
for (i=3;i<=n+1;++i)
{
x=i&1;
for (j=1;j<=i;++j)
f[x][j]=(f[!x][i-j+1]+f[x][j-1])%p;
}
ans=(f[x][n+1]*2)%p;
printf("%d\n",ans);
}
然而我们在很多神犇的博客里看到的题解是这样的(by hzwer学长)
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define ll long long
#define inf 1000000000
using namespace std;
inline ll read()
{
ll x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
int n,p;
int f[2][4505];
int main()
{
n=read();p=read();
f[1][1]=1;
for(int i=2;i<=n;i++)
for(int j=1;j<=n;j++)
{
int x=i&1;
f[x][j]=f[x][j-1]+f[x^1][i-j];
f[x][j]%=p;
}
printf("%d",n==1?1%p:f[n&1][n]*2%p);
return 0;
}
其实道理也很简单,黄学长把所有的i都向左平移了一位。。。