最小割的必须边
一定在最小割中的边、扩大容量后能增大最大流的边,Poj3204:
① 满流;② 残余网络中S能到入点、出点能到T。
从S开始DFS、T开始反向DFS,标记到达的点,然后枚举满流边即可。
最小割的可行边
被某一种最小割的方案包含的边,AHOI2009:
① 满流;② 删掉之后在残余网络中找不到u到v的路径。
在残余网络中tarjan求SCC,(u,v)两点在同一SCC中说明残余网络中存在u到v路径。
#include<cstdio>
#include<iostream>
#include<cstring>
#include<vector>
#include<queue>
#include<stack>
using namespace std;
const int maxn = 10005;
const int inf = 0x3f3f3f3f;
struct edge
{
int u,v,cap,flow;
};
vector<int> G[maxn];
vector<edge> f;
void add_edge(int u,int v,int cap)
{
f.push_back((edge){u,v,cap,0});
f.push_back((edge){v,u,0,0});
int m=f.size();
G[u].push_back(m-2);
G[v].push_back(m-1);
}
stack<int> S;
int d[maxn],cur[maxn],ans[maxn][2],x,y,z,n,m,s,t;
int pre[maxn],low[maxn],scc[maxn],dfn,cnt;
bool vis[maxn],map[505][505];
bool bfs()
{
memset(vis,false,sizeof(vis));
queue<int> Q;
Q.push(s);
vis[s]=true;d[s]=0;
while(!Q.empty())
{
int u=Q.front();
Q.pop();
for(int i=0;i<G[u].size();i++)
{
edge& e=f[G[u][i]];
if(e.cap>e.flow&&!vis[e.v])
{
d[e.v]=d[u]+1;
vis[e.v]=true;
Q.push(e.v);
}
}
}
return vis[t];
}
int dfs(int u,int a)
{
int flow=0,h;
if(u==t||a==0) return a;
for(int &i=cur[u];i<G[u].size();i++)
{
edge& e=f[G[u][i]];
if(d[e.v]==d[u]+1&&(h=dfs(e.v,min(e.cap-e.flow,a)))>0)
{
flow+=h;
f[G[u][i]^1].flow-=h;
e.flow+=h;
a-=h;
if(a==0) break;
}
}
return flow;
}
int Dinic(int s,int t)
{
int flow=0;
while(bfs()) {
memset(cur,0,sizeof(cur));
flow+=dfs(s,inf);
}
return flow;
}
void dfs(int u)
{
pre[u]=low[u]=++dfn;S.push(u);
for(int i=0;i<G[u].size();i++)
{
edge &e=f[G[u][i]];
if(e.cap>e.flow)
{
int v=e.v;
if(!pre[v]) dfs(v),low[u]=min(low[u],low[v]);
else if(!scc[v]) low[u]=min(low[u],pre[v]);
}
}
if(pre[u]==low[u]){
cnt++;
for(;;)
{
int x=S.top();
S.pop();scc[x]=cnt;
if(x==u) break;
}
}
}
int main()
{
scanf("%d%d%d%d",&n,&m,&s,&t);
for(int i=1;i<=m;i++)
{
scanf("%d%d%d",&x,&y,&z);
add_edge(x,y,z);
}
Dinic(s,t);
for(int i=1;i<=n;i++) if(!pre[i]) dfs(i);
for(int i=0;i<f.size();i+=2)
{
if(f[i].cap>f[i].flow) puts("0 0");
else
{
if(scc[f[i].u]!=scc[f[i].v]) printf("1 ");
else printf("0 ");
if(scc[f[i].u]==scc[s]&&scc[f[i].v]==scc[t]) puts("1");
else puts("0");
}
}
return 0;
}