题解 \(by\;zj\varphi\)

一道凸包题

对于每个导弹,它的飞行时间就是 \(tim=\frac{A}{a_i}+\frac{B}{b_i}\) 我们设 \(x=\frac{1}{a_i},y=\frac{1}{b_i}\) 那么 \(tim=Ax+By\)

化简后 \(y=-\frac{A}{B}x+\frac{tim}{B}\) 我们要让斜率最小,那么维护一个下凸包,但斜率只能是负的,所以我们只要左下凸包

所以对于 \(A,B\) 排序,按 \(A\) 加入点,最后记得维护左下凸包,不要把超过范围的也算进去

Code
  1. #include<bits/stdc++.h>
  2. #define ri register signed
  3. #define p(i) ++i
  4. using namespace std;
  5. namespace IO{
  6. char buf[1<<21],*p1=buf,*p2=buf;
  7. #define gc() p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++
  8. template<typename T>inline void read(T &x) {
  9. ri f=1;x=0;register char ch=gc();
  10. while(ch<'0'||ch>'9') {if (ch=='-') f=0;ch=gc();}
  11. while(ch>='0'&&ch<='9') {x=(x<<1)+(x<<3)+(ch^48);ch=gc();}
  12. x=f?x:-x;
  13. }
  14. }
  15. using IO::read;
  16. namespace nanfeng{
  17. #define FI FILE *IN
  18. #define FO FILE *OUT
  19. template<typename T>inline T cmax(T x,T y) {return x>y?x:y;}
  20. template<typename T>inline T cmin(T x,T y) {return x>y?y:x;}
  21. typedef double db;
  22. static const int N=3e5+7;
  23. int st[N],tp,pre[N],vis[N],n,mxa,mxb;
  24. db tk[N];
  25. struct node{int a,b,id;}mis[N];
  26. inline int operator<(const node &n1,const node &n2) {return n1.a==n2.a?n1.b>n2.b:n1.a>n2.a;}//每个相同的 a 只有最大的 b 有用
  27. inline db slope(node n1,node n2) {return ((db)n1.a*n2.a*(n2.b-n1.b))/((db)n1.b*n2.b*(n2.a-n1.a));}
  28. inline int main() {
  29. // FI=freopen("nanfeng.in","r",stdin);
  30. // FO=freopen("nanfeng.out","w",stdout);
  31. read(n);
  32. for (ri i(1);i<=n;p(i)) read(mis[i].a),read(mis[i].b),mis[i].id=i;
  33. sort(mis+1,mis+n+1);
  34. for (ri i(1);i<=n;p(i)) if(mxb<mis[i].b) mxb=mis[i].b,mxa=mis[i].a;
  35. st[p(tp)]=1;
  36. for (ri i(2);mxa<=mis[i].a&&i<=n;p(i)) {//mxa 防止右凸包也算进去
  37. if (mis[i].a==mis[st[tp]].a) {
  38. if (mis[i].b==mis[st[tp]].b)
  39. pre[mis[i].id]=pre[mis[st[tp]].id],pre[mis[st[tp]].id]=mis[i].id;
  40. continue;
  41. }
  42. while(tp>1&&slope(mis[st[tp]],mis[i])<tk[st[tp]]) --tp;
  43. tk[i]=slope(mis[st[tp]],mis[i]),st[p(tp)]=i;
  44. }
  45. while(tp)
  46. for (ri i(mis[st[tp--]].id);i;i=pre[i]) vis[i]=1;
  47. for (ri i(1);i<=n;p(i)) if (vis[i]) printf("%d ",i);
  48. return 0;
  49. }
  50. }
  51. int main() {return nanfeng::main();}

