poj 1265 皮克定理＋多边形面积＋边上点的数量的求取

最新推荐文章于 2021-01-26 21:35:24 发布

原创最新推荐文章于 2021-01-26 21:35:24 发布 · 690 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#C++ #计算几何

本文介绍了计算机几何领域的两个核心概念：皮克定理用于计算多边形面积，并通过整点计数实现；以及如何利用向量叉积计算多边形面积。文章提供了具体的算法实现示例。

计算机几何主要是结论：

１．皮克定理：　设sum_e是边上整点的个数，可以通过gcd求取（不算起点）的整点数，设sum_i为内部的多边形内部的整点，

多边形的面积就是　area = sum_i + sum_e/2 - 1

2 . 多边形面积：可以将多边形拆分为多个三角形，三角形可以利用向量的叉积求取面积，即x1y2 -x2y1的一半

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#define MAX 107

using namespace std;

struct Point 
{
    int x,y;
}p[MAX];

int gcd ( int a , int b )
{
    return b == 0 ? a : gcd ( b , a%b );
}

int det ( Point p1 , Point p2 )
{
    return p1.x * p2.y - p1.y * p2.x;
}

int main ( )
{
    int t,n;
    scanf ( "%d" , &t );
    int c = 1;
    while ( t-- )
    {
        scanf ( "%d" , &n );
        int tx,ty;
        int sum_i = 0 , sum_e = 0 ;
        int area = 0;
        p[0].x = p[0].y =  0; 
        for ( int i = 1 ; i <= n ; i ++ )
        {
            scanf ( "%d%d" , &tx , &ty );
            sum_e += gcd ( abs(tx) , abs(ty) );
            p[i].x = p[i-1].x + tx;
            p[i].y = p[i-1].y + ty;
            area += det ( p[i-1] , p[i] );
        }
        if ( area < 0 ) area = -area;
        printf ( "Scenario #%d:\n", c++ );
        printf ( "%d %d %.1f\n\n",(area-sum_e+2)/2 , sum_e , 0.5*area ) ;
    }
}