bzoj3190 [JLOI2013]レーシングハーフプレーンクロッシング

アイデアは凸包を見つけるのに少し似ています

最初に初期位置に従ってソートし、次に同じ位置を勾配の昇順でソートします

次に、直線を追加して激しく見つけます（2つのポイントを使用できます）

各直線の選択した間隔を保存すると、有効な勾配の単一の増加のため、それがカバーされている場合、それが最後の間隔になります

複数の直線が区間の端点で交差し、直線と一致することに注意してください。後者の点をカバーする必要があるため、前の直線を直接追加してから、この直線を最後に最も高い直線として追加します。最後の間隔

コード：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
using namespace std;
long double X;
double lin;
int i,j,n,top,sta[10004],cnt,ans[10005];
struct la
{long double k,y,bj;
int id;
}a[10005];
vector<int>v[10004];
bool cmp(la a,la b)
{
    if(a.y==b.y)return a.k<b.k;
    return a.y>b.y;
}
int main()
{	
    scanf("%d",&n);
    for(i=1;i<=n;i++)
    scanf("%lf",&lin),a[i].y=lin,a[i].id=i;
    for(i=1;i<=n;i++)
    scanf("%lf",&lin),a[i].k=lin;
    sort(a+1,a+1+n,cmp);
    sta[++top]=1;a[1].bj=1999999999999999999;
    for(i=2;i<=n;i++)
    {
        for(j=1;j<=top;j++)
        {
            if(a[i].k==a[sta[j]].k)break;
         //	cout<<a[sta[j]].y-a[i].y<<" "<<a[i].k-a[sta[j]].k<<endl;   
        X=(a[sta[j]].y-a[i].y)/(a[i].k-a[sta[j]].k);
       // cout<<" "<<a[sta[j]].y-a[i].y<<":"<<a[i].k-a[sta[j]].k<<"=="<<X<<endl;
        if(X<0)break;
        if(X+0.000000000001>a[sta[j]].bj&&X-0.000000000001<a[sta[j]].bj)//右边界重合 
        {
        v[sta[j]].push_back(a[sta[j+1]].id);//原来右边界的线，和这条线一样 
		  top=j+1;
            sta[top]=i;    a[i].bj=1999999999999999999;   	
		break;
		}
        else
        if(X<a[sta[j]].bj)//提前重合，右边的点无意义 
        {
        top=j+1;
            sta[top]=i;
            a[sta[j]].bj=X;
            a[i].bj=1999999999999999999;
            break;		
        }	
        }		
    }
    for(i=1;i<=top;i++)
    {
        int o=sta[i];
        ans[++cnt]=a[o].id;
        for(j=0;j<v[o].size();j++)
        ans[++cnt]=v[o][j];
    }
    sort(ans+1,ans+1+cnt);
    printf("%d\n",cnt);
    for(i=1;i<cnt;i++)
    printf("%d ",ans[i]);
    printf("%d",ans[cnt]);
}