Codeforces Round #250 (Div. 1) D. The Child and Sequence 分块

版权声明:存在错误或者不清楚的地方还望指出 https://blog.csdn.net/hypHuangYanPing/article/details/83272421
/*
Codeforces Round #250 (Div. 1) D. The Child and Sequence  分块
链接:https://codeforces.com/contest/438/problem/D
题意:
op 1 : 区间加
op 2 : 区间每个数取模
op 3 : 单点修改;

分块
操作1 : 区间加 完整块O(1),非完整块暴力累加即可;
操作2 :  取模 对于每个数 每次取模只可能减小 且每次取余减小的程度是log级别的,想法类似于开方;
         对于完整块,维护每块的最大值,看所给模数与块内最大值比是否需要进行更新操作;
         更新操作:维护块最大值,和,单点元素的值;
	 对于非完整块 直接取余,再更新当前块的和 和 最大值;
操作3 :  单点修改 , 更新最大值 和 和 即可;

**********tricks*********
(1%3+2%3)==?(1+2)%3;  艹........
输入挂加上 貌似第12组数据输入量较大; 莽夫......
times:1996ms
*/

#include<bits/stdc++.h>
#define ll long long
using namespace std;

/**********************************************Head-----Template****************************************/
bool Finish_read;
template<class T>inline void read(T &x){Finish_read=0;x=0;int f=1;char ch=getchar();while(!isdigit(ch)){if(ch=='-')f=-1;if(ch==EOF)return;ch=getchar();}while(isdigit(ch))x=x*10+ch-'0',ch=getchar();x*=f;Finish_read=1;}
template<class T>inline void print(T x){if(x/10!=0)print(x/10);putchar(x%10+'0');}
template<class T>inline void writeln(T x){if(x<0)putchar('-');x=abs(x);print(x);putchar('\n');}
template<class T>inline void write(T x){if(x<0)putchar('-');x=abs(x);print(x);}
ll gcd(ll a,ll b){return b==0?a:gcd(b,a%b);}
ll lcm(ll a,ll b){ll gg=gcd(a,b);a/=gg;if(a<=LLONG_MAX/b) return a*b;return LLONG_MAX;}
/********************************Head----Temlate**********************************************/

const int maxn=1e5+7;
int pos[maxn],blo,n,m;
ll sum[maxn];
int mx[maxn],s[maxn];

void push_down(int x){
	mx[x]=0;sum[x]=0;
	for(int i=(x-1)*blo+1;i<=min(n,x*blo);i++) mx[x]=max(mx[x],s[i]),sum[x]+=s[i];
}

void Update(int x,int mod){
	mx[x]=0;sum[x]=0;
	for(int i=(x-1)*blo+1;i<=min(blo*x,n);i++) s[i]%=mod,mx[x]=max(mx[x],s[i]),sum[x]+=s[i];
}

void update(int l,int r,int mod){
	for(int i=pos[l]+1;i<=pos[r]-1;i++) if(mod<=mx[i]) Update(i,mod);

	if(pos[l]!=pos[r]){
		for(int i=l;i<=pos[l]*blo;i++) s[i]%=mod;
		for(int i=(pos[r]-1)*blo+1;i<=r;i++) s[i]%=mod;
		push_down(pos[l]),push_down(pos[r]);
	}
	else {
		for(int i=l;i<=r;i++) s[i]%=mod;
		push_down(pos[l]);
	}
}

ll query(int l,int r){
	ll ans=0;
	for(int i=pos[l]+1;i<=pos[r]-1;i++) ans+=sum[i];
	if(pos[l]!=pos[r]){
		for(int i=l;i<=pos[l]*blo;i++) ans+=s[i];
		for(int i=(pos[r]-1)*blo+1;i<=r;i++) ans+=s[i];
	}
	else for(int i=l;i<=r;i++) ans+=s[i];
	return ans;
}

void change(int id,int x){
	s[id]=x;
	push_down(pos[id]);
}

int main (){
	scanf("%d %d",&n,&m);blo=(int)sqrt(n*1.0);
	for(int i=1;i<=n;i++) read(s[i]);
	for(int i=1;i<=n;i++) pos[i]=(i-1)/blo+1,mx[i]=0,sum[i]=0;
	for(int i=1;i<=pos[n];i++)
		for(int j=(i-1)*blo+1;j<=min(n,i*blo);j++)
			mx[i]=max(mx[i],s[j]),sum[i]+=s[j]*1ll;
	while(m--){
		int op;read(op);
		if(op==1) {
			int l,r;read(l),read(r);
			writeln(query(l,r));
		}
		else if(op==2){
			int l,r,x; read(l),read(r),read(x);;
			update(l,r,x);
		}
		else if(op==3){
			int k,x; read(k),read(x);
			change(k,x);
		}
	}
    return 0;
}

猜你喜欢

转载自blog.csdn.net/hypHuangYanPing/article/details/83272421