题解 \(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
#include<bits/stdc++.h>
#define ri register signed
#define p(i) ++i
using namespace std;
namespace IO{
char buf[1<<21],*p1=buf,*p2=buf;
#define gc() p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++
template<typename T>inline void read(T &x) {
ri f=1;x=0;register char ch=gc();
while(ch<'0'||ch>'9') {if (ch=='-') f=0;ch=gc();}
while(ch>='0'&&ch<='9') {x=(x<<1)+(x<<3)+(ch^48);ch=gc();}
x=f?x:-x;
}
}
using IO::read;
namespace nanfeng{
#define FI FILE *IN
#define FO FILE *OUT
template<typename T>inline T cmax(T x,T y) {return x>y?x:y;}
template<typename T>inline T cmin(T x,T y) {return x>y?y:x;}
typedef double db;
static const int N=3e5+7;
int st[N],tp,pre[N],vis[N],n,mxa,mxb;
db tk[N];
struct node{int a,b,id;}mis[N];
inline int operator<(const node &n1,const node &n2) {return n1.a==n2.a?n1.b>n2.b:n1.a>n2.a;}//每个相同的 a 只有最大的 b 有用
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));}
inline int main() {
// FI=freopen("nanfeng.in","r",stdin);
// FO=freopen("nanfeng.out","w",stdout);
read(n);
for (ri i(1);i<=n;p(i)) read(mis[i].a),read(mis[i].b),mis[i].id=i;
sort(mis+1,mis+n+1);
for (ri i(1);i<=n;p(i)) if(mxb<mis[i].b) mxb=mis[i].b,mxa=mis[i].a;
st[p(tp)]=1;
for (ri i(2);mxa<=mis[i].a&&i<=n;p(i)) {//mxa 防止右凸包也算进去
if (mis[i].a==mis[st[tp]].a) {
if (mis[i].b==mis[st[tp]].b)
pre[mis[i].id]=pre[mis[st[tp]].id],pre[mis[st[tp]].id]=mis[i].id;
continue;
}
while(tp>1&&slope(mis[st[tp]],mis[i])<tk[st[tp]]) --tp;
tk[i]=slope(mis[st[tp]],mis[i]),st[p(tp)]=i;
}
while(tp)
for (ri i(mis[st[tp--]].id);i;i=pre[i]) vis[i]=1;
for (ri i(1);i<=n;p(i)) if (vis[i]) printf("%d ",i);
return 0;
}
}
int main() {return nanfeng::main();}

最新文章

  1. thinkphp疑难解决4
  2. 第一届山东省ACM——Phone Number(java)
  3. FizzBuzzWhizz游戏的高效解法
  4. Java之数组了解
  5. atitit.sql server2008导出导入数据库大的表格文件... oracle mysql
  6. 通达信:显示K线图日期
  7. 修改placeholder颜色
  8. POJ 2728 Desert King
  9. Java_Activiti5_菜鸟也来学Activiti5工作流_之与Spring集成(三)
  10. &lt;a&gt;标签中href=&quot;javascript:;&quot;
  11. HDU 1084 - ACM
  12. c# 课堂总结7--函数
  13. HoloLens开发手记 - 构建2D应用 Building 2D apps
  14. Android服务器——使用TomCat实现软件的版本检测,升级,以及下载更新进度!
  15. C++设计模式——命令模式
  16. mybatis(3)---传参数的方法
  17. 用delete和trancate删除表记录的区别
  18. js检测当前设备是移动端还是PC端
  19. 为什么call比apply快
  20. android开发(37) android使用android_serialport_api 操作串口,解决权限问题

热门文章

  1. ZYNQ 中PS端GPIO EMIO使用
  2. ESP32使用SPIFFS文件系统笔记
  3. C++泛型编程-举例
  4. 痞子衡嵌入式:了解i.MXRTxxx系列ROM中灵活的串行NOR Flash启动硬复位引脚选择
  5. .Net Core with 微服务 - Polly 服务降级熔断
  6. Mycat读写分离的简单实现
  7. Spring总结之事务
  8. canvas实现任意正多边形的移动(点、线、面)
  9. Linux虚拟机扩展磁盘
  10. P4774-屠龙勇士-扩展中国剩余定理