青蛙的烦恼 (frog.pas/c/cpp)

 青蛙的烦恼 (frog.pas/c/cpp) 

 
池塘中有 n片荷叶恰好围成了一个凸多边形，有一只小青蛙恰好站在1号荷叶上，小青
蛙想通过最短的路程遍历所有的荷叶（经过一个荷叶一次且仅一次），小青蛙可以从一片荷
叶上跳到另外任意一片荷叶上。
 
输入数据(frog.in)
 
  第一行为整数n，荷叶的数量。
接下来 n 行，每行两个实数，为n 个多边形的顶点坐标，按照顺时针方向给出。保证不
会爆double。
 
输出数据(frog.out)：
遍历所有荷叶最短路程，请保留3位小数。
 
样例
输入:frog.in
4
50.0 1.0
5.0 1.0
0.0 0.0
45.0 0.0
 
输出:frog.out
50.211
 
数据范围：
对于所有数据,0<n<=720

=============================

写在前面

值得鼓励的是第二次写出了区间动规的转移方程

---------------------------------------------------------------

------------------------------------------------------------------------

具体方程看代码...

=============================

var
  n:longint;
  x,y:array[1..1440]of double;
  dis:array[1..1440,1..1440]of double;
  f:array[1..721,1..721,0..1]of double;
  //ans:double;
procedure init;
begin
  assign(input,'frog.in');
  assign(output,'frog.out');
  reset(input); rewrite(output);
end;

procedure terminate;
begin
  close(input); close(output);
  halt;
end;

function min(a,b:double):double;
begin
  if a>b then exit(b);
  exit(a);
end;

procedure main;
var
  i,j:longint;
  l,s:longint;
begin
  readln(n);
  for i:=1 to n do
    begin
      readln(x[i],y[i]);
      f[i,1,0]:=0;
      f[i,1,1]:=0;
    end;
    
  for i:=1 to 2*n do
    for j:=1 to 2*n do
      dis[i,j]:=sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
      
  for l:=2 to n do
    for s:=1 to n do
      begin
        f[s,l,0]:=min(dis[s,s+1]+f[s+1,l-1,0],
                      dis[s,s+l-1]+f[s+1,l-1,1]);
        f[s,l,1]:=min(dis[s+l-1,s+l-2]+f[s,l-1,1],
                      dis[s,s+l-1]+f[s,l-1,0]);
      end;
  writeln(f[1,n,0]:0:3);
end;

begin
  init;
  main;
  terminate;
end.