「雅礼集训 2017 Day1」市场

「雅礼集训 2017 Day1」市场

传送门

Loj

题解

如果没有(2)操作,那就是一个线段树区间修改模板.

你考虑除法每一次修改至少(/2),那么最多修改( ext{log})次就会变成(1).

所以只要暴力把(2)操作( ext{pushdown})即可.

复杂度(O(nlog^2 n))

代码

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<iostream>
using namespace std;
#define ll long long
#define REP(a,b,c) for(int a=b;a<=c;a++)
#define re register
#define file(a) freopen(a".in","r",stdin);freopen(a".out","w",stdout)
typedef pair<int,int> pii;
#define mp make_pair
inline int gi()
{
	int f=1,sum=0;char ch=getchar();
	while(ch>'9' || ch<'0'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0' && ch<='9'){sum=(sum<<3)+(sum<<1)+ch-'0';ch=getchar();}
	return f*sum;
}
const int N=200010;
int n,a[N];
int tag[N<<2],mn[N<<2],mx[N<<2];
ll s[N<<2];
void pushup(int o){s[o]=s[o<<1]+s[o<<1|1];mn[o]=min(mn[o<<1],mn[o<<1|1]);mx[o]=max(mx[o<<1],mx[o<<1|1]);}
void puttag(int o,int v,int len)
{
	s[o]+=1ll*v*len;
	mn[o]+=v;mx[o]+=v;
	tag[o]+=v;
}
void pushdown(int o,int l,int r)
{
	if(!tag[o])return;int mid=(l+r)>>1;
	puttag(o<<1,tag[o],mid-l+1);puttag(o<<1|1,tag[o],r-mid);
	tag[o]=0;
}
void build(int o,int l,int r)
{
	tag[o]=0;
	if(l==r){s[o]=mn[o]=mx[o]=a[l];return;}
	int mid=(l+r)>>1;build(o<<1,l,mid);build(o<<1|1,mid+1,r);
	pushup(o);
}
void modify_add(int o,int l,int r,int posl,int posr,int v)
{
	if(posl<=l&&r<=posr){puttag(o,v,r-l+1);return;}
	int mid=(l+r)>>1;pushdown(o,l,r);
	if(posl<=mid)modify_add(o<<1,l,mid,posl,posr,v);
	if(mid<posr)modify_add(o<<1|1,mid+1,r,posl,posr,v);
	pushup(o);
}
int get(int x,int d){return x>=0?x/d:(x-d+1)/d;}
void modify_div(int o,int l,int r,int posl,int posr,int v)
{
	if(posl<=l&&r<=posr)
		if(mn[o]-get(mn[o],v)==mx[o]-get(mx[o],v)){puttag(o,get(mn[o],v)-mn[o],r-l+1);return;}
	int mid=(l+r)>>1;pushdown(o,l,r);
	if(posl<=mid)modify_div(o<<1,l,mid,posl,posr,v);
	if(mid<posr)modify_div(o<<1|1,mid+1,r,posl,posr,v);
	pushup(o);
}
int query_mn(int o,int l,int r,int posl,int posr)
{
	if(posl<=l&&r<=posr)return mn[o];
	int mid=(l+r)>>1,mn=2e9;pushdown(o,l,r);
	if(posl<=mid)mn=min(mn,query_mn(o<<1,l,mid,posl,posr));
	if(mid<posr)mn=min(mn,query_mn(o<<1|1,mid+1,r,posl,posr));
	return mn;
}
ll query(int o,int l,int r,int posl,int posr)
{
	if(posl<=l&&r<=posr)return s[o];
	int mid=(l+r)>>1;ll sum=0;pushdown(o,l,r);
	if(posl<=mid)sum+=query(o<<1,l,mid,posl,posr);
	if(mid<posr)sum+=query(o<<1|1,mid+1,r,posl,posr);
	return sum;
}
int main()
{
	n=gi();int Q=gi();
	for(int i=1;i<=n;i++)a[i]=gi();
	build(1,1,n);
	while(Q--)
	{
		int opt=gi(),l=gi()+1,r=gi()+1;
		if(opt==1)modify_add(1,1,n,l,r,gi());
		if(opt==2)modify_div(1,1,n,l,r,gi());
		if(opt==3)printf("%d
",query_mn(1,1,n,l,r));
		if(opt==4)printf("%lld
",query(1,1,n,l,r));
	}
	return 0;
}