POJ 2391 Ombrophobic Bovines 网络流 建模
【题目大意】
给定一个无向图,点i处有Ai头牛,点i处的牛棚能容纳Bi头牛,求一个最短时间T使得在T时间内所有的牛都能进到某一牛棚里去。(1 <= N <= 200, 1 <= M <= 1500, 0 <= Ai <= 1000, 0 <= Bi <= 1000, 1 <= Dij <= 1,000,000,000)
一开始想拆点建图,0到x集合为汇,值为各个区域的牛数量, Y到终点连边,值为各个区域的容量,然后就是看怎么连x和y了
我一开始把可以连接的X和Y连起来,把可以互达的点在Y集合点那里连边,这样很麻烦,先跑一遍floyd把点到点的最短路求出来,然后直接X和Y集合可达即相连就行
二分结果,再建图,把在mid以内的X点对Y点连起来,跑最大流 判断结果即可
注意要用long long
一开始还没看清题意,一条路上可以同时走无数的牛,我一开始以为只能走一头,还敲MCMF去了。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
#define LL long long
#define INF 1LL<<60
using namespace std;
int f,p;
const int maxn=500;
struct Edge
{
int from,to,cap,flow;
};
struct Dinic
{
vector<Edge>edges;
vector<int> G[maxn];
int vis[maxn];
int cur[maxn];
int d[maxn];
void init(int n)
{
edges.clear();
for (int i=0; i<=n; i++)
{
G[i].clear();
}
}
void addedge(int from,int to,int cap)
{
int m;
edges.push_back((Edge)
{
from,to,cap,0
});
edges.push_back((Edge)
{
to,from,0,0
});
m=edges.size();
G[from].push_back(m-2);
G[to].push_back(m-1);
}
bool bfs(int s,int t)
{
memset(vis,0,sizeof vis);
queue<int> q;
q.push(s);
d[s]=0;
vis[s]=1;
while (!q.empty())
{
int u=q.front();
q.pop();
for (int i=0; i<G[u].size(); i++)
{
Edge& e=edges[G[u][i]];
if (!vis[e.to] && e.cap>e.flow)
{
vis[e.to]=1;
d[e.to]=d[u]+1;
q.push(e.to);
}
}
}
return vis[t];
}
int dfs(int x,int a,int t)
{
if (x==t || a==0) return a;
int flow=0,f;
for (int& i=cur[x]; i<G[x].size(); i++)
{
Edge& e=edges[G[x][i]];
if (d[x]+1==d[e.to] && (f=dfs(e.to,min(a,e.cap-e.flow),t))>0)
{
e.flow+=f;
edges[G[x][i]^1].flow-=f;
flow+=f;
a-=f;
if (a==0) break;
}
}
return flow;
}
int maxflow(int s,int t)
{
int flow=0;
while (bfs(s,t))
{ memset(cur,0,sizeof cur);
flow+=dfs(s,100000000,t);
}
return flow;
}
} mcmf;
int A[210],B[210];
LL path[210][210];
LL N;
void floyd()
{
for (int i=1; i<=f; i++)
{
for (int j=1; j<=f; j++)
{
for (int k=1; k<=f; k++)
{
if (j==k) continue;
path[j][k]=min(path[j][k],path[j][i]+path[i][k]);
N=max(N,path[j][k]);
}
}
}
}
int main()
{
//freopen("POJ_2391.in","r",stdin);
int a,b;
LL c;
while (scanf("%d%d",&f,&p)!=EOF)
{
int cur=0;
for (int i=1; i<=f; i++)
{
scanf("%d%d",&A[i],&B[i]);
cur+=A[i];
for(int j=1; j<=f; j++) path[i][j]=INF;
}
for (int i=1; i<=p; i++)
{
scanf("%d%d%lld",&a,&b,&c);
path[a][b]=min(path[a][b],c);
path[b][a]=min(path[b][a],c);
}
N=0;
floyd();
LL l,r,mid;
l=1,r=N;
//cout<<l<<" "<<r<<endl;
LL ans=-1;
while(l<r)
{
mcmf.init(2*f+10);
for (int i=1; i<=f; i++)
{
mcmf.addedge(0,i,A[i]);
mcmf.addedge(i,i+f,1<<30);
}
for (int i=1; i<=f; i++)
{
mcmf.addedge(i+f,2*f+5,B[i]);
}
mid=(r+l)>>1;
for (int i=1;i<=f;i++){
for (int j=1;j<=f;j++){
if (path[i][j]>mid || i==j) continue;
mcmf.addedge(i,f+j,1<<30);
}
}
int res=mcmf.maxflow(0,2*f+5);
if (res>=cur){
//cout<<res<<" "<<cur<<endl;
//cout<<mid<<endl;
ans=mid;
r=mid;
}
else{
l=mid+1;
} }
printf("%lld\n",ans);
}
return 0;
}
POJ 2391 Ombrophobic Bovines 网络流 建模的更多相关文章
- poj 2391 Ombrophobic Bovines, 最大流, 拆点, 二分, dinic, isap
poj 2391 Ombrophobic Bovines, 最大流, 拆点, 二分 dinic /* * Author: yew1eb * Created Time: 2014年10月31日 星期五 ...
- POJ 2391 Ombrophobic Bovines (Floyd + Dinic +二分)
Ombrophobic Bovines Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 11651 Accepted: 2 ...
- POJ 2391 Ombrophobic Bovines
Ombrophobic Bovines Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 18623 Accepted: 4 ...
- poj 2391 Ombrophobic Bovines(最大流+floyd+二分)
Ombrophobic Bovines Time Limit: 1000MSMemory Limit: 65536K Total Submissions: 14519Accepted: 3170 De ...
- POJ 2391 Ombrophobic Bovines【二分 网络流】
题目大意:F个草场,P条道路(无向),每个草场初始有几头牛,还有庇护所,庇护所有个容量,每条道路走完都有时间,问所有奶牛都到庇护所最大时间最小是多少? 思路:和POJ2112一样的思路,二分以后构建网 ...
- POJ 2391 Ombrophobic Bovines(二分+拆点+最大流)
http://poj.org/problem?id=2391 题意: 给定一个无向图,点i处有Ai头牛,点i处的牛棚能容纳Bi头牛,求一个最短时间T,使得在T时间内所有的牛都能进到某一牛棚里去. 思路 ...
- POJ 2391 Ombrophobic Bovines ( 经典最大流 && Floyd && 二分 && 拆点建图)
题意 : 给出一些牛棚,每个牛棚都原本都有一些牛但是每个牛棚可以容纳的牛都是有限的,现在给出一些路与路的花费和牛棚拥有的牛和可以容纳牛的数量,要求最短能在多少时间内使得每头牛都有安身的牛棚.( 这里注 ...
- POJ 2391 Ombrophobic Bovines ★(Floyd+二分+拆点+最大流)
[题意]有n块草地,一些奶牛在草地上吃草,草地间有m条路,一些草地上有避雨点,每个避雨点能容纳的奶牛是有限的,给出通过每条路的时间,问最少需要多少时间能让所有奶牛进入一个避雨点. 和POJ2112很类 ...
- POJ 2391.Ombrophobic Bovines (最大流)
实际上是求最短的避雨时间. 首先将每个点拆成两个,一个连接源点,一个连接汇点,连接源点的点的容量为当前单的奶牛数,连接汇点的点为能容纳的奶牛数. floyd求任意两点互相到达的最短时间,二分最长时间, ...
随机推荐
- ISR high memory参数
1.通过 show process memory 获取的数据参数解释: 来自 <http://blog.router-switch.com/2013/12/show-processes-memo ...
- kafka 分区
1. 默认存储 /tmp/kafka-logs test-0 0:分区数
- 19市赛 树状数组 第k个小的糖果
int find_kth(int k) { , cnt = , i; ; i >= ; i--)/ { ans += ( << i); if (ans >= maxn|| cn ...
- iOS 根据域名查询 IP 地址
在 iOS 开发中,如果需要知道网站的 IP 地址: #include <netdb.h> #include <arpa/inet.h> NSString *webSiteSt ...
- Educational Codeforces Round 69 (Rated for Div. 2)D(DP,思维)
#include<bits/stdc++.h>using namespace std;int a[300007];long long sum[300007],tmp[300007],mx[ ...
- 自定义更改 Xcode 新建 .h/.m 文件头部注释说明(文件名、工程名、作者、公司、版权等)信息
使用 Xcode 新建工程文件时,或默认生成一套注释说明信息在 .h/.m 文件的头部,说明了创建这个文件的名称.工程名.日期.作者.公司.版权等信息 // // ___FILENAME___ // ...
- [原]JointJS流程图
最近项目上需要用流程图来做问题定界分析,之前有同事用jsPlumb做过,但是阅读代码后觉得比较麻烦,所以自己又找了一圈,找到一个叫Dagre-D3的开源类库,画出来的效果如下图,Dagre-D3最大的 ...
- Maven项目-jar包冲突:MyServlet is not a Servlet解决方法
异常显示:运行时jar包冲突异常 问题所在: jar包冲突,自己导入的jar包与tomcat中的jar包相同. 解决方法: 对(servlet,jsp)jar包设置作用域. 再次启动项目,问题已解决.
- 图片的onload事件与better-scroll结合[ 当fastclick插件和better-scroll发生冲突导致点击事件失效时,可以给需要点击的元素加一个class="needsclick"]
<!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...
- gets和scanf区别
scanf 和 gets 读取字符串 深入了解scanf()/getchar()和gets()等函数 scanf与gets函数读取字符串的区别 今天看到一段话,大致是说gets比scanf()快,有点 ...