POJ 2455Secret Milking Machine(二分+网络流之最大流)
生活随笔
收集整理的這篇文章主要介紹了
POJ 2455Secret Milking Machine(二分+网络流之最大流)
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
題目地址:POJ2455
手殘真浪費時間啊。。又拖到了今天才找出了錯誤。。每晚兩道題不知不覺又變回了每晚一道題、。。sad。。
第一次在isap中忘記調用bfs,第二次則是遍歷的時候竟然是從1開始遍歷的。。。sad。。。
這題思路倒是非常easy,就是有一個比較坑的地方,就是這里的重邊要當兩條邊來用,曾經受最短路什么的影響,直接把慢的刪了,非常明顯不正確。。。這里的兩條重邊都是能夠走的。
建圖思路是,將1當作源點,將n當作匯點。右邊的地方就連邊,注意是無向邊。最后推斷最大流是否等于道路條數。二分找最小值。
代碼例如以下:
#include <iostream> #include <stdio.h> #include <string.h> #include <stdlib.h> #include <math.h> #include <ctype.h> #include <queue> #include <map> #include <algorithm>using namespace std; const int INF=1e9; int head[500], s, t, nv, cnt, sum; int num[500], d[500], pre[500], cur[500], q[500]; struct node {int u, v, cap, next; }edge[5000000]; struct N {int u, v, w; }bian[5000000]; void add(int u, int v, int cap) {edge[cnt].v=v;edge[cnt].cap=cap;edge[cnt].next=head[u];head[u]=cnt++;edge[cnt].v=u;edge[cnt].cap=cap;edge[cnt].next=head[v];head[v]=cnt++; } void bfs() {memset(num,0,sizeof(num));memset(d,-1,sizeof(d));int f1=0, f2=0, i;q[f1++]=t;d[t]=0;num[0]=1;while(f1>=f2){int u=q[f2++];for(i=head[u];i!=-1;i=edge[i].next){int v=edge[i].v;if(d[v]==-1){d[v]=d[u]+1;num[d[v]]++;q[f1++]=v;}}} } int isap() {memcpy(cur,head,sizeof(cur));int flow=0, u=pre[s]=s, i;bfs();while(d[s]<nv){if(u==t){int f=INF, pos;for(i=s;i!=t;i=edge[cur[i]].v){if(f>edge[cur[i]].cap){f=edge[cur[i]].cap;pos=i;}}for(i=s;i!=t;i=edge[cur[i]].v){edge[cur[i]].cap-=f;edge[cur[i]^1].cap+=f;}flow+=f;if(flow>=sum){return flow;}u=pos;}for(i=cur[u];i!=-1;i=edge[i].next){if(d[edge[i].v]+1==d[u]&&edge[i].cap){break;}}if(i!=-1){cur[u]=i;pre[edge[i].v]=u;u=edge[i].v;}else{if(--num[d[u]]==0) break;int mind=nv;for(i=head[u];i!=-1;i=edge[i].next){if(mind>d[edge[i].v]&&edge[i].cap){mind=d[edge[i].v];cur[u]=i;}}d[u]=mind+1;num[d[u]]++;u=pre[u];}}return flow; } int main() {int n, p, a, b, c, i, j, top;while(scanf("%d%d%d",&n,&p,&sum)!=EOF){top=0;while(p--){scanf("%d%d%d",&a,&b,&c);bian[top].u=a;bian[top].v=b;bian[top++].w=c;}int high=1000001, mid, low=1, ans, x;while(low<=high){mid=(high+low)/2;s=1;t=n;nv=t+1;memset(head,-1,sizeof(head));cnt=0;for(i=0;i<top;i++){if(bian[i].w<=mid){add(bian[i].u,bian[i].v,1);}}x=isap();//printf("%d %d %d\n",x, mid,top);if(x>=sum){ans=mid;high=mid-1;}else{low=mid+1;}}printf("%d\n",ans);}return 0; }
與50位技術專家面對面20年技術見證,附贈技術全景圖
總結
以上是生活随笔為你收集整理的POJ 2455Secret Milking Machine(二分+网络流之最大流)的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: list根据对象进行排序
- 下一篇: 三个打印函数printf()/sprin