BZOJ4592: [Shoi2015]脑洞治疗仪

BZOJ4592

（吐槽一下$dsy$，我=这个数据=了一天我$xxx！！$）

一开始看到这个题。诶，区间连续最大和诶！把脑洞的值为1，正常的赋为$-INF$
但是操作1怎么办。。
发现就是先统计出$l0到r0$区间内$-INF$的个数$tot$，然后区间赋为1.
然后统计$l1到r1$区间内1的个数$Tot$，如果$tot>Tot$就直接区间赋值了嘛。
其他情况就在$l1到r1$内二分一下然后区间赋值。
复杂度$O(nlog^2_n)$
打得轻松加愉快！$wori$怎么T了！$lz$明明最极限的数据都$2s$卡过去了,$EXM！！？$
然后就是$DT$的=数据时间。。
测出来还是$T$了一组$QAQ$，妈蛋被卡常了$QnQ$。哎两个$log$线段树还打得丑没得法。。
发现这个题区间最大连续和比较特殊。。只可能两个值。。改了一下写法卡卡常数终于$18s$混过去了。

【代码】

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#define N 200005
#define Mod 1000000007
#define INF 200001
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const ull base=31;

ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

int n,m;

class Seg_Tree{
    public:
        int l,r,tag,len,sum1,lsum,rsum,msum;
}e[N<<2]; 

void Build(int p,int l,int r)
{
    e[p].l=l,e[p].r=r,e[p].len=r-l+1;
    if(l==r) return;
    int mid=l+r>>1;
    Build(p<<1,l,mid);Build(p<<1|1,mid+1,r);
}

void pushup(Seg_Tree &p,Seg_Tree L,Seg_Tree R)
{
    p.sum1=L.sum1+R.sum1;
    p.lsum=L.lsum;p.rsum=R.rsum;
    if(L.sum1==L.len) p.lsum=max(p.lsum,L.sum1+R.lsum);
    if(R.sum1==R.len) p.rsum=max(p.rsum,R.sum1+L.rsum);
    p.msum=max(L.rsum+R.lsum,max(L.msum,R.msum));
}

void pushdown(int k)
{
    Seg_Tree &p=e[k];Seg_Tree &L=e[k<<1];Seg_Tree &R=e[k<<1|1];
    int t=p.tag;p.tag=0,L.tag=R.tag=t;
    L.sum1=L.msum=L.lsum=L.rsum=(t==1)?L.len:0;
    R.sum1=R.msum=R.lsum=R.rsum=(t==1)?R.len:0;
}

void Update(int p,int x,int y,int z)
{
    if(x>y) return;
    int l=e[p].l,r=e[p].r,mid=l+r>>1;
    if(l==x&&y==r) {
        e[p].tag=z;
        e[p].sum1=e[p].lsum=e[p].rsum=e[p].msum=(z==1)?e[p].len:0;
        return;
    }
    if(e[p].tag) pushdown(p);
    if(y<=mid) Update(p<<1,x,y,z);
    else if(x>mid) Update(p<<1|1,x,y,z);
    else Update(p<<1,x,mid,z),Update(p<<1|1,mid+1,y,z);
    pushup(e[p],e[p<<1],e[p<<1|1]);
}

void Solve1()
{
    static int x,y;x=read(),y=read();
    if(x>y) swap(x,y);
    Update(1,x,y,1);
}

int Find1(int p,int x,int y)
{
    int l=e[p].l,r=e[p].r,mid=l+r>>1;
    if(l==x&&y==r) return e[p].sum1;
    if(e[p].tag) pushdown(p);
    if(y<=mid) return Find1(p<<1,x,y);
    if(x>mid) return Find1(p<<1|1,x,y);
    return Find1(p<<1,x,mid)+Find1(p<<1|1,mid+1,y);
}

int Find_Pos(int L,int R,int x)
{
    int rtn=0,st=L;
    while(L<=R)
    {
        int mid=L+R>>1;
        if(Find1(1,st,mid)<=x) rtn=mid,L=mid+1;
        else R=mid-1;
    }
    return rtn;
}

void Solve2()
{
    static int x,y,xx,yy;
    x=read(),y=read(),xx=read(),yy=read();
    if(x>y) swap(x,y);if(xx>yy) swap(xx,yy);
    int tot=y-x+1-Find1(1,x,y);
    Update(1,x,y,1);
    int Tot=Find1(1,xx,yy);
    if(tot<Tot) yy=Find_Pos(xx,yy,tot);
    Update(1,xx,yy,-INF);
}

Seg_Tree Query(int p,int x,int y)
{
    int l=e[p].l,r=e[p].r,mid=l+r>>1;
    if(l==x&&y==r) return e[p];
    if(e[p].tag) pushdown(p);
    if(y<=mid) return Query(p<<1,x,y);
    if(x>mid) return Query(p<<1|1,x,y);
    Seg_Tree L=Query(p<<1,x,mid),R=Query(p<<1|1,mid+1,y);
    Seg_Tree rtn;pushup(rtn,L,R);rtn.l=L.l,rtn.r=R.r;
    return rtn;
}

void Solve3()
{
    static int x,y;x=read(),y=read();
    if(x>y) swap(x,y);
    Seg_Tree rtn=Query(1,x,y);
    printf("%d\n",rtn.msum);
}

int main()
{
    n=read();m=read();
    Build(1,1,n);
    for(int i=1;i<=m;i++)
    {
        static int f;f=read();
        if(f==0) Solve1();
        else if(f==1) Solve2();
        else Solve3();
    }
    return 0;
}