我正在研究一个非常大的数字(第100个数字)的斐波纳契算法。我需要让这个运行速度更快,但只需几秒钟我就没想完了。有没有办法让它更快?谢谢你的帮助。
#include <iostream>
using namespace std;
int main() {
string elem_major = "1";
string elem_minor = "0";
short elem_maj_int;
short elem_min_int;
short sum;
int length = 1;
int ten = 0;
int n;
cin >> n;
for (int i = 1; i < n; i++)
{
for (int j = 0; j < length; j++)
{
elem_maj_int = short(elem_major[j] - 48);
elem_min_int = short(elem_minor[j] - 48);
sum = elem_maj_int + elem_min_int + ten;
ten = 0;
if (sum > 9)
{
sum -= 10;
ten = 1;
if (elem_major[j + 1] == NULL)
{
elem_major += "0";
elem_minor += "0";
length++;
}
}
elem_major[j] = char(sum + 48);
elem_minor[j] = char(elem_maj_int + 48);
}
}
for (int i = length-1; i >= 0; i--)
{
cout << elem_major[i];
}
return 0;
}
答案 0 :(得分:6)
无论您对给定代码执行的优化程度如何,如果不更改基础算法,您只能稍微优化它。你的方法具有线性复杂性,对于大值,它会很快变慢。斐波那契数字的更快实现是通过在矩阵上进行矩阵exponentiation by squaring:
0 1
1 1
这种方法将具有对数复杂度,其中渐近更好。对此矩阵执行一些取幂,您会注意到n + 1
st Fibonacci数位于其右下角。
答案 1 :(得分:1)
我建议您使用类似cpp-bigint(http://sourceforge.net/projects/cpp-bigint/)的内容来表示您的大数字。 代码看起来像这样
#include <iostream>
#include "bigint.h"
using namespace std;
int main() {
BigInt::Rossi num1(0);
BigInt::Rossi num2(1);
BigInt::Rossi num_next(1);
int n = 100000;
for (int i = 0; i < n - 1; ++i)
{
num_next = num1 + num2;
num1 = std::move(num2);
num2 = std::move(num_next);
}
cout << num_next.toStrDec() << endl;
return 0;
}
我的机器上的快速基准测试:
time ./yourFib
real 0m8.310s
user 0m8.301s
sys 0m0.005s
time ./cppBigIntFib
real 0m2.004s
user 0m1.993s
sys 0m0.006s
答案 2 :(得分:0)
我会节省一些预先计算的点(特别是因为你正在寻找真正的大数字)
即说我保存了第500和第501个纤维数。那么如果有人问我600号纤维是什么?我会从502而不是从1开始计算。这样可以节省时间。
现在问题是你会保存多少分以及如何选择要保存的分数?
这个问题的答案完全取决于应用和可能的分布。