分块-bzoj2821: 作诗(Poetize)

http://hzwer.com/3663.html

我的程序好垃圾，常数好像很大的样子，本来是WA，然后改好之后交上去超时；

改来改去半天还是TLE，后来把分块长度变成sqrt(n/log2(n))，就过了；

k=sqrt(n/log(n))/log(2)

然而我自己并不会求这个；

唉数学差；

黄学长讲的很详细了；

因为是强制在线，我们先尽可能多的先预处理；

就是那个f[i][j];

因为我们的预处理，所以对于一个区间，我们需要暴力求解的数字最多不过两个块，区间里成块的都算好了放在f数组里面了；

对于不在整块里的数字，我们按黄学长说的就好啦；

但是统计次数的时候，我们用前坠和就炸；

所以我们要用前向星，排序一遍，把相同的数字按原始下标有序地排到一起；这样我们就可以通过两次二分来求出某数字在某区间出现的次数了；

代码思路是很清晰的，应该也好理解，就是跑得有点慢…..

#include<cstdio>//cfb
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
struct kuai{//就是分块
int ll,rr;//这个是区间范围
}c[1000];
struct cs{//a数组有两个用途，x,num是用来排序的；ll,rr,代表排序后数字i出现的范围
int x,num,ll,rr;
}a[100001];
int f[1000][1000],b[100001],aa[100001];//b是一个计数器，aa就是存放读入的数组；
int n,m,ll,x,y,ans,z,q;
inline int read()//黄学长的
{
int x=0,f=1;char ch=getchar();
while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
bool cmp(cs x,cs y){if(x.x!=y.x)return x.x<y.x;return x.num<y.num;}
void make(){
int l=1,k=sqrt(n/log(n))/log(2);
while(l<=n){//分块；
c[++ll].ll=l;
c[ll].rr=min(l+k-1,n);
l+=k;
}
for(int i=1;i<=ll;i++){
int sum=0;
for(int j=i;j<=ll;j++){
for(int k=c[j].ll;k<=c[j].rr;k++){
b[a[k].x]++;
if(b[a[k].x]>1&&(b[a[k].x]&1))sum--;
if(!(b[a[k].x]&1))sum++;
}
f[i][j]=sum;
}
for(int j=c[i].ll;j<=c[ll].rr;j++)b[a[j].x]=0;//比memset快好多；
}
sort(a+1,a+n+1,cmp);
for(int i=1;i<=n;i++){if(!a[a[i].x].ll)a[a[i].x].ll=i;a[a[i].x].rr=i;}
}
int max_min(int x,int y,int z){//在x,y区间里找到比z大的最小的数，用于区间左端；
int ans=y;
while(y>=x){
int mid=(y+x)>>1;
if(a[mid].num>=z){
ans=min(ans,mid);
y=mid-1;
}else x=mid+1;
}
return ans;
}
int min_max(int x,int y,int z){//在x,y区间里找到比z小的最大的数，用于区间左端；
int ans=x;
while(y>=x){
int mid=(y+x)>>1;
if(a[mid].num<=z){
ans=max(ans,mid);
x=mid+1;
}else y=mid-1;
}
return ans;
}
void check(int num,int x,int y,int l,int r){
x=max_min(a[num].ll,a[num].rr,x);
y=min_max(a[num].ll,a[num].rr,y);
l=max_min(a[num].ll,a[num].rr,l);
r=min_max(a[num].ll,a[num].rr,r);
x=y-x+1; y=r-l+1;//通过二分出来的下标来统计区间内数字出现次数
if(!y&&x&&!(x&1))ans++;//我一开始漏掉了Y=0的情况
if(y&&(y&1)&&!(x&1))ans++;
if(y&&!(y&1)&&(x&1))ans--;
}
void outit(int x,int y){
int first=0,last=0;
for(int i=1;i<=ll;i++)if(x<=c[i].ll&&c[i].rr<=y){if(!first)first=i;last=i;}
ans=f[first][last];//一大块的答案树算好了的
if(last){
for(int i=x;i<=c[first].ll-1;i++)if(!b[aa[i]]){check(aa[i],x,y,c[first].ll,c[last].rr),b[aa[i]]=1;}
for(int i=y;i>=c[last ].rr+1;i--)if(!b[aa[i]]){check(aa[i],x,y,c[first].ll,c[last].rr),b[aa[i]]=1;}
for(int i=x;i<=c[first].ll-1;i++)b[aa[i]]=0;
for(int i=y;i>=c[last ].rr+1;i--)b[aa[i]]=0;
}else{//找不到完整块的处理
for(int i=x;i<=y;i++)
if(!b[aa[i]]){
b[aa[i]]=1;
int l=max_min(a[aa[i]].ll,a[aa[i]].rr,x);
int r=min_max(a[aa[i]].ll,a[aa[i]].rr,y);
if(!((r-l+1)&1))ans++;
}
for(int i=x;i<=y;i++)b[aa[i]]=0;
}
printf("%d\n",ans);
}
int main()
{
n=read();q=read();m=read();
for(int i=1;i<=n;i++)a[i].x=read(),a[i].num=i,aa[i]=a[i].x;
make();
while(m--){
x=read(),y=read();
outit(min((x+ans)%n+1,(y+ans)%n+1),max((x+ans)%n+1,(y+ans)%n+1));//强制在线
}
}