链接:http://poj.org/problem?id=2112

题意:有k个挤奶器。编号1~k,c头牛,编号k+1~k+c,每一个挤奶器最多能给m头牛挤奶,给你一个k+c的邻接矩阵。要求每头牛都能挤奶而且要求c头牛须要走的全部路程中的最大路程最小,求这个最小的路。

思路:

1. 先用floyd处理出多源最短路

2. 用二分枚举答案的可能。初始上限应该为(200+30)*200。可是我这么开T了,可能由于代码太挫,改到1000,卡着时间过了,仅仅能说poj数据弱了。后来看别人的代码。和我的做法一样可是用了邻接表。就能设上限为40000了。在二分中:

(1)构造容量网络,以0点为源点。到每头牛的容量为1。以n+1点为汇点。每一个挤奶器到汇点的容量为m,当然反过来也能够,由于源点和汇点的流量是相等的(等于c)。对于每头牛和每一个挤奶器之间的距离,假设比枚举的距离还大。则容量为0,否则容量为1。

(2)Dinic找出网络最大流,非常明显最大流最大是c,当最大流是c的时候是一种答案。但不一定是最优。更新二分上限。假设最大流没达到c,则更新下限。

这是做完POJ2391知道了Dinic优化。改进后的写法。优化了三个地方:容量网络改为邻接表、Dinic优化、二分上限从floyd中返回,157MS ,我原以为Dinic优化应该是效率提高的主要原因。把Dinic优化去掉单纯用邻接表。266MS。原来邻接表才是这道题效率提高的主要原因,可是Dinic优化还是非常实用的。

假设floyd再加个剪枝,能够跑110MS

#include<cstring>
#include<string>
#include<fstream>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cctype>
#include<algorithm>
#include<queue>
#include<map>
#include<set>
#include<vector>
#include<stack>
#include<ctime>
#include<cstdlib>
#include<functional>
#include<cmath>
using namespace std;
#define PI acos(-1.0)
#define MAXN 200100
#define eps 1e-7
#define INF 0x7FFFFFFF
#define LLINF 0x7FFFFFFFFFFFFFFF
#define seed 131
#define mod 1000000007
#define ll long long
#define ull unsigned ll
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1 struct node{
int u,w,next;
}edge[MAXN];
int head[350],e[350][350];
int dist[350];
int cnt,n,m,k,c,src,sink;
void add_edge(int a,int b,int c){
edge[cnt].u = b;
edge[cnt].w = c;
edge[cnt].next = head[a];
head[a] = cnt++;
}
int floyd(){
int i,j,k;
int maxm = 0;
for(k=1;k<=n;k++){
for(i=1;i<=n;i++){
if(e[i][k]==INF) continue;
for(j=1;j<=n;j++){
if(e[k][j]!=INF&&e[i][k]+e[k][j]<e[i][j]){
e[i][j] = e[i][k] + e[k][j];
if(e[i][j]>maxm) maxm = e[i][j];
}
}
}
}
return maxm;
}
void build_graph(int minm){
int i,j;
memset(head,-1,sizeof(head));
cnt = 0;
for(i=1;i<=k;i++){
add_edge(i,n+1,m);
add_edge(n+1,i,0);
}
for(i=k+1;i<=n;i++){
add_edge(0,i,1);
add_edge(i,0,0);
}
for(i=k+1;i<=n;i++){
for(j=1;j<=k;j++){
if(e[i][j]<=minm){
add_edge(i,j,1);
add_edge(j,i,0);
}
}
}
}
int bfs(){
int i,j;
memset(dist,-1,sizeof(dist));
queue<int>q;
q.push(0);
dist[0] = 1;
while(!q.empty()){
int t = q.front();
q.pop();
for(i=head[t];i!=-1;i=edge[i].next){
if(dist[edge[i].u]==-1&&edge[i].w){
dist[edge[i].u] = dist[t] + 1;
q.push(edge[i].u);
}
}
}
if(dist[n+1]!=-1) return 1;
else return 0;
}
int dfs(int u,int delta){
int i,j;
int dd;
if(u==n+1) return delta;
for(i=head[u];i!=-1;i=edge[i].next){
if(dist[edge[i].u]==dist[u]+1&&edge[i].w&&(dd = dfs(edge[i].u,min(edge[i].w,delta)))){
edge[i].w -= dd;
edge[i^1].w += dd;
return dd;
}
}
dist[u] = -1;
return 0;
}
int main(){
int i,j;
while(scanf("%d%d%d",&k,&c,&m)!=EOF){
n = k + c;
for(i=1;i<=n;i++){
for(j=1;j<=n;j++){
scanf("%d",&e[i][j]);
if(e[i][j]==0) e[i][j] = INF;
}
}
int mid, l = 0,r=floyd();
int sum,temp;
while(l<r){
mid = (l+r)/2;
sum = 0;
build_graph(mid);
while(bfs()){
while(1){
temp = dfs(0,INF);
if(!temp) break;
sum += temp;
}
}
if(sum>=c) r = mid;
else l = mid + 1;
}
printf("%d\n",l);
}
return 0;
}

