UOJ#77 bzoj3218 a + b Problem 可持久化线段树优化建图+网络流

题目大意：
有n个点，如果第i点染黑得到bi的收益，染白得wi的收益。
每个点有属性ai，li，ri，pi。
如果第i个点染成黑色，并且存在一个点j满足
第j个点染白色,j< i 且 li<=aj<=ri
则i为奇怪的点，付出pi的代价。
求总收益最大。
n<=1e5

题目分析：
可以看出一个网络流的模型。
从源点向每个点连bi，每个点向汇点连wi。
对于第i个点，新建一个结点，从i向这个结点连pi，然后从这个点向所有能使i变成奇怪的点的结点连INF。
用所有b和w的和减去最大流即是答案。

但是n太大了，我们发现这样建边是n^2的。
所以需要优化一下。
我们对于所有的ai建一棵线段树，在叶子结点向对应的点连INF，其他结点向儿子连INF。
然后对于每个点，只要在线段树上查询l~r，然后连INF就行了。
但是由于要求 j < i ，所以把线段树可持久化一下。新建的结点向上一个版本的结点连INF即可。
然后点和边就都变成nlogn级别的了。

代码如下：

#include <cstdio>
#include <iostream>
#define MAXN 1000000000
#define MAXM 5200
#define N 520000
#define M 1200000
using namespace std;
const int INF=0x3f3f3f3f;
inline int Min(int x,int y) { return x<y?x:y; }
int n,S,T,ans;
int fir[N],nes[M],v[M],q[M],tot=1;
int dl[N],d[N],b[MAXM],w[MAXM];
struct segment{
    segment *ch[2];
    int id;
    void* operator new (size_t size,segment *x);
}*root[MAXM];
void* segment :: operator new(size_t size,segment *x)
{
    static segment *C=NULL,*mempool=NULL;
    if(C==mempool) mempool=(C=new segment[1<<15])+(1<<15);
    if(x==NULL) C->ch[0]=C->ch[1]=NULL;
    else C->ch[0]=x->ch[0],C->ch[1]=x->ch[1];
    C->id=T++;
    return C++;
}
void edge(int x,int y,int z)
{
    v[++tot]=y;
    q[tot]=z;
    nes[tot]=fir[x];
    fir[x]=tot;
}
void Edge(int x,int y,int z)
{
    edge(x,y,z);
    edge(y,x,0);
}
void update(segment *&c,int l,int r,int x,int y)
{
    int tmp=c?c->id:0;
    c=new (c) segment;
    if(tmp) Edge(c->id,tmp,INF);
    Edge(c->id,y,INF);
    if(l==r)
    {
        Edge(c->id,y,INF);
        return;
    }
    int mid=l+r>>1;
    if(x<=mid) update(c->ch[0],l,mid,x,y),Edge(c->id,c->ch[0]->id,INF);
    else update(c->ch[1],mid+1,r,x,y),Edge(c->id,c->ch[1]->id,INF);
}
void query(segment *c,int l,int r,int x,int y)
{
    if(c==NULL) return;
    if(x<=l && y>=r)
    {
        Edge(T,c->id,INF);
        return;
    }
    int mid=l+r>>1;
    if(x<=mid) query(c->ch[0],l,mid,x,y);
    if(y>mid) query(c->ch[1],mid+1,r,x,y);
}
bool bfs()
{
    static int c;
    int l=1,r=1;
    for(int i=S;i<=T;i++) d[i]=0;
    d[S]=1; dl[1]=S;
    while(l<=r)
    {
        c=dl[l++];
        for(int t=fir[c];t;t=nes[t])
        {
            if(!q[t] || d[v[t]]) continue;
            d[v[t]]=d[c]+1;
            dl[++r]=v[t];
            if(v[t]==T) return true;
        }
    }
    return false;
}
int dfs(int c,int flow)
{
    if(c==T || flow==0) return flow;
    int ans=0,tmp;
    for(int t=fir[c];t;t=nes[t])
    {
        if(!q[t] || d[v[t]]!=d[c]+1) continue;
        tmp=dfs(v[t],Min(q[t],flow));
        q[t]-=tmp; q[t^1]+=tmp;
        flow-=tmp; ans+=tmp;
        if(!flow) break;
    }
    if(!ans) d[c]=-1;
    return ans;
}
int dinic()
{
    int ans=0;
    while(bfs()) ans+=dfs(S,INF);
    return ans;
}
int main()
{
    scanf("%d",&n);
    T=n+1;
    for(int i=1,a,l,r,p;i<=n;i++)
    {
        scanf("%d%d%d%d%d%d",&a,&b[i],&w[i],&l,&r,&p);
        Edge(i,T,p);
        query(root[i-1],0,MAXN,l,r);
        T++;
        root[i]=root[i-1];
        update(root[i],0,MAXN,a,i);
    }
    for(int i=1;i<=n;i++)
    {
        Edge(S,i,b[i]);
        Edge(i,T,w[i]);
        ans+=b[i]+w[i];
    }
    printf("%d\n",ans-dinic());
    return 0;
}