bzoj 3035 二分答案+二分图最大匹配

　　大原来做的一道题，偷懒直接粘的原来的程序 

　　http://www.cnblogs.com/BLADEVIL/p/3433520.html

/**************************************************************
    Problem: 3035
    User: BLADEVIL
    Language: Pascal
    Result: Accepted
    Time:240 ms
    Memory:4216 kb
****************************************************************/
 
//By BLADEVIL
 
var
 
    n, m, t2, v                     :longint;
 
    t1                              :real;
 
    dis                             :array[0..100,0..100] of real;
 
    ans                             :real;
 
    pre, other, last                :array[0..200100] of longint;
 
    link                            :array[0..100] of longint;
 
    flag                            :array[0..100] of boolean;
 
    x1, x2, y1, y2                  :array[0..100] of longint;
 
    l                               :longint;
 
    len                             :array[0..200100] of real;
 
  
 
procedure init;
 
var
 
    i, j                            :longint;
 
begin
 
    read(n,m,t1,t2,v);
 
    for i:=1 to m do read(x2[i],y2[i]);
 
    for i:=1 to n do read(x1[i],y1[i]);
 
    t1:=t1/60;
 
    for i:=1 to n do
 
        for j:=1 to m do
 
            dis[i,j]:=sqrt((x1[i]-x2[j])*(x1[i]-x2[j])+(y1[i]-y2[j])*(y1[i]-y2[j]));
 
end;
 
  
 
procedure connect(x,y:longint; z:real);
 
begin
 
    inc(l);
 
    pre[l]:=last[x];
 
    last[x]:=l;
 
    other[l]:=y;
 
    len[l]:=z;
 
end;
 
  
 
function find(i:longint):boolean;
 
var
 
    q, p                            :longint;
 
begin
 
    q:=last[i];
 
    while q<>0 do
 
    begin
 
        p:=other[q];
 
        if (not flag[p]) then
 
        begin
 
            flag[p]:=true;
 
            if (link[p]=0) or (find(link[p])) then
 
            begin
 
                link[p]:=i;
 
                exit(true);
 
            end;
 
        end;
 
        q:=pre[q];
 
    end;
 
    exit(false);
 
end;
 
  
 
procedure judge(low,high:real);
 
var
 
    mid                             :real;
 
    tot                             :longint;
 
    i, j, ll                        :longint;
 
    k                               :longint;
 
    count                           :longint;
 
begin
 
    if high<low then exit;
 
    fillchar(last,sizeof(last),0);
 
    fillchar(link,sizeof(link),0);
 
    tot:=0;
 
    l:=1;
 
    mid:=(low+high)/2;
 
    k:=trunc(((mid-t1)/(t1+t2))+1);
 
    for i:=1 to n do
 
        for ll:=1 to k do
 
        begin
 
            inc(tot);
 
            for j:=1 to m do
 
                if (((mid-t1)-(ll-1)*(t1+t2))*v)>=dis[i,j]
 
                    then connect(tot,j,dis[i,j]/v+(ll-1)*(t1+t2)+t1);
 
        end;
 
    count:=0;
 
  
 
    for i:=1 to tot do
 
    begin
 
        fillchar(flag,sizeof(flag),false);
 
        if find(i) then inc(count);
 
    end;
 
  
 
    if (high-low<1e-8) then
 
        if count>=m then
 
        begin
 
            ans:=mid;
 
        end else exit;
 
    if count>=m then
 
    begin
 
        ans:=mid;
 
        judge(low,mid-1e-8);
 
    end else judge(mid+1e-8,high);
 
end;
 
  
 
procedure main;
 
begin
 
    judge(1,30000);
 
    writeln(ans:0:6);
 
end;
 
  
 
begin
 
    init;
 
    main;
 
end.