洛谷P4843 清理雪道【有上下界最小流】

最新推荐文章于 2026-04-19 02:42:20 发布

原创最新推荐文章于 2026-04-19 02:42:20 发布 · 226 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

网络流同时被 2 个专栏收录

22 篇文章

订阅专栏

最大流

9 篇文章

订阅专栏

本文介绍了一种基于网络流算法的滑雪场清理路径优化方案，通过建立有向无环图和应用最小流算法，实现最少飞行次数清理所有雪道的目标。

时空限制 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;
}