计蒜客--三角形内点
在一个平面坐标系中,我们可以选出三个不全在一条线上的点构成一个三角形。我们称一个在三角形内(不包含三角形的边上),横纵坐标皆为整数的点位这个三角形的内点。 对于一个由(0,0)、(n,m)、(p,0)作为顶点构成的三角形,请你设计程序求出他的内点数。
输入包括一行,包括三个用空格分隔的整数,分别为n,m,p(0 ≤ n < 32000,0 < m < 32000,0 < p < 32000)。
输出仅一个数,为这个三角形的内点的个数。
理解:
皮克定理是指一个计算点阵中顶点在格点上的多边形面积公式,该公式可以表示为2S=2a+b-2,其中a表示多边形内部的点数,b表示多边形边界上的点数,S表示多边形的面积。证明略。
由题目易得三角形面积为2S=p*m;
接下来就是计算三角形三条边上的格点。
x轴上的边上具有p+1个格点
然后对于其他两条边每条边上的格点分别是gcd(n,m)和gcd(abs(n-p),m);
以(0,0)--(n,m)边为例,如果我们求出gcd(n,m)=a
那么不难在[1,a]得出当且仅当、、... 、得到整数格点
因为此时(,)均为整数
由于两条边的交点(n,m)被重复计算所以边上的总格点数为b=p+gcd(n,m)+gcd(abs(n-p,m);
最后输出结果 ans=(2S+2-b)/2;即可
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<set>
#include<cmath>
#include<vector>
using namespace std;
typedef long long ll;
const int maxn=5*1e5+5;
const double eps=1e-8;
ll gcd(ll a,ll b){
return (b==0)?a:gcd(b,a%b);
}
int main()
{
ll n,m,p;
scanf("%lld%lld%lld",&n,&m,&p);
ll b=p+gcd(n,m)+gcd(abs(n-p),m);
//printf("%lld\n",b);
ll ans=p*m+2-b;
printf("%lld\n",ans/2);
return 0;
}