计算几何 pick定理

Pick定理：设以整数点为顶点的多边形的面积为S，多边形内部的整数点数为N，边界上的整数点数为L，则S=N+1/2*L-1

Pick定理主要是计算格点多边形（定点全是格点的不自交图形）P的面积与其边界和内部格点数之间的关系。

　　格点多边形的面积A(P)可以通过叉积计算出来，不过叉积计算出来的面积是实际面积的2倍；

　　边界上的格点B(P)可以通过计算相邻两点的横坐标之差与纵坐标之差的最大公约数的和得到；

　　内部的格点I(P)则通过公式得：I(P) = A(P)-B(P)/2+1计算出。

　　解释：

　　　a.关于边界格点计算两点横纵坐标之差就是以两个点构成的边做坐标轴，组成的三角形（或者线）的两个之角标求gcd

　　　b.格点多边形的面积是通过将多边形固定一个点，然后在遍历每两个点，三个点构成的三角形求面积。由于叉积可以为负，所以不必担心多加的三角形或者不在多边形内部的三角形，都会减去。

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
int L,N;
double S;
int n;
struct F
{
    int x,y;
}p[105];  //x为横坐标，y为纵坐标

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

void GetS()  //计算S
{
    for(int i=0;i<n;i++)
        S+=p[i].x*p[(i+1)%n].y-p[i].y*p[(i+1)%n].x;
    //S面积正负号代表方向
}

void OnEdge()
{
    for(int i=0; i<n; i++)
        L+=gcd(fabs(p[i].x-p[(i+1)%n].x),fabs(p[i].y-p[(i+1)%n].y));
}

void Inside()  //计算N
{
    N=(S-L+2)/2;
}
int main()
{
    cin>>n;
    for(int i=0;i<n;i++)
        cin>>p[i].x>>p[i].y;
    GetS(),OnEdge(),Inside();
    cout<<"S:"<<S*0.5<<endl;
    cout<<"OnEdge:"<<L<<endl;
    cout<<"InSide:"<<N<<endl;
}

 

poj 2954

http://poj.org/problem?id=2954

给出三角形的三个顶点左边，求三角形内部的点数

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
int S=0,L=0,N=0;
int n=3;
struct F
{
    int x,y;
}p[105];  //x为横坐标，y为纵坐标

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

void GetS()  //计算S
{
    for(int i=0;i<n;i++)
        S+=p[i].x*p[(i+1)%n].y-p[i].y*p[(i+1)%n].x;
    S=fabs(S);
}

void OnEdge()
{
    for(int i=0; i<n; i++)
        L+=gcd(fabs(p[i].x-p[(i+1)%n].x),fabs(p[i].y-p[(i+1)%n].y));
}

void Inside()  //计算N
{
    N=(S-L+2)/2;
}
int main()
{
    while(~scanf("%d %d %d %d %d %d", &p[0].x,&p[0].y,&p[1].x,&p[1].y,&p[2].x,&p[2].y))
    {
        S=0,N=0,L=0;
        if (p[0].x==0&&p[0].y==0&&p[1].x==0&&p[1].y==0&&p[2].x==0&&p[2].y==0)
            break;
        GetS();
        OnEdge();
        Inside();
        cout<<N<<endl;
    }
    return 0;
}