[APIO2014]回文串-（回文自动机裸题）

传送门
没看过题目的可以先看一下，题目要我们求的是回文串长度len*出现次数的最大值。
很符合回文自动机的性质的题目，可以用SAM+manacher做，也不难，但是对回文自动机来说就是裸题。
1，本质不同回文子串长度len数组。
2，出现次数可以用拓扑之后向前逆推求
（UPD：其实并不需要拓扑，因为回文数的节点本来有有先后顺序。之前做SAM习惯了）

总之就是很板子

#include <bits/stdc++.h>
using namespace std;
#define int long long
const int MAXN = 500000 , MAXK = 26;
char str[ MAXN + 5 ];
struct Palindrome_Automaton{
    int Size //总节点
    , Last , Root0 , Root1 , Trans[ MAXN + 5 ][ MAXK + 5 ] , Fail[ MAXN + 5 ];
    int Len[ MAXN + 5 ];int d[MAXN];queue<int> q;
    int weight[MAXN+5];
    Palindrome_Automaton( ) {
        Root0 = Size ++ , Root1 = Size; Last = Root1;
        Len[ Root0 ] = 0  , Fail[ Root0 ] = Root1;++d[Root1];
        Len[ Root1 ] = -1 , Fail[ Root1 ] = Root1;++d[Root1];
    }
    void Extend( int ch , int dex ) {
        int u = Last;
        while (str[ dex - Len[ u ] - 1 ] != str[ dex ] ) u = Fail[ u ]; //找到合格的后缀
        if( !Trans[ u ][ ch ] ) {         //无现成的边
            int Newnode = ++ Size , v = Fail[ u ];    //防止取掉整串
            Len[ Newnode ] = Len[ u ] + 2;
            while (str[ dex - Len[ v ] - 1 ] != str[ dex ] ) v = Fail[ v ];    //给他找一个fail指针
            Fail[ Newnode ] = Trans[ v ][ ch ] ,++d[Trans[ v ][ ch ]], Trans[ u ][ ch ] = Newnode;
        }
        Last = Trans[ u ][ ch ];
        weight[Last]++;
    }

    void Build( char *str ) {
        int len = strlen( str );
        for( int i = 0 ; i < len ; i ++ ) {
            Extend( str[ i ] - 'a' + 1 , i );
        }
    }
//
//    void getWeight()       ///多此一举的拓扑
//    {
//        int ma =0 ;
//        for (int k = 0; k <= Size; ++k) {
//            if (!d[k]) q.push(k);
//        }
//        while (!q.empty())
//        {
//            int u = q.front();q.pop();
//
//            int v = Fail[u];
//            weight[v]+=weight[u];
//            --d[v];
//            if (!d[v]) q.push(v);
//        }
//        for (int j = 2; j <= Size; ++j) {
//            ma = max(ma,Len[j]*weight[j]);
//        }
//        cout << ma << endl;
//    }
    void getWeight()
    {
        long long  ma =0 ;
        for (int k = Size; k >= 0; --k) {    ///向前逆推就好
            weight[Fail[k]] += weight[k];
        }
        for (int j = 0; j <= Size; ++j) {
            ma = max(ma,Len[j]*weight[j]);
        }
        cout << ma << endl;
    }

}PAM;

signed main( ) {
    scanf("%s", str );
    PAM.Build( str );
    PAM.getWeight();
    return 0;
}