考虑一段区间排序后是等差为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
*/