bzoj 4866: [Ynoi2017]由乃的商场之旅

设第i个字母的权值为1<<i,则一个可重集合可以重排为回文串,当且仅当这个集合的异或和x满足x==x&-x,用莫队维护区间内有多少对异或前缀和,异或后满足x==x&-x,这样端点移动的代价为字符集大小+1=27,因此时间复杂度为$O(27n\sqrt{m})$

#include<cstdio>
#include<cmath>
#include<algorithm>
char buf[],*ptr=buf-;
int _(){
int x=,c=*++ptr;
while(c<)c=*++ptr;
while(c>)x=x*+c-,c=*++ptr;
return x;
}
typedef unsigned int u32;
const int P=,N=;
int n,q,xs[N][];
int hx[P][],idp=;
int getid(int x){
int w=x%P;
while(hx[w][]){
if(hx[w][]==x)return hx[w][];
if((w+=)>=P)w-=P;
}
hx[w][]=x;
return hx[w][]=++idp;
}
u32 as[N],pos[N],B,ans=,t[N*];
struct Q{
int l,r,id;
}qs[N];
bool operator<(Q a,Q b){
if(pos[a.l]!=pos[b.l])return pos[a.l]<pos[b.l];
if(a.r!=b.r)return (a.r<b.r)^(pos[a.l]&);
return a.id<b.id;
}
void ins(int*x){
for(int i=;i<=;++i)ans+=t[x[i]];
++t[x[]];
}
void del(int*x){
--t[x[]];
for(int i=;i<=;++i)ans-=t[x[i]];
}
int main(){
fread(buf,,sizeof(buf),stdin)[buf]=;
n=_();q=_();
B=(n+)/sqrt(q+)+;
for(int i=;i<=n;++i)pos[i]=i/B;
while(*ptr<'a')++ptr;
for(int i=;i<=n;++i)xs[i][]=xs[i-][]^<<*ptr++-'a';
for(int i=;i<=n;++i){
for(int j=;j<;++j)xs[i][j]=xs[i][]^<<j;
}
for(int i=;i<=n;++i){
for(int j=;j<=;++j)xs[i][j]=getid(xs[i][j]);
}
for(int i=;i<q;++i){
qs[i].l=_()-;
qs[i].r=_();
qs[i].id=i;
}
std::sort(qs,qs+q);
int L=,R=;
for(int i=;i<q;++i){
int l=qs[i].l,r=qs[i].r;
while(L>l)ins(xs[--L]);
while(R<r)ins(xs[++R]);
while(L<l)del(xs[L++]);
while(R>r)del(xs[R--]);
as[qs[i].id]=ans;
}
for(int i=;i<q;++i)printf("%u\n",as[i]);
return ;
}
上一篇:Pig jline.Terminal错误


下一篇:pom.xml中添加远程仓库