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<'0'||c>'9')) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
int read()
{
    int x=0,f=1;char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
    while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
    return x*f;
}
int n,m,block,a[N],pre[N],suf[N];
unsigned int ans[N];
unsigned short cntpre[1<<26],cntsuf[1<<26];
char s[N];
struct data
{
    int k,l,r,i;
    bool operator <(const data&a) const
    {
        return k<a.k||k==a.k&&(k&1?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=2*sqrt(n);
    scanf("%s",s+1);
    for (int i=1;i<=n;i++) a[i]=1<<s[i]-'a';
    for (int i=1;i<=n;i++) pre[i]=pre[i-1]^a[i];
    for (int i=n;i>=1;i--) suf[i]=suf[i+1]^a[i];
    for (int i=1;i<=m;i++) q[i].l=read(),q[i].r=read(),q[i].k=q[i].l/block,q[i].i=i;
    sort(q+1,q+m+1);
    int l=1,r=0;
    for (int i=1;i<=m;i++)
    {
        ans[q[i].i]=ans[q[i-1].i];
        while (r<q[i].r)
        {
            ++r,cntpre[pre[r-1]]++,cntsuf[suf[r+1]]++;
            ans[q[i].i]+=cntpre[pre[r]];
            for (int j=0;j<26;j++) ans[q[i].i]+=cntpre[pre[r]^(1<<j)];
        }
        while (l>q[i].l)
        {
            --l,cntpre[pre[l-1]]++,cntsuf[suf[l+1]]++;
            ans[q[i].i]+=cntsuf[suf[l]];
            for (int j=0;j<26;j++) ans[q[i].i]+=cntsuf[suf[l]^(1<<j)];
        }
        while (r>q[i].r)
        {
            ans[q[i].i]-=cntpre[pre[r]];
            for (int j=0;j<26;j++) ans[q[i].i]-=cntpre[pre[r]^(1<<j)];
            cntpre[pre[r-1]]--,cntsuf[suf[r+1]]--,r--;
        }
        while (l<q[i].l)
        {
            ans[q[i].i]-=cntsuf[suf[l]];
            for (int j=0;j<26;j++) ans[q[i].i]-=cntsuf[suf[l]^(1<<j)];
            cntpre[pre[l-1]]--,cntsuf[suf[l+1]]--,l++;
        }
    }
    for (int i=1;i<=m;i++) printf("%u\n",ans[i]);
    return 0;
}