【BZOJ】【P3790】【神奇项链】【题解】【manacher+dp】

传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=3790

这题数据范围有问题,字符串长度我开到10^5才过,感谢lzr神犇的提醒

首先manacher,计算出所有的极长回文子串,问题转化为给定一些线段,用最少线段可重叠的覆盖整个区间,BIT优化dp即可

按右端点排序,f[i]=min{f[j],seg[j].sec>=seg[i].fst-1}+1

(好像倒着贪心可以做到线性?)

(哦还需要排序)

(基数排序……)

Code:

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<climits>
#define fst first
#define sec second
#define lw(x) (x&-x)
using namespace std;
const int maxn=100010;
typedef pair<int,int> pi;
char s[maxn],str[maxn<<1];
int n,m,f[maxn<<1],size;
pi seg[maxn<<1];
bool cmp(pi a,pi b){return a.sec!=b.sec?a.sec<b.sec:a.fst<b.fst;}
void get_seg(int l,int r){
	if(l>r)return;
	seg[++size]=pi(l,r);
}
int d[maxn];
int Qmin(int x){int o=d[x];while(x&&x<=n)o=min(o,d[x]),x+=lw(x);return o;}
void updata(int x,int f){while(x)d[x]=min(d[x],f),x-=lw(x);}
int main(){
	while(~scanf("%s",s+1)){
		memset(f,0,sizeof f);
		memset(d,0x7f,sizeof d);d[0]=0;
		memset(str,0,sizeof str);
		n=strlen(s+1);m=n*2+1;size=0;
		str[0]='+';str[1]='#';str[m+1]='-';
		for(int i=1;i<=n;i++)str[i<<1]=s[i],str[i<<1|1]='#';
		int k=1;f[1]=1;
		for(int i=2;i<=m;i++){
			if(f[k]+k>i)f[i]=min(f[k+k-1],f[k]+k-i);else f[i]=1;
			while(str[i-f[i]]==str[i+f[i]])f[i]++;
			get_seg((i-f[i]+2)/2,(i+f[i]-2)/2);
			if(f[i]+i>f[k]+k)k=i;
		}sort(seg+1,seg+1+size,cmp);
		memset(f,0,sizeof f);int ans=INT_MAX;
		for(int i=1;i<=size;i++){
			f[i]=Qmin(seg[i].fst-1)+1;
			if(seg[i].sec==n)ans=min(ans,f[i]);
			updata(seg[i].sec,f[i]);
		}cout<<ans-1<<endl;
	}
	return 0;
}