wiki 2490 导弹拦截塔

2013-09-23 21:16

二分答案+匈牙利判断

对于每一个时间，我们重新建一张二分图，由于每个塔可能打多次，所以要拆点，

对于每个拆的点的可行飞行距离为(mid-t1)-(ll-1)*(t1+t2)*v，其中mid为二分的答案

ll为将当前的点拆成第几个点（因为拆的点的时间是不一样的），然后依次判断该点和

入侵者的距离是否小于，是则加边。

建完图之后判断是否存在完美匹配，存在则向下二分，否则向上二分。

//吐槽下，靠靠，t1的单位是秒，t2的单位是分钟。。。

代码只是过了数据，好多地方都可以优化。

比较弱，二分没有标程写的好，标程可以直接得到最后答案。

//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.

转载于:https://www.cnblogs.com/BLADEVIL/p/3433520.html