滑雪场坐落在FJ省西北部的若干座山上。
从空中鸟瞰,滑雪场可以看作一个有向无环图,每条弧代表一个斜坡(即雪道),弧的方向代表斜坡下降的方向。
你的团队负责每周定时清理雪道。你们拥有一架直升飞机,每次飞行可以从总部带一个人降落到滑雪场的某个地点,然后再飞回总部。从降落的地点出发,这个人可以顺着斜坡向下滑行,并清理他所经过的雪道。
由于每次飞行的耗费是固定的,为了最小化耗费,你想知道如何用最少的飞行次数才能完成清理雪道的任务。
LYD:
类似 <有源汇上下界可行流> 的构图方法,但是不添加T到S的边,求一次超级源到超级汇的最大流。
加边(T,S,0,+∞),在上一步残量网络基础上再求一次超级源到超级汇的最大流。
流经T到S的边的流量就是最小流的值。
该算法的思想是在第一步中尽可能填充循环流,以减小最小流的代价。
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
#define inf 0x3f3f3f3f
#define N 155
#define M 400005
using namespace std;
int head[N],ver[M],nxt[M],f[M],tot,ch[N];
void add(int a,int b,int c)
{
tot++;nxt[tot]=head[a];head[a]=tot;ver[tot]=b;f[tot]=c;
tot++;nxt[tot]=head[b];head[b]=tot;ver[tot]=a;f[tot]=0;
return ;
}
queue<int>q;int S,T;
bool tell()
{
memset(ch,-1,sizeof(ch));
q.push(S);ch[S]=0;
while(!q.empty())
{
int tmp=q.front();q.pop();
for(int i=head[tmp];i;i=nxt[i])
{
if(f[i]&&ch[ver[i]]==-1)
{
ch[ver[i]]=ch[tmp]+1;
q.push(ver[i]);
}
}
}
return ch[T]!=-1;
}
int zeng(int a,int b)
{
if(a==T)return b;
int r=0;
for(int i=head[a];i&&b>r;i=nxt[i])
{
if(f[i]&&ch[ver[i]]==ch[a]+1)
{
int t=zeng(ver[i],min(f[i],b-r));
f[i]-=t;f[i^1]+=t;r+=t;
}
}
if(!r)ch[a]=-1;
return r;
}
int dinic()
{
int r=0,t;
while(tell())while(t=zeng(S,inf))r+=t;
return r;
}
int n;
int ru[N],chu[N];
int main()
{
scanf("%d",&n);tot=1;
int SS=n+1,TT=n+2;
for(int i=1;i<=n;i++)
{
int tmp,t;
scanf("%d",&tmp);
for(int j=1;j<=tmp;j++)
{
scanf("%d",&t);
add(i,t,inf);
chu[i]++;ru[t]++;
}
}
for(int i=1;i<=n;i++)
{
add(SS,i,inf);add(i,TT,inf);
}
S=n+3,T=n+4;
for(int i=1;i<=n;i++)
{
if(chu[i]>ru[i])
{
add(i,T,chu[i]-ru[i]);
}
else add(S,i,ru[i]-chu[i]);
}
dinic();
add(TT,SS,inf);
dinic();
printf("%d\n",f[tot]);
return 0;
}