poj 3068 "Shortest" pair of paths (最小费用最大流)

题目链接：哆啦A梦传送门

题意：有n个点，m条有权值的单向边，现在让你找出两条从0点到n-1点的路径，要求这两条路径的费用尽可能的少，并且这两条边没有公共点（也就是所给的边只能使用一次）。要是找不到两条边，输出：Not possible。反之输出总费用。

题解：我们先将题目给的数据点都先+1。

建模：

源点0到第一个点建边，权值为2，费用为0。

将题目给的数据边连边，权值为1，费用为该边的费用。

最后我们将n点与汇点建边，权值为2，费用为0。

最后我们跑个最小费用最大流就好了。

 

代码:

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>

using namespace std;

const int N=100;
const int INF=0x3f3f3f3f;

struct edge{
    int u,v,flow,cost,next;
    edge(){}
    edge(int _u,int _v,int _flow,int _cost,int _next){
        u=_u;v=_v;flow=_flow;cost=_cost;next=_next;
    }
}e[30000];

int head[N],cnt;
queue<int> que;

void init()
{
    memset(head,-1,sizeof(head));
    cnt=0;
}

void add(int x,int y,int z,int c)
{
    e[cnt]=edge(x,y,z,c,head[x]);
    head[x]=cnt++;
    e[cnt]=edge(y,x,0,-c,head[y]);
    head[y]=cnt++;
}

bool vis[N];
int dis[N],pre[N];

bool spfa(int s,int t)
{
    memset(dis,INF,sizeof(dis));
    memset(vis,0,sizeof(vis));
    memset(pre,-1,sizeof(pre));

    dis[s]=0;
    vis[s]=true;
    while(!que.empty()) que.pop();

    que.push(s);

    while(!que.empty())
    {
        int x=que.front();
        que.pop();
        vis[x]=false;

        for(int i=head[x];i!=-1;i=e[i].next){
            if(e[i].flow){
                int y=e[i].v;
                if(dis[x]+e[i].cost<dis[y]){
                    dis[y]=dis[x]+e[i].cost;
                    pre[y]=i;
                    if(!vis[y]){
                        que.push(y);
                        vis[y]=true;
                    }
                }
            }
        }

    }
     return pre[t]!=-1;

}

int costflow(int s,int t,int x)
{
    int ret=0;
    int max_flow=0;

    while(spfa(s,t))
    {
        int min_flow=INF;
        for(int i=pre[t];i!=-1;i=pre[e[i].u]){
            min_flow=min(e[i].flow,min_flow);
        }
        for(int i=pre[t];i!=-1;i=pre[e[i].u]){
            e[i].flow-=min_flow;
            e[i^1].flow+=min_flow;
        }
        ret+=dis[t]*min_flow;
        max_flow+=min_flow;
    }

    if(max_flow<x) return -1;///达不到满流
    return ret;
}

int  main()
{

    int n,m,k;
    int T=0;

    while(~scanf("%d%d",&n,&m))
    {
        if(n==0&&m==0) break;
        init();

        int s=0,t=n+1;
        add(s,1,2,0);
        add(n,t,2,0);

        int x,y,z;
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d%d",&x,&y,&z);
            x++;y++;
            add(x,y,1,z);
        }
        int ans=costflow(s,t,2);
        printf("Instance #%d: ",++T);
        if(ans==-1){
            puts("Not possible");
        }
        else printf("%d\n",ans);
    }
    return 0;
}