BZOJ3218 A+B problem

题目传送门

分析：

首先把公式换一下。。。

ans=sigma(i<=n){w[i]+b[i]} - sigma(i为白色){b[i]} -sigma(i为黑色){w[i]} - sigma(i奇怪){p[i]}

然后我们就需要后面那一块最小。。。

非黑即白是一对矛盾关系，奇怪的点也是一个矛盾关系

于是考虑最小割

首先每个点分别向S和T连边，容量分别是b[i]与w[i]

然后每个点 i 复制出一个 i' 连容量为p[i]的边

i' 再向所有矛盾的 j 连边，容量INF

demo，边数好像会炸

然后发现非法的 j 是区间

然后就线段树连边

编号还要比 i 小。。。

那就只有主席树了

这么强行地把主席树和网络流联系起来吗，无力吐槽2333

8是很懂，这个范围能过。。。

#include<cstdio>
#include<cmath>
#include<queue>
#include<vector>
#include<algorithm>
#include<iostream>
#include<cstring>
  
#define maxn 5005
#define maxm 1000005
#define INF 0x3f3f3f3f
#define opst(i) (((i-1)^1)+1)
 
using namespace std;
  
inline int getint()
{
    int num=0,flag=1;char c;
    while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
    while(c>='0'&&c<='9')num=num*10+c-48,c=getchar();
    return num*flag;
}
 
int n,S,T,N;
int fir[maxn*20],nxt[maxm],to[maxm],cap[maxm],cnt;
int h[20*maxn];
int A[maxn],L[maxn],R[maxn],B[maxn],W[maxn],P[maxn],Q[3*maxn];
struct node{int lc,rc;}t[maxn*20];
int rt[maxn],tot;
int ans,maxflow;
 
inline void newnode(int u,int v,int w)
{to[++cnt]=v,nxt[cnt]=fir[u],fir[u]=cnt,cap[cnt]=w;}
inline void Add(int u,int v,int w)
{newnode(u,v,w),newnode(v,u,0);}
 
inline int find(int num)
{
    int l=1,r=N;
    while(l<r)
    {
        int mid=(l+r)>>1;
        if(Q[mid]>=num)r=mid;
        else l=mid+1;
    }
    return l;
}
 
inline void Link(int now,int l,int r,int i)
{
    if(R[i]<l||L[i]>r)return;
    if(L[i]<=l&&r<=R[i]){Add(i+n,T+now,INF);return;}
    int mid=(l+r)>>1;
    if(t[now].lc)Link(t[now].lc,l,mid,i);
    if(t[now].rc)Link(t[now].rc,mid+1,r,i);
}
 
inline void insert(int &now,int pre,int l,int r,int i)
{
    if(!now||now==pre)now=++tot,t[now]=t[pre];
    if(l==r){Add(T+now,i,INF);if(pre)Add(T+now,T+pre,INF);return;}
    int mid=(l+r)>>1;
    if(A[i]<=mid)insert(t[now].lc,t[pre].lc,l,mid,i);
    else insert(t[now].rc,t[pre].rc,mid+1,r,i);
    if(t[now].lc)Add(T+now,T+t[now].lc,INF);
    if(t[now].rc)Add(T+now,T+t[now].rc,INF);
}
 
inline bool bfs()
{
    memset(h,-1,sizeof h);h[S]=0;
    queue<int>Q;Q.push(S);
    while(!Q.empty())
    {
        int u=Q.front();Q.pop();
        for(int i=fir[u];i;i=nxt[i])
            if(cap[i]&&!~h[to[i]])
                h[to[i]]=h[u]+1,Q.push(to[i]);
    }
    return ~h[T];
}
 
inline int dfs(int u,int flow)
{
    int used=0;
    if(u==T)return flow;
    for(int i=fir[u];i;i=nxt[i])
        if(cap[i]&&h[to[i]]==h[u]+1)
        {
            int w=flow-used;
            w=dfs(to[i],min(w,cap[i]));
            cap[i]-=w,cap[opst(i)]+=w,used+=w;
            if(used==flow)return flow;
        }
    if(!used)h[u]=-1;
    return used;
}
 
inline void dinic()
{while(bfs())maxflow+=dfs(S,INF);}
 
int main()
{
    n=getint();
    for(int i=1;i<=n;i++)
    {
        A[i]=getint(),ans+=B[i]=getint(),ans+=W[i]=getint();
        L[i]=getint(),R[i]=getint(),P[i]=getint();
        Q[++N]=A[i],Q[++N]=L[i],Q[++N]=R[i];
    }
    sort(Q+1,Q+N+1),N=unique(Q+1,Q+N+1)-(Q+1);
    T=2*n+1;
    for(int i=1;i<=n;i++)
    {
        A[i]=find(A[i]);
        L[i]=find(L[i]);
        R[i]=find(R[i]);
        Add(S,i,B[i]),Add(i,T,W[i]),Add(i,i+n,P[i]);
    }
    for(int i=1;i<=n;i++)
    {
        if(i>1)Link(rt[i-1],1,N,i);
        insert(rt[i],rt[i-1],1,N,i);
    }
    dinic();
    printf("%d
",ans-maxflow);
}