BZOJ 3218/UOJ #77 a + b Problem

VFK出的神题啊。。。
还是挺明显是最小割的吧。。
裸的最小割5000*5000会boom！
所以就要进行神奇的优化构图 神TM限制 j < i于是就要打可持久化

Po姐前面A+B的代码打的吼啊！滑稽
传送门 其实是图画的很好辣。。 讲的也挺清晰的
我就偷偷懒吧。。。

代码又打神TM短 会不会被吐槽。。这次代码并不快啊。。QAQ

#include<bits/stdc++.h>
using namespace std;
const int M=100005,N=5005,inf=1e9+7;
inline int read(){
    char ch=getchar(); int x=0,f=1;
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0'; ch=getchar();}
    return x*f;
}
struct E{int y,c,next;}e[M*10]; int len,first[M];
void ins(int x,int y,int c){
    e[++len]=(E){y,c,first[x]},first[x]=len;
    e[++len]=(E){x,0,first[y]},first[y]=len;
}
int h[M],t,tot; queue<int>q;
bool bfs(){
    for(int i=0;i<=tot;i++)h[i]=0;
    q.push(0); h[0]=1;
    while(!q.empty()){
        int x=q.front(); q.pop();
        for(int k=first[x];k;k=e[k].next){
            int y=e[k].y;
            if(!h[y] && e[k].c){h[y]=h[x]+1; q.push(y);}
        }
    }
    return h[t]>0;
}
int dfs(int x,int f){
    if(x==t)return f;
    int de=0;
    for(int k=first[x];k;k=e[k].next){
        int y=e[k].y;
        if(h[y]==h[x]+1){
            int sf=dfs(y,min(f-de,e[k].c));
            e[k].c-=sf,e[k^1].c+=sf,de+=sf;
            if(de==f)return f;
        }
    }
    if(!de)h[x]=0;
    return de;
}
int lc[M],rc[M],rt[M],a[N],L[N],R[N],Q[N*3];
void link(int x,int l,int r,int i)
{
    if(!x || L[i]>r || R[i]<l)return;
    if(L[i]<=l && r<=R[i])return ins(x,i*2+1,inf);
    int mid=(l+r)>>1;
    link(lc[x],l,mid,i),link(rc[x],mid+1,r,i);
}
void add(int &x,int u,int l,int r,int i){
    x=++tot;
    if(l==r){
        if(u)ins(u,x,inf);
        return ins(i*2,x,inf);
    }
    int mid=(l+r)>>1;
    if(a[i]<=mid)rc[x]=rc[u],add(lc[x],lc[u],l,mid,i);
    else lc[x]=lc[u],add(rc[x],rc[u],mid+1,r,i);
    if(lc[x])ins(lc[x],x,inf);
    if(rc[x])ins(rc[x],x,inf);
}
int main()
{
    t=1; int b,w,p,ans=0,n=read(),ql=0; len=1;
    for(int i=1;i<=n;i++){
        a[i]=read(),b=read(),w=read(),L[i]=read(),R[i]=read(),p=read();
        ins(0,i*2,w),ins(i*2,t,b),ins(i*2+1,i*2,p);
        ans+=b+w;
        Q[++ql]=L[i],Q[++ql]=R[i],Q[++ql]=a[i];
    }
    sort(Q+1,Q+1+ql);
    ql=unique(Q+1,Q+1+ql)-Q-1; tot=n*2+1;
    for(int i=1;i<=n;i++){
        a[i]=lower_bound(Q+1,Q+1+ql,a[i])-Q;
        L[i]=lower_bound(Q+1,Q+1+ql,L[i])-Q;
        R[i]=lower_bound(Q+1,Q+1+ql,R[i])-Q;
        if(i>1)link(rt[i-1],1,ql,i);
        add(rt[i],rt[i-1],1,ql,i);
    }
    while(bfs())ans-=dfs(0,inf);
    printf("%d\n",ans);
    return 0;
}