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. jmeter常用的内置变量

    1. vars   API:http://jmeter.apache.org/api/org/apache/jmeter/threads/JMeterVariables.html vars.get(& ...

  2. jmeter操作JDBC

    1. 依次添加计划.线程组.JDBC Connection Configuration.JDBC Request.HTTP请求.Debug Sampler.察看结果树 在计划中导入mysql的jdbc ...

  3. 第三十三篇 Python中关于OOP(面向对象)的常用术语

    面向对象的优点 从编程进化论可知,面向对象是一种更高等级的结构化编程方式,它的好处主要有两点: 1. 通过封装明确了内外,你做为类的缔造者,你就是女娲,女娲造物的逻辑别人无需知道,女娲想让你知道,你才 ...

  4. 基于语音转录的ted演讲推荐

    论文地址:https://arxiv.org/abs/1809.05350v1 二.  实现 我们从Kaggle[6]中获取了TED演讲数据集,其中包括2400个TED演讲的数据,包括标题.演讲者.标 ...

  5. LeetCode 61——旋转链表

    1. 题目 2. 解答 2.1. 方法一 将链表每个节点向右移动 1 个位置,其实就是让链表最后一个结点指向第一个结点. 因此,向右移动 k 个位置就重复上述过程 k 次即可. 然后,我们注意到,若链 ...

  6. PAT 1030 完美数列

    https://pintia.cn/problem-sets/994805260223102976/problems/994805291311284224 给定一个正整数数列,和正整数 p,设这个数列 ...

  7. mysql 数据包太小会引发错误信息

    Error querying database.  Cause: com.mysql.cj.jdbc.exceptions.PacketTooBigException: Packet for quer ...

  8. Leader Election

    Leader Election Zookeeper的基本操作 Zookeeper虽然是分布式系统,但它并不是为文件存储而设计的,Zookeeper里存储的一般是配置信息和源信息.实际上,Zookeep ...

  9. Bootstrap中的Affix插件

    我们为什么要用bootstrap?因为懒!哦....不,是因为方便,呃...意思差不多. 今天来说说Affix这个插件,它可以使导航栏固定,免去了自己手写的麻烦,用着非常方便,废话不多说,下面是用法. ...

  10. mac 倍速播放

    Mac的倍速播放,使用工具播放是quick time player,使用方式是,按住 “option” + “>>” 即可调整倍数: