泽泽在巴西

Description

　　泽泽帮助了英国某街道尽量减少酸雨的伤害，街道办主任非常感激他，就把他领到一扇门前，告诉他这扇门能通往好地方，具体好到什么程度要看泽泽人品。泽泽毫不犹豫地走了进去……
　　泽泽来到了足球王国——巴西。这可是个好地方，泽泽看来人品攒了不少了。这里大街小巷都在踢足球，其乐无穷。
突然，泽泽被一个人拎了起来，一看，是个足球流氓。他后面跟了一大群足球流氓，正虎视眈眈地看他。他们要求和泽泽比赛，输了就要揍他。
　　没办法，泽泽硬着头皮和足球流氓另外掳来的几个人一起组建了一只队伍，和足球流氓队比赛。
　　比赛开始，泽泽队率先发球。泽泽观察了四周，想怎么才能用最短的时间射门呢？
　　射门的时间为距离*2，而传球的时间是距离*1。所以泽泽想找一条用时最少的射门路径，来打败足球流氓。
　　足球流氓当然不会袖手旁观，他们会拦截。当泽泽队伍中的传球人、被传球人之间有某足球流氓并且他们在同一直线上时，传球不会成功，即不能这样传球。比如A（1，2）想传球给B（7，8），中间有个足球流氓C（3，4），则他们在同一直线，传球不成功。射门不受足球流氓影响。

Input

第1行3个整数x0，y0，n，m。x0，y0表示球门的坐标，n表示泽泽队伍有几个人，m表示足球流氓有几个人。
接下来的n行，分别有2个整数，表示泽泽球队的球员坐标。其中最前面的2个整数是泽泽的坐标。球一开始在泽泽脚下。
接下来的m行，分别有2个整数，表示足球流氓的球员坐标。
保证不会有2个人坐标相同。

Output

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

Sample Input

0 0 5 2
20 20
27 -14
0 16
-7 -9
23 38
22 24
3 0
Sample Output

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

数据范围
对于80%的数据，n<=10，m<=5
对于100%的数据，n<=300，m<=100
做法：最短路加计算机几何，先是用叉积公式判断是否在线段上，然后跑一遍最短路就好了。
代码如下：

var
 a,b:array[1..3000,1..3] of longint;
 x,y,ls,next,list,v:array[1..300000] of longint;
 w,d:array[1..300000]of real;
 i,j,n,m,q,p,r:longint;

function check(x1,y1,x2,y2,x,y:longint):boolean;
begin
 if (x>x1)and(x>x2)or(y>y1)and(y>y2)or(x<x1)and(x<x2)or(y<y1)and(y<y2)
        then exit(false)
        else exit(true);
end;

function cj(x1,y1,x2,y2,x,y:longint):longint;
var
 total:longint;
begin
 total:=(x1-x)*(y2-y)-(x2-x)*(y1-y);
 exit(total);
end;


procedure init;
var
 i,j,k:longint;
 f:boolean;
begin
 for i:=1 to n do
  for j:=i+1 to n+1 do
   begin
    f:=true;
    for k:=1 to m do
     if (cj(a[j,1],a[j,2],b[k,1],b[k,2],a[i,1],a[i,2])=0)and(check(a[i,1],a[i,2],a[j,1],a[j,2],b[k,1],b[k,2]))
      then begin
            f:=false;
            break;
           end
      else continue;
     if f then
          begin
           inc(r);
           x[r]:=a[i,3];
           y[r]:=a[j,3];
           if j<>n+1
            then w[r]:=sqrt(sqr(a[i,1]-a[j,1])+sqr(a[i,2]-a[j,2]))
            else w[r]:=2*sqrt(sqr(a[i,1]-a[j,1])+sqr(a[i,2]-a[j,2]));
           next[r]:=ls[x[r]];
           ls[x[r]]:=r;
           inc(r);
           x[r]:=y[r-1];
           y[r]:=x[r-1];
           w[r]:=w[r-1];
           next[r]:=ls[x[r]];
           ls[x[r]]:=r;
          end;
    end;
end;

procedure spfa;
var
 head,tail,t:longint;
begin
 head:=0; tail:=1;
 list[1]:=1;
 v[1]:=1;
 while head<>tail do
  begin
   inc(head);
   t:=ls[list[head]];
   while t>0 do
    begin
     if d[x[t]]+w[t]<d[y[t]]
      then begin
            d[y[t]]:=d[x[t]]+w[t];
            if v[y[t]]=0
             then begin
                   v[y[t]]:=1;
                   inc(tail);
                   list[tail]:=y[t];
                  end;
           end;
      t:=next[t];
     end;
   v[list[head]]:=0;
  end;
end;

begin
 readln(q,p,n,m);
 for i:=1 to n do
  begin
   read(a[i,1],a[i,2]);
   a[i,3]:=i;
  end;
 a[i+1,1]:=q; a[i+1,2]:=p; a[i+1,3]:=i+1;
 for i:=1 to m do
  readln(b[i,1],b[i,2]);
 init;
 for i:=1 to n+1 do
  d[i]:=maxlongint;
 d[1]:=0;
 spfa;
 writeln(d[n+1]:0:0);
end.