POJ-3304

我觉得这就是认真了就输了。。。
看见题里说了 You must assume that two floating point numbers a and b are equal if |a - b| < 10-8. 是吧。。。
就真写了个const double eps = 1e-8(一开始交的时候还是1e-10咧)
以及考虑eps的加减法和赋值
他说你就信?!
连样例都过不去信不信?!
所以说看到连样例都过不去,顿时吓尿了,对着那三个函数一通魔改,
样例终于过了。
那数据能过的了嘛?!(kuma)
于是乎不管那么多了,注释了add和dcmp,以及删除了相关的调用,又试了一发。。。
然后就特么过了!
精度真是坑人,
那么下面是代码:

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <climits>
#include <ctime>
#include <iostream>
#include <algorithm>
using namespace std;
const double eps = 1e-8;
struct P
{
    double x, y;
//    double add(double a, double b)
//    {
//        if(fabs(a + b) < fabs(a) + fabs(b)) return 0;
//        else return a + b;
//    }
//    int dcmp(double x)
//    {
//        if(fabs(x) < eps) return 0;
//        else return x < 0 ? -1 : 1;
//    }
    P() {}
    P(double x, double y) : x(x), y(y) {}
    P operator +(const P &p)
    {
        return P(x + p.x, y + p.y);
    }
    P operator -(const P &p)
    {
        return P(x - p.x, y - p.y);
    }
    P operator *(const double &d)
    {
        return P(x * d, y * d);
    }
    bool operator ==(const P &p)
    {
        return x == p.x && y == p.y;
    }
    double det(const P &p)
    {
        return x * p.y - y * p.x;
    }
    double dot(const P &p)
    {
        return x * p.x + y * p.y;
    }
    P operator =(const P &p)
    {
        x = p.x, y = p.y;
        return *this;
    }
};
P start[117], fine[117], all[217];

inline void init()
{
    memset(all, 0, sizeof(all));
    memset(start, 0, sizeof(start));
    memset(fine, 0, sizeof(fine));
return ;
}

inline bool on_point(P p1, P p2, P q1, P q2)
{
    return p1 == q1 || p1 == q2 || p2 == q1 || p2 == q2;
}

inline bool on_seg(P p1, P p2, P q)
{
    return (q - p1).det(p2 - q) == 0 && (q - p1).dot(p2 - q) >= 0;
}

inline bool on_inster(P p1, P p2, P q1, P q2)
{
    return ((q2 - q1).det(p1 - q1)) * ((p2 - q1).det(q2 - q1)) > 0;
}

int main()
{
    int t, n, i, j, k;
    while(scanf("%d", &t) != EOF)
    {
        while(t--)
        {
            scanf("%d", &n);
            init();
            for(i = 0; i < n; ++i)
            {
                scanf("%lf %lf %lf %lf", &start[i].x, &start[i].y, &fine[i].x, &fine[i].y);
                all[2 * i] = start[i], all[2 * i + 1] = fine[i];
            }
            if(n == 1 || n == 2)
            {
                puts("Yes!");
                continue;
            }
            bool flag;
            for(i = 0; i < 2 * n; ++i)
            {
                for(j = i + 1; j < 2 * n; ++j)
                {
                    flag = true;
                    P line_start = all[i], line_end = all[j];
                    for(k = 0; k < n; ++k)
                    {
                        if(on_point(start[k], fine[k], line_start, line_end)) continue;
                        if(on_seg(start[k], fine[k], line_start) || on_seg(start[k], fine[k], line_end) ||
                            on_seg(line_start, line_end, start[k]) || on_seg(line_start, line_end, fine[k])) continue;
                        if(on_inster(start[k], fine[k], line_start, line_end)) continue;
                        else
                        {
                            flag = false;
                            break;
                        }
                    }
                    if(flag == true) break;
                }
                if(flag == true) break;
            }
            if(flag == true) puts("Yes!");
            else puts("No!");
        }
    }
return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值