POJ 2954解题报告

来源：互联网发布：fanuc铣孔螺旋下刀编程编辑：程序博客网时间：2024/04/30 21:04

这道题自己想了个方法：遍历矩形区域(xmin, ymin, xmax, ymax)内的所有点，看是否落在两边之间。

然后就TLE了。

网上查解题报告。知道了pick定理:http://en.wikipedia.org/wiki/Pick's_theorem

//pick theorem
//A = i + b / 2 - 1

多边形的面积（取整）等于多边形内部的整数点个数(i) 加上边界整数点个数(b)的一半再减一。

证明方法见上述维基百科与http://jwilson.coe.uga.edu/emat6680fa05/schultz/6690/pick/pick_main.htm

我没有看啊。。。老了。。。

三角形的面积在知道三个顶点位置后可通过下面的公式得出：

//A = 0.5 * abs((x1 - x3)(y2 - y1) - (x1 - x2)(y3 - y1))

见http://en.wikipedia.org/wiki/Triangle

三角形一条边上的所有整数点（包括顶点）可以首先将这条边移到(0, 0)->(x, y)。这时，（x/gcd(x, y), y/gcd(x, y)）肯定在这条边上，并且是整数点，其余所有整数点的可以表示为k(x/gcd(x, y), y/gcd(x, y))。所以所有的整数点个数为gcd(x, y) + 1。即：

//b = gcd(x, y) + 1

这里解释得比较清楚:http://www.darkswordzone.com/?p=974

在上述A, b求出后（b = b1 + b1 + b3 - 3, 减去算了两边的三个顶点，因为顶点在题中一定是整数点），我们有：
//i = A + 1 - b / 2

代码如下：

#include <iostream>using namespace std;//pick theorem//A = i + b / 2 - 1//A = 0.5 * abs((x1 - x3)(y2 - y1) - (x1 - x2)(y3 - y1))//b = gcd(x, y) + 1//i = A + 1 - b / 2int gcd(int x, int y){if(x > y){if(y == 0)return x;elsereturn gcd(x - y, y);}else{if(x == 0)return y;elsereturn gcd(y - x, x);}}int main(){int x1, y1, x2, y2, x3, y3;while(true){cin>>x1>>y1>>x2>>y2>>x3>>y3;if(!x1 && !y1 && !x2 && !y2 && !x3 && !y3){return 0;}int A = 0.5 * abs((x1 - x3) * (y2 - y1) - (x1 - x2) * (y3 - y1));int b1 = gcd(abs(x2 - x1), abs(y2 - y1)) + 1;int b2 = gcd(abs(x3 - x1), abs(y3 - y1)) + 1;int b3 = gcd(abs(x3 - x2), abs(y3 - y2)) + 1;int i = A + 1 - (b1 + b2 + b3 - 3) / 2;cout<<i<<endl;}}