最大流,最小割
dinic
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#define ll long long
using namespace std;
const int Maxn=2e2+5,Maxm=(5e3+50)*2,INF=1e8;
struct edge{
int to,nxt;
ll f;
}e[Maxm];
int n,m,S,T,tot=1;
int head[Maxn],dep[Maxn],cur[Maxn];
inline void add(int u,int v,ll w)
{
e[++tot].to=v; e[tot].nxt=head[u]; head[u]=tot; e[tot].f=w;
e[++tot].to=u; e[tot].nxt=head[v]; head[v]=tot; e[tot].f=0;
}
inline bool bfs()
{
queue<int> q;
memset(dep,0,sizeof dep);
q.push(S);
dep[S]=1; cur[S]=head[S];
while(!q.empty())
{
int now=q.front();q.pop();
for(int i=head[now];i>1;i=e[i].nxt)
{
int nx=e[i].to;
if(dep[nx]||e[i].f==0) continue;
dep[nx]=dep[now]+1;
cur[nx]=head[nx];
if(nx==T) return 1;
q.push(nx);
}
}
return 0;
}
inline ll dfs(int x,ll limit)
{
if(x==T) return limit;
ll flow=0;
for(int i=cur[x];i>1&&flow<limit;i=e[i].nxt)
{
int nx=e[i].to;
if(dep[nx]==dep[x]+1&&e[i].f)
{
ll t=dfs(nx,min(limit-flow,e[i].f));
e[i].f-=t; e[i^1].f+=t;
flow+=t;
}
}
if(flow==0) dep[x]=0;
return flow;
}
inline ll dinic()
{
ll ret=0;
while(bfs())
{
ret+=dfs(S,INF);
}
return ret;
}
int main()
{
scanf("%d%d%d%d",&n,&m,&S,&T);
for(int i=1,a,b;i<=m;i++)
{
ll c;
scanf("%d%d%lld",&a,&b,&c);
add(a,b,c);
}
printf("%lld",dinic());
return 0;
}
费用流
EK
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
inline int read()
{
int x=0,f=1;
char ch=getchar();
while(ch<'0'||ch>'9'){ if(ch=='-') f=-1;ch=getchar(); }
while(ch>='0'&&ch<='9'){ x=(x<<1)+(x<<3)+(ch^48);ch=getchar(); }
return x*f;
}
const int Maxn=5e3+5,Maxm=5e4*2+500,INF=1e8;
struct edge{
int to,nxt,f,cost;
}e[Maxm];
queue<int> q;
int n,m,S,T,tot=1;
int head[Maxn],incf[Maxn],pre[Maxn],dis[Maxn];
bool vis[Maxn];
inline void add(int u,int v,int c,int d)
{
e[++tot].to=v; e[tot].nxt=head[u]; head[u]=tot; e[tot].f=c; e[tot].cost=d;
e[++tot].to=u; e[tot].nxt=head[v]; head[v]=tot; e[tot].f=0; e[tot].cost=-d;
}
inline bool spfa()
{
memset(dis,0x3f,sizeof dis);
memset(incf,0,sizeof incf);
q.push(S);
incf[S]=INF; dis[S]=0;
while(!q.empty())
{
int now=q.front();q.pop();
vis[now]=0;
for(int i=head[now];i>1;i=e[i].nxt)
{
int nx=e[i].to;
if(dis[nx]>dis[now]+e[i].cost&&e[i].f)
{
dis[nx]=dis[now]+e[i].cost;
incf[nx]=min(incf[now],e[i].f);
pre[nx]=i;
if(!vis[nx])
{
vis[nx]=1;
q.push(nx);
}
}
}
}
return incf[T]>0;
}
inline void EK(int& flow,int& cost)
{
flow=0; cost=0;
while(spfa())
{
int t=incf[T];
flow+=t; cost+=t*dis[T];
for(int i=T;i!=S;i=e[pre[i]^1].to)
{
e[pre[i]].f-=t;
e[pre[i]^1].f+=t;
}
}
}
int main()
{
n=read(); m=read(); S=read(); T=read();
for(int i=1,a,b,c,d;i<=m;i++)
{
a=read(); b=read(); c=read(); d=read();
add(a,b,c,d);
}
int flow,cost;
EK(flow,cost);
printf("%d %d",flow,cost);
return 0;
}
TIPS:
-
- 拆点之后的无向边加的时候,注意起点和终点的不同,不能直接加正向和反向流量一样的边
- 2.注意最小割点和最小割边的不同之处