题目传送门
。
解法:
跟水平可见直线很像啊。只不过这道题是从y轴开始的射线。
思路其实差不多。
单调栈维护下凸包。
先按照直线的k排序。
那么进来之后t其他直线。
当前直线为i。栈顶直线为x。
如果i和x的交点在y轴左侧。
那么说明在起点之前x才能超过i。
那么这是不合法的,所以把x踢掉。
然后剩下的就是水平可见直线了。
栈顶下一个元素为y。
如果i和x的交点在x和y交点的左边的话。那么x也要被踢掉
画个图一切都明白了。
最后栈里剩下的就是答案。
代码实现:
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cmath>
using namespace std;
struct node {int k,b,id;}a[51000];
int f[51000],top;const double esp=1e-8;
double jd(node a,node b) {return (double)(b.b-a.b)/(a.k-b.k);}
bool cmp(node n1,node n2) {if(n1.k!=n2.k)return n1.k<n2.k;else return n1.b<n2.b;}
bool cmp1(int n1,int n2) {return a[n1].id<a[n2].id;}int ans[11000];
int main() {
int n;scanf("%d",&n);
for(int i=1;i<=n;i++)scanf("%d",&a[i].b);
for(int i=1;i<=n;i++)scanf("%d",&a[i].k);
for(int i=1;i<=n;i++) a[i].id=i;
sort(a+1,a+1+n,cmp);top=0;
for(int i=1;i<=n;i++){
while(top>=1&&(a[i].b>a[f[top]].b ||
jd(a[i],a[f[top]])<jd(a[f[top]],a[f[top-1]]) ))top--;
f[++top]=i;
}
printf("%d\n",top);sort(f+1,f+1+top,cmp1);
printf("%d",a[f[1]].id);
for(int i=2;i<=top;i++)printf(" %d",a[f[i]].id);
return 0;
}