泽泽在巴西

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

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

代码

用计算几何算出哪些边是可以连的
再用spfa就好了

代码

const
  maxe=300000;
  maxv=10000000;

type
  arr=record
    x,y,next:longint;
    w:real;
  end;

var
  n,m,s,q:longint;
  ls:array[0..maxe] of longint;
  edge:array[0..maxe,1..2] of longint;
  a:array[0..maxv] of arr;
  f:array[0..maxe] of real;
  v:array[0..maxe] of longint;
  d:array[0..maxe] of longint;
  i,j,k:longint;
  ans:longint;
  x1,x2:longint;
  y1,y2:longint;

procedure spfa;
var
  i,j,k:longint;
  head,tail:longint;
begin
  fillchar(f,sizeof(f),$7f);
  head:=0;
  tail:=1;
  v[s]:=1;
  d[1]:=s;
  f[s]:=0;
  repeat
    head:=head+1;
    j:=ls[d[head]];
    while j<>0 do
      begin
        with a[j] do
          begin
            if f[x]+w<f[y]
              then
                begin
                  f[y]:=f[x]+w;
                  if v[y]=0
                    then
                      begin
                        tail:=tail+1;
                        d[tail]:=y;
                        v[y]:=1;
                      end;
                end;
            j:=next;
          end;
      end;
    v[d[head]]:=0;
  until head=tail;
end;

function pan(x1,y1,x2,y2,x3,y3:int64):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 add(x,y:longint;w:real);
begin
  q:=q+1;
  a[q].x:=x;
  a[q].y:=y;
  a[q].w:=w;
  a[q].next:=ls[x];
  ls[x]:=q;
end;

begin
  assign(input,'brazil.in');
  assign(output,'brazil.out');
  rewrite(output);
  reset(input);
  readln(edge[0,1],edge[0,2],n,m);
  q:=0;
  fillchar(ls,sizeof(ls),0);
  for i:=1 to n do
    readln(edge[i,1],edge[i,2]);
  for i:=n+1 to m+n do
    readln(edge[i,1],edge[i,2]);
  for i:=1 to n do
    begin
      for j:=i+1 to n do
        for k:=n+1 to m+n do
          if pan(edge[i,1],edge[i,2],edge[j,1],edge[j,2],edge[k,1],edge[k,2])
            then
              begin
                writeln(i);
                x1:=edge[i,1]; y1:=edge[i,2];
                x2:=edge[j,1]; y2:=edge[j,2];
                add(i,j,sqrt(sqr(x1-x2)+sqr(y1-y2)));
                add(j,i,sqrt(sqr(x1-x2)+sqr(y1-y2)));
                break;
              end;
    end;
  for i:=1 to n do
    begin
      x1:=edge[i,1]; y1:=edge[i,2];
      x2:=edge[0,1]; y2:=edge[0,2];
      add(i,0,sqrt(sqr(x1-x2)+sqr(y1-y2))*2);
      add(0,i,sqrt(sqr(x1-x2)+sqr(y1-y2))*2);
    end;
  s:=1;
  spfa;
  write(f[0]:0:0);
  close(input);
  close(output);
end.