时空限制 1000ms / 128MB
题目描述
滑雪场坐落在FJ省西北部的若干座山上。
从空中鸟瞰,滑雪场可以看作一个有向无环图,每条弧代表一个斜坡(即雪道),弧的方向代表斜坡下降的方向。
你的团队负责每周定时清理雪道。你们拥有一架直升飞机,每次飞行可以从总部带一个人降落到滑雪场的某个地点,然后再飞回总部。从降落的地点出发,这个人可以顺着斜坡向下滑行,并清理他所经过的雪道。
由于每次飞行的耗费是固定的,为了最小化耗费,你想知道如何用最少的飞行次数才能完成清理雪道的任务。
输入格式:
输入文件的第一行包含一个整数n (2 <= n <= 100) – 代表滑雪场的地点的数量。
接下来的n行,描述1~n号地点出发的斜坡,第i行的第一个数为mi (0 <= mi < n) ,后面共有mi个整数,由空格隔开,每个整数aij互不相同,代表从地点i下降到地点aij的斜坡。每个地点至少有一个斜坡与之相连。
输出格式:
输出文件的第一行是一个整数k – 直升飞机的最少飞行次数。
题目分析
算是有上下界最小流裸题吧
对于每条雪道,至少清理一次,可以反复经过,所以原图中每条边在网络中下界为1,上界为inf
建立虚拟源点s,向每个点连边,下界0,上界inf
每个点向虚拟汇点t连边,下界0,上界inf
分别表示每个点可以反复降落,也可以在这里结束清理
按有上下界的网络流建图即可
(因为一些弱智错误调了一万年)
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
using namespace std;
typedef long long lt;
int read()
{
int f=1,x=0;
char ss=getchar();
while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
return f*x;
}
const int inf=1e9;
const int maxn=50010;
int n;
struct node{int v,f,nxt;}E[maxn<<1];
int head[maxn],tot=1;
int deg[maxn],lev[maxn];
void add(int u,int v,int f)
{
E[++tot].nxt=head[u];
E[tot].v=v; E[tot].f=f;
head[u]=tot;
}
int bfs(int s,int t)
{
queue<int> q; q.push(s);
memset(lev,-1,sizeof(lev)); lev[s]=0;
while(!q.empty())
{
int u=q.front(); q.pop();
for(int i=head[u];i;i=E[i].nxt)
{
int v=E[i].v;
if(lev[v]==-1&&E[i].f)
{
lev[v]=lev[u]+1;
if(v==t) return 1;
q.push(v);
}
}
}
return 0;
}
int dfs(int u,int cap,int t)
{
if(u==t) return cap;
int flow=cap;
for(int i=head[u];i;i=E[i].nxt)
{
int v=E[i].v;
if(lev[v]==lev[u]+1&&flow&&E[i].f)
{
int f=dfs(v,min(flow,E[i].f),t);
E[i].f-=f; E[i^1].f+=f;
flow-=f;
}
}
return cap-flow;
}
int dicnic(int s,int t)
{
int maxf=0;
while(bfs(s,t)) maxf+=dfs(s,inf,t);
return maxf;
}
int main()
{
n=read();
for(int u=1;u<=n;++u)
{
int m=read();
while(m--)
{
int v=read();
add(u,v,inf); add(v,u,0);//这里inf-1还是inf
deg[u]--; deg[v]++;
}
}
int s=n+1,t=n+2;
for(int i=1;i<=n;++i)
{
add(s,i,inf); add(i,s,0);
add(i,t,inf); add(t,i,0);
}
int ss=n+3,tt=n+4;
for(int i=1;i<=n;++i)
{
if(deg[i]>0) add(ss,i,deg[i]),add(i,ss,0);
else if(deg[i]<0) add(i,tt,-deg[i]),add(tt,i,0);
}
add(t,s,inf); add(s,t,0);
dicnic(ss,tt);
int ans=E[tot].f;
head[t]=E[head[t]].nxt;
head[s]=E[head[s]].nxt;
printf("%d",ans-dicnic(t,s));
return 0;
}