USACO5.4 周游加拿大 Canada Tour

USACO5.4 周游加拿大 Canada Tour

题目

• 描述
  你赢得了一场航空公司举办的比赛，奖品是一张加拿大环游机票。旅行在这家航空公司开放的最西边的城市开始，然后一直自西向东旅行，直到你到达最东边的城市，再由东向西返回，直到你回到开始的城市。除了旅行开始的城市之外，每个城市只能访问一次，因为开始的城市必定要被访问两次（在旅行的开始和结束）。
  当然不允许使用其他公司的航线或者用其他的交通工具。
  给出这个航空公司开放的城市的列表，和两两城市之间的直达航线列表。找出能够访问尽可能多的城市的路线，这条路线必须满足上述条件，也就是从列表中的第一个城市开始旅行，访问到列表中最后一个城市之后再返回第一个城市。
• INPUT FORMAT
  Line 1: 航空公司开放的城市数 N 和将要列出的直达航线的数量 V。N 是一个不大于 100 的正整数。V 是任意的正整数。
  Lines 2..N+1: 每行包括一个航空公司开放的城市名称。城市名称按照自西向东排列。不会出现两个城市在同一条经线上的情况。每个城市的名称都 是一个字符串，最多15字节，由拉丁字母表上的字母组成；城市名称中没有空格。
  Lines N+2..N+2+V-1: 每行包括两个城市名称（由上面列表中的城市名称组成），用一个空格分开。这样就表示两个城市之间的直达双程航线。
• OUTPUT FORMAT
  Line 1: 按照最佳路线访问的不同城市的数量 M。如果无法找到路线，输出 1。
• SAMPLE INPUT (file tour.in)
  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
• SAMPLE OUTPUT (file tour.out)
  7
  也就是: Vancouver, Edmonton, Montreal, Halifax, Toronto, Winnipeg, Calgary, 和 Vancouver （回到开始城市，但是不算在不同城市之内）。

---

动态规划
把返回的路线反向，那么整条路线就变成了两条不相交的从起点到终点的路线。这样我们可以用DP解决。
状态设定：f[i，j] 为假定的甲乙两人，甲走到第i个城市，乙走到第j个城市时，两人走过的城市数目的和。
初始状态：f[1，1]=1。
状态转移方程：f[j，i]=f[i，j]=max{f[i，k]+1}(k到j存在飞机航线，以及f[i，k]>0，就是说存在f[i，k]的走法，1<=k

var n,m,i,j,ans,k:longint;
    c:char;
    s:string;
    city:array[1..100] of string;
    e:array[1..100,1..100] of boolean;
    f:array[1..100,1..100] of longint;

procedure getdata;
var i,j,x,y:longint;
begin
readln(n,m);
for i:=1 to n do readln(city[i]);
for i:=1 to m do
    begin
    s:='';
    read(c);
    while c<>' ' do
        begin
        s:=s+c;
        read(c);
        end;
    for j:=1 to n do
        if city[j]=s then begin x:=j;break;end;
    readln(s);
    for j:=1 to n do
        if city[j]=s then begin y:=j;break;end;
    e[x,y]:=true;
    e[y,x]:=true;
    end;
end;

begin
assign(input,'tour.in');reset(input);
assign(output,'tour.out');rewrite(output);
getdata;
f[1,1]:=1;
for i:=1 to n-1 do
    for j:=i+1 to n do
        begin
        for k:=1 to j-1 do
            if e[k,j] and (f[i,k]>0) and (f[i,k]+1>f[i,j]) then f[i,j]:=f[i,k]+1;
        f[j,i]:=f[i,j];
        end;
ans:=1;
for i:=1 to n-1 do
    if (e[i,n]) and (f[i,n]>ans) then ans:=f[i,n];
writeln(ans);
close(input);
close(output);
end.