NOIP 模拟 $18\; \rm 导弹袭击$的更多相关文章

  1. NOIP 模拟 $18\; \rm 老司机的狂欢$

    题解 \(by\;zj\varphi\) 一道很有趣的题,我用的动态开点线段树和倍增 首先对于第一问,不难想到要二分,二分时间,因为时间长一定不会比时间短能跑的人多 那么如何 check,先将所有老司 ...

  2. NOIP 模拟 $18\; \rm 炼金术士的疑惑$

    题解 \(by\;zj\varphi\) 高斯消元 根据高中化学知识,求解方程的就是一直方程凑出来的,焓值也一样 那么对于要求的方程和一直方程,我们做一次高斯消元,以每个物质为未知数,因为它保证有解, ...

  3. NOIP模拟 1

    NOIP模拟1,到现在时间已经比较长了.. 那天是6.14,今天7.18了 //然鹅我看着最前边缺失的模拟1,还是终于忍不住把它补上,为了保持顺序2345重新发布了一遍.. #   用  户  名   ...

  4. 20190902+0903合集-NOIP模拟

    一直没时间写QwQ 于是补一下. Day 1 晚饭吃的有点恶心…… $1s\,2s\,5s$ 还开 -O2 ?? 有点恐怖. T1 猛的一想: 把外面设成一个点, 向入口连一条权为排队时间的边 从出口 ...

  5. 2021.5.22 noip模拟1

    这场考试考得很烂 连暴力都没打好 只拿了25分,,,,,,,,好好总结 T1序列 A. 序列 题目描述 HZ每周一都要举行升旗仪式,国旗班会站成一整列整齐的向前行进. 郭神作为摄像师想要选取其中一段照 ...

  6. Noip模拟18 2021.7.17 (文化课专场)

    T1 导弹袭击(数学) 显然,我们要找到最优的A,B使得一组a,b优于其他组那么可以列出: $\frac{A}{a_i}+\frac{B}{b_i}<\frac{A}{a_j}+\frac{B} ...

  7. noip提高组1999 导弹拦截

    导弹拦截 背景 实中编程者联盟为了培养技术精湛的后备人才,必须从基础题开始训练. 描述 某国为了防御敌国的导弹袭击,研发出一种导弹拦截系统.但是这种导弹拦截系统有一个缺陷:虽然它的第一发炮弹能够到达任 ...

  8. 软件开发 [CJOJ 1101] [NOIP 模拟]

    Description 一个软件开发公司同时要开发两个软件,并且要同时交付给用户,现在公司为了尽快完成这一任务,将每个软件划分成m个模块,由公司里的技术人员分工完成,每个技术人员完成同一软件的不同模块 ...

  9. NOIP模拟

    1.要选一个{1,2,...n}的子集使得假如a和b在所选集合里且(a+b)/2∈{1,2,...n}那么(a+b)/2也在所选集合里 f[i]=2*f[i-1]-f[i-2]+g[i] g[n]:选 ...

随机推荐

  1. buu 新年快乐

    一.查壳 发现是upx的壳. 二.拖入ida,发现要先脱壳. 题外话.总结一下手动脱壳,esp定律: 1.先单步到只有esp红色时,右键数据窗口跟随. 2.到数据窗口后,左键硬件访问,byte和wor ...

  2. asp网页防止乱码

    <%@LANGUAGE="VBSCRIPT" CODEPAGE="65001"%> <%Session.CodePage=65001%> ...

  3. python 遍历文件夹中所有文件

    '''使用walk方法递归遍历目录文件,walk方法会返回一个三元组,分别是root.dirs和files. 其中root是当前正在遍历的目录路径:dirs是一个列表,包含当前正在遍历的目录下所有的子 ...

  4. RPC远程协议之Thrift入门

    在上一篇文章<RPC远程协议之原理分析>中,我介绍了RPC的工作原理及欲实现RPC框架功能应该做哪些事情,因为要做的事情太多,完全由开发人员研发实现,不是很现实,所以市面上出现了诸多RPC ...

  5. POJ3662

    poj3662 大意:n个点p条边的无向图,求在删去k条边后使1和n号点联通路径上的最长边最小值. 一开始理解错题意以为是分层图求最短路径,结果写完发现k太大了发现事情没有那么简单(讨厌英语题面!) ...

  6. Xshell记录Linux连接操作日志遇到的坑

    1.问题描述: 在Windows上,以前一直使用Secure CRT连接Linux主机进行远程操作,使用CRT的日志功能记录连接过程中的所有操作以及输出. 最近(2019-8-17)使用Xshell进 ...

  7. linux删除命令

    Linux下的命令,删除文件夹下的所有文件,而不删除文件夹本身rm -rf *

  8. Mysql分区表使用

    一.什么是分区表 分区表就是按一定规则将一张表分割成多个部分,达到和物理分表同样的效果,但操作起来更 简单,不同于物理分表那样使用时还需要指定使用的数据表.对于使用者来说和操作普通表无差别. 二.优势 ...

  9. C实现奇偶校验

    奇偶校验原理(来自百度百科):奇偶校验(Parity Check)是一种校验代码传输正确性的方法.根据被传输的一组二进制代码的数位中"1"的个数是奇数或偶数来进行校验.采用奇数的称 ...

  10. SignalR的客户端.NET Client介绍

    SignalR支持两种客户端:JavaScript Client和.NET Client.一个简单的SignalR例子中的SignalRDemo(点赞页面)就是JavaScript Client(HT ...