L3-014. 周游世界
周游世界是件浪漫事,但规划旅行路线就不一定了…… 全世界有成千上万条航线、铁路线、大巴线,令人眼花缭乱。所以旅行社会选择部分运输公司组成联盟,每家公司提供一条线路,然后帮助客户规划由联盟内企业支持的旅行路线。本题就要求你帮旅行社实现一个自动规划路线的程序,使得对任何给定的起点和终点,可以找出最顺畅的路线。所谓“最顺畅”,首先是指中途经停站最少;如果经停站一样多,则取需要换乘线路次数最少的路线。
输入格式:
输入在第一行给出一个正整数N(<= 100),即联盟公司的数量。接下来有N行,第i行(i=1, ..., N)描述了第i家公司所提供的线路。格式为:
M S[1] S[2] ... S[M]
其中M(<= 100)是经停站的数量,S[i](i=1, ..., M)是经停站的编号(由4位0-9的数字组成)。这里假设每条线路都是简单的一条可以双向运行的链路,并且输入保证是按照正确的经停顺序给出的 —— 也就是说,任意一对相邻的S[i]和S[i+1](i=1, ..., M-1)之间都不存在其他经停站点。我们称相邻站点之间的线路为一个运营区间,每个运营区间只承包给一家公司。环线是有可能存在的,但不会不经停任何中间站点就从出发地回到出发地。当然,不同公司的线路是可能在某些站点有交叉的,这些站点就是客户的换乘点,我们假设任意换乘点涉及的不同公司的线路都不超过5条。
在描述了联盟线路之后,题目将给出一个正整数K(<= 10),随后K行,每行给出一位客户的需求,即始发地的编号和目的地的编号,中间以一空格分隔。
输出格式:
处理每一位客户的需求。如果没有现成的线路可以使其到达目的地,就在一行中输出“Sorry, no line is available.”;如果目的地可达,则首先在一行中输出最顺畅路线的经停站数量(始发地和目的地不包括在内),然后按下列格式给出旅行路线:
Go by the line of company #X1 from S1 to S2. Go by the line of company #X2 from S2 to S3. ......
其中Xi是线路承包公司的编号,Si是经停站的编号。但必须只输出始发地、换乘点和目的地,不能输出中间的经停站。题目保证满足要求的路线是唯一的。
输入样例:4 7 1001 3212 1003 1204 1005 1306 7797 9 9988 2333 1204 2006 2005 2004 2003 2302 2001 13 3011 3812 3013 3001 1306 3003 2333 3066 3212 3008 2302 3010 3011 4 6666 8432 4011 1306 4 3011 3013 6666 2001 2004 3001 2222 6666输出样例:
2 Go by the line of company #3 from 3011 to 3013. 10 Go by the line of company #4 from 6666 to 1306. Go by the line of company #3 from 1306 to 2302. Go by the line of company #2 from 2302 to 2001. 6 Go by the line of company #2 from 2004 to 1204. Go by the line of company #1 from 1204 to 1306. Go by the line of company #3 from 1306 to 3001. Sorry, no line is available.
解题思路:最短路
#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <algorithm>
#include <stack>
#include <queue>
#include <climits>
#include <functional>
#include <vector>
#include <map>
#include <set>
using namespace std;
#define LL long long
const int INF=0x3f3f3f3f;
int a[15],b[1005];
int n;
int s[100005],e[1000005],nt[1000005],l[1000005],id[1000005];
int visit[100005],p[100005],sum[100005],pp[100005],dis[100005];
struct node
{
int id;
int dis,f;
friend bool operator <(node a,node b)
{
return a.dis>b.dis;
}
};
void Dijkstra(int ss)
{
memset(visit,0,sizeof visit);
memset(dis,INF,sizeof dis);
memset(p,-1,sizeof p);
memset(pp,-1,sizeof pp);
memset(sum,INF,sizeof sum);
dis[ss]=0,sum[ss]=0;
node pre,nt1;
pre.id=ss,pre.f=-1,pre.dis=0;
priority_queue<node>q;
q.push(pre);
while(!q.empty())
{
pre=q.top();
q.pop();
visit[pre.id]=1;
for(int i=s[pre.id];~i;i=nt[i])
{
int ee=e[i];
if(visit[ee]) continue;
if(dis[ee]>dis[pre.id]+l[i])
{
dis[ee]=dis[pre.id]+l[i];
if(pre.f!=id[i]&&pre.f!=-1) sum[ee]=sum[pre.id]+1;
else sum[ee]=sum[pre.id];
p[ee]=pre.id;
pp[ee]=id[i];
nt1.dis=dis[ee];
nt1.f=id[i];
nt1.id=ee;
q.push(nt1);
}
else
{
if(pre.f!=id[i]&&pre.f!=-1&&sum[ee]>sum[pre.id]+1)
{
sum[ee]=sum[pre.id]+1;
p[ee]=pre.id;
pp[ee]=id[i];
}
}
}
}
}
int main()
{
scanf("%d",&n);
int cnt=1;
memset(s,-1,sizeof s);
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
for(int j=1;j<=a[i];j++) scanf("%d",&b[j]);
for(int j=1;j<=a[i];j++)
{
for(int k=1;k<=a[i];k++)
{
if(b[j]!=b[k])
nt[cnt]=s[b[j]],s[b[j]]=cnt,l[cnt]=abs(j-k),id[cnt]=i,e[cnt++]=b[k];
}
}
}
int k;
scanf("%d",&k);
while(k--)
{
int ss,ee;
scanf("%d%d",&ss,&ee);
Dijkstra(ss);
if(dis[ee]==INF) {printf("Sorry, no line is available.\n");continue;}
printf("%d\n",dis[ee]);
int kk=ee;
stack<int>q;
stack<int>qq;
while(kk!=ss) q.push(kk),qq.push(pp[kk]),kk=p[kk];
kk=ss;
while(!q.empty())
{
printf("Go by the line of company #%d from %04d to %04d.\n",qq.top(),kk,q.top());
kk=q.top();q.pop();qq.pop();
}
}
return 0;
}