约数 AcWing 872. 最大公约数

约数 AcWing 872. 最大公约数

原题链接

AcWing 872. 最大公约数

算法标签

数学知识 约数 最大公约数 欧几里得算法

思路

求两个正整数 a 和 b 的 最大公约数 d
若b非零
则有 gcd(a,b) = gcd(b,a%b)
证明：
设a%b = a - kb 其中k = a/b(向下取整)
若d是(a,b)的公约数 则知 d|a 且 d|b 则易知 d|a-kb 故d也是(b,a%b) 的公约数
若d是(b,a%b)的公约数 则知 d|b 且 d|a-kb 则 d|a-kb+k*b = d|a 故而d|b 故而 d也是(a,b)的公约数
因此(a,b)的公约数集合和(b,a%b)的公约数集合相同 所以他们的最大公约数也相同 证毕
若b为零
0与任何数的最大公约数为该数本身， 返回a

代码

#include<bits/stdc++.h>
#define int long long
#define rep(i, a, b) for(int i=a;i<b;++i)
#define Rep(i, a, b) for(int i=a;i>b;--i)
#define x first
#define y second 
using namespace std;
const int mod = 1e9+7;
unordered_map<int, int> ump;
inline int read(){
   int s=0,w=1;
   char ch=getchar();
   while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
   while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();
   return s*w;
}
void put(int x) {
    if(x<0) putchar('-'),x=-x;
    if(x>=10) put(x/10);
    putchar(x%10^48);
}
signed main(){
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	int n=read();
	while(n--){
	    int a=read(), b=read();
	    printf("%lld\n", __gcd(a, b));
	}
	return 0;
}

最大公约数模板

int gcd(int a, int b){
    return b?gcd(b,a%b):a;
}

返回两数最大公约数相关函数

 __gcd(a, b)

原创不易
转载请标明出处
如果对你有所帮助 别忘啦点赞支持哈