[BZOJ]3790 神奇项链 Manacher+树状数组

3790: 神奇项链

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 548  Solved: 286
[Submit][Status][Discuss]

Description

母亲节就要到了，小 H 准备送给她一个特殊的项链。这个项链可以看作一个用小写字

母组成的字符串，每个小写字母表示一种颜色。为了制作这个项链，小 H 购买了两个机器。第一个机器可以生成所有形式的回文串，第二个机器可以把两个回文串连接起来，而且第二个机器还有一个特殊的性质：假如一个字符串的后缀和一个字符串的前缀是完全相同的，那么可以将这个重复部分重叠。例如：aba和aca连接起来，可以生成串abaaca或 abaca。现在给出目标项链的样式，询问你需要使用第二个机器多少次才能生成这个特殊的项链。 

Input

输入数据有多行，每行一个字符串，表示目标项链的样式。 

Output

多行，每行一个答案表示最少需要使用第二个机器的次数。 

Sample Input

abcdcba
abacada
abcdef

Sample Output

0
2
5

HINT

每个测试数据，输入不超过 5行 

每行的字符串长度小于等于 50000 

Source

[Submit][Status][Discuss]

HOME Back

题解

        这道题当我们用Manacher求出极大回文子串之后，问题，也就变成了询问最小线段数覆盖全区间，注意这里线段可以有交集.我们考虑dp贪心的思想，每个点的最优策略用值域树状数组来维护，f[i]表示1到i这个区间最少用多少线段覆盖，f[i]=min{f[j]} +1，其中f[j]的右端点要大于等于i的左端点，这样才能覆盖.

#include<stdio.h>#include<cstring>#include<algorithm>using namespace std;const int maxn=100010;const int inf=1000000010;char ss[maxn],s[maxn*2];int c[maxn],pal[maxn*2],m,n,ans,num;struct query{    int l,r;}q[maxn];inline bool cmp(query a,query b){return a.r<b.r;}inline void init(){    num=0;n=strlen(ss+1);m=2*n+1;ans=inf;    for(int i=1;i<=n;i++) s[i*2]=ss[i],s[i*2+1]='#';    s[0]='+',s[1]='#',s[m+1]='-';}inline void add(int l,int r){    l=l/2+1,r=r/2-1;    if(l>r) return;    q[++num].l=l,q[num].r=r;}inline void manacher(){    int mx=0,id=0;    for(int i=1;i<=m;i++){       if(mx>i) pal[i]=min(mx-i,pal[2*id-i]);       else pal[i]=1;       while(s[i-pal[i]]==s[i+pal[i]]) pal[i]++;       add(i-pal[i],i+pal[i]);   if(i+pal[i]>mx) mx=i+pal[i],id=i;}}inline void modify(int i,int x){    if(!i) return;    for(;i;i-=i&-i) c[i]=min(c[i],x);}inline int query(int i){    if(!i) return 0;    int tmp=inf;    for(;i<=n;i+=i&-i) tmp=min(tmp,c[i]);    return tmp;}int main(){    while(scanf("%s",ss+1)!=EOF){       init();       manacher();       for(int i=1;i<=n;i++) c[i]=inf;       for(int i=1;i<=num;i++){           int x=query(q[i].l-1)+1;           modify(q[i].r,x);           if(q[i].r==n) ans=min(ans,x);   }   printf("%d\n",ans-1);}}