【省选专题一】图论 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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值