计算几何之皮克定理的应用+线段上整数点个数+求三角形面积

最新推荐文章于 2023-09-17 17:01:09 发布

原创最新推荐文章于 2023-09-17 17:01:09 发布 · 1k 阅读

1 ·

CC 4.0 BY-SA版权

本文为博主原创文章，未经博主允许也可以转载。

计算几何专栏收录该内容

11 篇文章

订阅专栏

本文介绍了一种使用皮克定理计算由整数坐标定义的多边形面积的方法。通过分析多边形内部和边界上的整数点数，结合叉积原理求得面积，并给出具体实现代码。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

GYM-G
思路：皮克定理： s = in + on / 2 - 1
其中 s为多边形面积， in为在多边形内部的点数，on为在多边形边界上的点数。都为整数点数。
根据叉积求面积。
on可以通过 gcd( a.x - b.x , a.y - b.y )得到，一个线段上的整数点数（gcd( a.x - b.x , a.y , b.y )）。
Code：

#include <bits/stdc++.h>
#define LL long long
using namespace std;
const int AX = 1e5 + 666 ; 
struct Node{
	LL x , y ;
	double operator ^ ( const Node &b)const{  
		return x * b.y - y * b.x ;  
	}  
}p[AX];
int n ;
int gcd( int a , int b ){
	return !b ? a : gcd( b , a % b ) ;
}
double getArea(){
	double res = 0.0 ;
	for(int i = 0; i < n ; i++)
		res += ( p[i] ^ p[(i+1)%n] ) / 2 ;
	return fabs(res);
}
int main(){
	scanf("%d",&n) ;
	for( int i = 0 ; i < n ; i++ ){
		scanf("%lld%lld",&p[i].x,&p[i].y) ;
	}
	double s = getArea() ; 
	LL on = 0LL ;
	for( int i = 0 ; i < n ; i++ ){
		on += gcd( abs( p[i].x - p[(i+1)%n].x ) , abs( p[i].y - p[(i+1)%n].y ) ) ;
	}
	LL in = (LL)( 2 * s + 2 - on ) / 2LL ; // s = in + on / 2 - 1.0 ;
	printf("%lld\n",in);
	return 0 ; 
}