《算法笔记》读书记录DAY_9

CHAPTER_5  数学问题入门

5.1简单数学

有一类问题，他们不涉及很深的算法，但和数学息息相关。应对这类问题一般只要掌握基本的数理逻辑即可。

题目：

给定一个各位数字不完全相同的四位正整数N，先把N的四位数字按非递增排序得到数字B，再按非递减排序得到数字A，用B-A得到一个新数字。再将新数字做同样操作，一直重复这样做，很快就会停在一个神奇的数字6174。这个数叫Kaprekar常数，被称作数字黑洞。

现请编程演示其运算过程。

输入格式：

输入一个正整数N（N属于区间(0,10000)）  

输出格式：

如果N的四位数字全部相等，则在一行内输出“N-N=0000”；否则将计算的每一步在一行内输出，知道6174作为差值出现

输入样例： 

6767

输出样例：

7766-6677=1089

9810-0189=9621

9621-1269=8352

8532-2358=6174

参考代码： 

#include<stdio.h>
#include<algorithm>

using namespace std;

bool cmp(int a, int b) {
	if(a>b)
		return 1;
	else
		return 0;
}

int main() {
	int n,maxN,minN;     //n存放输入数字，maxN存放排序后的递减序列 ，minN存放排序后的递增序列  
	char num[5];         //num数组用于存放由整数转换而来的字符串 
	scanf("%d",&n);
	do{
		sprintf(num,"%d",n);      //将n转换成字符串便于操作
		sort(num,num+4);          //升序排列
		sscanf(num,"%d",&minN);    //升序结果转化成数字给minN
		sort(num,num+4,cmp);      //降序排列
		sscanf(num,"%d",&maxN);    //降序结果转化成数字给maxN
		n=maxN-minN;
		printf("%04d - %04d = %04d\n",maxN,minN,n);	   //控制格式输出每次运算结果	
	}while(n&&n!=6174);
	return 0;
}

5.2.1最大公约数

题目：

输入两个正整数，求其最大公约数

输入格式：

测试数据有多组，每组输入两个正整数

输出格式：

对于每组输入，请输出其最大公约数

样例输入： 

49 14 

25 15

样例输出： 

7 

5

思路： 

求最大公约数最常用的是欧几里得算法（辗转相除法），算法基于如下定理：

设a、b为正整数(假设b<a)，gcd()表示a与b的最大公约数，那么有gcd(a,b)=gcd(b,a%b)。

令c=a%b，此时一定有c<b。那么有gcd(a,b)=gcd(b,c)=gcd(c,b%c)。如此地推下去最终会得到gcd(a,b)=gcd(b,c)=gcd(c,b%c)=...=gcd(x,0)。此时x为最大公约数。

读者不妨用上面的样例代入计算，体验这个过程。

参考代码：

#include <iostream>
using namespace std;

int gcd(int a, int b) {
	if (a < b) {
		swap(a, b);
	}
	int c;
	while (b) {
		c = a % b;
		a = b;
		b = c;
	}
	return a;
}

int main() {
	int m, n;
	while (scanf("%d%d", &m, &n) != EOF) {
		cout<<gcd(m, n)<<endl;
	}
	return 0;
}

PS：上面为手动实现最大公约数函数gcd()。实际上在C++库<algorithm>中有函数_gcd()用于求最大公约数，但该函数并不是标准库函数，在某些平台可能不支持。

5.2.2最小公倍数

两个数的最小公倍数一般用lcm(a,b)表示。假设a与b的最大公约数为d，最小公倍数为k。那么有k=a*b/d。因此求最小公倍数的问题可以转换为求最大公约数。

PS：由于a*b在实际运算中可能会溢出，因此编程时最好的写法是a/d*b。

5.3分数的表示与运算

分数的表示：

在处理数学问题中我们经常要用到分数。我们将分数写成假分数的形式。

struct Fraction {

        int up,down;        //up为分子，down为分母

};

分数的化简：

化简步骤分为以下三步

（1）如果分母down为负数，那么令up和down都变为相反数(保证分母恒为正，负号加在分子上)；

（2）如果分子up为0，那么令down为1；

（3）约分，求出p=gcd(|up|,|down|)，然后令up和down都除以p。

Fraction reduction(Fraction result) {
	if(result.down<0) {
		result.up=-result.up;
		result.down=-result.down;
	}
	if(result.up==0) {
		result.down=1;
	}
	else {
		int p=gcd(abs(result.up),abs(result.down));
		result.up/=p;
		result.down/=p;
	}
	return result;
}

分数的加法：

计算公式

Fraction add(Fraction f1,Fraction f2) {
	Fraction result;
	result.up=f1.up*f2.down+f1.down*f2.up;
	result.down=f1.down*f2.down;
	result reduction(result);
}

分数的减法：

计算公式

Fraction minu(Fraction f1,Fraction f2) {
	Fraction result;
	result.up=f1.up*f2.down-f1.down*f2.up;
	result.down=f1.down*f2.down;
	result reduction(result);
}

分数的乘法：

计算公式

Fraction multi(Fraction f1,Fraction f2) {
	Fraction result;
	result.up=f1.up*f2.up;
	result.down=f1.down*f2.up;
	result reduction(result);
}

分数的除法：

计算公式

Fraction divide(Fraction f1,Fraction f2) {
	Fraction result;
	result.up=f1.up*f2.down;
	result.down=f1.down*f2.up;
	result reduction(result);
}