蓝桥杯2013年第四届真题-黄金连分数

题目

黄金分割数0.61803… 是个无理数，这个常数十分重要，在许多工程问题中会出现。有时需要把这个数字求得很精确。

对于某些精密工程，常数的精度很重要。也许你听说过哈勃太空望远镜，它首次升空后就发现了一处人工加工错误，对那样一个庞然大物，其实只是镜面加工时有比头发丝还细许多倍的一处错误而已，却使它成了“近视眼”！！

言归正传，我们如何求得黄金分割数的尽可能精确的值呢？有许多方法。

比较简单的一种是用连分数：

              1
黄金数 = ---------------------
                    1
         1 + -----------------
                      1
             1 + -------------
                        1
                 1 + ---------
                      1 + ...

题解

高精度

---

高精度模板【全】

---

首先要看出规律，精度由低到高，黄金数可以表示为：$1$、$\frac{1}{2}$、$\frac{2}{3}$、$\frac{3}{5}$、$\frac{5}{8}$、……

斐波那契数列，分子为Fib(n-1)，分母为Fib(n)。

---

如果要精确到小数点后100位，显然得用高精度才行。

我们思考一下如何完成这个计算工作，整理出要用到哪些高精度运算。

1. 斐波那契数列的计算需要采用高精度加法。当分子和分母的位数足够多时才能保证小数部分的位数超过100位。

2. 分子除以分母需要采用高精度除法（本质是高精度减法）。注意这里的高精度除法不是“高精度除以低精度”而是“高精度除以高精度”，所以需要使用减法来模拟。

3. 不进行真分数的除法，进行假分数的除法。如果分子比分母小，那么我们除法计算会很难进行，不如给分子后面多加105个0，这样就可以保证分子比分母大，而且商为至少存在101位的整数。

4. 得到101位的商后，观察一下第101位需不需要四舍五入到第100位，这里可以不采用高精度加法实现加一的操作。

5. 输入答案的时候不要忘记加上前面的“0.”。

---

如果还是写不出来，这里给出一个整体思路模板，大家只需要填写一些关键步骤就行了。

string add (string s1, string s2) { // 计算 s1 + s2 的值并返回

	return res;
}

string div (string s1, string s2) { // 计算 s1 / s2 的商并返回

	return res;
}

string Fib (int n) { // 计算 Fibonacci(n) // 这个如何定义可以根据自己对题目的理解
	
}

main() {
	// 得到Fib(1000) 和 Fib(999) // 足够大是为了保证小数点后的精度
	// res = div (Fib(999) * 10^105, Fib(1000))
	// cout << res[0~100] // 前100位
	// 观察四舍五入情况，并加上"0."后填入
}

代码

#include<bits/stdc++.h>
using namespace std;
const int N = 1e3+10;

struct strings { // 保存黄金数的分子和分母 
	string zi, mu;
	strings (string s1, string s2) {
		zi = s1, mu = s2;
	}
};

string add (string s1, string s2) {
	string s = "";
	int a[N], b[N], c[N];
	
	memset (a, 0, sizeof a);
	memset (b, 0, sizeof b);
	memset (c, 0, sizeof c);
	reverse (s1.begin (), s1.end ());
	reverse (s2.begin (), s2.end ());
	
	int n = s1.size (), m = s2.size ();
	int len = max(n, m);
	
	for (int i = 0;i < n;i ++) a[i] = s1[i] - '0';
	for (int i = 0;i < m;i ++) b[i] = s2[i] - '0';
	
	for (int i = 0;i < len;i ++) {
		c[i] += a[i] + b[i];
		c[i + 1] += c[i] / 10;
		c[i] %= 10;
	}
	
	while (c[len]) len ++;
	for (int i = len-1;i >= 0;i --)
		s += c[i] + '0';
	return s;
}

strings Fib (int n) { // Fib(1) 返回1，Fib(2) 返回1/2，Fib(3) 返回2/3，…… 
	if (n == 0) return strings("1", "1");
	
	string s1 = "1";
	string s2 = "1";
	
	while (n --) {
		string t = s2;
		s2 = add (s1, s2);
		s1 = t;
	}
	
	return strings (s1, s2);
}

bool great (string s1, string s2) { // 判断s1是否大于s2 
	if (s1.size () == s2.size ()) return s1 >= s2;
	if (s1.size () > s2.size ()) return true;
	return false;
}

string sub (string s1, string s2) {
	int a[N], b[N], c[N];
	string res = "";
	
	memset (a, 0, sizeof a);
	memset (b, 0, sizeof b);
	memset (c, 0, sizeof c);
	reverse (s1.begin (), s1.end ());
	reverse (s2.begin (), s2.end ());
	
	int n = s1.size (), m = s2.size ();
	int len = max (n, m);
	
	for (int i = 0;i < n;i ++) a[i] = s1[i] - '0';
	for (int i = 0;i < m;i ++) b[i] = s2[i] - '0';
	
	for (int i = 0;i < len;i ++) {
		c[i] += a[i] - b[i];
		if (c[i] < 0) c[i] += 10, c[i+1] --;
	}
	
	while (!c[len-1]) len --;
	for (int i = len-1;i >= 0;i --)
		res += c[i] + '0';
	return res;
}

string div (string s1, string s2) {
	string res = "", r = "";
	int n = s1.size (), len = 0, c[N];
	memset (c, 0, sizeof c);
	for (int i = 0;i < n;i ++) {
		int ans = 0;
		r += s1[i];
		while (great(r, s2))
			r = sub (r, s2),
			ans ++;
		c[len ++] = ans;
	}
	
	int i = 0; // 这里没有控制的特别细节，比如存在0的情况等，只要是因为本题不会出现特殊的情况 
	while (i < len - 1 && !c[i]) i ++;
	for (;i < len;i ++) res += c[i] + '0';
	return res;
}

int main()
{
	strings s = Fib (1000);
	cout << s.zi << endl << s.mu << endl << endl;
	
	for (int i = 0;i < 105;i ++) s.zi += "0"; // 添加105个0相当于乘以10^105 
	
	string ans = div (s.zi, s.mu);
	cout << ans << ' ' << ans.size() << endl << endl;
	
	string res = "";
	for (int i = 0;i < 100;i ++) res += ans[i];
	if (ans[100] - '0' >= 5) res = add (res, "1");
	res = "0." + res;
	cout << res << endl; 
	
	string right_ans = "0.6180339887498948482045868343656381177203091798057628621354486227052604628189024497072072041893911375";
	if (res == right_ans) cout << "yes" << endl;
	return 0;
}