c++ 基础数学问题

最大公约数

辗转相除法

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

更简洁写法

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

最小公倍数

a * b == 最大公约数 * 最小公倍数
因a*b可能溢出，采用 (a / d * b)，d是公因子可整除

分数基本运算

分数的表示

struct Fraction{
	int up, down;
}

分数的化简

Fraction reduction(Fraction res) {
	// (1)若分母为负数，令分子up和分母down为相反数
	if (res.down < 0) {
		res.up = -res.up;
		res.down = -res.down;
	}
	// (2)若分子up为0，令分母down为1
	if (res.up == 0)
		res.down = 1;
	// (3)求出分子绝对值与分母绝对值的最大公约数d，然后领分子分母同除以d
	else {
		int d = gcd(abs(res.up), abs(res.down));
		res.up /= d;
		res.down /= d;
	}
	return res;
}

分数的四则运算

同小学，Fraction res储存结果即可，都是整数运算用不上double

• 分数的除法，额外判断f2.up是否为0，若是则输出"Error"

分数的输出

void outputRes(Fraction r) {
	// (1)先化简
	r = reduction(r);
	// (2)整数
	if (r.down == 1) 
		printf("%lld", r.up);
	// (3)假分数，带分数形式输出【注意用绝对值比较up,down】
	else if (abs(r.up) > r.down)
		printf("%d %d/%d", r.up / r.down, abs(r.up) % r.down, r.down);
	// (4)真分数
	else
		printf("%d/%d", r.up, r.down);
}

• 因分数乘除易溢出，建议用long long

素数

素数基本判别

• 开根号写法【更安全】
  sqrt(1.0 * n)
• i平方写法【需用long long i, 仍不如法一安全，建议用法一】
  i * i < n

// 法一 sqrt
bool isPrime(int n) {
	if (n <= 1) return false;
	int sqr = (int)sqrt(1.0 * n);
	for (int i = 2; i < sqr; i++)
		if (n % i == 0)
			return false;
	return true;
}

素数表与素数筛法

质因子分解

考虑到2x3x5x7x11x13x17x19x23x29 已超过int范围
因而对int型，res数组只需开到10

#include<iostream>
#include<cmath>
using namespace std;
const int maxn = 100010;

bool isPrime(int n) { //判断n是否为素数
	if (n <= 1) return false;
	int sqr = (int)sqrt(1.0 * n);
	for (int i = 2; i <= sqr; i++)
		if (n % i == 0)
			return false;
	return true;
}

int prime[maxn], pNum = 0;
void Find_Prime() { //求素数表
	for (int i = 1; i < maxn; i++) {
		if (isPrime(i))
			prime[pNum++] = i;
	}
}

struct factor {
	int x, cnt; //x为质因子，cnt为其个数
}res[10];

int main() {
	//先打素数表
	Find_Prime(); //此句请务必记得写

	int n, num = 0; cin >> n; //num为n的不同质因子的个数
	if (n == 1) cout << "1=1"; //特判1的情况
	else {
		int sqr = (int)sqrt(1.0 * n); //n的根号
		//枚举根号n以内的质因子
		for (int i = 0; i < pNum && prime[i] <= sqr; i++) {
			if (n % prime[i] == 0) { //如果prime[i]是n的因子
				res[num].x = prime[i]; //记录该因子
				res[num].cnt = 0;
				while (n % prime[i] == 0) { //计算出质因子prime[i]的个数
					res[num].cnt++;
					n /= prime[i];
				}
				num++; //不同质因子个数加1
			}
			if (n == 1) break; //及时推出循环，节省点时间
		}
		if (n != 1) { //如果无法被根号n以内的质因子除尽
			res[num].x = n; //那么一定有一个大于根号n的质因子（超出素数表）
			res[num++].cnt = 1;
		}

		//按格式输出结果
		cout << n << '=';
		for (int i = 0; i < num; i++) {
			if (i > 0) cout << '*';
			cout << res[i].x;
			if (res[i].cnt > 1) {
				cout << '^' << res[i].cnt;
			}
		}
	}
	return 0;
}

参考资料：《算法笔记》