1. 程式人生 > >最大流Edmonds-Karp模板

最大流Edmonds-Karp模板

const int maxn=1e4+5;
struct Edge{
    int from,to,cap,flow;
    Edge(int u,int v,int c,int f):from(u),to(v),cap(c),flow(f){} 
};
struct EdmodsKarp{
    int n,m;
    vector<Edge> edges; 
    vector<int>G[maxn];//鄰接表,G[i][j]表示節點i的第j條邊在edges陣列中的需要
    int a[maxn];//當起點到i的可改進量
    int p[maxn];//最短路樹上p的入弧編號
void init(int n){ for(int i=0;i<n;i++) G[i].clear(); edges.clear(); } void AddEdge(int from,int to,int cap){ 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
); } int Maxflow(int s,int t){ int flow=0; for(;;){ //開始尋找增廣路 memset(a,0,sizeof(a)); queue<int> Q; Q.push(s); a[s]=inf; while(!Q.empty()){ int x=Q.front();Q.pop(); for(int
i=0;i<G[x].size();i++){ Edge &e=edges[G[x][i]]; if(!a[e.to]&&e.cap>e.flow){ p[e.to]=G[x][i]; a[e.to]=min(a[x],e.cap-e.flow); Q.push(e.to); } } if(a[t]) break;//找到了一條增廣路 } if(!a[t]) break;//不存在增廣路,當前流就是最大流 for(int u=t;u!=s;u=edges[p[u]].from){ edges[p[u]].flow+=a[t]; edges[p[u]^1].flow-=a[t]; } flow+=a[t]; } return flow; } };