HDU 5306 Gorgeous Sequence（线段树）

Description

一个长度为$n$的序列$a_{1},...,a_{n}$，有三种操作：

$0\ x\ y\ t$:对$\forall x\leq i\leq y$，令$a_{i}=min(a_{i},t)$

$1\ x\ y$:输出$a_{i}$的最大值$(x\leq i\leq y)$

$2\ x\ y$:输出$a_{i}$的和$(x\leq i\leq y)$

2.

3.

Input

第一行输入一整数$T$表示用例组数，每组用例首先输入两个整数$n$和$m$分别表示序列长度和操作数，之后输入$n$个整数$a_{1},...,a_{n}$表示该序列，最后$m$行每行一个操作$(T=100,\sum n\leq 10^6,\sum m\leq 10^6,1\leq x\leq y\leq n,0\leq a_{i},t<2^{31})$

Output

对于2,3操作，输出结果

Sample Input

1
5 5
1 2 3 4 5
1 1 5
2 1 5
0 3 5 3
1 1 5
2 1 5

Sample Output

5
15
3
12

Solution

线段树，一个节点t维护四个值，Sum[t]表示区间和，Max[t]表示区间最值，Tag[t]表示区间关于更新操作的lazy标记，Num[t]表示区间不大于Max[t]的值的数量，每次更新操作，假设把某段区间中的数与v取最小值，如果Max[t]<=v则不用更新，否则把该区间中所有大于v的子区间的Sum,Max,Tag,Num均变成0，这样就消除了该区间中大于v（更新后变成v）的数对所维护值的影响并知道了大于v的数的个数，之后O(1)的维护该区间的四个值即可，总时间复杂度O(nlogn)

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
#define maxn 1000005
#define ls (t<<1)
#define rs ((t<<1)|1)
ll Sum[maxn<<2];
int Num[maxn<<2],Max[maxn<<2],Tag[maxn<<2];
int T,n,m;
char buf[40*1024000+5],*ch;
void read(int &x)
{
    for(++ch;*ch<=32;ch++);
    for(x=0;'0'<=*ch;ch++)x=10*x+*ch-'0';
}
void push_up(int t)
{
    Sum[t]=Sum[ls]+Sum[rs];
    Max[t]=max(Max[ls],Max[rs]);
    Num[t]=Num[ls]+Num[rs];
}
void build(int l,int r,int t)
{
    Sum[t]=Num[t]=Max[t]=Tag[t]=0;
    if(l==r)
    {
        int a;
        //scanf("%d",&a);
        read(a);
        Sum[t]=Max[t]=Tag[t]=a,Num[t]=1;
        return ;
    }
    int mid=(l+r)/2;
    build(l,mid,ls),build(mid+1,r,rs);
    push_up(t);
}
void modify(int l,int r,int t,int v)
{
    if(Tag[t]&&Tag[t]<=v)return ;
    Tag[t]=v;
    if(Num[t]<r-l+1)
    {
        Max[t]=v;
        Sum[t]+=1ll*(r-l+1-Num[t])*v;
        Num[t]=r-l+1;
    }
}
void push_down(int l,int r,int t)
{
    if(!Tag[t])return ;
    int mid=(l+r)/2;
    modify(l,mid,ls,Tag[t]),modify(mid+1,r,rs,Tag[t]);
}
void clear(int l,int r,int t,int v)
{
    if(Max[t]<=v)return ;
    Tag[t]=0;
    if(l==r)
    {
        Sum[t]=Max[t]=Num[t]=0;
        return ;
    }
    int mid=(l+r)/2;
    clear(l,mid,ls,v),clear(mid+1,r,rs,v);
    push_up(t);
}
void update(int L,int R,int l,int r,int t,int v)
{
    if(Max[t]<=v)return ;
    if(L<=l&&r<=R)
    {
        clear(l,r,t,v);
        modify(l,r,t,v);
        return ;
    }
    push_down(l,r,t);
    int mid=(l+r)/2;
    if(L<=mid)update(L,R,l,mid,ls,v);
    if(R>mid)update(L,R,mid+1,r,rs,v);
    push_up(t);
}
ll query_sum(int L,int R,int l,int r,int t)
{
    if(L<=l&&r<=R)return Sum[t];
    push_down(l,r,t);
    ll ans=0;
    int mid=(l+r)/2;
    if(L<=mid)ans+=query_sum(L,R,l,mid,ls);
    if(R>mid)ans+=query_sum(L,R,mid+1,r,rs);
    return ans;
}
int query_max(int L,int R,int l,int r,int t)
{
    if(L<=l&&r<=R)return Max[t];
    int ans=0;
    push_down(l,r,t);
    int mid=(l+r)/2;
    if(L<=mid)ans=max(ans,query_max(L,R,l,mid,ls));
    if(R>mid)ans=max(ans,query_max(L,R,mid+1,r,rs));
    return ans;
}
int main()
{
    ch=buf-1;
    fread(buf,1,1000*35*1024,stdin);
    //scanf("%d",&T);
    read(T);
    while(T--)
    {
        //scanf("%d%d",&n,&m);
        read(n),read(m);
        build(1,n,1);
        while(m--)
        {
            int op,x,y,t;
            //scanf("%d%d%d",&op,&x,&y);
            read(op),read(x),read(y);
            if(op==0)
            {
                //scanf("%d",&t);
                read(t);
                update(x,y,1,n,1,t);
            }
            else if(op==1)printf("%d\n",query_max(x,y,1,n,1));
            else printf("%I64d\n",query_sum(x,y,1,n,1));
        }
    }
    return 0;
}