BZOJ4373 算术天才⑨与等差数列

考虑一段区间排序后是等差为k的数列的条件

首先特判k=0，如果区间最大值与最小值相等则合法

否则首先区间不能有重复，记录每个点的上一个与他权值相等的点的位置fro，区间取fro最大值之后判断是否在区间左边

然后设区间长度为len，则区间最大值必须等于区间最小值+(len-1)*k

满足上边两个条件之后，只要满足区间内每个点的权值模k都相等即可

考虑求关于x的同余方程a同余b模x的解，则

a可以表示为a'x+x'，b可以表示为b'x+x'

a'x+x'=b'x+x'(mod x)

(a'-b')x=0(mod x)

(a'-b')x=a-b，所以x为方程的解的充要条件是x是a-b的约数

设原数列为a，搞一个新数列c，c[i]=|a[i+1]-a[i]|

则我们只需要判断k是否为l到r-1的c值的最大公约数的约数即可

上边的所有东西都可以用线段树维护

维护前一个相等的数可以用set

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<map>
#include<set>
#include<bitset>
#include<queue>
#include<stack>
using namespace std;
#define MAXN 300010
#define MAXM 1010
#define INF 1000000000
#define MOD 1000000007
#define eps 1e-8
#define ll long long
int gcd(int x,int y){
	return !y?x:gcd(y,x%y);
}
struct data{
	int mx;
	int mn;
	int mxl;
	int v;
	friend data operator +(data x,data y){
		data z;
		z.mxl=max(x.mxl,y.mxl);
		z.mx=max(x.mx,y.mx);
		z.mn=min(x.mn,y.mn);
		z.v=gcd(x.v,y.v);
		return z;
	}
};
int n,m;
int a[MAXN];
map<int,int>h;
int fro[MAXN];
data v[MAXN<<2];
set<int>s[MAXN<<1];
int tot;
int la;
inline void ud(int x){
	v[x]=v[x<<1]+v[x<<1|1];
}
void build(int x,int y,int z){
	if(y==z){
		v[x].mx=v[x].mn=a[y];
		v[x].v=abs(a[y+1]-a[y]);
		v[x].mxl=fro[y];
		return ;
	}
	int mid=y+z>>1;
	build(x<<1,y,mid);
	build(x<<1|1,mid+1,z);
	ud(x);
}
void change(int x,int y,int z,int p){
	if(y==z){
		v[x].mx=v[x].mn=a[y];
		v[x].v=abs(a[y+1]-a[y]);
		v[x].mxl=fro[y];
		return ;
	}
	int mid=y+z>>1;
	if(p<=mid){
		change(x<<1,y,mid,p);
	}else{
		change(x<<1|1,mid+1,z,p);
	}
	ud(x);
}
data ask(int x,int y,int z,int l,int r){
	if(y==l&&z==r){
		return v[x];
	}
	int mid=y+z>>1;
	if(r<=mid){
		return ask(x<<1,y,mid,l,r);
	}else if(l>mid){
		return ask(x<<1|1,mid+1,z,l,r);
	}else{
		return ask(x<<1,y,mid,l,mid)+ask(x<<1|1,mid+1,z,mid+1,r);
	}
}
int main(){
	int o;
	int i,x,y,z;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;i++){
		scanf("%d",&a[i]);
		if(!h[a[i]]){
			h[a[i]]=++tot;
			s[tot].insert(0);
			s[tot].insert(n+1);
		}
		int t=h[a[i]];
		s[t].insert(i);
		fro[i]=*--(s[t].insert(i).first);
	}
	a[n+1]=a[n];
	build(1,1,n);
	while(m--){
		scanf("%d",&o);
		if(o==1){
			scanf("%d%d",&x,&y);
			x^=la;
			y^=la;
			int t=h[a[x]];
			int xx=*s[t].upper_bound(x);
			if(xx!=n+1){
				fro[xx]=fro[x];
				change(1,1,n,xx);
			}
			s[t].erase(x);
			if(!h[y]){
				h[y]=++tot;
				s[tot].insert(0);
				s[tot].insert(n+1);
			}
			t=h[y];
			s[t].insert(x);
			xx=*s[t].upper_bound(x);
			fro[x]=*--s[t].lower_bound(x);
			if(xx!=n+1){
				fro[xx]=x;
				change(1,1,n,xx);
			}
			a[x]=y;
			change(1,1,n,x);
			if(x!=1){
				change(1,1,n,x-1);
			}
		}
		if(o==2){
			scanf("%d%d%d",&x,&y,&z);
			x^=la;
			y^=la;
			z^=la;
			if(x==y){
				la++;
				printf("Yes\n");
				continue ;
			}
			data t=ask(1,1,n,x,y-1);
			data tt=t+ask(1,1,n,y,y);
			if(z==0&&tt.mx==tt.mn){
				la++;
				printf("Yes\n");
				continue ;
			}
			if(tt.mxl<x&&tt.mx==tt.mn+(y-x)*z){
				if(!(t.v%z)){
					la++;
					printf("Yes\n");
				}else{
					printf("No\n");
				}
			}else{
				printf("No\n");
			}
		}
	}
	return 0;
}

/*
5 3
1 3 2 5 6
2 1 5 1
1 5 4
2 1 5 1

*/