给定一张航空图,图中顶点代表城市,边代表两个城市间的直通航线。现要求找出一条满足下述限制条件的且途经城市最多的旅行路线。
从最西端城市出发,单向从西向东途经若干城市到达最东端城市,然后再单向从东向西飞回起点(可途经若干城市)。
除起点城市外,任何城市只能访问一次。
对于给定的航空图,试设计一个算法找出一条满足要求的最佳航空旅行路线。
输入格式
第一行有两个正整数 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;
}