【网络流24题-洛谷-P2770】航空路线问题（最长不相交路径、最小费用最大流）

题目描述

给定一张航空图，图中顶点代表城市，边代表 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 ;
}