uvaLive 6748 2D-Solar System

题目连接：https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=4760

题目描述：有n 个圆在一条直线上，直线无限长，每个圆有一个速度，可能往左，也可能往右，保证所有的圆在一开始的时候都相离。问能不能所有的圆永远都不相撞，如果相撞，输出相撞的时间。每个圆有3个量，半径 r，速度 a， 起始位置 b。

解题思路：

错解（我一开始的思路）：按照起始位置排序，考虑相邻的两个是否相撞，如果相撞，则更新答案（其实是有问题的，第一个圆可以跳过第二个圆与第三个圆相撞！）

正解：思路有点类似斜率优化dp。要维护一个next数组，next[i]表示i 左边第一个比i 的半径大的圆的位置。同样，是要先按照起始位置排序！！

复杂度是O(n)的，证明类似斜率优化dp的证明（每个位置只会入队一次，出队一次）。。。

ps：我实在太弱了。。。

详见代码。。

//#pragma comment(linker,"/STACK:102400000,102400000")
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<string>
#define ll long long
#define db double
#define PB push_back
using namespace std;

const int N = 50005;
const db INF = 1e50;
const db eps = 1e-7;

ll getll()
{
    char c=getchar();
    while(!(c>='0'&&c<='9')) c=getchar();
    ll res=c-'0';
    while((c=getchar())>='0'&&c<='9') res=res*10+c-'0';
    return res;
}

struct node
{
    ll r,a,b;
    void input()
    {
        scanf("%lld%lld%lld",&r,&a,&b);
    }
    bool operator < (const node &t) const
    {
        return b<t.b;
    }
} p[N];
int next[N];

int main()
{
#ifdef PKWV
    freopen("in.in","r",stdin);
#endif // PKWV
    int n;
    while(scanf("%d",&n)&&n)
    {
        for(int i=0;i<n;i++) p[i].input();
        sort(p,p+n);
        db ans=INF;
        for(int i=0;i<n;i++)
        {
            int j=i-1;
            while(j>=0)
            {
                db t=p[i].b-p[j].b-sqrt((db)((p[i].r+p[j].r)*(p[i].r+p[j].r)-abs(p[i].r-p[j].r)*abs(p[i].r-p[j].r)));
                t/=(db)(p[j].a-p[i].a);
                if(t>eps) ans=min(ans,t);
                if(p[j].r>p[i].r) break;
                j=next[j];
            }
            next[i]=j;
        }
        if(fabs(ans-INF)<1e-7)
        {
            printf("Collision-Free System\n");
        }else printf("%.2f\n",ans);
    }
    return 0;
}