hdu 2438

最新推荐文章于 2021-12-01 20:53:09 发布

Dorkdomain

最新推荐文章于 2021-12-01 20:53:09 发布

阅读量380

点赞数

本文链接：https://blog.csdn.net/Dorkdomain/article/details/45172495

版权

在网上看到两种方法始终不能理解另一种方法还好另外一种方法正确性显而易见

车转弯的时候车有段与地面的夹角角度是从0度变化到90度的。也就是转弯的时候需要一个最大的宽度才能过去。

否则就卡在那里了。这个宽度PH是先增加后减少的。是个凸型函数，因此是三分求的极值。

直线y的斜率为tan(θ)，还经过点（0， Lsin(θ)+D/cos(θ))因此得到y的直线方程。

y=xtan(θ)+Lsin(θ)+D/cos(θ)

求的PH就是当y=X(汽车当前在的街道的宽度)时，解出的x的值的绝对值。

-x=|x|=(lsinθ+w/cosθ-x)/tanθ; (l==L, D==w)

三分什么的都和计算几何有关系一直没时间去学这个东西一点也不会唉

#include <iostream>
#include <algorithm>
#include <queue>
#include <cstring>
#include <cstdio>
#include <vector>
#include <string>
#include <iterator>
#include <cmath>
#include <deque>
#include <stack>
#include <cctype>
#include <iomanip>
#include <stack>
using namespace std;

typedef long long ll;
typedef long double ld;

const int N = 1100;
const int INF = 0xfffffff;
const double EPS = 1e-8;
const ll MOD = 1e9 + 7;
const ld PI = acos(-1);

#define INFL 0x7fffffffffffffffLL
#define met(a, b) memset(a, b, sizeof(a))
#define rep(c, a, b) for (int c = a; c < b; c++)
#define nre(c, a, b) for (int c = a; c > b; c--)
#define put(a) cout << setiosflags(ios::fixed) << setprecision(a)

double x, y, l, w;
double cot(double x);

int main ()
{

    while (cin >> x >> y >> l >> w)
    {
        double l = 0, r = PI / 2, lm, rm;
        while (r - l > EPS)
        {
            lm = (l + r) / 2;
            rm = (lm + r) / 2;
            if (cot(lm) < cot(rm)) l = lm;
            else r = rm;
        }
        if (cot(lm) > y || x < w || y < w)
            cout << "no" << endl;
        else
            cout << "yes" << endl;
    }
    return 0;
}

double cot(double a)
{
    return ((l * sin(a) + w / cos(a) - x) / tan(a));
}