題目描述
給定一張航空圖，圖中頂點代表城市，邊代表 2 城市間的直通航線。現要求找出一條滿足下述限制條件的且途經城市最多的旅行路線。
(1)從最西端城市出發，單向從西向東途經若干城市到達最東端城市，然後再單向從東向西飛回起點(可途經若干城市)。
(2)除起點城市外，任何城市只能訪問 1 次。
對於給定的航空圖，試設計一個算法找出一條滿足要求的最佳航空旅行路線。
輸入格式
第 1 行有 2 個正整數 N 和 V，N 表示城市數，N<100，V 表示直飛航線數。
接下來的 N 行中每一行是一個城市名，可乘飛機訪問這些城市。城市名出現的順序是從西向東。也就是說，設 i,j 是城市表列中城市出現的順序，當 i>j 時，表示城市 i 在城市 j 的東邊，而且不會有 2 個城市在同一條經線上。城市名是一個長度不超過15 的字符串，串中的字符可以是字母或阿拉伯數字。例如，AGR34 或 BEL4。
再接下來的 V 行中，每行有 2 個城市名，中間用空格隔開，如 city1 city2 表示 city1到 city2 有一條直通航線，從 city2 到 city1 也有一條直通航線。
輸出格式
件第 1 行是旅行路線中所訪問的城市總數 M。 接下來的 M＋1 行是旅行路線的城市名，每行寫 1 個城市名。首先是出發城市名，然後按訪問順序列出其它城市名。 注意，最後 1 行（終點城市）的城市名必然是出發城市名。如果問題無解，則輸出“No Solution!”。
輸入輸出樣例
輸入 #1複製
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
輸出 #1複製
7
Vancouver
Edmonton
Montreal
Halifax
Toronto
Winnipeg
Calgary
Vancouver
說明/提示
感謝 @FlierKing 提供spj
思路：
拆點，將兩點之間的邊的花費設置爲-1,則要直接跑最小費用最大流就可以。需要注意的是起點和終點有連邊，這時最大流爲1的情況。之後做題需要 特殊考慮只有一個點，兩個點的情況。
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<string>
#include<vector>
#include<unordered_map>
#define mod (1000000007)
using namespace std;
typedef long long ll;
const int MAX = 2e5 + 5;
const int inf = 0x3f3f3f3f;
struct node {
int v,f,w,nxt;
} e[50005<<2];
int n,m,N;
int r[MAX];
int head[MAX],d[MAX],vis[MAX],tot=1,p[MAX];
void add(int u,int v,int f,int cost=0) {
e[++tot].v = v;e[tot].f = f;e[tot].w = cost;e[tot].nxt = head[u];head[u] = tot;
e[++tot].v = u;e[tot].f = 0; e[tot].w = -cost;e[tot].nxt = head[v];head[v] = tot;
}
bool bfs(int s,int t) {
for(int i = 0; i<=N; i++) d[i]=inf,vis[i]=0;
d[s]=0;
queue<int>q;
q.push(s);
while(!q.empty()) {
int u=q.front();
q.pop();
vis[u]=0;
for(int i=head[u]; ~i; i=e[i].nxt) {
int j=e[i].v;
if(e[i].f&&d[j]>d[u]+e[i].w) {
d[j]=d[u]+e[i].w;
p[j]=i;
if(!vis[j])vis[j]=1,q.push(j);
}
}
}
return d[t]<inf;
}
ll MCMF(int s,int t,ll &flow) {
ll ans=0;
while(bfs(s,t)) {
int x=t,f=inf;
while(x!=s) {
f = min(f,e[p[x]].f),x=e[p[x]^1].v;
}
flow += f;
ans+=1LL*d[t]*f;
x=t;
while(x!=s) {
e[p[x]].f-=f,e[p[x]^1].f+=f;
x=e[p[x]^1].v;
}
}
return ans;
}
unordered_map<string,int> mp;
unordered_map<int,string> name;
void bfs2(int st,int &ct,int a[]){
queue<int> q;
q.push(st);
while(!q.empty()){
int nn=q.front();q.pop();
a[ct++]=nn;
for(int i=head[nn];i!=-1;i=e[i].nxt){
if(e[i].f==0&&i<(i^1))//反向邊
q.push(e[i].v);
}
}
}
char s[1000],ss[1000];
int ans[1000],ans2[1000];
int main() {
scanf("%d%d",&n,&m);
for(int i=0;i<=2*n+1;i++) head[i]=-1;
for(int i=1;i<=n;i++){
scanf("%s",s);
mp[s]=i;name[i]=s;
if(i!=n&&i!=1)
add(i,i+n,1,-1);
else add(i,i+n,2,-1);
}
int st=1,ed=n+n,ff=0;N=ed;
for(int i=1;i<=m;i++){
scanf("%s%s",s,ss);
int t1=mp[s];
int t2=mp[ss];
if(t1>t2) swap(t1,t2);
if(t1==1&&t2==n)ff=1;
add(t1+n,t2,1,0);
}
ll fl=0,an=0;
an=MCMF(st,ed,fl);
if(fl<1)puts("No Solution!");
else if(fl==1){
if(ff){//如果起點到終點有直接的路徑，即使只有1條也是滿足條件的
cout<<2<<endl;
cout<<name[1]<<endl;
cout<<name[n]<<endl;
cout<<name[1]<<endl;
}
else puts("No Solution!");
}
else{
printf("%lld\n",-an-2);
queue<int> q;
for(int i=head[1+n];i!=-1;i=e[i].nxt){
if(e[i].f==0){
q.push(e[i].v);
}
}
int tp1=1,tp2=1;
ans[0]=1;
ans2[0]=1;
bfs2(q.front(),tp1,ans);q.pop();
bfs2(q.front(),tp2,ans2);q.pop();
for(int i=0;i<tp1;i++)
if(ans[i]<=n) //printf("%d\n",ans[i]);
cout<<name[ans[i]]<<endl;
for(int i=tp2-1;i>=0;i--)
if(ans2[i]<n) //printf("%d\n",ans2[i]);
cout<<name[ans2[i]]<<endl;
}
return 0 ;
}