【JZOJ 5295】【清华集训2017模拟】Create

Description

Solution

题目有一个比较显然的性质：因为我们不停的把连续的一段变成一个数的话，那么全局不同的数会越来越少；换句话说，就是如果我们在区间修改时，对于一个数都相同的区间我们才能修改，这样均摊下来复杂度并不是很高，

所以，现在的问题就变成了：如何快速求一段数相同的区间内对答案的贡献，

设这个区间的值为c，那么只有v小于等于c的三元组才有可能有贡献，
而贡献为多少，则是这个区间被这些可能的三元组覆盖的总面积（这么写大家看得懂吧。。。），

这个是标准的主席树嘛，

这样的复杂度为$O(n\log(n)^2)$，可能过不了，
我们发现，在主席树上遍历的过程和在线段树上的很像，考虑两个同时进行，这样既可。

Code

#include <cstdio>
#include <cstdlib>
#include <algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
#define min(q,w) ((q)>(w)?(w):(q))
#define max(q,w) ((q)<(w)?(w):(q))
#define CNT(L,R,E) (finds(1,n,root[fef[E]],L,R))
#define NW(e1,e) if(!(e)||(e1)==(e))b[(e)=++b0]=b[e1];
using namespace std;
typedef long long LL;
const int N=100500;
int read(int &n)
{
    char ch=' ';int q=0,w=1;
    for(;(ch!='-')&&((ch<'0')||(ch>'9'));ch=getchar());
    if(ch=='-')w=-1,ch=getchar();
    for(;ch>='0' && ch<='9';ch=getchar())q=q*10+ch-48;n=q*w;return n;
}
int m,n;
LL ans;
int a[N];
struct qqww
{
    int v,la,s;
}b1[N*4];
int root[N],b0;
struct qwqw
{
    int l,r,v;
    int la;
}b[N*70];
int sc[N][3],p[N],pn;
int fef[N];
bool PX(int q,int w){return sc[q][2]<sc[w][2];}
void build(int l,int r,int e1,int e,int l1,int r1)
{
    if(l1<=l&&r<=r1)
    {
        b[e].la++;
        b[e].v+=(r-l+1);
        return;
    }
    b[e].v+=(r1-l1+1);
    int t=(l+r)>>1;
    if(r1<=t)
    {
        NW(b[e1].l,b[e].l);
        build(l,t,b[e1].l,b[e].l,l1,r1);
    }
    else if(t<l1)
    {
        NW(b[e1].r,b[e].r);
        build(t+1,r,b[e1].r,b[e].r,l1,r1);
    }
    else 
    {
        NW(b[e1].l,b[e].l);NW(b[e1].r,b[e].r);
        build(l,t,b[e1].l,b[e].l,l1,t),build(t+1,r,b[e1].r,b[e].r,t+1,r1);
    }
}
int finds(int l,int r,int e,int l1,int r1)
{
    if(!e)return 0;
    if(l1==l&&r1==r)return b[e].v;
    int t=(l+r)>>1;
    int ans=(r1-l1+1)*b[e].la;
    if(r1<=t)ans+=finds(l,t,b[e].l,l1,r1);
    else if(t<l1)ans+=finds(t+1,r,b[e].r,l1,r1);
    else ans+=finds(l,t,b[e].l,l1,t)+finds(t+1,r,b[e].r,t+1,r1);
    return ans;
}
void merge(int l,int r,int e)
{
    if(b1[e*2].v==b1[e*2+1].v&&b1[e*2].v!=-1)
    {
        b1[e].v=b1[e*2].v;
        int t=(l+r)>>1;
        if(b1[e*2].s==-1)b1[e*2].s=CNT(l,t,b1[e].v);
        if(b1[e*2+1].s==-1)b1[e*2+1].s=CNT(t+1,r,b1[e].v);
        b1[e].s=b1[e*2].s+b1[e*2+1].s;
    }else b1[e].v=-1,b1[e].s=-1;
}
void build(int l,int r,int e)
{
    if(l==r)
    {
        b1[e].v=a[l];
        ans+=(b1[e].s=CNT(l,r,a[l]));
        return;
    }
    int t=(l+r)>>1;
    build(l,t,e*2);
    build(t+1,r,e*2+1);
    merge(l,r,e);
}
inline void doit(int l,int r,int e)
{
    if(!b1[e].la)return;
    b1[e].v=b1[e].la;
    b1[e].s=-1;
    if(l!=r)b1[e*2].la=b1[e].la,b1[e*2+1].la=b1[e].la;
    b1[e].la=0;
}
void change(int l,int r,int e,int l1,int r1,int l2)
{
    int t=(l+r)>>1;
    if(l!=r)doit(l,t,e*2),doit(t+1,r,e*2+1);
    if(l1<=l&&r<=r1)
    {
        if(b1[e].v!=-1)
        {
            if(b1[e].s==-1)b1[e].s=CNT(l,r,b1[e].v);
            ans-=b1[e].s;
            b1[e].la=l2;
            doit(l,r,e);
            b1[e].s=CNT(l,r,b1[e].v);
            ans+=b1[e].s;
            return;
        }
    }
    if(r1<=t)change(l,t,e*2,l1,r1,l2);
    else if(t<l1)change(t+1,r,e*2+1,l1,r1,l2);
    else change(l,t,e*2,l1,t,l2),change(t+1,r,e*2+1,t+1,r1,l2);
    merge(l,r,e);
}
int main()
{
    freopen("create.in","r",stdin);
    freopen("create.out","w",stdout);
    int q,w,e,_;
    read(n),read(m),read(_);
    fo(i,1,n)read(a[i]);
    fo(i,1,m)
    {
        p[i]=i;
        fo(j,0,2)read(sc[i][j]);
    }
    sort(p+1,p+1+m,PX);
    sc[0][2]=-2333;
    w=0;
    fo(i,1,m)
    {
        w+=(sc[q][2]!=sc[p[i]][2]);
        q=p[i];
        NW(root[w-1],root[w])
        build(1,n,root[w-1],root[w],sc[q][0],sc[q][1]);
        p[w]=sc[q][2];
    }
    pn=w;
    fo(i,1,pn)fef[p[i]]=i;
    fo(i,1,n+1)if(!fef[i])fef[i]=fef[i-1];
    ans=0;
    build(1,n,1);
    printf("%d\n",ans);
    fo(i,1,_)
    {
        read(q),read(w),read(e);
        q^=ans;w^=ans;e^=ans;
        change(1,n,1,q,w,e);
        printf("%lld\n",ans);
    }
    fclose(stdout);
    fclose(stdin);
    return 0;
}