C/C++ ACM GCD LCM 快速幂算法

GCD&LCM

原理

• GCD 辗转相除法
• LCM 利用公式 lcm(a,b)=a*b/gcd(a,b)

GCD&LCM模板

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

int lcm(int a,int b)
{
	return a/gcd(a,b)*b;
}

快速幂（取余）

模运算性质

• (a+b)%m=(a%m+b%m)%m
• (a*b)%m=(a%m*b%m)%m

原理

• 把 a^b 分解成 aⁿ¹ x aⁿ² x aⁿ³ …
• 当 b 为偶数次时可将 a^b 分解为 (a²)^b/2
• 当 b 为奇数次时可将 a^b 分解为 a^b-1 x a

快速幂模板

计算 a^b % c

#include <bits/stdc++.h>
using namespace std;

int main()
{
	long long a,b,c,ans;
	while(cin>>a>>b>>c)
	{
		ans=1;
		while(b>0)
		{
			if(b&1)ans=ans*a%c; //当剩余的b次幂是奇数次，将当前的a乘入答案使b剩余偶数次
			a=a*a%c; //分解偶数次
			b>>=1; //相当于b=b/2 将b右移一位
		}
		cout<<ans<<endl;
	}
	return 0;
}