3218: a + b Problem

题意自己看吧。。
图这个大就不贴了

感想

好神的题啊。。
第一次做这种题。。
感觉好劲啊。。
网络流构图原来还可以这样

solution

首先，这题我们先把他转为最小割
st–>x 权为b，代表黑色
x—>ed 权为w 代表白色
然后怎么处理这个p呢？
直接x–>y p（当然x，y合法）
当然不是，因为你这样会有一个问题，就是多个的话你的p会算多次，也就是x连了很多个y啦
这当然不是我们想要的
那怎么办呢？
对于这种多算问题，拆点总是一个很好的选择
考虑拆点x—>x+n
那么我们将x+n—>y就好了
写成代码张这样

init(st,u,b[u]);init(u,ed,w[u]);init(u,u+n,p[u]);
for (int i=1;i<u;i++)
    if (a[i]>=l[u]&&a[i]<=r[u])
        init(u+n,i,MAX);

但是你会发现一个问题，就是你在uoj上面只有60分。。
为什么？
因为这样的话很容易卡啊
我可以随便就让你一个点连了之前所有点
于是你的时空复杂度就是O(n^2)的了

怎么破？

我们观察到，要是一个点能使得另外一个点变成奇怪点，a值总在一个区间里面
那我们就像，可以不可以把这个区间做成一个点，然后让在这个范围里面的a连过来
这不就是线段树嘛！
每次就相当于让线段树相应区间的点连过来
然后包含a的，a就连边过去
这不就是一个理想的方案了嘛！

等等，我们好像没有考虑只受前面的点的影响！！
但这没关系，我们开一个可持久化线段树就好了啊！！

这样点数虽然变大了 2*n—->nlogn
但是我们的边数大大变小了 n^2——>nlogn

这样就可以跑过去了

#include<cstdio>
#include<iostream>
#include<queue>
#include<algorithm>
#include<cstring>
using namespace std;
const int N=5010;
const int MAX=1<<30;
int n;
int a[N],b[N],w[N],l[N],r[N],p[N];
struct qq{int x,y,z,last;}s[N*210];
int num,last[N*40];
void init (int x,int y,int z)
{
    if (x==0||y==0) return ;
    num++;
    s[num].x=x;s[num].y=y;s[num].z=z;
    s[num].last=last[x];
    last[x]=num;
    swap(x,y);z=0;
    num++;
    s[num].x=x;s[num].y=y;s[num].z=z;
    s[num].last=last[x];
    last[x]=num;
}
int st,ed;
int h[N*40];
bool Bt ()
{
    memset(h,-1,sizeof(h));h[st]=1;
    queue<int> q;
    q.push(st);
    while (!q.empty())
    {
        int x=q.front();q.pop();
        for (int u=last[x];u!=-1;u=s[u].last)
        {
            int y=s[u].y;
            if (s[u].z>0&&h[y]==-1)
            {
                h[y]=h[x]+1;
                q.push(y);
            }
        }
    }
    return h[ed]!=-1;
}
int mymin (int x,int y){return x<y?x:y;}
int dfs (int x,int f)
{
    if (x==ed) return f;
    int s1=0;
    for (int u=last[x];u!=-1;u=s[u].last)
    {
        int y=s[u].y;
        if (s[u].z>0&&h[y]==h[x]+1&&s1<f)
        {
            int lalal=dfs(y,mymin(s[u].z,f-s1));
            s1+=lalal;
            s[u].z-=lalal;
            s[u^1].z+=lalal;
        }
    }
    if (s1==0) h[x]=-1;
    return s1;
}
int s1[N*40],s2[N*40];
int root[N],lalal;
void add (int l,int r,int x,int y,int a,int id)
{
    init(y,x,MAX);
    init(y,id,MAX);
    if (l==r) return ;
    int mid=(l+r)>>1;
    if (a<=mid)
    {
        s1[y]=++lalal;s2[y]=s2[x];
        add(l,mid,s1[x],s1[y],a,id);
    }
    else
    {
        s2[y]=++lalal;s1[y]=s1[x];
        add(mid+1,r,s2[x],s2[y],a,id);
    }
}
void add2 (int l,int r,int now,int L,int R,int id)
{
    if (now==0) return ;
    if (l==L&&r==R)
    {
        init(id,now,MAX);
        return ;
    }
    int mid=(l+r)>>1;
    if (R<=mid) add2(l,mid,s1[now],L,R,id);
    else if (L>mid) add2(mid+1,r,s2[now],L,R,id);
    else add2(l,mid,s1[now],L,mid,id),add2(mid+1,r,s2[now],mid+1,R,id);
}
int main()
{
    int ans=0;
    num=1;memset(last,-1,sizeof(last));
    scanf("%d",&n);
    st=2*n+1,ed=st+1;lalal=ed;
    for (int u=1;u<=n;u++)
    {
        scanf("%d%d%d%d%d%d",&a[u],&b[u],&w[u],&l[u],&r[u],&p[u]);
        ans=ans+w[u]+b[u];
        init(st,u,b[u]);init(u,ed,w[u]);init(u,u+n,p[u]);
        root[u]=++lalal;
        add(0,1e9,root[u-1],root[u],a[u],u);
        add2(0,1e9,root[u],l[u],r[u],u+n);
    }
    while (Bt()==true) ans=ans-dfs(st,MAX);
    printf("%d\n",ans);
    return 0;
}