其他分享
首页 > 其他分享> > 21.10.3 T4

21.10.3 T4

作者:互联网

题意

给定一个长度为\(n\)的区间\(a_i\)与\(m\)个询问,每次询问给出\(l,r\),求

\[\sum_{i=l}^r\sum_{j=l}^r(max_{k=i}^ja_k)\times(min_{k=i}^{j}a_k) \]

\(n,m\leq10^5\)

sol

考虑分治。

对于每个区间我们只维护所有跨过了\(mid\)的区间\([i,j]\)对答案的,对于其他多于此的询问的其他区间我们继续下放处理。

对于一个区间,我们枚举右端点,根据\([l,mid]\)中每个点与右端点间\(min,max\)分别是在\([l,mid]\)或是\([mid+1,r]\)上分成四种情况,即\(min/max\)在或不在\([l,mid]\)区间内。容易发现这样的划分方式将\([l,mid]\)划分成了至多四段(可能有重),则该段答案即为四种情况之和。

对于这四种情况,我们建四棵线段树,维护区间和,支持区间加,来维护其系数。(其实显然可以树状数组,考场上脑子抽了)

具体的,我们预处理出整个大区间\([1,n]\)分治的答案,为方便下传询问,将询问打包成个\(vector\),一起下放即可。

但是这码\(T\)了(((

所以还是用树状数组吧\(QAQ\)

#include <bits/stdc++.h>
using namespace std;
#define mod 1000000007
#define N 100005
#define M N<<3
#define lid (id<<1)
#define rid (id<<1|1)
#define mid ((l+r)>>1)
int s[N],maxx[N],minn[N],ans[N],f[M],n,m;
inline int read(){
    int x=0,f=1;char c=getchar();
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9')x=(x<<1)+(x<<3)+(c&15),c=getchar();
    return x*f;
}
inline void add1(int &x,int y){x+=y;if(x>=mod) x-=mod;}
inline int add2(int x,int y){x+=y;if(x>=mod) x-=mod;return x;}
struct seg_tr{
	int lz[M],sum[M],siz[M],num[N];
	void build(int id,int l,int r){
		sum[id]=lz[id]=0;
		if(l==r) siz[id]=num[l];
		else{
			build(lid,l,mid),build(rid,mid+1,r);
			siz[id]=add2(siz[lid],siz[rid]);
		}
	}
	inline void pushdown(int id){
		lz,add1(lz[rid],lz[id]);
		add1(sum[lid],1LL*lz[id]*siz[lid]%mod);
		add1(sum[rid],1LL*lz[id]*siz[rid]%mod);
		lz[id]=0;
	}
	inline void add(int id,int l,int r,int x,int y,int z){
		if(l==x&&r==y) add1(lz[id],z),add1(sum[id],1LL*z*siz[id]%mod);
		else{
			if(lz[id]) pushdown(id);
			if(y<=mid) add(lid,l,mid,x,y,z);
			else if(x>mid) add(rid,mid+1,r,x,y,z);
			else add(lid,l,mid,x,mid,z),add(rid,mid+1,r,mid+1,y,z);
			sum[id]=add2(sum[lid],sum[rid]);
		}
	}
	inline int qry(int id,int l,int r,int x,int y){
		if(l==x&&r==y) return sum[id];
		if(lz[id]) pushdown(id);
		if(y<=mid) return qry(lid,l,mid,x,y);
		if(x>mid) return qry(rid,mid+1,r,x,y);
		return add2(qry(lid,l,mid,x,mid),qry(rid,mid+1,r,mid+1,y));
	}
}tr[5];

struct query{int l,r,id;};
vector<query >qry[M];
inline bool cmpr(query x,query y){return x.r<y.r;}
inline void solve(int id,int l,int r){
	int cnt=qry[id].size();
	if(l==r){
		f[id]=1LL*s[l]*s[l]%mod;
		for(int i=0;i<cnt;i++) add1(ans[qry[id][i].id],1LL*s[l]*s[l]%mod);
		return;
	}
	sort(qry[id].begin(),qry[id].end(),cmpr);
	int lsw=0;while(lsw<cnt&&qry[id][lsw].r<=mid) lsw++;
	maxx[mid]=minn[mid]=s[mid];
	for(int i=mid-1;i>=l;--i){
		maxx[i]=max(maxx[i+1],s[i]);
		minn[i]=min(minn[i+1],s[i]);
	}
	for(int i=l;i<=mid;i++){
		tr[0].num[i]=1;
		tr[1].num[i]=maxx[i];
		tr[2].num[i]=minn[i];
		tr[3].num[i]=1LL*maxx[i]*minn[i]%mod;
	}
	int ma,mi;
	ma=mi=s[mid+1];
	int x,y,z;x=y=z=mid+1;
	for(int i=0;i<4;i++) tr[i].build(1,l,mid);
	for(int i=mid+1;i<=r;i++){
		ma=max(ma,s[i]),mi=min(mi,s[i]);
		
		while(l<x&&minn[x-1]>=mi&&maxx[x-1]<=ma)x--;
		while(l<y&&minn[y-1]>=mi)y--;
		while(l<z&&maxx[z-1]<=ma)z--;
		
		if(x<=mid)tr[0].add(1,l,mid,x,mid,1LL*ma*mi%mod);
		if(y<x){
			tr[1].add(1,l,mid,y,x-1,mi);
			if(l<y)tr[3].add(1,l,mid,l,y-1,1);
		}
		if(z<x){
			tr[2].add(1,l,mid,z,x-1,ma);
			if(l<z)tr[3].add(1,l,mid,l,z-1,1);
		}
		while(lsw<cnt&&qry[id][lsw].r==i){
			if(qry[id][lsw].l>mid||qry[id][lsw].l==l&&qry[id][lsw].r==r){lsw++;continue;}
			for(int t=0;t<4;t++) add1(ans[qry[id][lsw].id],tr[t].qry(1,l,mid,qry[id][lsw].l,mid));
			lsw++;
		}
	}
	for(int i=0;i<4;i++) add1(f[id],tr[i].qry(1,l,mid,l,mid));
	for(int i=0;i<cnt;i++){
		if(qry[id][i].l==l&&qry[id][i].r==r) continue;
		if(qry[id][i].r<=mid) qry[lid].push_back(qry[id][i]);
		else if(qry[id][i].l>mid) qry[rid].push_back(qry[id][i]);
		else{
			qry[lid].push_back((query){qry[id][i].l,mid,qry[id][i].id});
			qry[rid].push_back((query){mid+1,qry[id][i].r,qry[id][i].id});
		}
	}
	solve(lid,l,mid),solve(rid,mid+1,r);
	add1(f[id],add2(f[lid],f[rid]));
	for(int i=0;i<cnt;i++)if(qry[id][i].l==l&&qry[id][i].r==r) add1(ans[qry[id][i].id],f[id]);
}
inline void file(){
	freopen("sequence.in","r",stdin);
	freopen("sequence.out","w",stdout);
}
int main(void){
	//file();
	n=read(),m=read();
	for(int i=1;i<=n;i++) s[i]=read();
	for(int i=1;i<=m;i++){
		qry[1].push_back((query){read(),read(),i});
	}
	solve(1,1,n);
	for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
	return 0;
}

``

标签:int,T4,mid,21.10,qry,rid,lz,id
来源: https://www.cnblogs.com/zyw4347/p/15364879.html