bzoj4866 [Ynoi2017]由乃的商场之旅 莫队

Description

---

由乃有一天去参加一个商场举办的游戏。商场派了一些球王排成一行。每个人面前有几堆球。说来也巧，由乃和你
一样，觉得这游戏很无聊，于是决定换一个商场。另一个商场是Deus的，他看到由乃来了，于是想出了一个更有趣
的游戏：写数据结构题这个题是这样的：
我们的情人，不过是随便借个名字，用幻想吹出来的肥皂泡，把信拿去吧，你可以使假戏成真。我本来是无病呻吟
，漫无目的的吐露爱情—现在这些漂泊不定的鸟儿有地方栖息了，你可以从信里看出来。拿去吧—由于不是出自
真心，话就说得格外动听，拿去吧，就这么办吧…由于世界会在7月20日完结，作为救世主，间宫卓司要在19日让
所有人回归天空现在已经是19日傍晚，大家集合在C栋的天台上，一共n个人在他们面前，便是终之空，那终结的天
空
回归天空是一件庄重的事情，所以卓司决定让大家分批次进行，给每个人给了一个小写字母’a’->'z’作为编号一个
区间的人如果满足他们的编号重排之后可以成为一个回文串，则他们可以一起回归天空，即这个区间可以回归天空
由于卓司是一个喜欢妄想的人，他妄想了m个区间，每次他想知道每个区间中有多少个子区间可以回归天空因为世
界末日要来了，所以卓司的信徒很多
由乃天天做数据结构已经快恶心死了，于是让您帮她做当然，您天天做数据结构题，肯定也觉得恶心死了，不过可
爱即正义！所以这个题还是要做的~

n,m<=60000

Solution

---

这种不强制在线的区间询问看起来就像莫队╮(╯▽╰)╭

注意到一个区间会产生1的贡献当且仅当所有字母出现偶数次或只有一个出现奇数次，那么就可以给每一个字符赋一个权值，我们只需要判断是否存在异或和为0或字符的权值就可以了
然后就是非常科学的复杂度了，但是喜闻乐见要卡常，于是我弃了也懒得拍了。。。

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#define rep(i,st,ed) for (register int i=st;i<=ed;++i)

const int N=80005;
const int B=1005;

struct Q {int l,r,id;} q[N];

int a[N],v[N],prt[N],ans;
int rec[33554435*2],sum,suf;

char str[N];

int read() {
	int x=0,v=1; char ch=getchar();
	for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
	for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
	return x*v;
}

void ins(int x,int opt) {
	int w=v[a[x]]; sum^=w;
	if (opt==1) {
		rec[suf]++; suf^=w;
		rep(i,0,25) ans+=rec[suf^v[i]];
		ans+=rec[suf];
	} else {
		rec[suf]++; suf^=sum;
		rep(i,0,25) ans+=rec[suf^v[i]];
		ans+=rec[suf]; rec[suf]++;
		suf^=sum; rec[suf]--;
	}
}

void del(int x,int opt) {
	int w=v[a[x]];
	if (opt) {
		rep(i,0,25) ans-=rec[v[i]^suf];
		ans-=rec[suf]; suf^=w; rec[suf]--;
	} else {
		suf^=sum;
		rec[suf]--; rec[suf^sum]++;
		rep(i,0,25) ans-=rec[v[i]^suf];
		ans-=rec[suf];
		rec[suf^sum]--; suf^=sum;
	}
	sum^=w;
}

void solve(int m) {
	for (int i=1,l=1,r=0;i<=m;++i) {
		for (;r<q[i].r;) ins(++r,1);
		for (;l>q[i].l;) ins(--l,0);
		for (;r>q[i].r;) del(r--,1);
		for (;l<q[i].l;) del(l++,0);
		prt[q[i].id]=ans;
	}
}

bool cmp(Q a,Q b) {
	int al=a.l/B,bl=b.l/B;
	return (al==bl)?(a.r<b.r):(al<bl);
}

int main(void) {
	rep(i,0,25) v[i]=1<<i;
	int n=read(),m=read();
	scanf("%s",str+1);
	rep(i,1,n) a[i]=str[i]-'a';
	rep(i,1,m) {
		q[i].l=read(),q[i].r=read();
		q[i].id=i;
	}
	std:: sort(q+1,q+m+1,cmp);
	solve(m);
	rep(i,1,m) printf("%d\n", prt[i]);
	return 0;
}