#6122. 「网络流 24 题」航空路线问题 两条不相交路径

题目链接：https://loj.ac/problem/6122

 

题意：

     给定一张航空图，图中顶点代表城市，边代表两个城市间的直通航线。现要求找出一条满足下述限制条件的且途经城市最多的旅行路线。

1. 从最西端城市出发，单向从西向东途经若干城市到达最东端城市，然后再单向从东向西飞回起点（可途经若干城市）。
2. 除起点城市外，任何城市只能访问一次。  

     对于给定的航空图，试设计一个算法找出一条满足要求的最佳航空旅行路线。

做法：

     其实题目的意思还是挺明显的就是去找从点1到点n的两条不相交路径，然后走一个环就好了，但是乍一想还真的是想不到网络流上去。。要不是因为这个是网络流24题的话。。如果想到网络流的话，应该也不是很难了，就是很普遍的拆点问题（防止路径相交），然后就去做dfs就好了。

---

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=400000;
const int maxm=1000000;
const int inf=0x3f3f3f3f;
int dis[maxn];
int vis[maxn],pre[maxn];
int head[maxn],cnt;
int n,m,sp,tp;
map<string,int> mp;
map<int,string> mps;
int ans=0;
struct node{
    int to,cap,cost,next,flag;
}e[maxm];
void add(int from,int to,int cap,int cost){

    //printf("from is %d, to is %d,cap is %d cost is %d\n",from,to,cap,cost);
    e[cnt].to=to; e[cnt].cap=cap;
    e[cnt].cost=cost; e[cnt].next=head[from];
    e[cnt].flag=1;head[from]=cnt++;

    e[cnt].to=from; e[cnt].cap=0;
    e[cnt].cost=-cost; e[cnt].next=head[to];
    e[cnt].flag=0;head[to]=cnt++;
}
bool spfa(int s,int t,int &flow,int &cost){
    queue<int> q;
    memset(dis,inf,sizeof(dis));
    memset(vis,0,sizeof(vis));
    memset(pre,-1,sizeof(pre));
    dis[s]=0;  q.push(s);
    vis[s]=1;
    int d=inf;
    while(!q.empty()){
        int u=q.front();
        q.pop();
        vis[u]=0;
        for(int i=head[u];~i;i=e[i].next){
            int v=e[i].to;
            if(e[i].cap>0&&dis[v]>dis[u]+e[i].cost){
                dis[v]=dis[u]+e[i].cost;
                pre[v]=i;
                if(!vis[v]){
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
    if(dis[t]==inf){
        return false;
    }
    for(int i=pre[t];~i;i=pre[e[i^1].to]){
        d=min(d,e[i].cap);
    }
    for(int i=pre[t];~i;i=pre[e[i^1].to]){
        e[i].cap-=d;
        e[i^1].cap+=d;
        cost+=e[i].cost*d;
    }
    flow+=d;
    return true;
}
int mcmf(int s,int t){
    int flow=0,cost=0;
    while(spfa(s,t,flow,cost)){
        //cout<<flow<<" "<<cost<<endl;
    }
    if(flow!=2) return -1;
    return cost;
}
int id;
int gain(string s,int ju){
    if(!mp[s]) {
        mp[s]=++id;
        mps[id]=s;
    }
    if(ju) return mp[s]+n;
    else return mp[s];
}
void dfs(int u,int f){
    if(u==n) return ;
    for(int i=head[u];~i;i=e[i].next){
        if(e[i].flag==0) continue;
        int v=e[i].to;
        if(v==f) continue;
        if(e[i].cap==0) {
            e[i].flag=0;
            if(v<=n) cout<<mps[v]<<endl;
            dfs(v,u);
            break;
        }
    }
}
stack<string> S;
void dfs1(int u,int f){
    for(int i=head[u];~i;i=e[i].next){
        if(e[i].flag==0) continue;
        int v=e[i].to;
        if(v==f) continue;
        if(e[i].cap==0) {
            e[i].flag=0;
            if(v<n) S.push(mps[v]);
            dfs1(v,u);
            break;
        }
    }
}
int main(){
    memset(head,-1,sizeof(head));
    ans=0;
    scanf("%d%d",&n,&m);
    sp=0;tp=201;
    add(sp,1,2,-1);
    add(n,tp,2,-1);
    for(int i=1;i<=n;i++){
        string s;
        cin>>s;
        gain(s,0);
        if(i!=1&&i!=n) add(gain(s,0),gain(s,1),1,-1);
    }
    for(int i=0;i<m;i++){
        string s1,s2;
        cin>>s1>>s2;
        int id1=gain(s1,0),id2=gain(s2,0);
        if(id1>id2) swap(id1,id2);
        if(id1!=1) id1+=n;
        if(id1==1&&id2==n){
            add(id1,id2,2,0);
        }
        else add(id1,id2,1,0);
    }
    ans=mcmf(sp,tp);
    if(ans==-1) printf("No Solution!\n");
    else {
        printf("%d\n",-(ans+2));
        cout<<mps[1]<<endl;
        dfs(1,-1);
        dfs1(1,-1);
        while(!S.empty()){
            string u = S.top();
            S.pop();
            cout<<u<<endl;
        }
        cout<<mps[1]<<endl;
    }

    return 0;
}