我决定将Haskell中Fibonacci程序的标准严格尾递归版本的性能与用C语言编写的性能进行比较,使用GMP进行比较,其结果是大到适合单词(在Haskell中我使用)多精度Integer
类型)。我将省略Haskell程序,因为这是关于C和GMP的问题。 C实现是这样的:
#include <stdio.h>
#include <stdlib.h>
#include <gmp.h>
void fib(unsigned int n){
mpz_t a, b, t;
mpz_init_set_ui(a, 0);
mpz_init_set_ui(b, 1);
mpz_init(t);
for(; n > 1; n --){
mpz_add(t, a, b);
mpz_set(a, b);
mpz_set(b, t);
}
//mpz_out_str(stdout, 10, b);
}
int main(int argc, char **argv){
unsigned long n, f;
if(argc != 2){
printf("Usage: fibc <number>\n");
return 1;
}
fib(atol(argv[1]));
return 0;
}
请注意,我注释掉了输出值的行,这大概花了一秒钟(我在Haskell版本中保留了这种行为)。
结果是:
time ./fibhs 1000000
./fibhs 1000000 5.77s user 0.05s system 99% cpu 5.831 total
time ./fibc 1000000
./fibc 1000000 11.19s user 0.00s system 100% cpu 11.194 total
我认为我必须使用GMP错误。任何人都可以在C代码中看到任何性能改进的可能性吗?
答案 0 :(得分:2)
打乒乓球。你有两个变量a和b以及一个临时t。您将结果添加到t中,然后将b复制到a和t复制到b。而是在将b添加到a并将b添加到b之间切换。最终结果是a或b,取决于n是奇数还是偶数。