poj 3070 Fibonacci(矩阵运算)

最新推荐文章于 2023-01-07 14:41:09 发布

acmj1991

最新推荐文章于 2023-01-07 14:41:09 发布

阅读量597

点赞数

分类专栏：矩阵文章标签： struct ini c

本文链接：https://blog.csdn.net/acmj1991/article/details/6943661

版权

矩阵专栏收录该内容

3 篇文章 0 订阅

订阅专栏

题意：求第n个Fibonacci数对10000取余

思路：由于给定的n非常大，所以不能够打表，而矩阵能计算两个数的加减

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;

#define A 10000
int n=2;
struct node{
	int ke[3][3];
};
node a,su,unit;

void print(node t)
{
	for(int i=0;i<n;i++){
		for(int j=0;j<n;j++)
			printf("%d ",t.ke[i][j]);
		printf("\n");
	}
}
node mul(node b,node c)
{
	node temp;
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
		{
			temp.ke[i][j]=0;
			for(int h=0;h<n;h++)
				temp.ke[i][j]+=b.ke[i][h]*c.ke[h][j];
			temp.ke[i][j]%=A;
		}
	return temp;
}
node get_sum(int k)
{
	node p=a,q=unit;
	while(k!=1)
	{
		if(k&1)
		{
			k--;
			q=mul(q,p);
		}
		else
		{
			k/=2;
			p=mul(p,p);
		}
	}
	return mul(p,q);
}

int main()
{
	int n;
	a.ke[0][0]=0,a.ke[0][1]=1;
	a.ke[1][0]=1,a.ke[1][1]=1;
	unit.ke[0][0]=1,unit.ke[1][1]=1;
	su.ke[0][0]=1,su.ke[0][1]=1;
	while(scanf("%d",&n)&&n!=-1)
	{
		if(n==0)printf("0\n");
		else if(n<=2)printf("1\n");
		else{
			node p=mul(su,get_sum(n-2));
			printf("%d\n",p.ke[0][1]%A);
		}
	}
}