[NOI 2016] 优秀的拆分

传送门

暴力做法：

作为NOI 2016的第一题，本题暴力分非常良心，足足有95之多。
前19个点 N<=2000 $N^2$的算法是可以通过的。
首先我们 $N^2$处理出 f[i] 表示在以i结尾处之前满足 AA的子串个数
然后在 $N^2$处理处 g[i] 表示在于i结尾处之前满足 AABB的子串个数
刚开始求区间hash值也用了$N^2$的方法，前19个点用了2000ms+
后来想起来求区间hash值只需要O(N)处理+O(1)查询就可以了 前19个点用了 400ms+
刚开始写的时候一直错，因为我得到区间哈希值的时候int撑不下了，用了unsigned long long才OK

#include <cstdio>
#include <iostream>
#include <cstring>
#define ull unsigned long long
using namespace std;
const int mod1=100271,mod2=500179;
char s[3000];
int f[3000],g[3000];
ull hash1[3000],hash2[3000];
ull base1[3000],base2[3000];
void make_hash(int n)
{
    base1[0]=base2[0]=1;
    for(int i=1;i<=n;i++)
    {
        base1[i]=base1[i-1]*28%mod1;
        base2[i]=base2[i-1]*28%mod2;
        hash1[i]=(hash1[i-1]*28%mod1+(s[i]-'a'+1))%mod1;
        hash2[i]=(hash2[i-1]*28%mod2+(s[i]-'a'+1))%mod2;
    }
}
ull gethash1(int l,int r)
{
    return (hash1[r]-(hash1[l-1]*base1[r-l+1])%mod1+mod1)%mod1;
}
ull gethash2(int l,int r)
{
    return (hash2[r]-(hash2[l-1]*base2[r-l+1])%mod2+mod2)%mod2;
}
void work()
{
    memset(hash1,0,sizeof(hash1)),memset(hash2,0,sizeof(hash2));
    memset(f,0,sizeof(f)),memset(g,0,sizeof(g));
    scanf("%s",(s+1));
    int len=strlen(s+1);
    make_hash(len);
    for(int i=1;i<=len;i++)
    {
        int j;
        i%2==1?j=2:j=1;
        for(;j<i;j+=2)
        {
            int x=i-j+1;
            x=j+x/2-1;
            if(gethash1(j,x)==gethash1(x+1,i)&&gethash2(j,x)==gethash2(x+1,i)) f[i]++;
        }
    }
    int ans=0;
    for(int i=3;i<=len;i++)
    {
        int v=f[i-1],x;
        for(int j=i+1;j<=len;j+=2)
        {
            x=j-i+1;x=i+x/2-1;
            if(gethash1(i,x)==gethash1(x+1,j)&&gethash2(i,x)==gethash2(x+1,j)) g[j]+=v;
        }
    }
    for(int i=4;i<=len;i++) ans+=g[i];
    printf("%d\n",ans);
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
     work();
}

正解：我这只蒟蒻暂时还不会

好像得用后缀的数据结构，然而我还没学。
日后补坑~
不过95分已经非常良心了，233