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.

wiki 2490 导弹拦截塔的更多相关文章

  1. Codevs2490 导弹防御塔

    2490 导弹防御塔 2490 导弹防御塔 时间限制: 1 s 空间限制: 64000 KB 题目等级 : 大师 Master         题目描述 Description Freda的城堡—— ...

  2. wiki oi 1044 拦截导弹

    题目描述 Description 某国为了防御敌国的导弹袭击,发展出一种导弹拦截系统.但是这种导弹拦截系统有一个缺陷:虽然它的第一发炮弹能够到达任意的高度,但是以后每一发炮弹都不能高于前一发的高度.某 ...

  3. HDU-1257 导弹拦截系统 http://acm.hdu.edu.cn/showproblem.php?pid=1257

    Problem Description 某国为了防御敌国的导弹袭击,发展出一种导弹拦截系统.但是这种导弹拦截系统有一个缺陷:虽然它的第一发炮弹能够到达任意的高度,但是以后每一发炮弹都不能超过前一发的高 ...

  4. AC日记——导弹拦截 洛谷 P1020 (dp+模拟)

    题目描述 某国为了防御敌国的导弹袭击,发展出一种导弹拦截系统.但是这种导弹拦截系统有一个缺陷:虽然它的第一发炮弹能够到达任意的高度,但是以后每一发炮弹都不能高于前一发的高度.某天,雷达捕捉到敌国的导弹 ...

  5. 【codevs1044】导弹拦截问题与Dilworth定理

    题目描述 Description 某国为了防御敌国的导弹袭击,发展出一种导弹拦截系统.但是这种导弹拦截系统有一个缺陷:虽然它的第一发炮弹能够到达任意的高度,但是以后每一发炮弹都不能高于前一发的高度.某 ...

  6. TYVJ P1020 导弹拦截 Label:水

    题目描述 某国为了防御敌国的导弹袭击,发展出一种导弹拦截系统.但是这种导弹拦截系统有一个缺陷:虽然它的第一发炮弹能够到达任意的高度,但是以后每一发炮弹都不能高于前一发的高度.某天,雷达捕捉到敌国的导弹 ...

  7. bzoj2044: 三维导弹拦截

    Description 一场战争正在A国与B国之间如火如荼的展开. B国凭借其强大的经济实力开发出了无数的远程攻击导弹,B国的领导人希望,通过这些导弹直接毁灭A国的指挥部,从而取得战斗的胜利!当然,A ...

  8. bzoj3035: 导弹防御塔

    Description Freda的城堡——“Freda,城堡外发现了一些入侵者!”“喵...刚刚探究完了城堡建设的方案数,我要歇一会儿嘛lala~”“可是入侵者已经接近城堡了呀!”“别担心,rain ...

  9. nyoj 79 导弹拦截

    点击打开链接 拦截导弹 时间限制:3000 ms  |  内存限制:65535 KB 难度:3 描述 某国为了防御敌国的导弹袭击,发展中一种导弹拦截系统.但是这种导弹拦截系统有一个缺陷:虽然它的第一发 ...

随机推荐

  1. 问题:调用 ASP.Net Core WebAPI的HTTP POST方法时,从 [FromBody] 中读取的 MongoDB GeoJsonObjectModel成员总是null

    问题描述: POST/PUT to ASP.Net Core with [FromBody] to a MongoDB GeoJsonObjectModel member is always null ...

  2. html5特效库

    Swiper是纯javascript打造的滑动特效插件,面向手机.平板电脑等移动终端.能实现触屏焦点图.触屏Tab切换.触屏多图切换等常用效果. delaunay.js是一款能在图片中形成无数个三角形 ...

  3. 你可能会用到的"奇技赢巧"

    工作中偶尔会遇到一些不常见的问题,但是解决起来又极其麻烦,通常要找很多资料才能搞定,这里我总结了近段时间遇到的一些问题,可能会对你有帮助,高手勿喷. 1.关于iPhone最下面会弹出奇怪框框的问题 就 ...

  4. 第一篇 Charles的配置及相关使用

    // Charles Proxy License // 适用于Charles任意版本的注册码,谁还会想要使用破解版呢. // Charles 4.2目前是最新版,可用.   Registered Na ...

  5. Tomcat突然用开发工具启动不起来,只报了个红色的警告,没有其他任何异常

    碰到这个问题是,是因为我的catalina.bat文件做了配置修改,导致与工具这边的启动设置起了冲突 下面这个是我在Catalina.bat中新增的配置,删掉这个就可以了 set JAVA_OPTS= ...

  6. LeetCode 4——两个排序数组中的中位数

    1. 题目 2. 解答 2.1. 方法一 由于两个数组都是排好序的,因此首先可以想到的思路就是利用归并排序把两个数组合并成一个有序的长数组,然后直接取出中位数即可. class Solution: d ...

  7. ISAP 最大流 最小割 模板

    虽然这道题用最小割没有做出来,但是这个板子还是很棒: #include<stdio.h> #include<math.h> #include<string.h> # ...

  8. Java 实现一个带提醒的定时器

    定时闹钟预览版EXE下载链接:https://files.cnblogs.com/files/rekent/ReadytoRelax_jar.zip 功能说明: 实现了一个休息提醒器,用户首先设定一个 ...

  9. lintcode-47-主元素 II

    47-主元素 II 给定一个整型数组,找到主元素,它在数组中的出现次数严格大于数组元素个数的三分之一. 注意事项 数组中只有唯一的主元素 样例 给出数组[1,2,1,2,1,3,3] 返回 1 挑战 ...

  10. 算法(7)Majority Element II

    题目:找出数组中出现次数大于n/3次的数字 思路:摩尔投票法.所有的帖子中都说:先遍历一遍数组找到备选元素,然后再遍历一遍数组考察下这个元素是否是真的超过n/3,然后就直接上代码,但是现在的问题是:我 ...