Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 625 Accepted Submission(s): 367
#include<iostream> #include<cstdio> #include<cstring> #include <algorithm> #include <math.h> #include <queue> using namespace std; const int N = 1005; const int INF = 999999999; struct Edge{ int v,w,next; }edge[N*N]; int head[N]; int level[N]; int tot,max_increase; void init() { memset(head,-1,sizeof(head)); tot=0; } void addEdge(int u,int v,int w,int &k) { edge[k].v = v,edge[k].w=w,edge[k].next=head[u],head[u]=k++; edge[k].v = u,edge[k].w=0,edge[k].next=head[v],head[v]=k++; } int BFS(int src,int des) { queue<int>q; memset(level,0,sizeof(level)); level[src]=1; q.push(src); while(!q.empty()) { int u = q.front(); q.pop(); if(u==des) return 1; for(int k = head[u]; k!=-1; k=edge[k].next) { int v = edge[k].v; int w = edge[k].w; if(level[v]==0&&w!=0) { level[v]=level[u]+1; q.push(v); } } } return -1; } int dfs(int u,int des,int increaseRoad){ if(u==des||increaseRoad==0) { return increaseRoad; } int ret=0; for(int k=head[u];k!=-1;k=edge[k].next){ int v = edge[k].v,w=edge[k].w; if(level[v]==level[u]+1&&w!=0){ int MIN = min(increaseRoad-ret,w); w = dfs(v,des,MIN); if(w > 0) { edge[k].w -=w; edge[k^1].w+=w; ret+=w; if(ret==increaseRoad){ return ret; } } else level[v] = -1; if(increaseRoad==0) break; } } if(ret==0) level[u]=-1; return ret; } int Dinic(int src,int des) { int ans = 0; while(BFS(src,des)!=-1) ans+=dfs(src,des,INF); return ans; } int d,n,m,src,des; bool vis[N]; void dfs1(int u,int ans){ if(u==des){ max_increase = max(max_increase,ans); return ; } for(int k=head[u];k!=-1;k=edge[k].next){ int v = edge[k].v,w = edge[k].w; if(!vis[v]){ vis[v] = true; dfs1(v,min(w,ans)); ///最大流量由最小容量边决定 vis[v] = false; } } } int main() { int tcase; scanf("%d",&tcase); while(tcase--){ init(); max_increase = -1; scanf("%d%d%d%d%d",&d,&n,&m,&src,&des); for(int i=1;i<=m;i++){ int u,v,w; scanf("%d%d%d",&u,&v,&w); addEdge(u,v,w,tot); } memset(vis,false,sizeof(vis)); dfs1(src,99999999); ///deal int max_flow = Dinic(src,des); printf("%d %.3lf\n",d,max_flow*1.0/max_increase); } return 0; }
原文:http://www.cnblogs.com/liyinggang/p/5723249.html