uoj #46[清华集训2014]玄学

uoj

因为询问是关于一段连续区间内的操作的,所以对操作构建线段树,这里每个点维护若干个不交的区间,每个区间\((l,r,a,b)\)表示区间\([l,r]\)内的数要变成\(ax+b\)

每次把新操作加入线段树中下一个叶子,然后如果某个节点里所有操作都加进去了,就条到父亲,把两个儿子的信息合并到父亲上.这里合并就是把两个区间集合合并成一个,例如两个区间\([a,c]\)和\([b,d](a\le b\le c\le d)\)会合并成\([a,b),[b,c),[c,d]\).合并出来的区间如果属于两个集合中元素的交,那么对应的\(ax+b\)就是\(a_r(a_l+b_l)+b_r=a_ra_lx+a_rb_l+b_r\).

查询时依次找到对应的线段树节点,然后在区间集合上二分出\(k\)元素所在的区间,把\(ans\)更新成\(a*ans+b\)即可

#include<bits/stdc++.h>
#define LL long long
#define uLL unsigned long long
#define db double

using namespace std;
const int N=6e5+10;
LL rd()
{
    LL x=0,w=1;char ch=0;
    while(ch<'0'||ch>'9'){if(ch=='-') w=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
    return x*w;
}
int n,mod,m,q,tt,ans;
struct node
{
    int l,r,a,b;
    node(){}
    node(int nl,int nr,int na,int nb){l=nl,r=nr,a=na,b=nb;}
    node operator + (const node &bb) const {return node(max(l,bb.l),min(r,bb.r),1ll*a*bb.a%mod,(1ll*b*bb.a+bb.b)%mod);}
}a[N*40];
#define mid ((l+r)>>1)
int va[N],aq[N<<2][2];
void modif(int o,int l,int r,node x)
{
    if(l==r){a[++tt]=x,aq[o][0]=aq[o][1]=tt;return;}
    if(m<=mid) modif(o<<1,l,mid,x);
    else modif(o<<1|1,mid+1,r,x);
    if(m<r) return;
    int l1=aq[o<<1][0],r1=aq[o<<1][1],l2=aq[o<<1|1][0],r2=aq[o<<1|1][1];
    aq[o][0]=tt+1;
    while(l1<=r1||l2<=r2)
    {
        if(l2>r2||(l1<=r1&&a[l1].r<a[l2].l)) a[++tt]=a[l1];
        else if(l1>r1||(l2<=r2&&a[l2].r<a[l1].l)) a[++tt]=a[l2];
        else
        {
            if(a[l1].l<a[l2].l) a[++tt]=node(a[l1].l,a[l2].l-1,a[l1].a,a[l1].b);
            else if(a[l1].l>a[l2].l) a[++tt]=node(a[l2].l,a[l1].l-1,a[l2].a,a[l2].b);
            if(tt>aq[o][0]) a[tt].l=max(a[tt].l,a[tt-1].r+1);
            if(a[tt].l>a[tt].r) --tt;
            a[++tt]=a[l1]+a[l2];
        }
        if(tt!=aq[o][0]) a[tt].l=max(a[tt].l,a[tt-1].r+1);
        if(a[tt].l>a[tt].r) --tt;
        if(l1<=r1&&l2<=r2&&a[l1].r==a[l2].r) ++l1,++l2;
        else if(l2>r2||(l1<=r1&&a[l1].r<a[l2].r)) ++l1;
        else ++l2;
    }
    aq[o][1]=tt;
}
void wk(int o,int x)
{
    int ll=aq[o][0],rr=aq[o][1],z=0;
    while(ll<=rr)
    {
        int md=(ll+rr)>>1;
        if(a[md].l<=x) z=md,ll=md+1;
        else rr=md-1;
    }
    if(a[z].r>=x) ans=(1ll*ans*a[z].a+a[z].b)%mod;
}
void quer(int o,int l,int r,int ll,int rr,int x)
{
    if(ll<=l&&r<=rr){wk(o,x);return;}
    if(ll<=mid) quer(o<<1,l,mid,ll,rr,x);
    if(rr>mid) quer(o<<1|1,mid+1,r,ll,rr,x);
}

int main()
{
    int ty=rd()&1;
    n=rd(),mod=rd();
    for(int i=1;i<=n;++i) va[i]=rd();
    q=rd();
    for(int i=1;i<=q;++i)
    {
        int op=rd(),l=rd()^(ans*ty),r=rd()^(ans*ty);
        if(op==1)
        {
            int aa=rd(),bb=rd();
            ++m;
            modif(1,1,q,node(l,r,aa,bb));
        }
        else
        {
            int x=rd()^(ans*ty);
            ans=va[x];
            quer(1,1,q,l,r,x);
            printf("%d\n",ans);
        }
    }
    return 0;
}

原文地址：https://www.cnblogs.com/smyjr/p/11508151.html