Hoj 1188 Orchard Trees

最新推荐文章于 2016-11-23 22:41:24 发布

niuox

最新推荐文章于 2016-11-23 22:41:24 发布

阅读量748

点赞数

分类专栏： Algorithm_计算几何

本文链接：https://blog.csdn.net/niuox/article/details/8593144

版权

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

15 篇文章 0 订阅

订阅专栏

题目：http://acm.hit.edu.cn/hoj/problem/view?id=1188

本题核心在于如何判定一个点P是否在三角形ABC内（包括在边上）。

方法是：看PAB、PAC、PBC面积之和是否等于ABC的面积即可。

求三角形面积可以用叉积来求，注意X的范围在[1,99],不包括0和100,Y一样。

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <vector>
using namespace std;

struct Point
{
    double x;
    double y;
};
Point p[3];

//求面积绝对值
double cross(Point O,Point A,Point B)
{
    return fabs((A.x - O.x) *  (B.y - O.y) - (B.x - O.x) * (A.y - O.y))/2;
}

bool judge()
{
    for(int i=0; i<3; i++)
    {
        if(p[i].x!=0 || p[i].y!=0)
        {
            return true;
        }
    }
    return false;
}

int main()
{
#ifndef ONLINE_JUDGE
    freopen("in.txt","r",stdin);
#endif
    while(scanf(" %lf %lf %lf %lf %lf %lf",&p[0].x,&p[0].y,&p[1].x,&p[1].y,&p[2].x,&p[2].y) && judge())
    {
        double maxX = 0,minX = 150;
        double maxY = 0,minY = 150;
        for(int i=0; i<3; i++)
        {
            if(p[i].x>maxX)
                maxX = p[i].x;
            if(p[i].y>maxY)
                maxY = p[i].y;
            if(p[i].x<minX)
                minX = p[i].x;
            if(p[i].y<minY)
                minY = p[i].y;
        }
        //求三角形的面积
        double area = cross(p[0],p[1],p[2]);
        int num = 0;
        for(int i=ceil(minX); i<=floor(maxX); i++)
        {
            if(i <= 0 || i>=100)
            {
                continue;
            }
            for(int j=ceil(minY); j<=floor(maxY); j++)
            {
                if(j <= 0 || j>=100)
                {
                    continue;
                }
                double s = 0;
                Point temp;
                temp.x = i;
                temp.y = j;
                s += cross(temp,p[0],p[1]);
                s += cross(temp,p[1],p[2]);
                s += cross(temp,p[2],p[0]);
                if(fabs(area-s)<=1e-8)
                {
                    num++;
                }
            }
        }
        printf("%4d\n",num);
    }
}