hihocoder1398 网络流五之最大权闭合子图
虽然我自己现在总结不好最大权闭合子图。但也算稍稍理解辣。
网络流起步ing~~~(～￣▽￣)～
#include<iostream>
#include<cstdio>
#include<queue>
#include<algorithm>
using namespace std;
struct node
{
int point;
int nxt;
int value;
};
node line[500000];
int head[5000],tail=-1;
int n,m;
void add(int x,int y,int z)
{
line[++tail].point=y;
line[tail].value=z;
line[tail].nxt=head[x];
head[x]=tail;
line[++tail].point=x;
line[tail].value=0;
line[tail].nxt=head[y];
head[y]=tail;
}
int a[5000],b[5000];
int dep[50000],cur[50000];
bool BFS(int begin,int end)
{
for(int i=begin;i<=end;i++)
{
dep[i]=0;
cur[i]=head[i];
}
queue<int>q;
q.push(begin);
dep[begin]=1;
while(!q.empty())
{
int pas=q.front();
q.pop();
for(int i=head[pas];i!=-1;i=line[i].nxt)
if(!dep[line[i].point]&&line[i].value)
{
dep[line[i].point]=dep[pas]+1;
q.push(line[i].point);
}
}
if(dep[end])
return true;
return false;
}
int dfs(int now,int aim,int limte)
{
if(now==aim||!limte)
return limte;
int flow=0,f;
for(int i=cur[now];i!=-1;i=line[i].nxt)
{
cur[now]=i;
if(dep[line[i].point]==dep[now]+1&&(f=dfs(line[i].point,aim,min(limte,line[i].value))))
{
flow+=f;
limte-=f;
line[i].value-=f;
line[i^1].value+=f;
if(!limte)
break;
}
}
return flow;
}
int dinic(int begin,int end)
{
int res=0;
while(BFS(begin,end))
res+=dfs(begin,end,0x7fffffff);
return res;
}
int main()
{
int ans=0;
scanf("%d%d",&n,&m);
for(int i=0;i<=n+m+1;i++)
head[i]=-1;
for(int i=1;i<=m;i++)
scanf("%d",&b[i]);
int num,dat;
for(int i=1;i<=n;i++)
{
scanf("%d%d",&a[i],&num);
ans+=a[i];
for(int j=1;j<=num;j++)
{
scanf("%d",&dat);
add(i,n+dat,0x7fffffff);
}
}
for(int i=1;i<=n;i++)
add(0,i,a[i]);
for(int i=n+1;i<=n+m;i++)
add(i,n+m+1,b[i-n]);
ans-=dinic(0,n+m+1);
printf("%d",ans);
}
我的代码，让你感到真真的实惠石灰（逃）
