POJ 1265 Area pick定理 + 多边形面积求法

点击打开链接

题目大意：

给一个平面上的简单多边形，求边上的点，多边形内的点，多边形面积。

解题思路：
这个题用了很多知识点：

1、以格子点为顶点的线段，覆盖的点的个数为GCD(dx,dy)，其中，dxdy分别为线段横向占的点数和纵向占的点数。如果dx或dy为0，则覆盖的点数为dy或dx。
2、Pick公式：平面上以格子点为顶点的简单多边形的面积=边上的点数/2+内部的点数 - 1。
3、任意一个多边形的面积等于按顺序求相邻两个点与原点组成的向量的叉积之和。

设Ω是m边形(如下图)，顶点沿边界正向排列，，坐标依次为

建立Ω的多边形区域向量图。

由图知坐标原点与多边形任意相邻的两个顶点构成一个三角形，而三角形的面积可由三个顶点构成的两个平面向量的外积求得。

任意多边形的面积公式

多边形计算公式的计算和原点的选取没有关系，通常可以选点（0,0）或者多边形的第一个点（这个时候比较直观了，看起来就是把多边形分成一个个三角形和加起来，读者自己可以画个图）就可以了。

用G++ 交WA 用C++ 交AC了

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <cstring>
using namespace std;

int GCD(int a, int b) {
    while (b) {
        int r = b;
        b = a % b;
        a = r;
    }
    return a;
}

int main()
{
    int t, cas, n, i, x1, y1, x2, y2, dx, dy, on, in, area;

    scanf ("%d", &cas);
    for (t = 1; t <= cas; t++) {
        scanf ("%d", &n);
        x1 = y1 = on = area = 0;
        for (i = 1; i <= n; i++) {
            scanf ("%d%d", &dx, &dy);
            x2 = x1 + dx;
            y2 = y1 + dy;
            on += GCD(abs(dx), abs(dy));
            area += x1 * y2 - x2 * y1;
            x1 = x2;
            y1 = y2;
        }
        if(area<0)area=-area;
        in = (area + 2 - on) / 2;
        printf ("Scenario #%d:\n%d %d %.1f\n\n", t, in, on, area/2.0);
    }
    return 0;
}