BZOJ 1200: [HNOI2005]木梳

虐哭了QAQ

首先简单分析一下会发现将最终的序列游程编码之后不会出现连续的长度为3的递增或递减数列,然后就弃疗了。

于是只好去翻题解……卧槽这贪心是什么鬼。

看了半天证明感觉有点方,于是决定把代码敲(抄)一遍就算了。

其实我觉得既然复杂度是n*c^2*d^2,不如把c和d设大一点,反正也证不出来。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int N=100000+5;
const ll inf=1LL<<60;
int b[N][5][3],a[N];
ll f[N][5][3][2];
int main(){
	//freopen("a.in","r",stdin);
	int n;scanf("%d",&n);
	for(int i=1;i<=n;i++)scanf("%d",&a[i]);
	memset(f,122,sizeof(f));
	for(int j=0;j<=2;j++)
	if(1+j<=n)
	for(int k=-1;k<=1;k++)
	if(a[1+j]+k<=a[1]){
		b[1][j+2][k+1]=a[1+j]+k;
		f[1][j+2][k+1][0]=a[1]-b[1][j+2][k+1];
		f[1][j+2][k+1][1]=a[1]-b[1][j+2][k+1];
	}
	for(int i=2;i<=n;i++)
	for(int j=-2;j<=2;j++)
	if(i+j<=n&&i+j>0)
	for(int k=-1;k<=1;k++)
	if(a[i+j]+k<=a[i]){
		b[i][j+2][k+1]=a[i+j]+k;
		int d=a[i]-b[i][j+2][k+1];
		for(int j1=-2;j1<=2;j1++)
		for(int k1=-1;k1<=1;k1++)
		if(b[i-1][j1+2][k1+1]>b[i][j+2][k+1])
		f[i][j+2][k+1][0]=min(f[i][j+2][k+1][0],f[i-1][j1+2][k1+1][1]+d);
		else if(b[i-1][j1+2][k1+1]<b[i][j+2][k+1])
		f[i][j+2][k+1][1]=min(f[i][j+2][k+1][1],f[i-1][j1+2][k1+1][0]+d);
		else{
			f[i][j+2][k+1][0]=min(f[i][j+2][k+1][0],f[i-1][j1+2][k1+1][0]+d);
			f[i][j+2][k+1][1]=min(f[i][j+2][k+1][1],f[i-1][j1+2][k1+1][1]+d);
		}
	}
	ll ans=inf;
	for(int j=-2;j<=2;j++)
	for(int k=-1;k<=1;k++)
	ans=min(ans,min(f[n][j+2][k+1][0],f[n][j+2][k+1][1]));
	printf("%lld\n",ans);
	return 0;
}


你可能感兴趣的:(BZOJ 1200: [HNOI2005]木梳)