luogu P7446 [Ynoi2007] rfplca

题面传送门
震惊!Ynoi居然不卡常。
首先你看到这个东西开了\(2.5s\)而且数据范围只有\(4\times 10^5\)就大概知道这个是个根号了。
然后又发现这个有很强的弹飞绵羊既视感。
于是考虑分块。
考虑维护和弹飞绵羊一样的一个东西\(B_i\)表示\(i\)跳出当前块到哪里。
然后零散块暴力重构即可。
整块不太好做。
可以发现这个东西只增不减就想到五彩斑斓的世界那道题。
然后发现在进行块长次操作后一定会很好维护。
所以记录一下能不能修改即可。
时间复杂度\(O((n+q)\sqrt n)\)
code:

#include <vector>
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<cmath>
#include<algorithm>
#include<bitset>
#include<set>
#include<map>
#define I inline
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define abs(x) ((x)>0?(x):-(x))
#define l(x) x<<1
#define r(x) x<<1|1
#define re register
#define ll long long
#define db double
#define N 500000
#define M 2000
#define eps (1e-5)
#define mod (1<<31)
#define U unsigned int
using namespace std;
int n,m,A[N+5],B[N+5],L[N+5],R[N+5],k,F[N+5],G[N+5],x,y,z,op,nx,ny,Lastans;
I int GetA(int x){return max(A[x]-G[x/k],0);}
I int GetB(int x){return max(B[x]-G[x/k],0);}
I void Make(int x){for(re int i=L[x];i<=R[x];i++) B[i]=(A[i]<L[x]?A[i]:B[A[i]]);}
I void swap(int &x,int &y){x^=y^=x^=y;}
int main(){
	freopen("1.in","r",stdin);
	re int i,j;scanf("%d%d",&n,&m);n--;k=sqrt(n);for(i=1;i<=n;i++) scanf("%d",&A[i]),A[i]--;
	for(i=0;i<=n/k;i++)L[i]=max(i*k,1),R[i]=min(i*k+k-1,n);
	for(i=0;i<=n/k;i++){
		Make(i);for(j=L[i];j<=R[i];j++)A[j]>=L[i]&&(F[i]=1);
	}
	while(m--){
		scanf("%d%d%d",&op,&x,&y);x^=Lastans;y^=Lastans;x>y&&(swap(x,y),0);x--;y--;if(op==1){
			scanf("%d",&z);z^=Lastans;nx=x/k;ny=y/k;if(nx==ny){for(i=x;i<=y;i++) A[i]=max(A[i]-z,0);Make(nx);continue;}
			for(i=x;i<=R[nx];i++) A[i]=max(A[i]-z,0);Make(nx);
			for(i=L[ny];i<=y;i++) A[i]=max(A[i]-z,0);Make(ny);
			for(i=nx+1;i<ny;i++){
				if(!F[i]) {G[i]+=z;G[i]=min(G[i],2*n);continue;}F[i]=0;
				for(j=L[i];j<=R[i];j++) A[j]=max(A[j]-z,0),A[j]>=L[i]&&(F[i]=1);Make(i);
			} 
		}
		else{
			while(x^y)x>y&&(swap(x,y),0),nx=x/k,ny=y/k,nx^ny?(y=GetB(y)):(GetB(x)^GetB(y)?(x=GetB(x),y=GetB(y)):(y=GetA(y)));
			printf("%d\n",Lastans=x+1);
		}
	}
}
上一篇:图论问题复习总结


下一篇:[一本通学习笔记] 广度优先搜索和优化