泽泽在巴西（jzoj普及组模拟第三题）（超难）

题目：

题目描述

泽泽帮助了英国某街道尽量减少酸雨的伤害，街道办主任非常感激他，就把他领到一扇门前，告诉他这扇门能通往好地方，具体好到什么程度要看泽泽人品。泽泽毫不犹豫地走了进去……

泽泽来到了足球王国——巴西。这可是个好地方，泽泽看来人品攒了不少了。这里大街小巷都在踢足球，其乐无穷。

突然，泽泽被一个人拎了起来，一看，是个足球流氓。他后面跟了一大群足球流氓，正虎视眈眈地看他。他们要求和泽泽比赛，输了就要揍他。

没办法，泽泽硬着头皮和足球流氓另外掳来的几个人一起组建了一只队伍，和足球流氓队比赛。

比赛开始，泽泽队率先发球。泽泽观察了四周，想怎么才能用最短的时间射门呢？

射门的时间为距离*2，而传球的时间是距离*1。所以泽泽想找一条用时最少的射门路径，来打败足球流氓。

足球流氓当然不会袖手旁观，他们会拦截。当泽泽队伍中的传球人、被传球人之间有某足球流氓并且他们在同一直线上时，传球不会成功，即不能这样传球。比如A（1，2）想传球给B（7，8），中间有个足球流氓C（3，4），则他们在同一直线，传球不成功。射门不受足球流氓影响。

输入

第1行4个整数x0，y0，n，m。x0，y0表示球门的坐标，n表示泽泽队伍有几个人，m表示足球流氓有几个人。

接下来的n行，分别有2个整数，表示泽泽球队的球员坐标。其中最前面的2个整数是泽泽的坐标。球一开始在泽泽脚下。

接下来的m行，分别有2个整数，表示足球流氓的球员坐标。

保证不会有2个人坐标相同。

输出

输出一个整数，是最短时间四舍五入取整后的结果。

样例输入

0 0 5 2

20 20

27 -14

0 16

-7 -9

23 38

22 24

3 0

样例输出

52

数据范围限制

提示

【样例说明】

泽泽（20，20）传给3号队员（0，16），3号队员再射门（0，0），总共用了52。

【限制】

对于80%的数据，n<=10，m<=5

对于100%的数据，n<=300，m<=100

作者思路：floyd算法。
这道题难点是预处理，需要求出斜率，用勾股求长度，再用Floyd。

花了我一个半小时啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。

代码：

var a:array[0..401,0..401] of longint;
    f:array[0..401,0..401] of real;
    n,m:longint;
function min(x,y:real):real;
begin
  if x<y then exit(x) else exit(y);
end;
function bool(x1,y1,x2,y2,x3,y3:longint):boolean;
begin
  if (x1-x3)*(y2-y3)-(x2-x3)*(y1-y3)<>0 then exit(true);
  if (sqrt(sqr(x1-x2)+sqr(y1-y2))<sqrt(sqr(x1-x3)+sqr(y1-y3)))
  or(sqrt(sqr(x1-x2)+sqr(y1-y2))<sqrt(sqr(x2-x3)+sqr(y2-y3)))
  then exit(true);
  exit(false);
end;
procedure init;
var i,j,k:longint;
begin
  read(a[0,1],a[0,2],n,m);
  fillchar(f,sizeof(f),$7f {div 10});
  for i:=1 to n do read(a[i,1],a[i,2]);
  for i:=1+n to n+m do read(a[i,1],a[i,2]);
  for i:=1 to n do
    for j:=i+1 to n do
      for k:=n+1 to m+n do
        if (i<>j)and(i<>k)and(j<>k)and
        bool(a[i,1],a[i,2],a[j,1],a[j,2],a[k,1],a[k,2]) then
        begin
          f[i,j]:=sqrt(sqr(a[i,1]-a[j,1])+sqr(a[i,2]-a[j,2]));
          f[j,i]:=f[i,j];
        end;
  for i:=1 to n do
  begin
    f[i,0]:=sqrt(sqr(a[i,1]-a[0,1])+sqr(a[i,2]-a[0,2]))*2;
    f[0,i]:=f[i,0];
  end;
end;
procedure floyd;
var i,j,k:longint;
begin
  for i:=0 to n do
    for j:=0 to n do
      for k:=0 to n do
        f[j,k]:=min(f[j,k],f[j,i]+f[i,k]);
  write(f[0,1]:0:0);
end;
begin
  {assign(input,'brazil.in');
  assign(output,'brazil.out');
  reset(input);
  rewrite(output);     }
  init;
  floyd;
  {close(input);
  close(output);      }
end.