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 ...
随机推荐
- python中reversed()函数的用法
reversed(a) # 返回一个反转的迭代器 举例 a = [1, 2, 3, 4, 5] b = reversed(a) # b是一个迭代器 # print(list(b)) = [5, 4, ...
- ArrayList,LinkedList和String
import java.util.ArrayList; public class Demo{ public static void main(String[] args) throws Excepti ...
- java基础阶段几个必会面试题
摘自:https://www.cnblogs.com/zn19961006/p/11869182.html java基础阶段几个必会面试题 目录 1.说出你对面向对象的理解 在我理解,面向对象是向现实 ...
- react 脚手架 及路由和 redux
前提是我们需要下载 nodejs 使用 npm 下载 react 的脚手架,react-router-dom,redux 全局下载 react 的脚手架:npm i create-react-app ...
- vuex 讲解
vuex 状态的管理状态,它采用集中式存储管理应用的所有组件的状态,尤其是在中大型项目,则是很好的开发利器 vuex 的流程图 vuex 的优势: 1. vuex 的存储状态,响应式的 2. 他是所有 ...
- pandas.Series.value_counts
pandas.Series.value_counts Series.value_counts(normalize=False, sort=True, ascending=False, bins=Non ...
- el-mement表单校验-校验失败时自动聚焦到失败的input框
思路:调用input的focus()事件,聚焦到失败的input框上, 问题:如何获取失败的input框,结合element的校验方法,查看可以提供入口的地方 调用this.$refs[formNam ...
- ArrayList && HashMap扩容策略
ArrayList扩容策略:默认10 扩容时是base + base/2, 即10 15 22 33 49...扩容时不安全:grow方法扩容时,赋值 elementData = Arrays.cop ...
- ubuntu最近升级到最新的linux内核后,网络无法使用怎么办?
答:进入旧的内核中编译需要的网卡模块 1. 启动旧的内核进入系统 2. 安装新内核源码 3. 找出当前的网卡型号 4. 尝试卸载某个与网卡相关的内核模块,观察是否影响当前网卡的使用,如果有影响,那么便 ...
- Spring Bean学习创建及使用<一>
转自:http://glzaction.iteye.com/blog/1299441 Spring中依赖注入有三种注入方式: 一.构造器注入: 二.设值注入(setter方式注入): 三.Feild方 ...