[分块] 分块练习三题

分块嘛…
冯巨道德午餐讲过的

唯一不同的是可能需要st[N],ed[N]来记录一下范围

例 教主的魔法

板

#include<bits/stdc++.h>
using namespace std;
#define in Read()
int in{
	int i=0,f=1;char ch=0;
	while(!isdigit(ch)&&ch!='-') ch=getchar();
	if(ch=='-') ch=getchar(),f=-1;
	while(isdigit(ch)) i=(i<<1)+(i<<3)+ch-48,ch=getchar();
	return i*f;
}

const int N=1e6+5;
int n,q,a[N],b[N],add[N];
int siz,m,st[N],ed[N];

void init(){
	n=in,q=in;
	for(int i=1;i<=n;++i) a[i]=b[i]=in;
	siz=sqrt(n);
	m=n/siz+1;
	for(int i=1;i<=m;++i){
		st[i]=siz*(i-1)+1;
		ed[i]=siz*i;
	}
	ed[m]=n;
	for(int i=1;i<=m;++i)
		sort(b+st[i],b+ed[i]+1);
	return;
}

void frag_modify(int l,int r,int num,int c){
	for(int i=l;i<=r;++i) a[i]+=c;
	for(int i=st[num];i<=ed[num];++i) b[i]=a[i];
	sort(b+st[num],b+ed[num]+1);
	return;
}

void blok_modify(int blol,int blor,int c){
	for(int i=blol;i<=blor;++i) add[i]+=c;
	return;
}

void modify(int l,int r,int c){
	int blol=l/siz+1,blor=r/siz+1;
	if(blol==blor) frag_modify(l,r,blol,c);
	else if(blor-blol==1) frag_modify(l,ed[blol],blol,c),frag_modify(st[blor],r,blor,c);
	else frag_modify(l,ed[blol],blol,c),frag_modify(st[blor],r,blor,c),blok_modify(blol+1,blor-1,c);
	return;
}

int frag_inquir(int l,int r,int num,int c){
	int cnt=0;
	for(int i=l;i<=r;++i) if(a[i]+add[num]>=c) ++cnt;
//	printf("[%d,%d] %d\n",l,r,num);
	return cnt;
}

int blok_inquir(int blol,int blor,int c){
	int cnt=0;
	for(int i=blol;i<=blor;++i){
		int pos=lower_bound(b+st[i],b+ed[i]+1,c-add[i])-b;
		cnt+=ed[i]-pos+1;
	}
	return cnt;
}

void inquir(int l,int r,int c){
	int res=0;
	int blol=l/siz+1,blor=r/siz+1;
	if(blol==blor) res+=frag_inquir(l,r,blol,c);
	else if(blor-blol==1) res+=(frag_inquir(l,ed[blol],blol,c)+frag_inquir(st[blor],r,blor,c));
	else res+=frag_inquir(l,ed[blol],blol,c)+frag_inquir(st[blor],r,blor,c)+blok_inquir(blol+1,blor-1,c);
//	printf("%d %d %d\n",frag_inquir(l,ed[blol],blol,c),frag_inquir(st[blor],r,blor,c),blok_inquir(blol+1,blor-1,c));
	printf("%d\n",res);
	return;
}

int main(){
	init();
	
//	for(int i=1;i<=m;++i) printf("[%d,%d] ",st[i],ed[i]);
	
	for(int i=1;i<=q;++i){
		char opt[5];int a,b,c;
		scanf("%s",opt);a=in,b=in,c=in;
		if(opt[0]=='M') modify(a,b,c);
		if(opt[0]=='A') inquir(a,b,c);
	}
	return 0;
}

例 Argestes and Sequence

板

#include<bits/stdc++.h>
using namespace std;
#define in Read()
#define int long long
int in{
	int i=0,f=1;char ch=0;
	while(!isdigit(ch)&&ch!='-') ch=getchar();
	if(ch=='-') ch=getchar(),f=-1;
	while(isdigit(ch)) i=(i<<1)+(i<<3)+ch-48,ch=getchar();
	return i*f;
}

const int N=1e5+5;
int b[1000][13][10],a[N];
int n,m,T,siz,tot,st[N],ed[N];

int qpw(int a,int b){
	int res=1;
	while(b){
		if(b&1) res*=a;
		a*=a;
		b>>=1;
	}
	return res;
}