之前的写法,1875MS擦边过了。

#include<cstring>
#include<string>
#include<fstream>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cctype>
#include<algorithm>
#include<queue>
#include<map>
#include<set>
#include<vector>
#include<stack>
#include<ctime>
#include<cstdlib>
#include<functional>
#include<cmath>
using namespace std;
#define PI acos(-1.0)
#define MAXN 50100
#define eps 1e-7
#define INF 0x7FFFFFFF
#define seed 131
#define mod 1000000007
#define ll long long
#define ull unsigned ll
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1 int edge[300][300],customer[300][300];
int vis[300],dist[300][300];
int n,m,k,c;
void floyd(){
int i,j,k;
for(k=1;k<=n;k++){
for(i=1;i<=n;i++){
for(j=1;j<=n;j++){
if(edge[i][k]!=INF&&edge[k][j]!=INF&&edge[i][k]+edge[k][j]<edge[i][j])
edge[i][j] = edge[i][k] + edge[k][j];
}
}
}
}
void build_graph(int minm){
int i,j;
memset(customer,0,sizeof(customer));
for(i=1;i<=k;i++) customer[i][n+1] = m;
for(i=k+1;i<=n;i++) customer[0][i] = 1;
for(i=k+1;i<=n;i++){
for(j=1;j<=k;j++){
if(edge[i][j]<=minm) customer[i][j] = 1;
}
}
}
int bfs(){
int i,j;
memset(vis,0,sizeof(vis));
memset(dist,0,sizeof(dist));
queue<int>q;
q.push(0);
vis[0] = 1;
while(!q.empty()){
int t = q.front();
q.pop();
for(i=0;i<=n+1;i++){
if(!vis[i]&&customer[t][i]){
vis[i] = 1;
dist[t][i] = 1;
q.push(i);
}
}
}
if(vis[n+1]) return 1;
else return 0;
}
int dfs(int u,int delta){
int i,j,s;
if(u==n+1) return delta;
s = delta;
for(i=0;i<=n+1;i++){
if(dist[u][i]){
int dd = dfs(i,min(customer[u][i],delta));
customer[u][i] -= dd;
customer[i][u] += dd;
delta -= dd;
}
}
return s - delta;
}
int main(){
int i,j;
while(scanf("%d%d%d",&k,&c,&m)!=EOF){
n = k + c;
for(i=1;i<=n;i++){
for(j=1;j<=n;j++){
scanf("%d",&edge[i][j]);
if(edge[i][j]==0) edge[i][j] = INF;
}
}
floyd();
int mid, l = 0, r = 10000;
int sum;
while(l<r){
mid = (l+r)/2;
sum = 0;
build_graph(mid);
while(bfs()) sum += dfs(0,INF);
if(sum==c) r = mid;
else l = mid + 1;
}
printf("%d\n",l);
}
return 0;
}

