皮克定理，多边形面积以及线段上整点个数

皮克定理，多边形面积以及拓展欧几里得算法的奇特应用

皮克定理

皮克定理给出了顶点在格点上的多边形的面积公式:

引自百度百科的图¹:

$$\begin{array}{r}2S=2a+b-2\end{array}$$
其中$a$表示多边形内部的点数，$b$表示多边形边界上的整数点个数，$S$表示多边形的面积。

多边形面积

多边形可以划分成多个三角形，其面积就等于三角形面积之和。而三角形面积可以通过叉乘来求。

非常自然的，我们可以从多边形任意一个顶点$p_0$出发，连一条到其他每一个点的连线。

引自大佬博客的图²
通过对每一个三角形叉乘求面积，就可以求出多边形面积

简化：
事实上，$p_0$点大可不必选择多边形顶点。对于平面上任意一点，与$p_1，p_2，...，p_n$所得三角形的面积就等于多边形面积(证明略，一种可行的证明方法是设出$p_0$点，列出叉乘的表达式，最后发现展开之后与$p_0$无关)。因此，方便起见，我们可以选择原点作为$p_0$。

代码：

scanf("%d", &n); // n表示点的个数
for (int i = 0; i < n; i++)
{
	scanf("%d%d", &x[i], &y[i]);
}
long long s = 0;
for (int i = 1; i <= n; i++) { 
	// 向量p0->px就是点px的坐标
	// 通过取模最后算p0和p(n-1)的叉乘
	s += x[i % n] * y[i - 1] - x[i - 1] * y[i % n];
}
printf("%lld", s);

线段上整数点个数

对于一条两顶点在整数点上的线段，顶点为$(x_1,y_1)，(x_2,y_2)$。想要求在这条线段上的整数点的个数。

结论是：整数点的个数$s=gcd(|x_1-x_2|,|y_1-y_2|)+1$
$x_2=\Delta x=pd，y_2=\Delta y=qd$

直观理解：
观察方程$g=gcd(|x_1-x_2|,|y_1-y_2|)，g$是$|x_1-x_2|,|y_1-y_2|$的公因子，所以
$$\{\begin{array}{l}\text{g=p∣x1−x2∣}\\ \text{g=q∣y1−y2∣}\end{array}$$
因此,$|x_1-x_2|,|y_1-y_2|$可分为$g$个长为$p$的份,$g$的长度为$q$的份$(g\in Z)$,也即整个线段被分成$g$份,所以整数点有$g+1$个

证明：
为简化过程，不妨以线段的一个顶点重新建立坐标系
则可以设$\Delta x=x_2,\Delta y=y_2,d=gcd(\Delta x,\Delta y)$,再设$p,q$，与$x,y,d$之间有如下关系：
$$\{\begin{array}{l}{x}_2=\Delta x=pd\\ y_2=\Delta y=qd\end{array}$$

设点$(x_1+a,y_1+h)$是线段上一点，由三角形相似得：
$\frac{a}{\Delta x}=\frac{h}{\Delta y}=>a=h\frac{\Delta x}{\Delta y}=h\frac{pd}{qd}=h\frac{p}{q}$

点在线段上，所以h的取值范围为$[0,y_2]$（别忘了我们已经重新建立了坐标系）。我们可以枚举$h$为$[0,y_2]$中每一个整数，如果此时$a$也是整数，说明该点是一个整数点

而$a$是整数时，$p$是$q$的因子。因此，整数点的个数$sum$可以表示成（记号$[命题]$表示命题为真时该表达式值为1，否则值为0。$[q|ip]$即表示$q|ip$为真时值为1，否则值为0，下式即是对$i$从$0$到$y_2$的情况求和）：

$$\begin{array}{r}sum=\sum _{i=0}^{y_2}[q|ip]\end{array}$$

又$gcd(p,q)=1,pq$互质，所以$q|ip⟺q|i$，又因为$y_2=qd$，上式可化简为：

$$\begin{array}{rl}sum& =\sum _{i=0}^{y_2}[q|i]\\ & =\sum _{i=0}^{qd}[q|i]\end{array}$$

继续化简(当$i$可以为$q,2q,...,dq$时$[q|i]$为真，共d个)，
$$\begin{array}{rl}sum& =\sum _{i=0}^d\\ & =d+1\end{array}$$

例题 GCPC2017G Water Testing

最后给一道例题

You just bought a large piece of agricultural land, but you noticed that – according to regulations you have to test the ground water at specific points on your property once a year. Luckily the description of these points is rather simple. The whole country has been mapped using a Cartesian Coordinate System (where (0; 0) is the location of the Greenwich Observatory). The corners of all land properties are located at integer coordinates according to this coordinate system. Test points for ground water have to be erected on every point inside a property whose coordinates are integers.

Input The input consists of:
one line with a single integer $n(3\le n\le 100000)$, the number of corner points of your property;
n lines each containing two integers x and y $(-10^6\le x,y\le 10^6)$, thecoordinates of each corner.
The corners are ordered as they appear on the border of your property and the polygon described by the points does not intersect itself.

Output The number of points with integer coordinates that are strictly inside your property. Sample Input 1 Sample Output
4
0 0
0 10
10 10
10 0

AC代码:

#include <bits/stdc++.h>
using namespace std;

long long x[1000000], y[1000000];

int main () {
	int n; cin >> n;
	for (int i = 0; i < n; i++) cin >> x[i] >> y[i];

	long long a = 0, b = 0;
	for (int i = 1; i <= n; i++) {
		a += x[i%n]*y[i-1] - x[i-1]*y[i%n];
		b += __gcd(abs(x[i%n]-x[i-1]),abs(y[i%n]-y[i-1]));
	}
	
	cout << (abs(a)-b+2)/2 << endl;	
}

---

1. https://baike.baidu.com/item/%E7%9A%AE%E5%85%8B%E5%AE%9A%E7%90%86/7554085?fr=aladdin ↩︎

2. https://www.cnblogs.com/xiexinxinlove/p/3708147.html ↩︎