bzoj3218 a+b Problem（最小割+主席树优化建边）

由于6.22博主要学测，大半时间学文化课，近期刷题量&写题解的数量会急剧下降。

这题出得挺经典的，首先一眼最小割，考虑朴素的做法：与S联通表示白色，与T联通表示黑色，S向i连流量为w[i]的边，i向T连流量为b[i]的边，然后i'向i连容量为p[i]的边，所有满足条件的j向i'连一条容量为无穷大的边（只要满足其一就要割掉）。然后边数显然不合法，一眼线段树优化，然而发现线段树无法连边，考虑主席树连边，对原序列建主席树，每个点的区间对应主席树上的O(logn)个节点，然后连接容量为无穷大的边即可。

#include<bits/stdc++.h>
using namespace std;
const int N=2e5+,M=1e6+,inf=1e9;
struct edge{int x,id;}e[N];
int n,m,ecnt,tot,T,ans,mx;
int L[N],R[N],val[N],seq[N],hd[N],v[M],w[M],nxt[M],q[N],lv[N],lc[N],rc[N],root[N],cnt[N];
void adde(int x,int y,int z)
{
    v[++ecnt]=y,w[ecnt]=z,nxt[ecnt]=hd[x],hd[x]=ecnt;
    v[++ecnt]=x,w[ecnt]=,nxt[ecnt]=hd[y],hd[y]=ecnt;
}
int findl(int x)
{
    int l=,r=mx,ans=mx;
    while(l<=r)
    {
        int mid=l+r>>;
        if(seq[mid]>=x)ans=mid,r=mid-;else l=mid+;
    }
    return ans;
}
int findr(int x)
{
    int l=,r=mx,ans=;
    while(l<=r)
    {
        int mid=l+r>>;
        if(seq[mid]<=x)ans=mid,l=mid+;else r=mid-;
    }
    return ans;
}
int update(int prt,int l,int r,int i)
{
    int rt=++tot;
    if(l==r)lc[rt]=rc[rt]=,cnt[rt]=cnt[prt]+;
    else{
        int mid=l+r>>;
        if(val[i]<=mid)lc[rt]=update(lc[prt],l,mid,i),rc[rt]=rc[prt];
        else rc[rt]=update(rc[prt],mid+,r,i),lc[rt]=lc[prt];
        cnt[rt]=cnt[lc[rt]]+cnt[rc[rt]];
    }
    adde(i,*n++rt,inf);
    adde(*n++prt,*n++rt,inf);
    return rt;
}
void query(int rt,int l,int r,int i)
{
    if(!cnt[rt])return;
    if(L[i]<=l&&r<=R[i]){adde(*n++rt,n+i,inf);return;}
    int mid=l+r>>;
    if(L[i]<=mid)query(lc[rt],l,mid,i);
    if(R[i]>mid)query(rc[rt],mid+,r,i);
}
bool bfs()
{
    memset(lv,-,sizeof lv);
    int qs=,qe=;
    q[]=lv[]=;
    while(qs<qe)
    {
        int u=q[qs++];
        for(int i=hd[u];i;i=nxt[i])if(w[i]&&lv[v[i]]==-)lv[v[i]]=lv[u]+,q[qe++]=v[i];
    }
    if(lv[T]==-)return ;
    return ;
}
int dfs(int u,int low)
{
    if(u==T||!low)return low;
    int sum=;
    for(int i=hd[u];i;i=nxt[i])
    if(w[i]&&lv[v[i]]==lv[u]+)
    {
        int tmp=dfs(v[i],min(low,w[i]));
        w[i]-=tmp,w[i^]+=tmp,low-=tmp,sum+=tmp;
        if(!low)return sum;
    }
    if(low)lv[u]=-;
    return sum;
}
bool cmp(edge x,edge y){return x.x<y.x;}
int main()
{
    scanf("%d",&n);
    ecnt=,T=1e5;
    for(int i=,b,w,p;i<=n;i++)
    {
        scanf("%d%d%d%d%d%d",&e[i].x,&b,&w,&L[i],&R[i],&p);
        e[i].id=i,ans+=b+w,adde(,i,w),adde(i,T,b),adde(n+i,i,p);
    }
    sort(e+,e+n+,cmp);
    mx=val[e[].id]=,seq[]=e[].x;
    for(int i=;i<=n;i++)
    if(e[i].x==e[i-].x)val[e[i].id]=mx;else val[e[i].id]=++mx,seq[mx]=e[i].x;
    for(int i=;i<=n;i++)
    {
        L[i]=findl(L[i]),R[i]=findr(R[i]);
        root[i]=update(root[i-],,mx,i),query(root[i-],,mx,i);
    }
    while(bfs())ans-=dfs(,inf);
    printf("%d",ans);
}