一些数学问题——公约数公倍数、质数判断、质因数分解

最大公约数与最小公倍数

最大公约数

辗转相除法：

int gcd(int a,int b){
	if(b==0) return a;
	else return gcd(b,a%b);
}

最小公倍数

在求出a b的最大公约数d之后，可以马上得到a b的最小公倍数是ab/d

素数

素数的判断

判断一个数字从2到其根号是否能被它整除

bool isPrime(int n) {
	if (n <= 1)return false;//注意1是不算质数的
	int sqr = (int)sqrt(1.0*n);
	for (int i = 2; i <= sqr; i++) {
		if (n%i == 0)return false;
	}
	return true;
}

可以基于这个配合枚举列出1-100的所有素数：

#include<iostream>
using namespace std;

const int maxn = 101;
int prime[maxn], pNum = 0;//素数表和素数个数
//bool p[101] = { 0 }; 

bool isPrime(int n) {
	if (n <= 1)return false;//注意1不是素数
	int sqr = sqrt(n);
	for (int i = 2; i <= sqr; i++) {
		if (n%i == 0)return false;
	}
	return true;
}

void find_prime() {
	for (int i = 1; i < maxn; i++) {
		if (isPrime(i)) {
			prime[pNum++] = i;
			//p[i] = true;
		}
	}
}

int main() {
	find_prime();
	for (int i = 0; i < pNum; i++) {
		cout << prime[i] << " ";
	}
	return 0;
}

时间复杂度为O(n根号n）在10^5之内可以接收。更高效的算法：筛法 O(nloglogn)
例子：

筛选法

#include<iostream>
using namespace std;

const int maxn = 101;
int prime[maxn], pNum = 0;//素数表和素数个数
bool p[maxn];//默认是素数


void find_prime() {
	for (int i = 2; i < maxn; i++) {
		if (p[i] == true) {
			//是素数
			prime[pNum++] = i;
			for (int j = i + i; j < maxn; j += i) {
				p[j] = false;//筛掉i的倍数
			}
		}
	}
}

int main() {
	fill(p, p + maxn, true);
	find_prime();
	for (int i = 0; i < pNum; i++) {
		cout << prime[i] << " ";
	}
	return 0;
}

质因子分解

将一个整数n写成一个或者多个质数的乘积

struct factor{
	int x,cnt;
}fac[10];

为什么十个就够了 因为2 3 5 7 11 13 17 19 23 29这些数乘起来就已经超过int范围了。

结论：对一个正整数n，如果存在1和本身之外的因子，那么一定是在sqrt(n)的左右成对出现。如果它在[2,n]内存在质因子，要么这些质因子全部小于sqrt(n)，要么只存在一个大于sqrt(n)的质因子，其余质因子全部小于等于sqrt(n)

if(n%prime[i]==0){
	//如果prime[i]是n的因子
	fac[num].x=prime[i];
	fac[num].cnt=0;
	while(n%prime[i]==0){
		fac[num].cnt++;
		n/=prime[i];
	}
	num++;
}

if(n!=1){
 //如果无法被根号n以内的质因子除尽
 fac[num].x=n;//一定有一个大于根号n的质因子
 fac[num++].cnt=1;
}

综合应用：给一个长整型（64位）求其分解的因式

#include<vector>
#include<iostream>
using namespace std;

const int maxn = 2001;//判断上限
int prime[1000];
int primeNum = 0;
bool p[maxn];

struct factor {
	int x, cnt;
};
vector<factor> facList;

//筛选法求质数表,注意筛选法是每个数都要看
void getPrime() {//生成2000以内质数
	fill(p, p + maxn, true);//默认都是质数
	for (int i = 2; i < maxn; i++) {
		if (p[i] == true) {
			prime[primeNum++] = i;
			for (int j = i + i; j < maxn; j += i) {
				p[j] = false;//i的倍数都不是质数
			}
		}
	}
}

//求因式
void get_prime_re(long long n) {
	factor tempFac;
	if (n <= 2) {
		tempFac.x = n;
		tempFac.cnt = 1;
		facList.push_back(tempFac);
		return;
	}
	
	for (int i = 0; i < primeNum; i++) {
		if (n%prime[i] == 0) {//如果可以呗prime[i]整除
			tempFac.x = prime[i];//指定质因子种类
			tempFac.cnt = 0;//初始化个数
			while (n%prime[i] == 0) {//算这个质因子的个数
				tempFac.cnt++;
				n /= prime[i];
			}
			facList.push_back(tempFac);//加到总表
		}
		if (n == 0)break;
		
	}
}

int main() {

	//输入
	long long inputNum;
	cin >> inputNum;

	getPrime();
	get_prime_re(inputNum);
	cout << inputNum << "=";
	for (int i = 0; i < facList.size(); i++) {
		if (facList[i].cnt != 0) {
			if (facList[i].cnt == 1) {
				cout << facList[i].x;
			}
			else {
				cout << facList[i].x << "^" << facList[i].cnt;
			}
			if (i != facList.size() - 1)cout << "*";
		}
	}
	return 0;
}