欧几里得算法-最大公约数，最小公倍数；大整数运算

最大公约数：欧几里得算法

原理：a，b的公约数就是b，a%b的公约数（a>b)
辗转相除法。

int gcd(int a, int b)
{
	if (b == 0) return a;
	else return gcd(b, a % b);
}
/* 如果b是0    那!b表示 非0
   如果b不是0  那!b表示为0 */
int gcd(int a, int b)
{
	return !b a:gcd(b,a%b);
}
int main()
{
	int a, b;
	while (cin>>a>>b && a!=EOF)
	{
		a > b ? cout << gcd(a, b) : cout << gcd(b, a);
		cout<<endl;
	}
}

最小公倍数：基于最大公约数。

a/gcd(a,b) *b就是最小公倍数。
质因子分解：最多存在一个大于sqrt（n）的质因子。
for (int i = 0; i < num && prime[i]<=sqrt(n); i++)

大整数运算

超过数据类型存储的运算
输入用字符串读入，然后转存为数组。

减法 注意最后1234-1231 = 0003 要去0 减len

bign sub(bign a, bign b)
{
	bign c;
	int sum = 0, in = 0;
	for (int i = 0; i < a.len || i < b.len; i++)
	{
		if (a.d[i] < b.d[i]) {
			--a.d[i + 1];
			a.d[i] += 10;
		}
		c.d[c.len++] = a.d[i] - b.d[i];
	}
	/* 要去除高位的0 1234-1231 = 0003  c.len至少=1*/
	while (c.len - 1 > 0 && c.d[c.len - 1] == 0)--c.len;
	return c;
}

除法：可能有余数，可以记录下来。如果需要最后的余数的话
同样前面可能都是0

bign divide(bign a, int b, int &r)//可能会产生的余数
{
	bign c;
	int sum;
	c.len = a.len;//最多位数
	for (int i = a.len - 1; i >= 0; i--)
	{
		r = r * 10 + a.d[i];
		if (r < b) c.d[i] = 0;
		else {
			c.d[i] = r / b;
			r = r % b;
		}
	}
	while (c.len - 1 > 0 && c.d[c.len - 1] == 0)--c.len;
	return c;
}