poj3070Fibonacci

最新推荐文章于 2019-03-26 00:15:50 发布

mazheng1989

最新推荐文章于 2019-03-26 00:15:50 发布

阅读量414

点赞数

CC 4.0 BY-SA版权

分类专栏： ACM 文章标签： matrix output c

本文链接：https://blog.youkuaiyun.com/mazheng1989/article/details/7014641

ACM 专栏收录该内容

44 篇文章

订阅专栏

本文介绍了一种利用快速幂和矩阵乘法计算特定数学问题的方法。通过将问题转化为矩阵乘法，并采用快速幂技巧，实现了高效求解。代码示例展示了如何通过递归分解和二进制拆分来优化计算过程。

第一时间想到把n二进制拆分。

0ms通过的代码：

 #include<stdio.h>

struct Matrix22
{
	int d11;
	int d12;
	int d21;
	int d22;
   Matrix22 operator =(Matrix22 a)
	{
		d11=a.d11;
		d12=a.d12;
		d21=a.d21;
		d22=a.d22;
	}
    Matrix22 operator *(Matrix22 b)
	{
		Matrix22 c;
		c.d11=(d11*b.d11+d12*b.d21)%10000;
		c.d12=(d11*b.d12+d12*b.d22)%10000;
		c.d21=(d21*b.d11+d22*b.d21)%10000;
		c.d22=(d21*b.d12+d22*b.d22)%10000;
		return c;
	}	
}Matrix[32];
int main()
{

	freopen("input.txt","r",stdin);
	Matrix[0].d11=1;
	Matrix[0].d12=1;
	Matrix[0].d21=1;
	Matrix[0].d22=0;
	for(int i=1;i<32;i++)
	{
		Matrix[i]=Matrix[i-1]*Matrix[i-1];
	}
	
	int n;
	while(scanf("%d",&n)&&n!=-1)
	{
		Matrix22 output;
		output.d11=1;
		output.d12=0;
		output.d21=0;
		output.d22=1;
		int i=0;
		while(n>0)
		{
			if(n&1)
			{
			    output=output*Matrix[i];
			}
			i++;
			n=n>>1;
		}
		printf("%d\n",output.d12%10000);
	}
}