码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj 3035 pascal

时间:2014-07-12 13:32:17      阅读:220      评论:0      收藏:0      [点我收藏+]

标签:cti   for   io   代码   算法   line   

这道题也是一个二分图匹配的题,也需要用匈牙利算法,但需要小小的改进

 

中间需要一个二分搜答案的过程,这也是其中的一个小小的变型吧

 

需要注意一下读入,他是有分钟有秒,所以不要被坑哦

 

最后不说废话,直接上代码

 

/**************************************************************
    Problem: 3035
    User: Davidliu
    Language: Pascal
    Result: Accepted
    Time:236 ms
    Memory:4216 kb
****************************************************************/
 
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.

bzoj 3035 pascal,布布扣,bubuko.com

bzoj 3035 pascal

标签:cti   for   io   代码   算法   line   

原文地址:http://www.cnblogs.com/Davidlxy/p/3836082.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!