1797: [Ahoi2009]Mincut 最小割

来源:互联网 发布:如何利用淘宝赚钱 编辑:程序博客网 时间:2024/06/03 06:13

题目链接

题目大意:给出一个有向图,每条边有流量,给出源点汇点s、t。对于每条边,询问:(1)是否存在一个最小割包含该边(2)是否所有的最小割都包含该边

题解:套路题……知道结论就好了

我的收获:套路++

#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>#include <stack>using namespace std;const int N=4005;const int M=120005;const int INF=1e9;int n,m,st,ed;int t,head[N],last[N];int num[N],d[N];int scnt,col[N];int tim,dfn[N],low[N];bool Exit,ins[N];stack<int> s;struct edge{int fro,to,c,nex;}e[M*2];void add(int u,int v,int w){e[t]=(edge){u,v,w,head[u]};last[u]=head[u]=t++;}void insert(int x,int y,int z){add(x,y,z),add(y,x,0);};int dfs(int x,int in){    if(x==ed) return in;    int ans=0,f;    for(int i=last[x];i!=-1;last[x]=i=e[i].nex)    {        int v=e[i].to;        if(e[i].c&&d[v]==d[x]-1){            f=dfs(v,min(in-ans,e[i].c));            ans+=f;e[i].c-=f;e[i^1].c+=f;            if(Exit||ans==in) return ans;        }    }    if(--num[d[x]]==0) Exit=1;//GAP优化    d[x]++,num[d[x]]++,last[x]=head[x];    return ans;}void tarjan(int x){    int now=0;    dfn[x]=low[x]=++tim;    s.push(x);ins[x]=1;    for(int i=head[x];i!=-1;i=e[i].nex){        int v=e[i].to;        if(!e[i].c) continue;        if(!dfn[v]) tarjan(v),low[x]=min(low[x],low[v]);        else if(ins[v]) low[x]=min(low[x],dfn[v]);    }    if(low[x]==dfn[x]){        scnt++;        while(x!=now){            now=s.top();s.pop();            ins[now]=0;col[now]=scnt;        }    }}void Isap(){    Exit=0;    while(!Exit) dfs(st,INF);}void work(){    Isap();    for(int i=1;i<=n;i++) if(!dfn[i]) tarjan(i);    for(int i=0;i<t;i+=2){        putchar((!e[i].c&&col[e[i].fro]!=col[e[i].to])?'1':'0');        putchar(' ');        putchar((!e[i].c&&col[st]==col[e[i].fro]&&col[ed]==col[e[i].to])?'1':'0');        putchar('\n');    } }void init(){    memset(head,-1,sizeof(head));    memset(last,-1,sizeof(last));    cin>>n>>m>>st>>ed;num[0]=n;    for(int x,y,z,i=1;i<=m;i++) scanf("%d%d%d",&x,&y,&z),insert(x,y,z);}int main(){    init();    work();    return 0;}