2018.09.25 poj3070 Fibonacci（矩阵快速幂）

SC.ldxcaicai

于 2018-09-25 09:47:41 发布

阅读量151

点赞数 2

CC 4.0 BY-SA版权

分类专栏： # 矩阵快速幂文章标签：数论

本文链接：https://blog.youkuaiyun.com/dreaming__ldx/article/details/82834905

数学专题同时被 2 个专栏收录

72 篇文章

订阅专栏

矩阵快速幂

15 篇文章

订阅专栏

本文分享了在POJ平台上解决矩阵快速幂题目的一次实战经验。通过代码示例详细展示了如何利用矩阵快速幂方法高效求解特定类型的问题，特别强调了在实际编程中调整通用代码库的重要性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门
矩阵快速幂板题，写一道来练练手。
这一次在poj做题总算没忘了改万能库。
代码：

#include<iostream>
#include<cstdio>
#define mod 10000
#define A a[0][0]
#define B a[0][1]
#define C a[1][0]
#define D a[1][1]
using namespace std;
int n;
struct Matrix{int a[2][2];Matrix(){A=0,B=C=D=1;}};
inline Matrix operator*(Matrix a,Matrix b){
	Matrix ret;
	ret.A=(a.A*b.A+a.B*b.C)%mod;
	ret.B=(a.A*b.B+a.B*b.D)%mod;
	ret.C=(a.C*b.A+a.D*b.C)%mod;
	ret.D=(a.C*b.B+a.D*b.D)%mod;
	return ret;
}
inline int ksm(){
	if(n==0)return 0;
	if(n==1)return 1;
	n-=2; 
	Matrix x,ret;
	while(n){
		if(n&1)ret=ret*x;
		x=x*x,n>>=1;
	}
	return ret.D;
}
int main(){
	while(scanf("%d",&n)&&~n)printf("%d\n",ksm());
	return 0;
}