C语言名题100的小程序（2.10）(用矩阵算)

最新推荐文章于 2023-08-03 22:47:56 发布

nk_ysg

最新推荐文章于 2023-08-03 22:47:56 发布

阅读量203

点赞数

CC 4.0 BY-SA版权

文章标签： c 语言 matrix include

本文链接：https://blog.youkuaiyun.com/nk_ysg/article/details/4134872

本文介绍了一种利用矩阵快速幂算法高效计算斐波那契数的方法。通过递归实现矩阵乘法，有效降低了计算复杂度。适用于需要大量斐波那契数计算的应用场景。

#include <stdio.h> void Matrix_Mutiply(int n, long a, long b,long c, long d, long *a1, long *b1, long *c1, long *d1) { long tempa,tempb, tempc,tempd; if (1 == n) { *a1 = a; *b1 = b; *c1 = c; *d1 = d; return; } else if ( (n & 0x01) == 0) { Matrix_Mutiply(n/2,a*a+b*c,a*b+b*d,a*c+c*d,b*c+d*d, a1,b1,c1,d1); return; } else { Matrix_Mutiply(n/2,a*a+b*c,a*b+b*d,a*c+c*d,b*c+d*d, &tempa,&tempb,&tempc,&tempd); *a1 = a * tempa + b * tempc; *b1 = a * tempb + b * tempd; *c1 = c * tempa + d * tempc; *d1 = c * tempb + d * tempd; return; } } long fib(int n) { long a, b,c, d; if ( n <= 2) return 1; Matrix_Mutiply(n-2,1L,1L,1L,0L,&a,&b,&c,&d); return a + b; } int main(void) { int iCnt; for (iCnt = 2; iCnt <= 30; ++iCnt) printf("%ld/n",fibonacci(iCnt)); return 0; }