网络流24题1 飞行员配对方案问题 ssl 2601

17 篇文章 0 订阅
14 篇文章 0 订阅

这里写图片描述

分析

用网络流写二分图匹配,建一个点s,连接二分图的一边,一个点t,连接二分图的另一边。整个图的每一条边的流量都是1。
从s到t的最大流就是ans。

网络流的方法:
Dinic
这个算法的效率是相当高的,算法流程如下:
1、每次以源点为起始点bfs,求出每个点的编号d[i],d[i]表示从源点到i点通过至少几条残余流量大于0的边,能够到达i点。
2、只有那些满足d[u]+1=d[v]的边(u,v)才被视为存在,然后在这里面不断DFS找增广路并增广(其实随便走都可以),如果没有增广路了,那么返回步骤1,如果BFS不到汇点,证明算法结束。

code

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<string>
#include<algorithm>

using namespace std;

struct arr{
    int x,y,w;
    int op;
    int next;
}edge[1000000];
int ls[4000];
int v[4000];
int w[4000];
int f[4000];
int queue[5000];
int dis[4000];
int s,t;
int n,m,nn;
int min1;
int ans;

void er(int ww)
{
    w[nn]=ww;
}

void add(int x,int y,int w)
{
    nn++;
    edge[nn].x=x;
    edge[nn].y=y;
    edge[nn].w=w;
    edge[nn].op=nn+1;
    edge[nn].next=ls[x];
    ls[x]=nn;
    v[x]=nn;
    er(w);

    nn++;
    edge[nn].x=y;
    edge[nn].y=x;
    edge[nn].w=w;
    edge[nn].op=nn-1;
    edge[nn].next=ls[y];
    ls[y]=nn;
    v[y]=nn;
    er(0);
}

bool bfs()
{
    memset(dis,-1,sizeof(dis));
    memset(queue,0,sizeof(queue));
    int head,tail;
    head=0; tail=1;
    queue[head]=0;
    dis[0]=0;
    do
    {
        head++;
        int i=ls[queue[head]];
        while (i!=0)
        {
            if ((w[i]>0)&&(dis[edge[i].y]==-1))
            {
                dis[edge[i].y]=dis[edge[i].x]+1;
                tail++;
                queue[tail]=edge[i].y;
                if (edge[i].y==t)  return true;
            }
            i=edge[i].next;
        }
    }while (head<tail);
    return false;
}

bool find(int x,int num)
{
    if (x==t) return true;
    if (x!=s) min1=min(min1,w[num]);
    for (int i=v[x];i!=0;i=edge[i].next)
    {
        int x=edge[i].x;
        int y=edge[i].y;
        if ((dis[x]+1==dis[y])&&(w[i]!=0)&&(find(y,i)))
        {
            f[x]=y;
            w[i]-=min1;
            w[edge[i].op]+=min1;
            return true;
        }
    }
    return false;
}

void dinic()
{
    ans=0;
    while (bfs())
    {
        min1=2000000000;
        find(s,0);
        ans+=min1;
    }
}

int main()
{
    scanf("%d%d",&n,&m);
    if ((n==-1)&&(m==-1)) 
        return 0;
    int i=0;
    int x,y;
    while (1==1)
    {
        i++;
        scanf("%d%d",&x,&y);
        if ((x==-1)&&(y==-1))
          break;
        add(x,y,1);
    }
    s=0;
    t=n+m+1;
    for (int i=1;i<=n;i++)
        add(0,i,1);
    for (int i=n+1;i<=n+m;i++)
        add(i,t,1);
    dinic();
    if (ans==0)
        printf("No Solution!");
    else
        printf("%d\n",ans); 
    return 0;
}
网络算法可以用来解决最大最小割问题,而飞行员配对问题可以转化成最大问题。下面是一个基于 Python 的网络算法求解飞行员配对方案问题的示例代码: ```python from collections import defaultdict class Graph: def __init__(self): self.graph = defaultdict(dict) def add_edge(self, u, v, w): self.graph[u][v] = w def bfs(self, s, t, parent): visited = set() queue = [s] visited.add(s) while queue: u = queue.pop(0) for v in self.graph[u]: if v not in visited and self.graph[u][v] > 0: queue.append(v) visited.add(v) parent[v] = u return True if t in visited else False def max_flow(self, source, sink): parent = {} max_flow = 0 while self.bfs(source, sink, parent): path_flow = float('inf') s = sink while s != source: path_flow = min(path_flow, self.graph[parent[s]][s]) s = parent[s] max_flow += path_flow v = sink while v != source: u = parent[v] self.graph[u][v] -= path_flow self.graph[v][u] += path_flow v = parent[v] return max_flow def pilot_pairing(pilots, planes): num_pilots = len(pilots) num_planes = len(planes) graph = Graph() source = 's' sink = 't' for i in range(num_pilots): graph.add_edge(source, 'P{}'.format(i), 1) for j in range(num_planes): if pilots[i][0] <= planes[j][0] and pilots[i][1] >= planes[j][1]: graph.add_edge('P{}'.format(i), 'A{}'.format(j), 1) graph.add_edge('A{}'.format(i), sink, 1) return graph.max_flow(source, sink) # 示例用法 pilots = [(1, 5), (2, 4), (3, 6), (2, 5)] planes = [(1, 4), (2, 5), (3, 6)] print(pilot_pairing(pilots, planes)) # 输出为 3 ``` 在上面的示例代码中,我们定义了一个 `Graph` 类来表示图,其中 `add_edge` 方法用于添加边,`bfs` 方法用于执行广度优先搜索,`max_flow` 方法用于计算最大。 在 `pilot_pairing` 函数中,我们首先创建了一个 `Graph` 对象,然后为源点和汇点分别添加一条边,并遍历所有的飞行员和飞机,如果某个飞行员可以驾驶某个飞机,则在他们之间添加一条边。最后调用 `max_flow` 方法计算最大。 在 `max_flow` 方法中,我们首先执行广度优先搜索来寻找一条增广路径,然后计算路径上的最小剩余容量,更新路径上的边的量,并更新最大的值。重复执行这个过程,直到没有增广路径为止。最后返回最大的值。 在上面的示例中,我们输入了 4 个飞行员和 3 个飞机,输出得到了最大为 3,说明有 3 对飞行员和飞机可以配对成功。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值