Our Journey of Dalian Ends && Our Journey of Xian Ends 最小费用最大流

2017 ACM-ICPC 亚洲区（乌鲁木齐赛区）网络赛 Our Journey of Dalian Ends

题意：要求先从大连到上海，再从上海打西安，中途会经过其他城市，每个城市只能去一次，出一次，给出航班信息，问最小花费。

每个城市只能去一次，出一次，那么很明显需要对每个城市拆点，就分成入点和出点，然后如果按照题意说法，把汇点连大连和上海，然后把上海和西安连汇点，那么很明显不对，

因为跑出来的可能是从上海直接到的上海，然后大连到西安，所以不能把跟汇点和源点相连的设在同一点，那么我们就可以汇点跟上海连，大连和西安跟源点连，或者反过来也可以。

#include<cstdio>
#include<cmath>
#include<queue>
#include<iostream>
#include<algorithm>
#include<tr1/unordered_map>
using namespace std;
const int N=2e4+11,M=1e6+11,inf=1e9+7;
struct Side{
    int v,ne,w,val;
}S[M<<1];
string s1[M],s2[M];
tr1::unordered_map<string,int> mmp;
int n,sn,sb,se,head[N],ww[N],vis[N],flow[N],lu[N],dis[N];
void init(){
    sn=0;
    sb=0;se=n*2+1;
    for(int i=sb;i<=se;i++) head[i]=-1;
}
void add(int u,int v,int w,int val){
    S[sn].w=w;S[sn].val=val;
    S[sn].v=v;S[sn].ne=head[u];
    head[u]=sn++;
}
void addE(int u,int v,int w,int val){
    add(u,v,w,val);add(v,u,0,-val); 
}
bool spfa(){
    queue<int> q;
    for(int i=sb;i<=se;i++){
        dis[i]=inf;
        vis[i]=0;
        flow[i]=inf;
        lu[i]=-1;
    }
    dis[sb]=0;
    vis[sb]=1;
    q.push(sb);
    int u,v;
    while(!q.empty()){
        u=q.front();q.pop();vis[u]=0;
        for(int i=head[u];~i;i=S[i].ne){
            v=S[i].v;
            if(S[i].w>0&&dis[v]>dis[u]+S[i].val){
                lu[v]=i;
                dis[v]=dis[u]+S[i].val;
                flow[v]=min(flow[u],S[i].w);
                if(!vis[v]){
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
    return dis[se]!=inf; 
}
int mfml(){
    int ans=0,ansc=0;
    while(spfa()){
        ans+=flow[se];
        ansc+=flow[se]*dis[se];
        for(int i=lu[se];~i;i=lu[S[i^1].v]){
            S[i].w-=flow[se];
            S[i^1].w+=flow[se];
        }
    }
    if(ans==2) return ansc;
    return -1;
}
int main(){
    int t,m;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&m);
        n=0;mmp.clear();
        for(int i=0;i<m;i++){
            cin>>s1[i]>>s2[i]>>ww[i];
            if(!mmp.count(s1[i])) mmp[s1[i]]=++n;
            if(!mmp.count(s2[i])) mmp[s2[i]]=++n;
        }
        init();
        for(int i=1;i<=n;i++) addE(i,i+n,1,0);
        addE(mmp["Shanghai"],mmp["Shanghai"]+n,1,0);
        addE(sb,mmp["Xian"],1,0);
        addE(sb,mmp["Dalian"],1,0);
        addE(mmp["Shanghai"]+n,se,2,0);
        for(int i=0;i<m;i++){
            addE(mmp[s1[i]]+n,mmp[s2[i]],1,ww[i]);
            addE(mmp[s2[i]]+n,mmp[s1[i]],1,ww[i]);
        }
        printf("%d
",mfml());
    }
    return 0;
} 

 ACM-ICPC 2017 Asia QingdaoOur Journey of Xian Ends

题意：跟上面类似，不过就是变成了，上海分成了两个之间有高速路的机场，虹桥跟浦东，然后这个高速路不计入那个进出里面。

明白题意之后，可以知道路线就西安到虹桥到青岛到浦东还有西安到浦东到虹桥到青岛到虹桥到浦东，这两条。

而不管是第一条还是第二条，如果只算航班信息的话，我们都可以统一出一个建图方法，源点连虹桥，流量为2，源点连浦东，流量为1，西安连汇点，流量为1，青岛连汇点流量为2，还有要注意青岛跟虹桥间的流量应该为2，应该它可以直接通过这个航班做一个轮回。

#include<cstdio>
#include<cmath>
#include<queue>
#include<iostream>
#include<algorithm>
#include<tr1/unordered_map>
using namespace std;
const int N=2e4+11,M=1e6+11,inf=1e9+7;
struct Side{
    int v,ne,w,val;
}S[M<<1];
string s1[M],s2[M];
tr1::unordered_map<string,int> mmp;
int n,sn,sb,se,head[N],ww[N],vis[N],flow[N],lu[N],dis[N];
void init(){
    sn=0;
    sb=0;se=n*2+1;
    for(int i=sb;i<=se;i++) head[i]=-1;
}
void add(int u,int v,int w,int val){
    S[sn].w=w;S[sn].val=val;
    S[sn].v=v;S[sn].ne=head[u];
    head[u]=sn++;
}
void addE(int u,int v,int w,int val){
    add(u,v,w,val);add(v,u,0,-val); 
}
bool spfa(){
    queue<int> q;
    for(int i=sb;i<=se;i++){
        dis[i]=inf;
        vis[i]=0;
        flow[i]=inf;
        lu[i]=-1;
    }
    dis[sb]=0;
    vis[sb]=1;
    q.push(sb);
    int u,v;
    while(!q.empty()){
        u=q.front();q.pop();vis[u]=0;
        for(int i=head[u];~i;i=S[i].ne){
            v=S[i].v;
            if(S[i].w>0&&dis[v]>dis[u]+S[i].val){
                lu[v]=i;
                dis[v]=dis[u]+S[i].val;
                flow[v]=min(flow[u],S[i].w);
                if(!vis[v]){
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
    return dis[se]!=inf; 
}
int mfml(){
    int ans=0,ansc=0;
    while(spfa()){
        ans+=flow[se];
        ansc+=flow[se]*dis[se];
        for(int i=lu[se];~i;i=lu[S[i^1].v]){
            S[i].w-=flow[se];
            S[i^1].w+=flow[se];
        }
    }
    if(ans==3) return ansc;
    return -1;
}
int main(){
    int t,m;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&m);
        n=0;mmp.clear();
        for(int i=0;i<m;i++){
            cin>>s1[i]>>s2[i]>>ww[i];
            if(!mmp.count(s1[i])) mmp[s1[i]]=++n;
            if(!mmp.count(s2[i])) mmp[s2[i]]=++n;
        }
        init();
        for(int i=1;i<=n;i++) addE(i,i+n,1,0);
        addE(mmp["Hongqiao"],mmp["Hongqiao"]+n,1,0);
        addE(mmp["Qingdao"],mmp["Qingdao"]+n,1,0);
        addE(sb,mmp["Pudong"],1,0);
        addE(sb,mmp["Hongqiao"],2,0);
        addE(mmp["Xian"]+n,se,1,0);
        addE(mmp["Qingdao"]+n,se,2,0);
        for(int i=0,w;i<m;i++){
            if(s1[i]=="Qingdao"&&s2[i]=="Hongqiao") w=2;
            else if(s2[i]=="Qingdao"&&s1[i]=="Hongqiao") w=2;
            else w=1;
            addE(mmp[s1[i]]+n,mmp[s2[i]],w,ww[i]);
            addE(mmp[s2[i]]+n,mmp[s1[i]],w,ww[i]);
            
        }
        printf("%d
",mfml());
    }
    return 0;
}