bzoj1797: [Ahoi2009]Mincut 最小割

首先一定跑一遍最小割，在残量网络上tarjan一遍
（1） 满足u和v不在同一个强连通分量且有流通过。
（2）满足（1）且u在s的强连通分量，v在t的强连通分量。

#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#include<queue>#include<stack>using namespace std;stack<int> s1;queue<int> q;int read(){    char ch=getchar();int f=0;    while(ch>'9'||ch<'0') ch=getchar();    while(ch>='0'&&ch<='9') {f=f*10+(ch^48);ch=getchar();}    return f;}int n,m,s,t,dis[4005],head[4005],dfn[4005],low[4005],pos[4005],size[4005];int tim,tot,cnt;long long ans;bool vis[4005];struct node{    int from;    int to;    int next;    int w;}edge[120005];void add(int u,int v,int w){    edge[tot].from=u;    edge[tot].to=v;    edge[tot].w=w;    edge[tot].next=head[u];    head[u]=tot++;}bool bfs(){    memset(dis,0,sizeof(dis));    dis[s]=1;    q.push(s);    while(!q.empty())    {        int x=q.front();        q.pop();        for(int i=head[x];i!=-1;i=edge[i].next)        {            if(!dis[edge[i].to]&&edge[i].w)            {                dis[edge[i].to]=dis[x]+1;                q.push(edge[i].to);            }        }    }    if(!dis[t])    return 0;    return 1;}int dfs(int x,int flow){    if(x==t||!flow)    return flow;    int ret=0;    for(int i=head[x];i!=-1;i=edge[i].next)    {        if(dis[edge[i].to]==dis[x]+1)        {            int f=dfs(edge[i].to,min(flow,edge[i].w));            edge[i].w-=f;edge[i^1].w+=f;            ret+=f;            flow-=f;            if(!flow)            break;        }    }    if(!ret)    dis[x]=-1;    return ret;}void dinic(){    while(bfs())    {        ans+=dfs(s,0x7f7f7f7f);    }}void tarjan(int x){    dfn[x]=low[x]=++tim;    s1.push(x);    vis[x]=1;    for(int i=head[x];i!=-1;i=edge[i].next)    {        if(edge[i^1].w)        {            if(!dfn[edge[i].to])            {                tarjan(edge[i].to);                 low[x]=min(low[x],low[edge[i].to]);            }            else if(vis[edge[i].to])            {                low[x]=min(low[x],dfn[edge[i].to]);            }         }    }    if(low[x]==dfn[x])    {        cnt++;        while(s1.top()!=x)        {            pos[s1.top()]=cnt;            size[cnt]++;            vis[s1.top()]=0;            s1.pop();        }         s1.pop();        pos[x]=cnt;        size[cnt]++;        vis[x]=0;    } }int main(){    //freopen("input.txt","r",stdin);    //freopen("output.txt","w",stdout);    memset(head,-1,sizeof(head));    n=read(),m=read(),s=read(),t=read();    int x,y,z;    for(int i=1;i<=m;i++)    {        x=read(),y=read(),z=read();        add(x,y,z);        add(y,x,0);    }    dinic();    for(int i=1;i<=n;i++)    {        if(!dfn[i])        tarjan(i);    }    for(int i=0;i<tot;i+=2)    {        if(pos[edge[i].from]!=pos[edge[i].to]&&!edge[i].w)        {            printf("1 ");            if(pos[s]==pos[edge[i].from]&&pos[t]==pos[edge[i].to])            {                puts("1");            }            else puts("0");         }        else        {            puts("0 0");        }    }    //for(int i=1;i<=n;i++)    //printf("%d ",pos[i]);}