题意自己看吧。。
图这个大就不贴了
感想
好神的题啊。。
第一次做这种题。。
感觉好劲啊。。
网络流构图原来还可以这样
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;
}