Gym 101873G - Water Testing - [皮克定理]
题意:
在点阵上,给出 N 个点的坐标(全部都是在格点上),将它们按顺序连接可以构成一个多边形,求该多边形内包含的格点的数目。
多边形的边不一定是水平或者竖直的,可以倾斜
题解:
首先由皮克定理
( S 是多边形面积,a 是多边形内部格点数目,b 是多边形边界上的格点数目)
可知我们需要求得:多边形的面积和边界上格点的数目
-
对于两端点 (x1,y1),(x2,y2) 都再格点上的一条线段,该线段上的格点数目为。
这很好理解,对于横坐标差值和纵坐标差值求得的最大公因数 g,相当于将横坐标差值分成 g 份,由于是整除,因此显然每份的左右端点都是整数,对于纵坐标也是同样的道理,由于是最大公因数,所以不可能再分更多份,因此 即求得两端点间最多能分成多少段由格点分割的线段,再加上 1 即整条线段上的格点数目。 -
对于格点按顺序给出的多边形,
设 且 O 为原点,则面积为
这个画个图模拟一下也非常容易理解。
代码:
#include<bits/stdc++.h> using namespace std; typedef long long ll; typedef pair<ll,ll> pll; const int maxn=1e5+10; int n; pll p[maxn]; inline ll gcd(ll m,ll n){return n?gcd(n,m%n):m;} int main() { cin>>n; for(int i=0;i<n;i++) scanf("%lld%lld",&p[i].first,&p[i].second); ll S2=0, b=0; for(int i=0;i<n;i++) { S2+=p[i].first*p[(i+1)%n].second-p[i].second*p[(i+1)%n].first; b+=gcd(abs(p[i].first-p[(i+1)%n].first),abs(p[i].second-p[(i+1)%n].second)); } cout<<(abs(S2)-b+2)/2<<endl; }
本文作者:kingwzun
本文链接:https://www.cnblogs.com/kingwz/p/16316664.html
版权声明:本作品采用知识共享署名-非商业性使用-禁止演绎 2.5 中国大陆许可协议进行许可。
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步