• 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
    ",ans[i]);
    	return 0;
    }
    
  • 相关阅读:
    PHP(Zend Studio)入门视频
    PHP开发工具
    ggplot2 pdf import in Adobe Illustrator missing font AdobePiStd
    Python逐行读取文件内容
    R语言合并data.frame
    shell删除最后一列、删除第一行、比较文件
    C语言的fopen函数(文件操作/读写)
    Android 必须知道2018年流行的框架库及开发语言,看这一篇就够了!
    谈谈你对Glide和Picasso他们的对比的优缺点
    Google推荐——Glide使用详解(图片加载框架)
  • 原文地址:https://www.cnblogs.com/lokiii/p/8626718.html
Copyright © 2020-2023  润新知