BZOJ4866 Ynoi2017由乃的商场之旅（莫队）

　　显然能重排为回文串相当于出现次数为奇数的字母不超过一个。考虑莫队，问题在于如何统计添加/删除一位的贡献。将各字母出现次数奇偶性看做二进制数，做一个前缀和一个后缀和。在右端添加一位时，更新区间的前缀、后缀集合，将后缀集合中合法的统计入答案即可。另三种更新类似。复杂度O(26*n√n)。

　　但是这个毒瘤题他不仅卡常还卡空间。捣鼓了半天最后欢快地弃疗了。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define ll long long
#define N 60010
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<''||c>'')) c=getchar();return c;}
int gcd(int n,int m){return m==?n:gcd(m,n%m);}
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
int n,m,block,a[N],pre[N],suf[N];
unsigned int ans[N];
unsigned short cntpre[<<],cntsuf[<<];
char s[N];
struct data
{
    int k,l,r,i;
    bool operator <(const data&a) const
    {
        return k<a.k||k==a.k&&(k&?r>a.r:r<a.r);
    }
}q[N];
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj4866.in","r",stdin);
    freopen("bzoj4866.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read(),m=read();block=*sqrt(n);
    scanf("%s",s+);
    for (int i=;i<=n;i++) a[i]=<<s[i]-'a';
    for (int i=;i<=n;i++) pre[i]=pre[i-]^a[i];
    for (int i=n;i>=;i--) suf[i]=suf[i+]^a[i];
    for (int i=;i<=m;i++) q[i].l=read(),q[i].r=read(),q[i].k=q[i].l/block,q[i].i=i;
    sort(q+,q+m+);
    int l=,r=;
    for (int i=;i<=m;i++)
    {
        ans[q[i].i]=ans[q[i-].i];
        while (r<q[i].r)
        {
            ++r,cntpre[pre[r-]]++,cntsuf[suf[r+]]++;
            ans[q[i].i]+=cntpre[pre[r]];
            for (int j=;j<;j++) ans[q[i].i]+=cntpre[pre[r]^(<<j)];
        }
        while (l>q[i].l)
        {
            --l,cntpre[pre[l-]]++,cntsuf[suf[l+]]++;
            ans[q[i].i]+=cntsuf[suf[l]];
            for (int j=;j<;j++) ans[q[i].i]+=cntsuf[suf[l]^(<<j)];
        }
        while (r>q[i].r)
        {
            ans[q[i].i]-=cntpre[pre[r]];
            for (int j=;j<;j++) ans[q[i].i]-=cntpre[pre[r]^(<<j)];
            cntpre[pre[r-]]--,cntsuf[suf[r+]]--,r--;
        }
        while (l<q[i].l)
        {
            ans[q[i].i]-=cntsuf[suf[l]];
            for (int j=;j<;j++) ans[q[i].i]-=cntsuf[suf[l]^(<<j)];
            cntpre[pre[l-]]--,cntsuf[suf[l+]]--,l++;
        }
    }
    for (int i=;i<=m;i++) printf("%u\n",ans[i]);
    return ;
}