POJ--2112--Optimal Milking【Floyd+Dinic+二分答案】的更多相关文章

  1. POJ 2112 Optimal Milking (Dinic + Floyd + 二分)

    Optimal Milking Time Limit: 2000MS   Memory Limit: 30000K Total Submissions: 19456   Accepted: 6947 ...

  2. POJ 2112 Optimal Milking (Floyd+二分+最大流)

    [题意]有K台挤奶机,C头奶牛,在奶牛和机器间有一组长度不同的路,每台机器每天最多能为M头奶牛挤奶.现在要寻找一个方案,安排每头奶牛到某台机器挤奶,使得C头奶牛中走过的路径长度的和的最大值最小. 挺好 ...

  3. POJ 2112 Optimal Milking【网络流+二分+最短路】

    求使所有牛都可以被挤牛奶的条件下牛走的最长距离. Floyd求出两两节点之间的最短路,然后二分距离. 构图: 将每一个milking machine与源点连接,边权为最大值m,每个cow与汇点连接,边 ...

  4. POJ 2112 Optimal Milking 最短路 二分构图 网络流

    题意:有C头奶牛,K个挤奶站,每个挤奶器最多服务M头奶牛,奶牛和奶牛.奶牛和挤奶站.挤奶站和挤奶站之间都存在一定的距离.现在问满足所有的奶牛都能够被挤奶器服务到的情况下,行走距离的最远的奶牛的至少要走 ...

  5. POJ 2112 Optimal Milking(最大流+二分)

    题目链接 测试dinic模版,不知道这个模版到底对不对,那个题用这份dinic就是过不了.加上优化就WA,不加优化TLE. #include <cstdio> #include <s ...

  6. POJ 2112 Optimal Milking (二分 + floyd + 网络流)

    POJ 2112 Optimal Milking 链接:http://poj.org/problem?id=2112 题意:农场主John 将他的K(1≤K≤30)个挤奶器运到牧场,在那里有C(1≤C ...

  7. Poj 2112 Optimal Milking (多重匹配+传递闭包+二分)

    题目链接: Poj 2112 Optimal Milking 题目描述: 有k个挤奶机,c头牛,每台挤奶机每天最多可以给m头奶牛挤奶.挤奶机编号从1到k,奶牛编号从k+1到k+c,给出(k+c)*(k ...

  8. POJ 2112 Optimal Milking (二分+最短路径+网络流)

    POJ  2112 Optimal Milking (二分+最短路径+网络流) Optimal Milking Time Limit: 2000MS   Memory Limit: 30000K To ...

  9. POJ 2112—— Optimal Milking——————【多重匹配、二分枚举答案、floyd预处理】

    Optimal Milking Time Limit:2000MS     Memory Limit:30000KB     64bit IO Format:%I64d & %I64u Sub ...

  10. POJ 2112 Optimal Milking ( 经典最大流 && Floyd && 二分 )

    题意 : 有 K 台挤奶机器,每台机器可以接受 M 头牛进行挤奶作业,总共有 C 头奶牛,机器编号为 1~K,奶牛编号为 K+1 ~ K+C ,然后给出奶牛和机器之间的距离矩阵,要求求出使得每头牛都能 ...

随机推荐

  1. Windows下Go语言 幽灵蛛的配置

    这里是环境变量 目录结构 在这其中,我主要需要观察src

  2. class com.sun.jersey.core.impl.provider.entity.XMLJAXBElementProvider$Text

    运行mapreduce遇到的错: Java.lang.ClassCastException: classcom.sun.jersey.core.impl.provider.entity.XMLJAXB ...

  3. Test Doubles - Fakes, Mocks and Stubs.

    https://dev.to/milipski/test-doubles---fakes-mocks-and-stubs This text was originally posted at Prag ...

  4. Constructing Roads --hdoj

    Constructing Roads Time Limit : 2000/1000ms (Java/Other)   Memory Limit : 65536/32768K (Java/Other) ...

  5. if,elif,else的关系 input print int的用法

    qian=input("找劳保网是什么网站?:")if qian=="zhaolaobaowang.com": print("正确")els ...

  6. 同一sql程序执行比数据库执行慢

    最近项目发现同一个sql在java端执行比在数据库执行慢很多,原因可能是程序的sql参数类型与数据库字段的类型不一致.

  7. WPF使用Winform PDFView控件

    最近开发wpf项目中有一个模块需要显示PDF文件内容.由于WPF本身没有PDF加载控件(似乎有收费的我查到过类似的资料.如果有新的pdf控件也请通知我一下谢谢). 项目使用之前也是从网上获取的资料,因 ...

  8. shopping car 1.0

    #!/usr/bin/env python# -*- coding: utf-8 -*-# @File : 20180510001.py# @Author: Anthony.waa# @Date : ...

  9. 高级I/O函数

    给套接口上的I/O设置超时 1.调用alarm,在调用超过指定时间时产生SIGALARM信号,这涉及到信号处理,而且可能和进程中其他的alarm冲突 2.使用select阻塞在等待I/O上,selec ...

  10. React 学习笔记:1-react 入门

    接下来的项目里有用到react,最近一段时间主要关注于react 的学习.大部门都是网上的资料,学习整理并记录,加深记忆. React 是Facebook推出的用来构建用户界面的JavaScript库 ...