51nod3149 Triangle

3149 Triangle

Alice不是个安分的学生。

Alice太聪明了，以至于每次上课数学老师要讲的东西她很快就懂了，之后就开始玩手机。

为了挫挫Alice的锐气，数学老师给她出了一道这样的题目：

在一个平面直角坐标系内，有一个以整点为顶点的三角形，问这个三角形内有多少个整点？

你会做吗？

输入

多组数据，每组数据输入一行六个整数x1 y1 x2 y2 x3 y3，分别代表三角形三个顶点的坐标
输入最后以0 0 0 0 0 0 结尾，对于这一行输入不用输出答案

输出

对于每组数据，输出一行表示答案

数据范围

对于30%的数据： 数据组数和坐标绝对值<=300
对于60%的数据：数据组数和坐标绝对值<=3000
对于100%的数据： 数据组数和坐标绝对值<=50000

输入样例

输入样例1：
0 0 1 0 0 1
输入样例2：
0 0 5 0 0 5
输入样例3：
0 0 1 0 0 1
0 0 5 0 0 5

输出样例

输出样例1：
0 
输出样例2：
6
输出样例3：
0
6

分析：

运用PICK定理，定义如下：

PICK定理：格子上的多边形面积=边界上格子点数/2+内部点数-1。

放代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <vector>
#include <queue>
#define INF 2139062143
#define MAX 0x7ffffffffffffff
#define del(a,b) memset(a,b,sizeof(a))
#define Rint register int
using namespace std;
typedef long long ll;
template<typename T>
inline void read(T&x)
{
    x=0;T k=1;char c=getchar();
    while(!isdigit(c)){if(c=='-')k=-1;c=getchar();}
    while(isdigit(c)){x=x*10+c-'0';c=getchar();}x*=k;
}
struct G{
    double x,y;
    G(double x=0.0,double y=0.0):x(x),y(y){}
};
G operator + (const G &a ,const G& b){return G(a.x+b.x,a.y+b.y);}
G operator - (const G &a ,const G& b){return G(a.x-b.x,a.y-b.y);}
double cross(G a,G b){
    return fabs(a.x*b.y-a.y*b.x);
}
struct node{
    int x,y;
};
node a[3];
int gcd(int a,int b){
    return b==0?a:gcd(b,a%b);
}
int get_node(node a,node b){
    double x=fabs((double)a.x-b.x),y=fabs((double)a.y-b.y);
    if(x==0&&y==0) return 0;
    return gcd(x,y)-1;
}
int main()
{
    while(1){
        bool f=0;
        for(int i=0;i<3;i++){
            read(a[i].x),read(a[i].y);
            if(a[i].x || a[i].y) f=1;
        }
        if(!f) break;
        G b=G(a[1].x-a[0].x,a[1].y-a[0].y),c=G(a[2].x-a[0].x,a[2].y-a[0].y);
        int s=int(cross(b,c)/2);
        int L=get_node(a[0],a[1])+get_node(a[0],a[2])+get_node(a[2],a[1])+3;
        printf("%d\n",s+1-L/2);
    }
    return 0;
}