poj 3070 Fibonacci（矩阵快速幂，简单）

题目

还是一道基础的矩阵快速幂。

具体的居者的幂公式我就不明示了。

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

int num,mod=;
struct matrix
{
    int a[][];
};

matrix multiply(matrix x,matrix y)//矩阵乘法
{
       matrix temp;
       for(int i=;i<num;i++)
       {
               for(int j=;j<num;j++)
               {
                       int ans=;
                       for(int k=;k<num;k++)
                       {
                               ans+=((x.a[i][k]*y.a[k][j])%mod);
                       }
                       temp.a[i][j]=ans%mod;
               }
       }
       return temp;
}

matrix calc(matrix origin,matrix answ,int n)//矩阵快速幂——answ*origin^n
{
     while(n)
     {
             if(n%==)
                    answ=multiply(origin,answ);
             origin=multiply(origin,origin);
             n/=;
     }
     return answ;
}

int main()
{

  // freopen("in.txt", "r", stdin);
  // freopen("out.txt", "w", stdout);
   int n;
    while(scanf("%d",&n)!=EOF)
    {
        if(n==-)break;
        matrix origin= {,,
                        ,};
        matrix answ={,,
                     ,};
        num=;
        if(n==||n==)printf("%d\n",n);
        else
        {
            answ=calc(origin,answ,n-);
            printf("%d\n",answ.a[][]);
        }
    }
    return ;
}