poj3429（有错）

最新推荐文章于 2019-04-17 23:23:03 发布

zhengnanlee

最新推荐文章于 2019-04-17 23:23:03 发布

阅读量829

点赞数

分类专栏：计算几何学文章标签： ACM题解报告计算几何学

本文链接：https://blog.csdn.net/zhengnanlee/article/details/11861731

版权

计算几何学专栏收录该内容

40 篇文章 0 订阅

订阅专栏

不知道为什么错。。

/*************************************************************************
> File Name: poj3429.cpp
> Author: zhengnanlee
> Mail: zhengnanlee@hotmail.com
> Created Time: 2013年09月21日星期六 09时10分12秒
************************************************************************/

#include <iostream>
#include <math.h>

using namespace std;

#define eps 1e-5
#define zero(x) (((x)>0?(x):-(x))<eps)
#define LL double

struct point{ LL x, y; }p[105];

int parallel(point u1, point u2, point v1, point v2)
{
    return zero((u1.x - u2.x)*(v1.y - v2.y) - (v1.x - v2.x)*(u1.y - u2.y));
}

point intersection(point u1, point u2, point v1, point v2)
{
    point ret = u1;
    LL t = ((u1.x - v1.x)*(v1.y - v2.y) - (u1.y - v1.y)*(v1.x - v2.x)) / ((u1.x - u2.x)*(v1.y - v2.y) - (u1.y - u2.y)*(v1.x - v2.x));
    ret.x += (u2.x - u1.x)*t;
    ret.y += (u2.y - u1.y)*t;
    return ret;
}

double distance(point a)
{
    return sqrt(a.x * a.x + a.y * a.y);
}

int main()
{
    int N;
    cin >> N;
    int countt = N;
    for (int i = 1; i <= N; i++)
    {
        cin >> p[i].x >> p[i].y;
    }
    int M;
    cin >> M;
    int cnt = 0;
    for (int i = 0; i < M; i++)
    {
        int a, b, c, d;
        cin >> a >> b >> c >> d;
        if (!parallel(p[a], p[b], p[c], p[d]))
        {
            point temp = intersection(p[a], p[b], p[c], p[d]);
            p[++countt] = temp;
            if (zero(distance(temp))) cnt++;
        }
    }
    cout << cnt << endl;
}