皮克定理

计蒜客--三角形内点

在一个平面坐标系中,我们可以选出三个不全在一条线上的点构成一个三角形。我们称一个在三角形内(不包含三角形的边上),横纵坐标皆为整数的点位这个三角形的内点。 对于一个由(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]得出当且仅当\frac{m}{n}*\frac{n}{a}\frac{m}{n}*\frac{2n}{a}、... 、\frac{m}{n}*\frac{a*n}{a}得到整数格点

因为此时(\frac{n}{a},\frac{m}{n}*\frac{n}{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;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值