(3070)POJ

Amon-

于 2015-02-14 18:44:34 发布

阅读量347

点赞数

CC 4.0 BY-SA版权

分类专栏： ACM-ICPC

本文链接：https://blog.youkuaiyun.com/wang57389675/article/details/43818733

ACM-ICPC 专栏收录该内容

166 篇文章

订阅专栏

本文介绍了矩阵快速幂运算的基本概念，并通过实例展示了如何使用C++实现矩阵快速幂运算，解决特定数学问题。

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

#include<iostream>
#include<cstdio>
#include<string.h>
#include<string>
#include<set>
#include<algorithm>
#include<cmath>

#define ll __int64
#define MAX 1000009
using namespace std;

int n;
const int mod = 10000;
struct Matrix
{
int m[12][12];
};

Matrix mul(Matrix a,Matrix b) //¾ØÕó³Ë·¨
{
int i,j,k;
Matrix c;
for (i = 0 ; i < n; i++)
for (j = 0; j < n; j++)
{
c.m[i][j] = 0;
for (k = 0; k < n; k++)
c.m[i][j] += (a.m[i][k] * b.m[k][j])%mod;
c.m[i][j] %= mod;
}
return c;
}

Matrix pow(Matrix a,Matrix b,int x)//µÝ¹éÇó¿ìËÙÃÝ
{
while (x)
{
if (x & 1)
b = mul(b,a);
x = x >> 1;
a = mul(a,a);
}
return b;
}

int main()
{
int k;
int t;
Matrix A,B;
while(~scanf("%d",&k))
{
if(k==-1)break;
n = 2;
A.m[0][0] = 1;
A.m[1][0] = 1;
A.m[0][1] = 1;
A.m[1][1] = 0;
Matrix C = pow(A,A,k);
int ans = 0;
printf("%d\n",C.m[1][1]);
}
return 0;
}