BZOJ1524: [POI2006]Pal

1524: [POI2006]Pal

Time Limit: 5 Sec  Memory Limit: 357 MB
Submit: 308  Solved: 101
[Submit][Status]

Description

给出n个回文串s1, s2, …, sn 求如下二元组(i, j)的个数 si + sj 仍然是回文串 规模 输入串总长不超过2M bytes

Input

The first line of input file contains the number of strings n. The following n lines describe each string: The i+1-th line contains the length of the i-th string li, then a single space and a string of li small letters of English alphabet. You can assume that the total length of all strings will not exceed 2,000,000. Two strings in different line may be the same.

Output

Print out only one integer, the number of palindromes

Sample Input

6
2 aa
3 aba
3 aaa
6 abaaba
5 aaaaa
4 abba

Sample Output

14

HINT

Source

题解：
说一下做这题的艰辛过程。。。
刚开始看见题画了画图发现好像短串必须是长串的前缀，然后就开开心心的打程序，然后就WA了。。。发现 a 和 aba  显然不能构成回文串。。。
然后又想  发现短串好像需要不重叠覆盖长串，然后又开始写程序，然后又开始WA，后来发现 aa 和 aaa 能构成回文串。。。
无奈之下请教vfleaking，然后发现了这样的算法：
字典序hash。
我们先把所有串插入一个trie树，然后统计每个串的hash，然后再枚举每个串，沿trie树向下走，枚举每一个前缀，判断他们俩连起来的字符串正着和反着是否一样，
直接hash判断即可。
代码写起来不容易，我的代码快垫底了。。。
代码：

#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 2000000+10

#define maxm 2000000

#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
#define base 131

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,tot,t[maxn][27],g[maxn],len[maxm];
ll f[maxn],h[maxn],ans,ha[maxm];
char s[maxn];
string st[maxm];

int main()

{

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

    freopen("output.txt","w",stdout);
    h[0]=1;
    for1(i,maxn)h[i]=h[i-1]*base;

    n=read();ans=-n;
    for1(k,n)
    {
        len[k]=read();
        scanf("%s",s+1);st[k]=s+1;
        int now=0;ll hash=0;
        for1(i,len[k])
         {
             int x=s[i]-'a'+1;
             if(!t[now][x])t[now][x]=++tot;
             now=t[now][x];
             hash=hash*base+x;
         }
        ha[k]=hash; 
        f[now]=k;g[now]++;
    }
    for1(k,n)
    {
        int now=0; 
        for0(i,len[k]-1) 
         {
             int x=st[k][i]-'a'+1;
             now=t[now][x];
             if(g[now]&&ha[f[now]]*h[len[k]]+ha[k]==ha[k]*h[i+1]+ha[f[now]])ans+=(ll)g[now]*2;
         }
    }
    printf("%lld\n",ans);

    return 0;

}