【省选专题一】图论 jzoj 1665. 【AHOI2009】最小割

Description

　　A,B两个国家正在交战，其中A国的物资运输网中有N个中转站，M条单向道路。设其中第i (1≤i≤M)条道路连接了vi,ui两个中转站，那么中转站vi可以通过该道路到达ui中转站，如果切断这条道路，需要代价ci。现在B国想找出一个路径切断方案，使中转站s不能到达中转站t，并且切断路径的代价之和最小。
　　小可可一眼就看出，这是一个求最小割的问题。但爱思考的小可可并不局限于此。现在他对每条单向道路提出两个问题：
　　问题一：是否存在一个最小代价路径切断方案，其中该道路被切断？
　　问题二：是否对任何一个最小代价路径切断方案，都有该道路被切断？
　　现在请你回答这两个问题。

Input

　　第一行有4个正整数，依次为N,M,s和t。第2行到第(M+1)行每行3个正整数v,u,c表示v中转站到u中转站之间有单向道路相连，单向道路的起点是v，终点是u,切断它的代价是c(1≤c≤100000)。
　　注意:两个中转站之间可能有多条道路直接相连。
　　同一行相邻两数之间可能有一个或多个空格。

Output

　　对每条单向边，按输入顺序，依次输出一行，包含两个非0即1的整数，分别表示对问题一和问题二的回答(其中输出1表示是，输出0表示否)。
　　同一行相邻两数之间用一个空格隔开，每行开头和末尾没有多余空格。

Sample Input

6 7 1 6
1 2 3
1 3 2
2 4 4
2 5 1
3 5 5
4 6 2
5 6 3

Sample Output

1 0
1 0
0 0
1 0
0 0
1 0
1 0

Data Constraint

Hint

【样例说明】
　　设第(i+1)行输入的边为i号边，那么{1,2},{6,7},{2,4,6}是仅有的三个最小代价切割方案。它们的并是{1,2,4,6,7}，交是 。
【数据规模和约定】
　测试数据规模如下表所示:
数据编号 1 2 3 4 5 6 7 8 9 10
N= 10 20 200 200 1000 1000 1000 2000 3000 4000
M= 50 200 2000 2000 20000 20000 40000 50000 60000 60000

分析：
首先对原图跑最大流，得到残余网络。
很容易得知若一条边u->v存在u->v的增广路，则这条边一定不是割边，反之亦然则可能是割边。
那么我们就把点分成了s集和t集还有一部分不确定。若一条边连接了s集和t集，则这条边一定是必然存在度最小割中的边。

那么我们就得到了这样一个算法：
把残余网络中每个SCC缩成一个点，这样图中剩下的边则必然是满流边，且图中的每一个割都对应了原图中的一个最小割。设点i所在的强连通分量为id[i].那么一条边可能存在于割集中当且仅当id[u]!=id[v]。一条边一定存在于割集中当且仅当id[u]==id[s]且id[v]==id[t]。
（by ymw%%%） dalao表示已经忘记证明了= =

代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;

const int N=4005;
const int M=60005;

int n,m,s,t,cnt,last[N],cur[N],dfn[N],low[N],tim,sta[N],top,dis[N],tot,id[N];
bool ins[N];
struct edge{int to,next,c;}e[M*2];
queue<int> que;

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*10+ch-'0';ch=getchar();}
    return x*f;
}

void addedge(int u,int v,int c)
{
    e[++cnt].to=v;e[cnt].c=c;e[cnt].next=last[u];last[u]=cnt;
    e[++cnt].to=u;e[cnt].c=0;e[cnt].next=last[v];last[v]=cnt;
}

bool bfs()
{
    for (int i=1;i<=n;i++) dis[i]=0;
    while (!que.empty()) que.pop();
    dis[s]=1;que.push(s);
    while (!que.empty())
    {
        int u=que.front();que.pop();
        for (int i=last[u];i;i=e[i].next)
            if (e[i].c&&!dis[e[i].to])
            {
                dis[e[i].to]=dis[u]+1;
                if (e[i].to==t) return 1;
                que.push(e[i].to);
            }
    }
    return 0;
}

int dfs(int x,int maxf)
{
    if (x==t|!maxf) return maxf;
    int ret=0;
    for (int &i=cur[x];i;i=e[i].next)
        if (e[i].c&&dis[e[i].to]==dis[x]+1)
        {
            int f=dfs(e[i].to,min(e[i].c,maxf-ret));
            e[i].c-=f;
            e[i^1].c+=f;
            ret+=f;
            if (maxf==ret) break;
        }
    return ret;
}

void dinic()
{
    while (bfs())
    {
        for (int i=1;i<=n;i++) cur[i]=last[i];
        dfs(s,1000000000);
    }
}

void tarjan(int x)
{
    dfn[x]=low[x]=++tim;
    ins[x]=1;sta[++top]=x;
    for (int i=last[x];i;i=e[i].next)
    {
        if (!e[i].c) continue;
        if (!dfn[e[i].to])
        {
            tarjan(e[i].to);
            low[x]=min(low[x],low[e[i].to]);
        }
        else if (ins[e[i].to]) low[x]=min(low[x],dfn[e[i].to]);
    }
    if (low[x]==dfn[x])
    {
        tot++;int y=sta[top--];id[y]=tot;ins[y]=0;
        while (y!=x) y=sta[top--],id[y]=tot,ins[y]=0;
    }
}

int main()
{
    n=read();m=read();s=read();t=read();cnt=1;
    for (int i=1;i<=m;i++)
    {
        int x=read(),y=read(),z=read();
        addedge(x,y,z);
    }
    dinic();
    for (int i=1;i<=n;i++)
        if (!dfn[i]) tarjan(i);
    for (int i=2;i<=cnt;i+=2)
        if (e[i].c) puts("0 0");
        else
        {
            if (id[e[i].to]!=id[e[i^1].to]) printf("1 ");
                else printf("0 ");
            if (id[e[i].to]==id[t]&&id[e[i^1].to]==id[s]) puts("1");
                else puts("0");
        }
    return 0;
}