P3501 [POI2010]ANT-Antisymmetry

P3501 [POI2010]ANT-Antisymmetry

二分+hash

注意：答案超出int范围

------------

先拿一个反对称串来做栗子：010101

我们可以发现 0101（左边右边各削掉1个）,01（左边右边各削掉2个）都是反对称串

多举几个例子，我们可以总结出一个性质：一个反对称串的所有同中心的子串都是反对称串

∴长为 n 的子串中的反对称子串数= n/2

------------
于是我们就可以设计算法了

每次枚举中心点，然后二分查找反对称串的最长长度。

对于反对称的问题，我们可以将主串正序逆序都计算一遍hash值。为了方便逆序的可以直接取反

每次判断时将子串取出，就可以达O(1)判断

复杂度O(nlogn)

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef unsigned long long ull;
inline int min(int &a,int &b) {return a<b ?a:b;}
const int base=;
char q[]; int n;
ull h1[],h2[],fac[],ans; //自然溢出hash
int main(){
    scanf("%d",&n); fac[]=;
    scanf("%s",q);
    for(int i=;i<=n;++i){
        h1[i]=h1[i-]*base+(q[i-]=='');
        fac[i]=fac[i-]*base; //通用取串方法：利用fac数组取出子串hash值，其中fac[i]=base^i
    }
    for(int i=n;i>=;--i) h2[i]=h2[i+]*base+(q[i-]==''); //逆序计算（直接取反）
    for(int i=;i<n;++i){
        int l=,r=min(i,n-i),mid;
        while(l<r){ //二分查找长度的一半
       　　 mid=l+((r-l)>>)+;
       　　 ull p1=h1[i+mid]-h1[i-mid]*fac[mid*];
      　　  ull p2=h2[i-mid+]-h2[i+mid+]*fac[mid*]; //取出子串
      　　  if(p1==p2) l=mid;
      　　  else r=mid-;
        }
        ans+=l;
    }cout<<ans;
    return ;
}