Jzoj5231 序列问题

最新推荐文章于 2019-12-24 17:32:26 发布

扩展的灰

最新推荐文章于 2019-12-24 17:32:26 发布

阅读量489

点赞数 1

分类专栏： OI 数据结构 ----前缀和求解策略 ----分治 --------cdq分治 Jzoj 文章标签：分治 cdq分治前缀和扩展的灰

本文链接：https://blog.csdn.net/JacaJava/article/details/78545228

版权

OI 同时被 3 个专栏收录

321 篇文章 1 订阅

订阅专栏

Jzoj

227 篇文章 3 订阅

订阅专栏

求解策略

153 篇文章 0 订阅

订阅专栏

这是一道使用CDQ分治算法解决的序列问题，要求求解序列中连续子序列的最大值与最小值乘积之和。通过预处理前缀最大值和最小值，然后二分找到满足条件的端点，利用指针维护单调性，最终实现O(n log n)的时间复杂度解决方案。

摘要由CSDN通过智能技术生成

给你序列A，定义f(l,r)=max{al,al+1,....,ar} , g(l,r)=min{al,al+1,....,ar}

求ΣΣf(i,j)*g(i,j) { i∈[1,n],j∈[i,n] }

哇第一道cdq分治的题诶，当时好像还不会这个玩意2333

我们cdq分治一下，先求出[1,mid]和[mid+1,r]的答案

让后我们考虑枚举j∈[mid+1,r]

对于每一个j，我们可以通过预处理出前缀min，max再二分出两个端点p1,p2满足

f(p1,j)=g(p1,j)=A[j]

f(p2,p1)=A[j]或者g(p2,p1)=A[j]

最后f(l,p2)和g(l,p2)与A[j]无关

每次统计答案可以分成这三个区间来维护

又发现，由于A[j]的单调，p2和p1肯定是单调递减的，可以用指针来维护，成功将复杂度降为O(n lg n)

#include<stdio.h>
#include<algorithm>
#define M 1000000007
#define L long long
using namespace std;
L mx[500010],mn[500010];
L sx[500010],s[500010];
L sm[500010],sn[500010];
int n;
L cdq(int l,int r){
	if(l==r) return s[l]*s[r]%M;
	int m=l+r>>1,pi=m,pj=m;
	L ans=(cdq(l,m)+cdq(m+1,r))%M;
	mx[m+1]=-M; mn[m+1]=M; 
	sm[m+1]=sn[m+1]=sx[m+1]=0;
	for(int i=m;i>=l;--i){
		mx[i]=max(mx[i+1],s[i]);
		mn[i]=min(mn[i+1],s[i]);
		sm[i]=(sm[i+1]+mx[i])%M;
		sn[i]=(sn[i+1]+mn[i])%M;
		sx[i]=(sx[i+1]+mx[i]*mn[i]%M)%M;
	}
	L rm=-M,rn=M;
	for(int j=m+1;j<=r;++j){
		rm=max(s[j],rm);
		rn=min(s[j],rn);
		while(pi>=l&&mx[pi]<rm) pi--;
		while(pj>=l&&mn[pj]>rn) pj--;
		if(pi>pj){
			ans=(ans+(m-pi)*rm%M*rn%M)%M;
			ans=(ans+rn%M*(sm[pj+1]-sm[pi+1])%M)%M;
			ans=(ans+sx[l]-sx[pj+1]+M)%M;
		} else {
			ans=(ans+(m-pj)*rm%M*rn%M)%M;
			ans=(ans+rm*(sn[pi+1]-sn[pj+1])%M)%M;
			ans=(ans+sx[l]-sx[pi+1]+M)%M;
		}
	}
	return ans;
}
int main(){
	freopen("seq.in","r",stdin);
	freopen("seq.out","w",stdout);
	scanf("%d",&n);
	for(int i=1;i<=n;++i) scanf("%lld",s+i);
	printf("%lld\n",(M+cdq(1,n))%M);
}