poj 2954

题目概述

直角坐标系中，给定一个三角形的三个顶点坐标，求其内部格点（横纵坐标皆为整数的点）数，顶点和在边上的格点不计

时限

1000ms/3000ms

输入

每行六个整数，分表描述三个顶点的坐标，输入以六个0结束

限制

所有数字绝对值<=15000

输出

每行一个数，为所求内部格点数

样例输入

0 0 1 0 0 1
0 0 5 0 0 5
0 0 0 0 0 0

样例输出

0
6

讨论

计算几何，pick公式应用，和刚刚做的poj 1265几乎一致（甚至更简单），pick公式为A=I+E/2-1，其中A为面积，I为内部格点数，E为边上格点数，先求出三角形的面积，然后利用gcd求出边上格点数，最后算出内部格点数

题解状态

168K，0MS，C++，1135B

题解代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define INF 0x3f3f3f3f  
#define MAXN 10003
#define memset0(a) memset(a,0,sizeof(a))
#define EPS 1e-6

int gcd(int a, int b)
{
    while (b) {
        int c = a%b;
        a = b;
        b = c;
    }
    return a;
}
int xp(int x1, int y1, int x2, int y2, int x3, int y3)
{
    return (x1 - x2)*(y3 - y2) - (y1 - y2)*(x3 - x2);
}
int fun(int x1, int y1, int x2, int y2, int x3, int y3)
{
    int A = abs(xp(x1, y1, x2, y2, x3, y3));//求出三角形面积
    int E = 0;//初始化边上格点数 下面计算每条边上的格点数
    if ((x1 - x2) && (y1 - y2))
        E += gcd(abs(x1 - x2), abs(y1 - y2));
    else
        E += abs(x1 - x2) + abs(y1 - y2);
    if ((x2 - x3) && (y2 - y3))
        E += gcd(abs(x2 - x3), abs(y2 - y3));
    else
        E += abs(x2 - x3) + abs(y2 - y3);
    if ((x3 - x1) && (y3 - y1))
        E += gcd(abs(x3 - x1), abs(y3 - y1));
    else
        E += abs(x3 - x1) + abs(y3 - y1);
    return (A + 2 - E) / 2;//利用pick公式得到内部格点数
}
int main(void)
{
    //freopen("vs_cin.txt", "r", stdin);
    //freopen("vs_cout.txt", "w", stdout);

    int x1, y1, x2, y2, x3, y3;//在math.h中有名为y1和y0的两个函数 故不能作为全局变量
    while (scanf("%d%d%d%d%d%d", &x1, &y1, &x2, &y2, &x3, &y3) && (x1 || y1 || x2 || y2 || x3 || y3))//input
        printf("%d\n", fun(x1, y1, x2, y2, x3, y3));//output
}

EOF