BZOJ3790: 神奇项链

3790: 神奇项链

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 79  Solved: 40
[Submit][Status]

Description

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

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

Input

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

Output

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

Sample Input

abcdcba
abacada
abcdef

Sample Output

0
2
5

HINT

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

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

题解：

我想呵呵。。。
我是这样想得，我们先manacher一遍，令f[i]表示打出1--i的最小花费
然后可以O（n）求出以每个点结尾的最长回文子串，不妨设最远扩展到x，
然后就可以用min（f[j]）（x-1<=j<=i-1） +1来更新f[i]
就可以线段树搞了。
代码：

#include<cstdio>

#include<cstdlib>

#include<cmath>

#include<cstring>

#include<algorithm>

#include<iostream>

#include<vector>

#include<map>

#include<set>

#include<queue>

#include<string>

#define inf 1000000000

#define maxn 500000+5

#define maxm 20000000+5

#define eps 1e-10

#define ll long long

#define pa pair<int,int>

#define for0(i,n) for(int i=0;i<=(n);i++)

#define for1(i,n) for(int i=1;i<=(n);i++)

#define for2(i,x,y) for(int i=(x);i<=(y);i++)

#define for3(i,x,y) for(int i=(x);i>=(y);i--)

#define mod 1000000007

using namespace std;

inline int read()

{

    int x=0,f=1;char ch=getchar();

    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}

    return x*f;

}
int n,m,a[maxn],f[maxn],p[maxn];
char s[maxn];
struct seg{int k,l,r,mi;}t[4*maxn];
inline void build(int k,int l,int r)
{
    t[k].l=l;t[k].r=r;int mid=(l+r)>>1;t[k].mi=inf;
    if(l==r){if(!l)t[k].mi=-1;return;}
    build(k<<1,l,mid);build(k<<1|1,mid+1,r);
}
inline void pushup(int k)
{
    t[k].mi=min(t[k<<1].mi,t[k<<1|1].mi);
}
inline void change(int k,int x,int y)
{
    int l=t[k].l,r=t[k].r,mid=(l+r)>>1;
    if(l==r){t[k].mi=min(t[k].mi,y);return;}
    if(x<=mid)change(k<<1,x,y);else change(k<<1|1,x,y);
    pushup(k);
}
inline int query(int k,int x,int y)
{
    int l=t[k].l,r=t[k].r,mid=(l+r)>>1;
    if(l==x&&r==y)return t[k].mi;
    if(y<=mid)return query(k<<1,x,y);
    else if(x>mid)return query(k<<1|1,x,y);
    else return min(query(k<<1,x,mid),query(k<<1|1,mid+1,y));
}

int main()

{

    freopen("input.txt","r",stdin);

    freopen("output.txt","w",stdout);

    while(scanf("%s",s)!=EOF)
    {
        n=strlen(s);
        build(1,0,n);
        for1(i,n)a[i<<1]=s[i-1];
        m=n;n=2*n+1;
        int id=0,mx=0;
        for1(i,n)
        {
            if(mx>i)p[i]=min(p[2*id-i],mx-i);
            while(i-p[i]-1>0&&i+p[i]+1<=n&&a[i-p[i]-1]==a[i+p[i]+1])p[i]++;
            if(i+p[i]>mx)mx=i+p[i],id=i;
        }
        int now=0;
        for1(i,n)if(i+p[i]>now)
        {
            for2(j,now+1,i+p[i])f[j]=2*i-j;
            now=i+p[i];
        }
        for1(i,m)change(1,i,query(1,max(0,(f[i<<1]>>1)-1),i-1)+1);
        printf("%d
",query(1,m,m));
        memset(a,0,sizeof(a));
        memset(s,0,sizeof(s));
        memset(p,0,sizeof(p));
        memset(f,0,sizeof(f));
    }

    return 0;

}  

刚开始装逼用gets结果WA了好几发感觉人生无望，这种题都A不了，就颓了一节课T_T

 UPD：orzky的题解http://blog.csdn.net/iamzky/article/details/41852799