hdu1005 矩阵快速幂

#include<iostream>

#include<cstdio>

#include<cstring>

#include<set>

#include<queue>

#include<vector>

#include<map>

using namespace std;

#define ll long long

int m[][],a,b,n,t[][];

int temp[][];

void power(int x)

{

    if(x==)

    {

        m[][]=m[][]=;

        m[][]=m[][]=;

        return ;

    }

    if(x==)

    {

        m[][]=a%;

        m[][]=b%;

        m[][]=;

        m[][]=;

        return ;

    }

    power(x>>);

    for(int i=;i<=;i++)

    {

        for(int j=;j<=;j++)

        {

            temp[i][j]=;

            for(int k=;k<=;k++)

            {

                temp[i][j]=(temp[i][j]+(m[i][k]*m[k][j])%)%;

            }

         }

    }

    if(x&)

    {

    for(int i=;i<=;i++)

    {

        for(int j=;j<=;j++)

        {

            m[i][j]=;

            for(int k=;k<=;k++)

            {

                m[i][j]=(m[i][j]+(temp[i][k]*t[k][j])%)%;

            }

         }

    }

}

    else

    {

            for(int i=;i<=;i++)

    {

        for(int j=;j<=;j++)

        {

            m[i][j]=temp[i][j];

         }

    }

    }

    for(int i=;i<=;i++)

    {

        for(int j=;j<=;j++)

        {

            m[i][j]=m[i][j]%;

        }

    }

    return ;

}

int main()

{

  while(scanf("%d%d%d",&a,&b,&n),a||b||n)

  {

        t[][]=a%;

        t[][]=b%;

        t[][]=;

        t[][]=;

      if(n==||n==)

      {

      printf("1\n");

      continue;

}

      else

      {

          power(n-);

          printf("%d\n",(m[][]+m[][])%);

      }

  }

  return ;

 } 

鸣谢https://blog.csdn.net/cambridgeacm/article/details/7703809