王茗仟 @ 2022-11-04 20:28:48
#define int long long
#define re register
#define in inline
#define max(a,b) ((a)>(b)?a:b)
#define min(a,b) ((a)<(b)?a:b)
using namespace std;
const int MAXN=1e4+10;
const int MAXM=2e5+19;
const int INF=0x3f3f3f3f3f;
in int read(){
re int x=0,f=0;re char c=getchar();
while(!isdigit(c)) f|=(c=='-'),c=getchar();
while(isdigit(c)) x=(x<<3)+(x<<1)+c-'0',c=getchar();
return f?-x:x;
}
int n,m,s,t,dis[MAXN],cur[MAXN];
struct edge{
int u,v,w;
int nx;
}e[MAXM];
int head[MAXN],tot;
in void add(int u,int v,int w){
e[tot].u=u;
e[tot].v=v;
e[tot].w=w;
e[tot].nx=head[u];
head[u]=tot++;
}
bool bfs(){
//广搜一下测一下剩下点的深度,顺便找出增广路
queue<int>q;
memset(dis,-1,sizeof(dis));
q.push(s);
dis[s]=0;
cur[s]=head[s];//cur是当前弧优化的宝贝
while(q.size()){
int u=q.front();
q.pop();
for(int i=head[u];~i;i=e[i].nx){
int v=e[i].v;
if(dis[v]==-1&&e[i].w){
dis[v]=dis[u]+1;
cur[v]=head[v];
if(v==t){
return 1;
}
q.push(v);
}
}
}
return 0;//不存在增广路了,可以结束le
}
int dfs(int u,int lim){
// 当前节点 能增加的最大流量
//深搜??也就是增广一下,给下流
if(u==t){
return lim;
}
int flow=0;//从u点分出去的
for(int i=cur[u];~i;i=e[i].nx){
cur[i]=i;//当前弧优化,很新对你来说,(走过就不再走了)
int v=e[i].v;
if(dis[v]==dis[u]+1&&e[i].w){
int minf=dfs(v,min(e[i].w,lim-flow));//求要给v多少
//limit-flow给这条边能分的最多有多少
e[i].w-=minf;
e[i^1].w+=minf;
flow+=minf;
if(flow==lim){
//全分完了直接结束吧
return flow;
}
}
}
return flow;
}
int dinic(){
int ans=0;
while(bfs()){
// cout<<ans<<endl;
ans+=dfs(s,INF);
}
return ans;
}
signed main(){
n=read();m=read();s=read();t=read();
memset(head,-1,sizeof(head));
for(re int i=1;i<=m;i++){
int u,v,w;
u=read();v=read();w=read();
add(u,v,w);add(v,u,0);
}
cout<<dinic()<<endl;
return 0;
}
by 王茗仟 @ 2022-11-04 20:40:45
@野生林登万
边计数是从0开始的,0也有一条边
by treeqwq @ 2022-11-04 20:40:50
@野生林登万 head初始化为-1了啊,为啥要改这里,而且用异或找反边必须这样写啊
by 王茗仟 @ 2022-11-04 20:42:45
@LSOT 我进行调试,发现他的bfs一直返回1,导致while一直进行,但是dfs返回一直是0,答案一直是20
by 野生林登万 @ 2022-11-04 20:44:45
这是咱的部分程序(cur默认为0):
int sta = cur[u];
if(!sta)sta = head[u];
for(int i = sta;i;i = e[i].nxt){
cur[u] = i;
const int v = e[i].to;
if(e[i].v <= 0 || val <= 0 || dep[v] != dep[u]+1)continue;
int push_push_uwu = min(val,e[i].v);
int actually_push_awa = dfs(v,push_push_uwu);
if(actually_push_awa <= 0)continue;
val -= actually_push_awa;
e[i].v -= actually_push_awa;
e[i^1].v += actually_push_awa;
ret += actually_push_awa;
}
可能写法有点不一样吧,重要的是for循环里面那部分,确实不知道为什么要这样取(反边为什么要在for循环里面取呢)。可不可尝试下从2开始取边
by 王茗仟 @ 2022-11-04 20:45:42
@野生林登万
这些都是一样的
by IceYukino @ 2022-11-04 20:47:53
@王茗仟 你当前弧优化出锅了
by treeqwq @ 2022-11-04 20:47:56
@王茗仟
《const int INF=0x3f3f3f3f3f;
》
楽
by 王茗仟 @ 2022-11-04 20:48:59
@LSOT
这没有问题的
by IceYukino @ 2022-11-04 20:49:06
@王茗仟 cur的初始化在if外面啊
by 王茗仟 @ 2022-11-04 20:49:26
@gqh_juruo
具体是怎样,我看不出来,而且我改了for的条件还是错误的