CF 782B 和 HDU 4717 经典移点问题 【 二分 和 三分 时间】

CF 782B 传送门

/** @Cain*/
const int maxn = 1e5+5;
const db eps = 1e-7;
int a[maxn],b[maxn];
int main()
{
    int n;
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
    }
    for(int i=1;i<=n;i++){
        scanf("%d",&b[i]);
    }
    db l = 0, r = inf , mid;
    while(r - l > eps){
        mid = (r+l)/2.0;
        db maxx = -inf, minn = inf;
        for(int i=1;i<=n;i++){
            db xl = a[i]*1.0 - b[i]*mid;
            db xr = a[i]*1.0 + b[i]*mid;
            maxx = max(maxx,xl);  //千万不要找错了,是找最小的进
 //行比较, 而不是最大的, 既是往左和右移动距离小的那两个.
 //如果直接找最大的两个点则有可能有交叉部分.
            minn = min(minn,xr);
        }
        if(fabs(maxx - minn) < eps)
            break;
        else if(maxx < minn) r = mid;
        else l = mid;
    }
    printf("%.12lf\n",mid);
}

HDU 4717 传送门
//思路：在平面中任意两个点的距离一定是先减小, 后增大的, 不信可以画出来看一看. 所以满足单峰函数, 所以就可以用三分来做.

/** @Cain*/
const int maxn = 3e2+5;
const db eps = 1e-6;
int x[maxn],y[maxn],vx[maxn],vy[maxn];
int n;
int cas = 1;
db cal(db t)
{
    db res = 0;
    for(int i=1;i<=n;i++){
        for(int j=i+1;j<=n;j++){
            db x1 = x[i]*1.0+vx[i]*t;
            db y1 = y[i]*1.0+vy[i]*t;
            db x2 = x[j]*1.0+vx[j]*t;
            db y2 = y[j]*1.0+vy[j]*t;
            db maxx = sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
            res = max(maxx,res);
        }
    }
    return res;
}
void solve()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%d%d%d%d",&x[i],&y[i],&vx[i],&vy[i]);
    }
    db l = 0, r = inf,len;
    db lm,rm;
    while(r - l > eps){
        len = (r-l)/3.0; 
        lm = l + len;
        rm = r - len;
        db dislm = cal(lm);
        db disrm = cal(rm);
        if(dislm < disrm) r = rm;
        else l = lm;
    }
    printf("Case #%d: %.2f %.2f\n",cas++,l,cal(l));
}