bzoj 4540: [Hnoi2016]序列【单调栈+线段树】

强烈安利：http://blog.csdn.net/qq_34637390/article/details/51313126

这篇讲标记讲的非常好，这个标记非常神奇……

首先last表示扫描到last这个点了，val[x]表示x到last中的最小值，sum[x]表示last分别等于1,2,3....last时的和（也就是所有val历史版本的和）

然后答案就是当前last扫到r，l到r之间所有sum的和（稍微画一下就理解了

然后考虑现在要挪动last，新的last的影响范围是他之前第一个比他小的位置（单调栈维护）到新last这一段，考虑对这一段的修改，就是把所有的val修改成a[last]，sum加上原来的val

这里用abcd做lazytag，每次修改这样搞一下：（以下为截图

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
const int N=100005;
int n,m,a[N],top,s[N];
long long ans[N];
struct wen
{
	int l,r,id;
}q[N];
bool cmp(const wen &a,const wen &b)
{
	return a.r<b.r;
}
struct lazy
{
	long long a,b,c,d;
	lazy(long long A=0,long long B=0,long long C=0,long long D=0)
	{
		a=A,b=B,c=C,d=D;
	}
	lazy operator + (const lazy &y) const
	{
		return lazy(a*y.a,y.b+b*y.a,c+y.c*a,d+y.d+b*y.c);
	}
};
struct qwe
{
	int l,r;
	long long val,sum;
	lazy lz;
}t[N<<2];
int read()
{
	int r=0,f=1;
	char p=getchar();
	while(p>'9'||p<'0')
	{
		if(p=='-')
			f=-1;
		p=getchar();
	}
	while(p>='0'&&p<='9')
	{
		r=r*10+p-48;
		p=getchar();
	}
	return r*f;
}
void build(int ro,int l,int r)
{
	t[ro].l=l,t[ro].r=r,t[ro].lz=lazy(1,0,0,0);
	if(l==r)
		return;
	int mid=(l+r)>>1;
	build(ro<<1,l,mid);
	build(ro<<1|1,mid+1,r);
}
void add(int ro,lazy p)
{
	t[ro].sum+=p.c*t[ro].val+p.d*(t[ro].r-t[ro].l+1);
	t[ro].val=p.a*t[ro].val+p.b*(t[ro].r-t[ro].l+1);
	t[ro].lz=t[ro].lz+p;
}
void pd(int ro)
{
	if(t[ro].lz.a!=1||t[ro].lz.b||t[ro].lz.c||t[ro].lz.d)
	{
		add(ro<<1,t[ro].lz);
		add(ro<<1|1,t[ro].lz);
		t[ro].lz=lazy(1,0,0,0);
	}
}
void update(int ro,int l,int r,lazy p)
{
	if(t[ro].l==l&&t[ro].r==r)
	{
		add(ro,p);
		return;
	}
	pd(ro);
	int mid=(t[ro].l+t[ro].r)>>1;
	if(r<=mid)
		update(ro<<1,l,r,p);
	else if(l>mid)
		update(ro<<1|1,l,r,p);
	else
	{
		update(ro<<1,l,mid,p);
		update(ro<<1|1,mid+1,r,p);
	}
	t[ro].sum=t[ro<<1].sum+t[ro<<1|1].sum;
	t[ro].val=t[ro<<1].val+t[ro<<1|1].val;
}
long long ques(int ro,int l,int r)
{
	if(t[ro].l==l&&t[ro].r==r)
		return t[ro].sum;
	pd(ro);
	int mid=(t[ro].l+t[ro].r)>>1;
	if(r<=mid)
		return ques(ro<<1,l,r);
	else if(l>mid)
		return ques(ro<<1|1,l,r);
	else
		return ques(ro<<1,l,mid)+ques(ro<<1|1,mid+1,r);
}
int main()
{
	n=read(),m=read();
	for(int i=1;i<=n;i++)
		a[i]=read();
	for(int i=1;i<=m;i++)
		q[i].l=read(),q[i].r=read(),q[i].id=i;
	sort(q+1,q+1+m,cmp);
	build(1,1,n);
	for(int i=1,j=1;i<=n;i++)
	{
		while(top&&a[s[top]]>=a[i])
			top--;
		update(1,s[top]+1,i,lazy(0,a[i],0,0));
		add(1,lazy(1,0,1,0));
		s[++top]=i;
		for(;j<=m&&q[j].r==i;j++)
			ans[q[j].id]=ques(1,q[j].l,q[j].r);
	}
	for(int i=1;i<=m;i++)
		printf("%lld\n",ans[i]);
	return 0;
}