masonpop @ 2023-05-23 17:09:08
感觉很离谱,自己照着网上一篇博客学的 dinic,感觉就算最坏
#include <bits/stdc++.h>
using namespace std;
#define Mem(a) memset((a),0,sizeof(a))
#define int long long
const int inf=1e15;
const int maxn=4e4+10;
int n,m,s,t,head[maxn],to[maxn],nxt[maxn],ver[maxn],tot;
int dep[maxn];
//flow:到当前点1能增广的最大流量
//dep:深度,用于分层
int ans;
queue<int> q;
inline void add(int x,int y,int z)
{
to[++tot]=y;
nxt[tot]=head[x];
head[x]=tot;
ver[tot]=z;
}
inline int rev(int x)//反边
{
if(x&1)return x+1;
return x-1;
}
int flag;
inline void bfs()
{
memset(dep,0,sizeof(dep));
while(!q.empty())q.pop();//清空
q.push(s);dep[s]=1;
while(!q.empty())
{
int x=q.front();q.pop();
for(int i=head[x];i;i=nxt[i])
{
int y=to[i];
if(dep[y] || !ver[i])continue;//访问过或者枯竭了
dep[y]=dep[x]+1;
q.push(y);
}
}
}
inline int dfs(int x,int flow)
{
if(x==t)
{
ans+=flow;
flag=1;
return flow;
}
for(int i=head[x];i;i=nxt[i])
{
int y=to[i];
if(dep[y]!=dep[x]+1 || !ver[i])continue;//层数不对或枯竭了
int delt=dfs(y,min(flow,ver[i]));
if(flag)
{
ver[i]-=delt;
ver[rev(i)]+=delt;
return delt;
}
}
return 0;
}
inline void dinic()
{
while(1)
{
bfs();
if(!dep[t])return;//到不了
while(1)
{
flag=0;
dfs(s,inf);
if(!flag)break;
}
}
}
signed main()
{
scanf("%lld%lld%lld%lld",&n,&m,&s,&t);
for(int i=1;i<=m;i++)
{
int x,y,z;
scanf("%lld%lld%lld",&x,&y,&z);
add(x,y,z);add(y,x,0);
}
dinic();
printf("%lld\n",ans);
return 0;
}
by masonpop @ 2023-05-23 18:33:04
问题解决了,是在 dfs 的过程中找到一条增广路之后不能立即停止,这一点和 FF 不同。
by Tibrella @ 2023-06-21 08:08:43
@masonpop 是当前弧优化的问题吧,我用找到一条增广路不立即停止的写法也会 T(
by masonpop @ 2023-06-21 08:11:35
@Tibrella 但是貌似我用多路增广不加当前弧优化也 A 了?
by Tibrella @ 2023-06-21 08:46:39
@masonpop 看看代码,感谢,我多路 T 了(
by masonpop @ 2023-06-21 08:56:58
代码。
@Tibrella
by Tibrella @ 2023-06-21 09:43:31
@masonpop 感谢,我加了一个剪枝之后过了