bzoj1797: [Ahoi2009]Mincut 最小割(最小割+强联通tarjan)

1797: [Ahoi2009]Mincut 最小割

题目：传送门 

---

题解：

　　 感觉是一道肥肠好的题目。

　　 第二问其实比第一问简单？

　　 用残余网络跑强联通，流量大于0才访问。

　　 那么如果两个点所属的联通分量分别处于st和ed，那一定会被割掉，那么第一问就也会是1

　　 但是第一问单独处理，就有点GG

　　 膜了一发题解，发现贼尼玛niu bi：

　　 还是利用联通分量，如果这条边满流，且连接的两个点所处的联通分量不同，就ok

　　 太菜了不会证明...大佬hzwer

 　　 

---

代码：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define inf 999999999
using namespace std;
typedef long long LL;
struct node
{
    int x,y,next,other;LL c;
}a[2110000];int len,last[1110000];
void ins(int x,int y,LL c)
{
    int k1,k2;
    k1=++len;
    a[len].x=x;a[len].y=y;a[len].c=c;
    a[len].next=last[x];last[x]=len;
    
    k2=++len;
    a[len].x=y;a[len].y=x;a[len].c=0;
    a[len].next=last[y];last[y]=len;
    
    a[k1].other=k2;
    a[k2].other=k1;
}
int n,m,st,ed,head,tail;
int list[1110000],h[1110000];
bool bt_h()
{
    memset(h,0,sizeof(h));h[st]=1;
    list[1]=st;head=1;tail=2;
    while(head!=tail)
    {
        int x=list[head];
        for(int k=last[x];k;k=a[k].next)
        {
            int y=a[k].y;
            if(h[y]==0 && a[k].c>0)
            {
                h[y]=h[x]+1;
                list[tail++]=y;
            }
        }
        head++;
    }
    if(h[ed]>0)return true;
    return false;
}
LL find_flow(int x,LL flow)
{
    if(x==ed)return flow;
    LL s=0,t;
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(h[y]==h[x]+1 && a[k].c>0 && s<flow)
        {
            s+=t=find_flow(y,min(a[k].c,flow-s));
            a[k].c-=t;a[a[k].other].c+=t;
        }
    }
    if(s==0)h[x]=0;
    return s;
}
int low[1110000],dfn[1110000],belong[1110000],sta[1110000];
int tp,id,cnt;
bool v[1110000];
void dfs(int x)
{
    low[x]=dfn[x]=++id;
    sta[++tp]=x;v[x]=true;
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(a[k].c!=0)
        {
            if(dfn[y]==-1)
            {
                dfs(y);
                low[x]=min(low[x],low[y]);
            }
            else
                if(v[y]==true)
                    low[x]=min(low[x],dfn[y]);
        }
    }
    if(low[x]==dfn[x])
    {
        int i;cnt++;
        do{
            i=sta[tp--];
            v[i]=false;
            belong[i]=cnt;
        }while(i!=x);
    }
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&st,&ed);
    len=0;memset(last,0,sizeof(last));
    for(int i=1;i<=m;i++)
    {
        int x,y;LL c;
        scanf("%d%d%lld",&x,&y,&c);
        ins(x,y,c);
    }
    LL ans=0;
    while(bt_h())ans+=find_flow(st,inf);
    memset(sta,0,sizeof(sta));
    memset(dfn,-1,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(v,false,sizeof(v));
    id=tp=cnt=0;
    for(int i=1;i<=n;i++)
        if(dfn[i]==-1)
            dfs(i);
    for(int i=1;i<len;i+=2)
    {
        int x=a[i].x,y=a[i].y;
        if(a[i].c!=0){printf("0 0
");continue;}
        if(belong[x]!=belong[y])printf("1 ");
        else printf("0 ");
        if((belong[x]==belong[st] && belong[y]==belong[ed]))
        printf("1
");
        else printf("0
");
    }
    return 0;
}