BZOJ 1797: [Ahoi2009]Mincut 最小割

题目在这里呀！

别话：对于凸包我表示暂时放一放（几何这个东西与我终生无缘），先做几道别的题吧！（虽然有点不甘心）

题意：给定一张图，询问图中的每条有向边是否能出现在某个最小割中，和是否一定属于所有的最小割。

题解：

当然要在残余网络上做啦，所以先是一遍最小割！然后进行强连通缩点（日常），S和T不在同一个强连通块中。

给定一条边u,v。如果u到v这条边没有满流，那么一定不属于任何最小割。

如果两个点不在同一个强连通块内，则可以属于某个最小割中。

如果belong[u]==belong[S] && belong[v]==belong[T]，那么u属于S的出流，v属于T的入流，则一定会割掉这条边，这条边就一定属于任何一最小割。

然后这道题就做完啦~

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <queue>
#define N 10000
#define M 300000
using namespace std;
const int inf=1e9;
struct edge{
	int from,to,next,cap;
}e[M+M];
int cnt,n,m,S,T,u[M],v[M],xx,dis[N],cur[N],head[N],ans,low[N],dfn[N],stack[N],belong[N],c[N],Time,t,tot;

inline void add_edge(int u,int v,int c){
	e[++cnt]=(edge){u,v,head[u],c};head[u]=cnt;
	e[++cnt]=(edge){v,u,head[v],0};head[v]=cnt;
}
bool bfs(){
	queue<int>q;
	for(int i=0;i<=n;i++) dis[i]=-1;
	q.push(S);dis[S]=0;
	while(!q.empty()){
		int now=q.front();q.pop();
		for(int i=head[now];i;i=e[i].next)
			if(e[i].cap && dis[e[i].to]==-1){
				dis[e[i].to]=dis[now]+1;
				q.push(e[i].to);
			}
	}
	return dis[T]!=-1;
}
int dfs(int x,int f){
	if(x==T) return f;
	int w,used=0;
	for(int i=cur[x];i;i=e[i].next)
		if(dis[e[i].to]==dis[x]+1){
			w=f-used;
			w=dfs(e[i].to,min(w,e[i].cap));
			e[i].cap-=w;e[i ^ 1].cap+=w;
			if(e[i].cap) cur[x]=i;
			used+=w;
			if(used==f) return f;
		}
	if(!used) dis[x]=-1;
	return used;
}
void dinic(){
	while(bfs()){
		for(int i=0;i<=n;i++) cur[i]=head[i];
		ans+=dfs(S,inf);
	}
}
void tarjan(int u){
	low[u]=dfn[u]=++Time;
	c[u]=1;
	stack[++tot]=u;
	for(int i=head[u];i;i=e[i].next){
		int v=e[i].to;
		if(!e[i].cap) continue;
		if(!dfn[v]){
			tarjan(v);
			low[u]=min(low[v],low[u]);
		}
		else if(c[v]) low[u]=min(low[u],dfn[v]);
	}
	if(low[u]==dfn[u]){
		int tmp=0;
		t++;
		while(tmp != u){
			tmp=stack[tot--];
			c[tmp]=0;
			belong[tmp]=t;
		}
	}
}
int main()
{
	scanf("%d%d%d%d",&n,&m,&S,&T);
	cnt=1;
	for(int i=1;i<=m;i++){
		scanf("%d%d%d",&u[i],&v[i],&xx);
		add_edge(u[i],v[i],xx);
	}
	dinic();
	for(int i=1;i<=n;i++) if(!dfn[i]) tarjan(i);
	for(int i=1;i<=m;i++){
		if(e[i+i].cap) puts("0 0");
		else{
			int ans1=0,ans2=0;
			if(belong[u[i]] != belong[v[i]]) ans1=1;
			if(belong[u[i]]==belong[S] && belong[v[i]]==belong[T]) ans2=1;
			printf("%d %d\n",ans1,ans2);
		}
	}
	return 0;
}