Codeforces932G Palindrome Partition【回文树DP（附border的若干证明）】

洛谷链接

转化为求偶回文串划分的过程就不细说了。

PS：下面的图截自skyh的PPT 水题再选讲，先Orz为敬。

弱周期引理

长度大于等于一半的串的匹配形成等差数列

根据2|u|>=|v|应用弱周期引理得到周期之间的整除关系从而推出等差。

长度大于等于|S|/2的border长度形成等差数列

回文串的回文后缀(border)长度可以表示成log个等差数列

本题做法

考虑这么做的正确性，（以下只讨论$len[fail[x]]\ge len[x]/2$的情况，若小于$1/2$显然构不成等差数列，直接继承是对的）。
我们继承$g[fail[x]]$，实际上是想继承 $i-d$ 位置的值，而$fail[x]$一定对应第 $i-d$ 位置吗？
如果在 $(i-d,i)$ 这个区间中出现了 $fail[x]$ 对应的串，由引理1，$fail[x]$ 这个串在 $x$ 中的匹配一定构成等差数列，而 $fail[x]$ 在 $x$ 的开头出现了，所以可以找到 $x$ 的一个比 $len[x]-len[fail[x]]$ 更小的周期，这与 $fail[x]$ 是最大border矛盾，所以继承不会出错。

而引理2,3则保证了本题的复杂度为$O(n\log n)$

Upd：因为等差数列是跳过去的，所以我们每次只更新了等差数列开头的 $g$，而没有更新中间的，于是我们要保证 $g[fail[x]]$ 在 $i-d$ 处一定是以等差数列的开头项出现的，这样值才是对的。而如果不是开头项，说明对于 $i$ 来说还有更长的回文后缀，与其最大矛盾。

顺带一提，为了方便转移，$g[x]$实际上是存的当前等差数列除了开头一项的贡献，即：

Code：

#include<bits/stdc++.h>
#define maxn 1000005
using namespace std;
const int mod = 1e9+7;
int n,f[maxn],g[maxn];
char s[maxn],s2[maxn];
namespace PAM{
	int fa[maxn]={1,1},ch[maxn][26],len[maxn]={0,-1},last,sz=1,anc[maxn],dif[maxn];
	void extend(int x){
		int p=last,c=s[x]-'a';
		while(s[x]!=s[x-len[p]-1]) p=fa[p];
		if(!ch[p][c]){
			int q=fa[p]; len[++sz]=len[p]+2;
			while(s[x]!=s[x-len[q]-1]) q=fa[q];
			fa[sz]=ch[q][c], dif[sz]=len[sz]-len[fa[sz]];
			anc[sz]= dif[sz]==dif[fa[sz]] ? anc[fa[sz]] : fa[sz];
			ch[p][c]=sz;
		}
		last=ch[p][c];
	}
}
using namespace PAM;
int main()
{
	scanf("%s",s+1),n=strlen(s+1);
	for(int i=1;i<=n;i+=2) s2[i]=s[(i+1)/2];
	for(int i=2;i<=n;i+=2) s2[i]=s[n-i/2+1];
	for(int i=1;i<=n;i++) s[i]=s2[i];
	f[0]=1;
	for(int i=1;i<=n;i++){
		extend(i);
		for(int x=last;x;x=anc[x]){
			g[x]=f[i-len[anc[x]]-dif[x]];
			if(anc[x]!=fa[x]) g[x]=(g[x]+g[fa[x]])%mod;
			if(!(i&1)) f[i]=(f[i]+g[x])%mod;
		}
	}
	printf("%d\n",f[n]);
}