[英]Fibonacci numbers with GMP lib
編輯:解決了。 ans var需要在每次迭代中設置為0.這是我忽略的一個錯誤。一個mpz_set_ui(ans,0); 在每個循環開始時解決了這個問題。感謝Oli Charlesworth的努力!
我使用GMP庫在C中編碼斐波那契函數,我使用矩陣乘法算法,但出於某些意想不到的原因,我得到了錯誤的結果(我堅信算法是對的)。它必須與Gmp有關結構和inits.Any想法?
/*********************************************************
Find n th Fibonacci number.Matrix algorithm.
+-- --+ +-- --+
|f(n+1) f(n) | | 1 1 | ^n
| | == | |
|f(n) f(n-1) | | 1 0 |
+-- --+ +-- --+
*********************************************************/
int
fastfib(mpz_t result,int n)
{
int cons[2][2]={{1,1},{1,0}}; //"constant" matrix
mpz_t mat[2][2]; //actual matrix holding fibonacci numbers
mpz_t ans; //holds f(n+1)
mpz_t return_val; //used for calculations
mpz_t temp; //used for calculations as well
//initialize (automatically set to zero)
mpz_init(ans);
mpz_init(return_val);
mpz_init(temp);
mpz_init(mat[0][0]);
mpz_init(mat[0][1]);
mpz_init(mat[1][0]);
mpz_init(mat[1][1]);
//start with n=1
mpz_set_ui(mat[0][0],1);
mpz_set_ui(mat[1][0],1);
mpz_set_ui(mat[0][1],1);
mpz_set_ui(mat[1][1],0);
//some trivial cases
if(n==0){
mpz_set_ui(result,0);
return 0;
}
if(n==1){
mpz_set_ui(result,1);
return 0;
}
if(n==2){
mpz_set_ui(result,1);
return 0;
}
n--;
while(n>1){
//fib[n+1]
//ans=mat[0][0]*cons[0][0]+mat[0][1]*cons[1][0];
mpz_set_ui(ans,0);
mpz_mul_ui(temp,mat[0][0],cons[0][0]);
mpz_add(ans,ans,temp);
mpz_mul_ui(temp,mat[0][1],cons[1][0]);
mpz_add(ans,ans,temp);
//update matrix
mpz_set(mat[1][1],mat[1][0]); //mat[1][1]=mat[1][0];
mpz_set(mat[1][0],mat[0][0]); //mat[0][1]=mat[1][0]=mat[0][0];
mpz_set(mat[0][1],mat[0][0]);
mpz_set(mat[0][0],ans); //mat[0][0]=ans;
n--;
}
//clear vars
mpz_clear(ans);
mpz_clear(return_val);
mpz_clear(temp);
mpz_set(result,mat[0][0]);
return 0;
}
一些調試信息
results for n=2 // expected
2 1 // 2 1
1 1 // 1 1
results for n=3
5 2 // 3 2
2 1 // 2 1
results for n=4
12 5 // 5 3
5 2 // 3 2
results for n=5
29 12 // 8 5
12 5 // 5 3
results for n=6
70 29 // 13 8
29 12 // 8 5
70
解決了。 ans var需要在每次迭代中設置為0.這是我忽略的一個錯誤。一個mpz_set_ui(ans,0); 在每個循環開始時解決了這個問題。感謝Oli Charlesworth的努力!
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.