"Shortest" pair of paths (费用流)

http://poj.org/problem?id=3068

题意:就是两个东西从S到T且走的路不能重复,图中边有运输费用,问你最少花多少钱能到T

题解:   这样想,题中给的边先按边权为费用,容量为1,然后建立超级源点和超级汇点,源点连题中出发的点容量为2,边权为0,

汇点连题中给的到的点容量为2,边权为0。这样跑最大流,即满足最大流的最小费用(抽象一点想,边边之间容量为1,一次只能过一个东西,但汇点的地方又是容量为2的,就可以满足两个东西不重复的走到T啦)

代码如下:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int maxn=1000000;
const int inf=0x3f3f3f3f;
int head[maxn];//邻接表 
int d[maxn];//最短距离 dis 
int s,e,no;
int vis[maxn];//标记数组 
int pre[maxn];//记录前驱结点 方便  “反悔” 
struct point{
    int u,v,flow,next,cost;
    point(){};
    point(int x,int y,int z,int w,int c):u(x),v(y),next(z),flow(w),cost(c){};
}p[maxn];
void add(int x,int y,int z,int c){
    p[no]=point(x,y,head[x],z,c);   
    head[x]=no++;
    p[no]=point(y,x,head[y],0,-c);  //加负边 
    head[y]=no++;
}
void init(){
    memset(head,-1,sizeof(head));
    no=0;
}
bool spfa(){
    int i,x,y;
    queue<int>q;
    memset(d,0x3f,sizeof(d));
    memset(vis,false,sizeof(vis));
    memset(pre,-1,sizeof(pre));
    d[s]=0;   
    vis[s]=true;  
    q.push(s);
    while(!q.empty()){
        x=q.front();  
        q.pop();    
        vis[x]=false;
        for(i=head[x];i!=-1;i=p[i].next){
            if(p[i].flow&&d[y=p[i].v]>d[x]+p[i].cost){
                d[y]=d[x]+p[i].cost;   
                pre[y]=i;
                if(vis[y])  
                    continue;
                vis[y]=true;  
                q.push(y);
            }  
        }  
    }
    return d[e]!=d[e+1];
}
int mcmf(){
    int mincost=0,maxflow=0,minflow,i;
    while(spfa()){
        minflow=inf;
        for(i =pre[e];i!=-1;i=pre[p[i].u])
            minflow=min(minflow,p[i].flow);
        for(i=pre[e];i!=-1;i=pre[p[i].u]){
            p[i].flow-=minflow;
            p[i^1].flow+=minflow;
        }
        mincost+=d[e]*minflow; 
        maxflow+=minflow;
    }
    if(maxflow==2)//判读是否满流 判断最大流等于2 则有两条路 
        return mincost;
    return -1;
}
int n,m;
int main(){
    int res=1;//输出 Case :的 
    while(scanf("%d%d",&n,&m),n||m){
        init();//初始化 
        s=0;//源点为0 
        e=n+1;//汇点为N+1 
        add(s,1,2,0);//源点向点1件容量为2,花费为0的边 
        add(n,e,2,0);//点N向汇点建容量为2,花费为0的边 
        while(m--){
            int u,v,c;
            scanf("%d%d%d",&u,&v,&c);
            add(u+1,v+1,1,c);//每条边的起点和终点建容量为1,花费为边权的的边 
          //  add(v,u,1,c)//无向图加这个  把无向图拆成两个方向不同的有向图 
        }
        int ans=mcmf();
        printf("Instance #%d: ",res++);
        if(ans==-1)//不满流 
            printf("Not possible\n");
        else//满流则输出最小花费 
            printf("%d\n",ans);
    }
    return 0;
}

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值