第04章_函数和递归
Example_0401_组合数.cpp
Example_0403_孪生素数.cpp
Example_0406_用递归法计算阶乘.cpp
个人认为,这章的内容比较简单,而且习题也比较少,稍要思考的有第一个和第二个例子。
其中,第一个例子中需要去思考如何使用数学的方法去优化程序,而第二个程序则是用到数学上的定理~
// Example_0401_组合数.cpp
/**
* 题目名称:组合数
* 题目描述:输入非负整数m和n,输出组合数C(n,m)=n!/( (n-m)!*m! ), 其中 m <= n <= 20
**/
/**
* 下面注释部分代码为经典做法,直接使用一个函数代替这个阶乘,但是,这个方法容易在中间计算时数字范围溢出。
* 为了解决这一问题,我们采用了另一种思路做法,以将这个溢出问题解决。
*
* #include <stdio.h>
*
* inline double f(double n)
* {
* double i, m = 1;
* for (i = 1; i <= n; i++){
* m *= i;
* }
* return m;
* }
*
* int main()
* {
* double m, n;
* scanf("%lf%lf", &m, &n);
* printf("%.0lf\n", f(n) / (f(m) * f(n - m) ));
* return 0;
* }
**/
#include <iostream>
using namespace std;
long long cnm(int n, int m)
{
long s = 1;
if(m > n / 2){ // 先做最小优化,C(m, n)的结果与 C(m, m - n)相同
m = n - m;
} // 经过上面的if处理后,m的值就一定比n-m要小了。
// 对于1, 2, 3, ..., m, m+1, ..., n-m, ..., n, 公式中的n! 与 m! (n-m)!会有重叠部分
// 我们可以将n!与(n-m)!重复部分去掉,那么i的初值就可以是i = n-m+1,
// 嵌入到这个循环中的while循环,用于除数部分
// 为了防止溢出,一边计算(n-m+1)!,一边除掉1, 2, 3...m能除的数。
int k = 1;
for(int i = n - m + 1; i<= n; ++i){
s *= (long long) i;
// cout << "i = " << i << " s = " << s << endl;
while(k <= m && s % k == 0){
s /= (long long) k;
// cout << "进入了while循环:" << " k = " << k << " s = " << s << endl;
++k;
}
}
return s;
}
int main()
{
long long m, n;
cin >> n >> m;
long long result = cnm(n, m);
cout << result << endl;
return 0;
}
/**
* 观察输出过程:
* 7 3
* i = 5 s = 5
* 进入了while循环: k = 1 s = 5
* i = 6 s = 30
* 进入了while循环: k = 2 s = 15
* 进入了while循环: k = 3 s = 5
* i = 7 s = 35
* 35
**/
// Example_0403_孪生素数.cpp
/**
* 题目名称:孪生素数
* 题目描述:
* 如果n和n+2都是素数,则称它们是孪生素数。输入m,输出两个数均不超过m的最大孪生素数。 5 <= m <= 10000.
* 样例输入: 20
* 样例输出: 17 18
* 样例输入: 1000
* 样例输出: 881 883
**/
#include <iostream>
#include <cmath>
#include <assert.h>
using namespace std;
int is_prime(int x)
{
int i, m;
assert( x >= 0); // 宏assert用于错误检测,当表达式为零时,宏写错误到STDERR并退出程序执行。
// 在程序调用的"调试"中一般会使用这个宏,不过在写一些工程时,一般不会用到这个。
if(x == 1) { // 1不是素数
return 0;
}
static_cast <int> (m = floor(sqrt(x) + 0.5)); // 这里前面提到过,使用+0.5是为了减少误差,防止0.9999不能为1的情况。
for( i = 2; i <= m; ++i){ // 这里,还将判断的界限设置为原来数的开方,减少判断的次数,这里不使用i*i<=x是为了防止内存溢出。
if(x % i == 0){
return 0; // 有余数则非素数
}
}
return 1; // 全过测试通过,为素数。
}
// 补充:定理,如果n不是素数,则n有满足1<d<=sqrt(n)的一个因子d
// 证明:如果n不是素数,则由定义n有一个因子d满足1<d<n,
// 如果d大于sqrt(n), 则n/d是满足1<n/d<=sqrt(n)的一个因子。
int main()
{
int m;
cin >> m;
for(int i = m ; i >= 5; --i){
if (is_prime(i) && is_prime(i - 2) ){
cout << i - 2 << " " << i << endl;
break;
}
}
return 0;
}
// Example_0406_用递归法计算阶乘.cpp
/**
* 题目名称:用递归法计算阶乘
* 题目描述:阶乘的函数f(n) = n!可以定义为:f(0) = 1, f(n) = f(n - 1) * n ( n >= 1)
**/
#include <stdio.h>
int f(int n)
{
return n == 0 ? 1 : f(n - 1) * n;
}
int main()
{
printf("%d\n", f(3));
return 0;
}