POJ - 2112 Optimal Milking (dijkstra + 二分 + 最大流Dinic)
题目分析
题意:在一个农场中有k台挤奶器和c只奶牛,每个挤奶器最多只能为m只奶牛挤奶,每个挤奶器和奶牛都视为一个点,将编号1~k记为挤奶器的位置,编号k+1~k+c记为奶牛的位置,奶牛只能在这k+c个位置之间移动,输入将给出每个位置和其余k+c个位置的之间道路距离,其中0代表无法到达
问让所有奶牛进行挤奶的情况下(也就是让每头奶牛都走到一个挤奶器的位置上去,而且这个挤奶器上的奶牛不得超过m个),求c只奶牛中走的最远的奶牛的最小移动总距离。
思路:首先思考到这题目要二分答案,因为移动最远的奶牛的移动总距离越大,就越有可能满足要求
然后我们用c次dijkstra求出每个奶牛到所有挤奶器的最近距离,再根据我们二分的奶牛移动最远距离,限制了每个奶牛可以到达的挤奶器
随后我们就要判断这种情况下是否所有奶牛都可以到达一个挤奶器,而且每个挤奶器上的奶牛不超过m个,此时就是一个明显的求最大流问题了,建图如下:
1)由源点向每个奶牛建一条容量为1的边
2)由每个挤奶器向汇点建一条容量为m的边
3)由每个奶牛向其可以到达的所有挤奶器建一条容量为inf的边
最后,我们跑出这个图的最大流,如果最大流等于c,说明当前的情况是满足的,需要尝试更小的最远的距离;如果最大流不等于c,说明当前的最远距离太小了,以至于有奶牛无法到达挤奶器,需要尝试更大的最远距离。
代码区
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#include<string>
#include<fstream>
#include<vector>
#include<stack>
#include <map>
#include <iomanip> #define bug cout << "**********" << endl
#define show(x, y) cout<<"["<<x<<","<<y<<"] "
#define LOCAL = 1;
using namespace std;
typedef long long ll;
const int inf = 0x3f3f3f3f;
const ll mod = 1e9 + ;
const int Max = 1e4 + ;
const int Max2 = 3e2+ ; struct Edge //构建的残余网络
{
int to, flow, next;
} edge[Max << ]; struct Edge2 //原图结点之间的关系
{
int to, dis;
}; int k, c, m, s, t;
vector<Edge2> edge_raw[Max2];
int dist[Max2][Max2]; //奶牛和机器的最近距离
bool vis[Max2]; int head[Max2], tot;
int dis[Max2], cur[Max2]; void init()
{
memset(dist, inf, sizeof(dist));
for (int i = ; i < Max2; i++)
edge_raw[i].clear();
} void add(int u, int v, int flow)
{
edge[tot].to = v;
edge[tot].flow = flow;
edge[tot].next = head[u];
head[u] = tot++;
} void dijkstra(int start)
{
memset(vis, , sizeof(vis));
priority_queue<pair<int, int> > q;
dist[start][start] = ;
q.push(make_pair(, start));
while (!q.empty())
{
int u = q.top().second;
q.pop();
if (vis[u])
continue;
vis[u] = true;
for (vector<Edge2>::iterator it = edge_raw[u].begin(); it != edge_raw[u].end(); ++it)
{
if (!vis[it->to] && dist[start][it->to] > dist[start][u] + it->dis)
{
dist[start][it->to] = dist[start][u] + it->dis;
q.push(make_pair(-dist[start][it->to], it->to));
}
}
}
} bool bfs() //判断图是否连通
{
queue<int> q;
memset(dis, -, sizeof(dis));
dis[s] = ;
q.push(s);
while (!q.empty())
{
int u = q.front();
q.pop();
for (int i = head[u]; i != -; i = edge[i].next)
{
int v = edge[i].to;
if (dis[v] == - && edge[i].flow > ) //可以借助边i到达新的结点
{
dis[v] = dis[u] + ; //求顶点到源点的距离编号
q.push(v);
}
}
}
return dis[t] != -; //确认是否连通
} int dfs(int u, int flow_in)
{
if (u == t)
return flow_in;
int flow_out = ; //记录这一点实际流出的流量
for (int i = cur[u]; i != -; i = edge[i].next)
{
cur[u] = i;
int v = edge[i].to;
if (dis[v] == dis[u] + && edge[i].flow > )
{
int flow_part = dfs(v, min(flow_in, edge[i].flow));
if (flow_part == )
continue; //无法形成增广路
flow_in -= flow_part; //流出了一部分,剩余可分配流入就减少了
flow_out += flow_part; //记录这一点最大的流出 edge[i].flow -= flow_part;
edge[i ^ ].flow += flow_part; //减少增广路上边的容量,增加其反向边的容量
if (flow_in == )
break;
}
}
return flow_out;
} int Dinic(int ans)
{
int sum = ;
while (bfs())
{
for (int i = ; i <= ans; i++)
cur[i] = head[i];
sum += dfs(s, inf);
}
return sum;
} int main()
{
#ifdef LOCAL
//freopen("input.txt", "r", stdin);
//freopen("output.txt", "w", stdout);
#endif
while (scanf("%d%d%d", &k, &c, &m) != EOF)
{
init();
for (int i = ; i <= k + c; i++)
{
for (int j = , way; j <= k + c; j++)
{
scanf("%d", &way);
if (way)
{
Edge2 e;
e.to = j;
e.dis = way;
edge_raw[i].push_back(e);
e.to = i;
edge_raw[j].push_back(e); //POJ不是C14以上的语言,所以只能写的麻烦一些了
}
}
}
for (int i = k + ; i <= k + c; i++) //求每个奶牛和所有挤奶器的最短距离
{
dijkstra(i);
} int l = , r = Max, len = ; //二分答案
s = k + c + , t = ; //源点,汇点
while (l <= r)
{
memset(head, -, sizeof(head));
tot = ;
int mid = (l + r) >> ;
for (int i = k + ; i <= k + c; i++)
{
for (int j = ; j <= k; j++)
{
if (dist[i][j] <= mid) //根据限制的最远距离,选取奶牛可以到达的挤奶器建边
{
add(i, j, inf);
add(j, i, );
}
}
add(s, i, );
add(i, s, ); //由源点到每头奶牛的边容量均为1
}
for (int i = ; i <= k; i++)
{
add(i, t, m);
add(t, i, ); //由机器到汇点的边容量为机器最大接纳奶牛数
} int max_flow = Dinic(k + c + );
if (max_flow == c) //满足要求,进一步缩小最远距离
{
r = mid - ;
len = mid;
}
else
{
l = mid + ;
}
}
printf("%d\n", len);
}
return ;
}
POJ - 2112 Optimal Milking (dijkstra + 二分 + 最大流Dinic)的更多相关文章
- POJ 2112 Optimal Milking(二分+最大流)
http://poj.org/problem?id=2112 题意: 现在有K台挤奶器和C头奶牛,奶牛和挤奶器之间有距离,每台挤奶器每天最多为M头奶挤奶,现在要安排路程,使得C头奶牛所走的路程中的最大 ...
- POJ 2112 Optimal Milking (Floyd+二分+最大流)
[题意]有K台挤奶机,C头奶牛,在奶牛和机器间有一组长度不同的路,每台机器每天最多能为M头奶牛挤奶.现在要寻找一个方案,安排每头奶牛到某台机器挤奶,使得C头奶牛中走过的路径长度的和的最大值最小. 挺好 ...
- POJ 2112 Optimal Milking (二分+最短路径+网络流)
POJ 2112 Optimal Milking (二分+最短路径+网络流) Optimal Milking Time Limit: 2000MS Memory Limit: 30000K To ...
- POJ 2112 Optimal Milking (二分 + floyd + 网络流)
POJ 2112 Optimal Milking 链接:http://poj.org/problem?id=2112 题意:农场主John 将他的K(1≤K≤30)个挤奶器运到牧场,在那里有C(1≤C ...
- POJ 2112: Optimal Milking【二分,网络流】
题目大意:K台挤奶机,C个奶牛,每台挤奶器可以供M头牛使用,给出奶牛和和机器间的距离矩阵,求所有奶牛走最大距离的最小值 思路:最大距离的最小值,明显提示二分,将最小距离二分之后问题转化成为:K台挤奶机 ...
- Poj 2112 Optimal Milking (多重匹配+传递闭包+二分)
题目链接: Poj 2112 Optimal Milking 题目描述: 有k个挤奶机,c头牛,每台挤奶机每天最多可以给m头奶牛挤奶.挤奶机编号从1到k,奶牛编号从k+1到k+c,给出(k+c)*(k ...
- POJ 2112 Optimal Milking (二分 + 最大流)
题目大意: 在一个农场里面,有k个挤奶机,编号分别是 1..k,有c头奶牛,编号分别是k+1 .. k+c,每个挤奶机一天最让可以挤m头奶牛的奶,奶牛和挤奶机之间用邻接矩阵给出距离.求让所有奶牛都挤到 ...
- POJ 2112—— Optimal Milking——————【多重匹配、二分枚举答案、floyd预处理】
Optimal Milking Time Limit:2000MS Memory Limit:30000KB 64bit IO Format:%I64d & %I64u Sub ...
- POJ 2112 Optimal Milking (Dinic + Floyd + 二分)
Optimal Milking Time Limit: 2000MS Memory Limit: 30000K Total Submissions: 19456 Accepted: 6947 ...
随机推荐
- JavaWeb_(Hibernate框架)Hibernate中事务
Hibernate中事务 事务的性质 事物的隔离级别 配置事务的隔离级别 事务的性质 原子性:原子,不可再分,一个操作不能分为更小的操作,要么全都执行,要么全不执行. 一致性:事务在完成时,必须使得所 ...
- Jmeter -- 循环控制器和线程并发(关注执行顺序)
测试计划中包含两个线程组,分别设置如下: 线程组1:线程数为2,循环次数为4 线程组2:线程数为4 执行顺序 监听器View Results in Table(用表格查看结果)中,ThreadName ...
- 识别C++编译器编译标准
cout << __cplusplus<<endl; C++03:__cplusplus = 199711L C++11:__cplusplus = 201103L
- IP输出 之 分片ip_fragment、ip_do_fragment
概述 ip_fragment函数用于判断是否进行分片,在没有设置DF标记的情况下进入分片,如果设置了DF标记,则继续判断,如果不允许DF分片或者收到的最大分片大于MTU大小,则回复ICMP,释放skb ...
- js设计模式-代理模式
1.什么是设计模式? 设计模式:在软件设计过程中常用的代码规范,针对特定的场景 2.应用场景: 麦当劳点餐 观察者模式 规定的代码格式 花店送花 :代理模式 真实对象(男同学)-----代理对 ...
- 离线安装nuget包EPPlus
1先去https://www.nuget.org/packages/EPPlus/4.1.0下载,epplus.4.1.0.nupkg 2找到本地文件位置:H:\DOWNLOAD\ 3在vs的程序包管 ...
- python - paramiko模块 修改logging提示为静默
import paramiko import logging log = paramiko.util.get_logger("paramiko") log.setLevel(log ...
- 能否保证service不被杀死?
Service设置成START_STICKY kill 后会被重启(等待5秒左右),重传Intent,保持与重启前一样 提升service优先级 在AndroidManifest.xml文件中对于in ...
- Activity缓存方法
有a.b两个Activity,当从a进入b之后一段时间,可能系统会把a回收,这时候按back,执行的不是a的onRestart而是onCreate方法,a被重新创建一次,这是a中的临时数据和状态可能就 ...
- Java生成三位随机数
转: [转]Java生成三位随机数 public class Test2 { public static void main(String [] srgs) { int i=(int)(Math.ra ...