void init(){
	
	memset(b,0,sizeof b);
	
	n=in,m=in;
	siz=sqrt(n);
	tot=n/siz+1;
	for(int i=1;i<=tot;++i){
		st[i]=siz*(i-1)+1;
		ed[i]=siz*i;
	}
	ed[tot]=n;
	
	int block=1;
	for(int i=1;i<=n;++i){
		int x=a[i]=in,pos=0;
		if(i>ed[block]) ++block;
		while(x){
			++b[block][++pos][x%10];
			x/=10;
		}
		for(int i=pos+1;i<=10;++i) ++b[block][i][0];
	}
	
	return;
}

void modify(int x,int y){
	int t=a[x],block=x/siz+(bool)(x%siz),pos=0;
	while(t){
		--b[block][++pos][t%10];
		t/=10;
	}
	for(int i=pos+1;i<=10;++i) --b[block][i][0];
	a[x]=y,pos=0;
	while(y){
		++b[block][++pos][y%10];
		y/=10;
	}
	for(int i=pos+1;i<=10;++i) ++b[block][i][0];
	return;
}

int frag(int l,int r,int d,int p){
	int cnt=0;
	for(int i=l;i<=r;++i){
		if((a[i]/qpw(10,d-1))%10==p) ++cnt;
	}
	return cnt;
}

int blok(int l,int r,int d,int p){
	int cnt=0;
	for(int i=l;i<=r;++i) cnt+=b[i][d][p];
	return cnt;
}

void query(int l,int r,int d,int p){
	int blol=l/siz+(bool)(l%siz),blor=r/siz+(bool)(r%siz),res;
	if(blor-blol<=1) res=frag(l,r,d,p);
	else res=frag(l,ed[blol],d,p)+frag(st[blor],r,d,p)+blok(blol+1,blor-1,d,p);
	printf("%d\n",res);
	return;
}

signed main(){
	T=in;
	while(T--){
		init();
		for(int i=1;i<=m;++i){
			char opt[3];
			scanf("%s",opt);
			if(opt[0]=='S'){
				int a=in,b=in;
				modify(a,b);
			}else if(opt[0]=='Q'){
				int a=in,b=in,c=in,d=in;
				query(a,b,c,d);
			}
		}
	}
	return 0;
} 

例 [HNOI2010]弹飞绵羊

想学LCT

考虑分块
单点查，单点改

考虑维护一个块：
记录块中点跳出块所需的次数step和跳出去的点to
然后就可以跳块了
预处理时倒着加

分块的核心也许就在于把单点暴力合并成块吧

吐槽
这题下标太坑了
而且下标错了就T掉

#include<bits/stdc++.h>
using namespace std;
#define in Read()
#define int long long
int in{
	int i=0,f=1;char ch=0;
	while(!isdigit(ch)&&ch!='-') ch=getchar();
	if(ch=='-') ch=getchar(),f=-1;
	while(isdigit(ch)) i=(i<<1)+(i<<3)+ch-48,ch=getchar();
	return i*f;
}

const int N=2e5+5;
int n,m,bl[N],siz,tot,k[N],st[N],to[N],step[N];

void init(){
	
	n=in;
	siz=sqrt(n);
	for(int i=1;i<=n;++i){
		k[i]=in;
		bl[i]=(i-1)/siz+1;
		if(bl[i]^bl[i-1]) st[bl[i]]=i;
	}
	st[bl[n]+1]=n+1;
	m=in;
	
	for(int i=n;i>=1;--i){
		if(i+k[i]>=st[bl[i]+1]){
			step[i]=1;
			to[i]=i+k[i];
		}else{
			step[i]=step[i+k[i]]+1;
			to[i]=to[i+k[i]];
		}
	}
	
	return;
	
}

void query(int p){
	int res=0;
	while(p<=n){
		res+=step[p];
		p=to[p];
	}
	printf("%lld\n",res);
	return;
}

void modify(int p,int t){
	k[p]=t;
	for(int i=st[bl[p]+1]-1;i>=st[bl[p]];--i){
		if(i+k[i]>=st[bl[p]+1]){
			step[i]=1;
			to[i]=i+k[i];
		}else{
			step[i]=step[i+k[i]]+1;
			to[i]=to[i+k[i]];
		}
	}
	return;
}

signed main(){
	init();
	for(int i=1;i<=m;++i){
		int opt=in;
		if(opt==1){
			query(in+1);
		}else{
			int a=in+1,b=in;
			modify(a,b);
		}
	}
	return 0;
}