[BZOJ3790]神奇项链-manachar

last update at 2017.11.6

1.修改了少量口胡的语句
2.删除了多嘴的内容
3.规范了板块和格式

---

说在前面

= =这博客写着写着电脑蓝屏了，还好写的不是很多，不伤心不伤心…

---

题目：

BZOJ3790传送门

题面

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

输入和输出

Input：
输入数据有多行，每行一个字符串，表示目标项链的样式。
Output：
多行，每行一个答案表示最少需要使用第二个机器的次数。

数据范围

每个测试数据，输入不超过 5行
每行的字符串长度小于等于 50000

---

解法

感觉好像很难的样子，然而想一想发现，第一个机器其实并没有什么卵用，重点在第二个机器。
如果我们能够算出这个字符串里所有的回文串（因为他们都可以由机器1生成出来），然后把他们当成线段处理（回文串的开始与末尾分别为左右端点），求出如何用最少的线段覆盖整个区域，就成了一个线段覆盖问题。

至于线段覆盖怎么算= =….
自己查去
（大概就是O（n）for一遍，每次选取已扩展区域中，右端点最右的线段，然后更新当前区域右端点）

---

贴代码：

解释一些东西：

p数组的含义是：

如果你不知道，快去学manachar！

lin数组的含义是：以i为起点的回文串长度，最大是lin[i]，相当于是存线段。
maxx变量的含义是：当前最右位置，用于更新区域最右位置i
manachar_()函数内，变量pos的含义是：当前回文串起点为pos位置处。

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;

char s[50005],str[100005];
int len,p[100005],mx,id,topp=0;
int lin[50005],maxx=0,ans=0;

void manachar_(){
    for(int i=0;i<=topp;i++){
        if(mx > i)
            p[i]=(p[2*id - i] < mx - i ? p[2*id - i] : mx - i);
        else
            p[i]=1;
        while(str[i+p[i]]==str[i-p[i]]) p[i]++;
        if(i+p[i]>mx){
            mx=i+p[i];
            id=i;
        }
    }
}

void dp(){
    int pos;
    for(int i=2;i<=len*2;i++){
            pos=i/2-p[i]/2+1;
            lin[pos]=max(lin[pos],p[i]-1);
    }
    for(int i=1;i<=len;){
        maxx=max(maxx,lin[i]+i-1);
        ans++;
        int tmp=maxx;
        for(int j=i+1;j<=tmp;j++)
            maxx=max(maxx,lin[j]+j-1);
        i=tmp+1;
    }
}

int main(){
    while(scanf("%s",s)!=EOF){
        p[0]=0;
        memset(lin,0,sizeof(lin));
        topp=0;mx=0;ans=0;maxx=0;
        str[topp]='$';
        str[++topp]='#';
        len=strlen(s);
        for(int i=0;i<len;i++){
            str[++topp]=s[i];
            str[++topp]='#';
        }
        manachar_();
        dp();
        printf("%d\n",ans-1);
    }
    return 0;
}