題目傳送門
。
解法:
跟水平可見直線很像啊。只不過這道題是從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;
}