洛谷P3376【模板】網絡最大流 Dinic模板
阿新 • • 發佈:2017-09-17
span -c -s blog name sca print 技術 pop
之前的Dinic模板照著劉汝佳寫的vector然後十分鬼畜跑得奇慢無比,雖然別人這樣寫也沒慢多少但是自己的就是令人捉急。
改成鄰接表之後快了三倍,雖然還是比較慢但是自己比較滿意了。雖然一開始ecnt從0開始WA了一發。。。
之前的碼風也十分鬼畜呀縮進只縮1、2格不懂自己怎麽想的。。
反正今天就安心劃劃水。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<queue>
#include <vector>
typedef long long LL;
const int maxn=20050,maxm=200050,INF=0x7f7f7f7f;
using namespace std;
int n,m,s,t,u,v,w,ecnt=1,fir[maxn],dis[maxn],cur[maxn],ans;
struct edge {
int from,to,cap,flow,nxt;
edge(){}
edge(int from,int to,int cap,int flow,int nxt):from(from),to(to),cap(cap),flow(flow),nxt(nxt){}
}e[maxm];
void add(int u,int v,int w) {
e[++ecnt]=edge(u,v,w,0,fir[u]);
e[++ecnt]=edge(v,u,0,0,fir[v]);
fir[u]=ecnt-1; fir[v]=ecnt;
}
void init() {
scanf("%d%d%d%d",&n,&m,&s,&t);
for(int i=1;i<=m;i++) {
scanf("%d%d%d",&u,&v,&w);
add(u,v,w);
}
}
queue <int>que;
int bfs(int s,int t) {
memset(dis,0,sizeof(dis));
dis[s]=1;
que.push(s);
while(!que.empty()) {
int x=que.front() ;que.pop();
for(int i=fir[x];i;i=e[i].nxt)
if(!dis[e[i].to]&&e[i].flow<e[i].cap) {
dis[e[i].to]=dis[x]+1;
que.push(e[i].to);
}
}
return dis[t];
}
int dfs(int x,int a) {
if(x==t||a==0) return a;
int fl,f=0;
for(int &i=cur[x];i;i=e[i].nxt) {
if(dis[e[i].to]==dis[x]+1&&(fl=dfs(e[i].to,min(a,e[i].cap-e[i].flow)))){
e[i].flow+=fl;
e[i^1].flow-=fl;
a-=fl;
f+=fl;
}
if(a==0) break;
}
return f;
}
int Dinic(int s,int t) {
int res=0;
while(bfs(s,t)) {
for(int i=1;i<=n;i++) cur[i]=fir[i];
res+=dfs(s,INF);
}
return res;
}
void work() {
ans=Dinic(s,t);
printf("%d\n",ans);
}
int main()
{
init();
work();
return 0;
}
View Code
洛谷P3376【模板】網絡最大流 Dinic模板