已经AC 触宝编程题 第二题 三分时间

#大概的思路就是 默认这个距离(最大两点的距离)函数是时间变量的一个凸函数(我猜的,具体证明不会),凸函数找最大值,用三分法就可以了。 然后就三分时间则可以求得答案


#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;

const int MX = 308;
const double eps = 0.00001;
struct V{
    double _x, _y;
    V(){}
    V(double x, double y):_x(x),_y(y){}
    V operator *(double mul) const{
        return V(mul*_x, mul*_y);
    }
    V operator +(const V & v1) const{
        return V(_x+v1._x, _y+v1._y);
    }
    double dis(const V & v1){
        double sum1 = _x - v1._x; sum1 *= sum1;
        double sum2 = _y - v1._y; sum2 *= sum2;
        return sqrt(sum1+sum2);
    }
    void pf(){
        printf("point :%lf %lf\n", _x, _y);
    }
};
int n;
V point[MX];
V speed[MX];
V now[MX];
double get_dis(double tim){
    for(int i = 0; i < n; i++){
        now[i] = point[i] + speed[i]*tim;
    }
    double tmax = 0.0;
    for(int i = 0; i < n; i++){
        for(int j = i+1; j < n; j++){
            tmax = max(tmax, now[i].dis(now[j]));
        }
    }
    return tmax;

}
int main(){
    while(cin>>n){
        for(int i = 0;i < n;i++){
            double x,y;scanf("%lf%lf",&x,&y);
            point[i] = V(x,y);
            double sx,sy;scanf("%lf%lf",&sx,&sy);
            speed[i] = V(sx,sy);
        }
        double res = get_dis(0.00);
        double flag_time = 0.00;
        double l = 0.00;
        double r = 1000000000.00;
        while(l < r - eps){
            double mid1 = (l+r)*0.5;
            double mid2 = (l+mid1)*0.5;
            double res1 = get_dis(mid1);
            double res2 = get_dis(mid2);
            if(res1 < res2){
                l = mid2;
                if(res1 < res){
                    flag_time = mid1;
                    res = res1;
                }
            }
            else{
                r = mid1;
                if(res2 < res){
                    flag_time = mid1;
                    res = res1;
                }
            }
        }
        printf("%.2lf %.2lf\n",flag_time, res);
    }
    return 0;
}

#C++工程师#
全部评论
好厉害啊~
点赞 回复 分享
发布于 2017-09-05 21:28
AC了吗
点赞 回复 分享
发布于 2017-09-05 21:05
这我就不得不膜拜了
点赞 回复 分享
发布于 2017-09-05 21:06
AC的代码
点赞 回复 分享
发布于 2017-09-05 21:26
COOL.
点赞 回复 分享
发布于 2017-09-05 21:26
大佬啊,我还用暴力做的…
点赞 回复 分享
发布于 2017-09-05 21:34
膜拜大佬!!!!!!!!!!!!
点赞 回复 分享
发布于 2017-09-05 21:38
整体应该是凸函数,而且没有局部解,相当于聚类只有一个解。
点赞 回复 分享
发布于 2017-09-05 21:42
大佬,有题目么
点赞 回复 分享
发布于 2017-09-05 22:03
总算有个人懂我的意思了......
点赞 回复 分享
发布于 2017-09-05 22:06

相关推荐

躺尸修仙中:因为很多92的也去卷中小厂,反正投递简历不要钱,面试不要钱,时间冲突就推,不冲突就面试积累经验
点赞 评论 收藏
分享
10-25 12:05
已编辑
湖南科技大学 Java
若梦难了:我有你这简历,已经大厂乱杀了
点赞 评论 收藏
分享
评论
点赞
6
分享
牛客网
牛客企业服务