用程序表示的辗转相除法——求最大公约数

线段上格点的个数

  给定平面上两个格点 P 1 = ( x 1 , y 1 ) P_1=(x_1,y_1) P1=(x1,y1) P 2 = ( x 2 , y 2 ) P_2=(x_2,y_2) P2=(x2,y2),线段P1P2上,除P1和P2以外一共有几个格点?

限制条件

  • − 1 0 9 ≤ x 1 , x 2 , y 1 , y 2 ≤ 1 0 9 -10^9 \le x_1,x_2,y_1,y_2 \le 10^9 109x1,x2,y1,y2109

输入

P1 = (1,11),P2  = (5,3)

输出

(2,9)、(3,7)、(4、5)

答案是 ∣ x 1 − x 2 ∣ |x_1-x_2| x1x2 ∣ y 1 − y 2 ∣ |y_1-y_2| y1y2的最大公约数-1(当做一个结论来记吧,我也不知道咋来的)
 
主题

计算最大公约数——辗转相除法

g c d ( a , b ) gcd(a,b) gcd(a,b)是计算自然数a和b的最大公约数的函数,a除以b得到的商和余数分别为p和q。因为 a = b × p − q a=b\times p-q a=b×pq,所以 g c d ( b , q ) gcd(b,q) gcd(b,q)既整除a又整除b(假设 g c d ( b . q ) gcd(b.q) gcd(b.q)求得的结果记作r,即b与q的最大公约数是r,即可表示为 b % r = = 0 , q % r = = 0 b\%r==0 , q\%r==0 b%r==0,q%r==0),又 a = b × p − q a=b\times p-q a=b×pq,所以 a ÷ r = = b × p − q r → b r × p − q r a\div r == \frac{b\times p - q}{r} \rightarrow \frac{b}{r}\times p-\frac{q}{r} a÷r==rb×pqrb×prq,所以得出 g c d ( b . q ) gcd(b.q) gcd(b.q)的结果既能整除a又能整除b,也就能整除 g c d ( a , b ) gcd(a,b) gcd(a,b)。反之,因为 q = a − b × p q = a-b\times p q=ab×p,同理可证 g c d ( a , b ) gcd(a,b) gcd(a,b)整除 g c d ( b , q ) gcd(b,q) gcd(b,q)。因此可以知道 g c d ( a , b ) = = g c d ( b , a % b ) gcd(a,b)==gcd(b,a\%b) gcd(a,b)==gcd(b,a%b)。不断这样操作下去,由于gcd的第二个参数总是不断缩小的,最终会得到 g c d ( a , b ) = = g c d ( c , 0 ) gcd(a,b)==gcd(c,0) gcd(a,b)==gcd(c,0)。0和c的最大公约数是c,所以 g c d ( c , 0 ) = = c gcd(c,0)==c gcd(c,0)==c,这样就计算出了 g c d ( a , b ) gcd(a,b) gcd(a,b)。程序实现如下:
 

int gcd(int a,int b)
{
	if(b==0) return a;
	return gcd(b,a%b);
}
// 或者是这样
int gcd(int a,int b)
{
	return b?gcd(b,a%b):a;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值