loj6122「网络流 24 题」航空路线问题

(http://www.elijahqi.win/2017/12/18/loj6122%E3%80%8C%E7%BD%91%E7%BB%9C%E6%B5%81-24-%E9%A2%98%E3%80%8D%E8%88%AA%E7%A9%BA%E8%B7%AF%E7%BA%BF%E9%97%AE%E9%A2%98/%20%E2%80%8E)
题目描述

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

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

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

第一行有两个正整数 NNN 和 VVV，NNN 表示城市数，VVV 表示直飞航线数。
接下来的 NNN 行中每一行是一个城市名，可乘飞机访问这些城市。城市名出现的顺序是从西向东。也就是说，设 i,ji,ji,j 是城市表列中城市出现的位置次序，当 i>ji>ji>j 时，表示 城市 iii 在城市 jjj 的东边，而且不会有两个城市在同一条经线上。城市名是一个长度不超过 151515 的字符串，串中的字符可以是大小写字母或阿拉伯数字。例如，AGR34\text{AGR34}AGR34 或 BEL4\text{BEL4}BEL4。
再接下来的 VVV 行中，每行有两个城市名，中间用空格隔开，如 city1 city2\text{city1 city2}city1 city2 表示 city1\text{city1}city1 到 city2\text{city2}city2 有一条直通航线，从 city2\text{city2}city2 到 city1\text{city1}city1 也有一条直通航线。
输出格式

输出最佳航空旅行路线。
第一行是旅行路线中所访问的城市总数 MMM。
接下来的 M+1M+1M+1 行是旅行路线的城市名，每行一个。首先是出发城市名，然后按访问顺序列出其它城市名。注意，最后一行（终点城市）的城市名必然是出发城市名。如果有多组最优解，输出任意一组均可；如果问题无解，则输出 No Solution!。
样例
样例输入

8 9
Vancouver
Yellowknife
Edmonton
Calgary
Winnipeg
Toronto
Montreal
Halifax
Vancouver Edmonton
Vancouver Calgary
Calgary Winnipeg
Winnipeg Toronto
Toronto Halifax
Montreal Halifax
Edmonton Montreal
Edmonton Yellowknife
Edmonton Calgary

样例输出

7
Vancouver
Edmonton
Montreal
Halifax
Toronto
Winnipeg
Calgary
Vancouver

数据范围与提示

对于所有数据，N<100N < 100N<100

题意要求相当于我需要求两条不重合的1~n的路径 使得他们经过的点最多

那么我可以把每个点拆成两个点中间容量为1 限制一下 只能使用一次 然后每个点的使用费用为1 然后那么我跑最大费用最大流 就可以求出了 最后我

两遍dfs把路径搜索出来即可 注意存在一些情况诸如只存在起点和终点直接可达 但是他们中间是不存在其他边可到达 的这时候我就需要判断一下最大流是不是1

如果是1 就直接输出起点和终点即可

#include<map>
#include<queue>
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#define N 220
#define inf 0x3f3f3f3f
using namespace std;
map<string,int> mp;bool flag[N];
int num=1,h[N],f[N];
struct node{
    int x,y,z,next,c;
}data[N*N];
inline void insert1(int x,int y,int z,int c){
    data[++num].y=y;data[num].z=z;data[num].next=h[x];h[x]=num;data[num].c=c;data[num].x=x;
    data[++num].y=x;data[num].z=0;data[num].next=h[y];h[y]=num;data[num].c=-c;data[num].x=y;
}
int T,pre[N],path[N],n,v;
inline bool spfa(){
    memset(f,128,sizeof(f));memset(flag,0,sizeof(flag));flag[1]=1;f[1]=0;queue<int>q;memset(pre,-1,sizeof(pre));q.push(1);
    while(!q.empty()){
        int x=q.front();q.pop();flag[x]=0;
        for (int i=h[x];i;i=data[i].next){
            int y=data[i].y,z=data[i].z,c=data[i].c;
            if (f[x]+c>f[y]&&z){
                f[y]=f[x]+c;pre[y]=x;path[y]=i;
                if (!flag[y]) flag[y]=1,q.push(y);
            }
        }
    }
    if (pre[T]==-1) return 0;else return 1;
}
string name[N];
int main(){
    freopen("2770.in","r",stdin);
//  freopen("2770.out","w",stdout);
    std::ios::sync_with_stdio(false);
    cin>>n>>v;string st1,st2;T=n<<1;
    for (int i=1;i<=n;++i){cin>>st1;mp[st1]=i;name[i]=st1;insert1(i,i+n,1,0);}
    for (int i=1;i<=v;++i){
        cin>>st1>>st2;int x=mp[st1],y=mp[st2];if (x>y) swap(x,y);insert1(x+n,y,1,1);
    //  printf("%d %d\n",x,y);
    }insert1(1,1+n,1,0);insert1(n,n<<1,1,0);int ans=0,ans1=0;
    //for (int i=2;i<=num;++i) printf("%d %d %d %d\n",data[i].x,data[i].y,data[i].z,data[i].c);
    while(spfa()){
        int minn=inf,now=T;
        while(now!=1) minn=min(minn,data[path[now]].z),now=pre[now];ans+=minn;now=T;
        while(now!=1) {ans1+=data[path[now]].c*minn;data[path[now]].z-=minn;data[path[now]^1].z+=minn;now=pre[now];}
    }bool visit[N];
    if(ans==1){cout<<2<<endl<<name[1]<<endl<<name[n]<<endl<<name[1];return 0;} 
    if (ans!=2){printf("No Solution!");return 0;}bool flag1=0,print=1;int x=1;cout<<ans1<<endl;
    while(1){
        flag1=0;if (print)cout<<name[x]<<endl;visit[x]=1;
        for (int i=h[x];i;i=data[i].next){
            if (data[i].z==0){
                x=data[i].y;if (x>n) print=0;else print=1;flag1=1;
            }
        }
        if (!flag1) break;
    }visit[1]=visit[1+n]=0;flag1=0,print=0;x=n+1;vector<string> aq;
    while(1){
        flag1=0;if (print) aq.push_back(name[x]);visit[x]=1;
        for (int i=h[x];i;i=data[i].next){
            if (data[i].z==0&&!visit[data[i].y]){
                x=data[i].y;if (x>n) print=0;else print=1;flag1=1;
            }
        }
        if (!flag1) break;
    }
    for (int i=aq.size()-1;i>=0;--i) cout<<aq[i]<<endl;cout<<name[1]<<endl;
//  for (int i=1;i<=n;++i) printf("%d %d %d\n",i,visit[i],visit[i+n]);
    return 0;
}