【BZOJ 1925】【SDOI 2010】地精部落

代码短的惨绝人寰，思维复杂度高的惨绝人寰。

首先很容易抽象出模型，求长度为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都向左平移了一位。。。