小萨的烦恼(代理解，待改进。。。)

小萨的烦恼

Time Limit:10000MS  Memory Limit:65536K
Total Submit:195 Accepted:66
Case Time Limit:1000MS

Description

　　圣诞节又要到了！小萨希望和自己心仪的MM一起出去度过一个浪漫的的圣诞节。他进行了详尽的准备，找到了N个很适合他们去约会的好地方，但无奈小萨和他的MM都是初三学校，晚上必须回学校上晚修，没有足够的时间让他们走遍每一个地方。迫于无奈，小萨选择了一个最合适的地方S。小萨打算和他的MM一起步行到那个地方，那样他们就有很多时间来聊天玩乐，但是无奈时间有限，小萨决定坐公共汽车过去。然而，小萨发觉，自己的钱只够买一张公共汽车车票了！小萨该怎么办呢？他找到了你……
【题目描述】
　　小萨将给出一张地图，它可以看做一个有N（N<=100）个结点的图。这张地图有N个公共汽车站，小萨只可以在这N个汽车站上公共汽车。有些车站之间存在一条双向通路，无论是公共汽车还是小萨，都只能走这些通路。若两个车站之间的距离为d，步行所需要的时间为2*d秒，坐公共汽车所需要的时间为d秒。小萨他们只有T秒的时间，他希望你求出他和他的MM出去约会后能否及时赶回来，若不能，则输出“You are day dreaming!”，否则输出他们所需要花费的最少时间。
注意：一张公共汽车车票可以使用两次（只可以搭乘同一辆公共汽车），即可以认为小萨去约会地点的时间和回来的时间相等。

Input

第一行有三个整数N、T、S
接下来是一个N*N的邻接矩阵。两两车站之间的距离不超过10^9。A[I,j]为0表示I和J车站不连通。

Output

若小萨他们不能在限定时间内赶回学校，那么输出“You are day dreaming!”（不包括引号）
否则，输出一个整数，为他们所需要花费的最少时间。

Sample Input

4 5 4
0 1 1 1
1 0 1 1
1 1 0 1
1 1 1 0

Sample Output

2

Hint

【数据范围】
20%的数据：N<=10
100%的数据：N<=70 T<=10^9

Source

fhn

• var
    i,j,k,n,m:longint;
    d:array[-1..100] of longint;
    v:array[-1..100] of boolean;
    a:array[-1..100,-1..100] of longint;
  
  procedure dij;
  var
    min,u,i:longint;
  begin
    repeat
      min:=maxlongint div 3;
      u:=0;
      for i:=1 to n do
        if (v[i])and(d[i]<min) then
        begin
          min:=d[i];
          u:=i;
        end;
      if u>0 then
      begin
        v[u]:=false;
        for i:=1 to n do
          if (v[i])and(d[u]+a[u,i]<d[i]) then
            d[i]:=d[u]+a[u,i];
      end;
    until u=0;
  end;
  
  procedure main;
  var
    t,s,i,j,k,min:longint;
  begin
    readln(n,t,s);
    for i:=1 to n do
      for j:=1 to n do
      begin
        read(a[i,j]);
        a[i,j]:=a[i,j]*2;
        if a[i,j]=0 then a[i,j]:=maxlongint div 3;
      end;
    min:=maxlongint div 3;
    for i:=1 to n-1 do
      for j:=i+1 to n do
      begin
        if a[i,j]=maxlongint div 3 then continue;
        a[i,j]:=a[i,j] div 2;
        a[j,i]:=a[j,i] div 2;
        fillchar(v,sizeof(v),true);
        for k:=2 to n do
          d[k]:=maxlongint div 3;
        dij;
        if d[s]<min then min:=d[s];
        a[i,j]:=a[i,j]*2;
        a[j,i]:=a[j,i]*2;
      end;
    if min*2<=t
      then writeln(min*2)
      else writeln('You are day dreaming!');
  end;
  
  begin
    main;
  end.