多校1005 HDU5785 Interesting (manacher)

 // 多校1005 HDU5785 Interesting
 // 题意：给你一个串，求相邻两个回文串左边端点*右边端点的和
 // 思路：马拉车算出最长回文半径，求一个前缀和，既得到每个点对答案的贡献。
 // ans=L[i]*R[i-1]
 // L[i] 以i开始的所有回文串结尾坐标的和
 // R[i] 以i结尾的所有回文串开始坐标的和
 // 这题我们关键是求L[] R[]
 // 另开两个数组a[] b[] 分别记录当前点对答案贡献 和每个点出现的次数

 // #pragma comment(linker, "/STACK:102c000000,102c000000")
 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <sstream>
 #include <string>
 #include <algorithm>
 #include <list>
 #include <map>
 #include <vector>
 #include <queue>
 #include <stack>
 #include <cmath>
 #include <cstdlib>
 // #include <conio.h>
 using namespace std;
 #define clc(a,b) memset(a,b,sizeof(a))
 #define inf 0x3f3f3f3f
 #define lson l,mid,rt<<1
 // #define rson mid+1,r,rt<<1|1
 const int N = 1e6+;
 const int MOD = 1e9+;
 #define LL long long
 #define LB long double
 // #define mi() (l+r)>>1
 double const pi = acos(-);
 const double eps = 1e-;
 void fre(){freopen("in.txt","r",stdin);}
 inline int read(){int x=,f=;char ch=getchar();while(ch>''||ch<'') {if(ch=='-') f=-;ch=getchar();}while(ch>=''&&ch<='') { x=x*+ch-'';ch=getchar();}return x*f;}

 char str[N];
 char s[N<<];
 int Len[N<<];
 void Manacher(int n){
     clc(s,'#');
     for(int i=;i<n;i++)
         s[i*+]=str[i];
         int len=n*+;s[]='$';
         int mx=,po=;
         for(int i=;i<=len;i++)  {
             if(mx>i)
                 Len[i]=min(Len[po*-i],mx-i);
             else Len[i]=;
             while(s[i+Len[i]]==s[i-Len[i]])  Len[i]++;
             if(i+Len[i]>mx) { mx=i+Len[i];po=i;}
         }
 }

 LL a[N<<],b[N<<];
 LL L[N],R[N];
 int main(){
     while(~scanf("%s",str)){
         int len=strlen(str);
         Manacher(len);
         len=(len<<)+;
         for(int i=;i<=len;i++) a[i]=,b[i]=;
         for(int i=len;i>=;i--){
             a[i-Len[i]+]+=i;
             a[i+]-=i;
             b[i-Len[i]+]++;
             b[i+]--;
         }
         LL a1=,a2=;
         for(int i=;i<=len;i++){
             a1+=a[i],a2+=b[i];
             if(i%==)
             L[i/]=(a1-a2*i/)%MOD;
         }

         for(int i=;i<=len;i++) a[i]=,b[i]=;
         a1=,a2=;
         for(int i=len;i>=;i--){
             a[i]+=i;
             a[i+Len[i]]-=i;
             b[i]++;
             b[i+Len[i]]--;
         }
         for(int i=;i<=len;i++){
             a1+=a[i],a2+=b[i];
             if(i%==)
             R[i/]=(a1-a2*i/)%MOD;
         }

         LL ans=;
         for(int i=;i<=len/;i++){
             ans=(ans+L[i]*R[i-]%MOD)%MOD;
         }
         printf("%I64d\n",ans);
     }
     return ;
 }