Reading comprehension HDU - 4990

思路：

哎，比较菜，不会找规律，就只好构造了一个奇数矩阵和偶数矩阵！  构造对应的base矩阵 就是 ji*ou 

判断如果偶数矩阵就再乘一下！

#include <iostream>
#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <queue>
using namespace std;
typedef long long ll;
ll n,m;
const int maxn=3;
struct Matrix//矩阵的类
{
  ll a[maxn][maxn];
  void init()           //将其初始化为单位矩阵
  {
    memset(a,0,sizeof(a));
    for(int i=1;i<maxn;i++)
      a[i][i]=1;
  }
}  ;

Matrix mul(Matrix a,Matrix b)  //(a*b)%mod  矩阵乘法
{
  Matrix ans;
  for(int i=1;i<maxn;i++)
    for(int j=1;j<maxn;j++)
    {
      ans.a[i][j]=0;
      for(int k=1;k<maxn;k++)
        ans.a[i][j]+=a.a[i][k]*b.a[k][j];
      ans.a[i][j]%=m;
    }
    return ans;
}

Matrix pow(Matrix a)
{
    Matrix res;
    for(int i=1;i<=2;i++)
    {
        for(int j=1;j<=2;j++)
        if(i==j)
            res.a[i][j]=1;
        else
            res.a[i][j]=0;
    }

    while(n)
    {
        if(n&1)
        {
            res=mul(a,res);
        }
        a=mul(a,a);
        n/=2;
    }
    return res;

}


int main()
{
    while(scanf("%lld%lld",&n,&m)!=EOF)
    {
        Matrix res,ji,ou,base,ans;
        ans.a[1][1]=1,ans.a[2][1]=1,ans.a[1][2]=ans.a[2][2]=0;
        ji.a[1][1]=2,ji.a[1][2]=1,ji.a[2][1]=0,ji.a[2][2]=1;
        ou.a[1][1]=2,ou.a[2][2]=1,ou.a[1][2]=ou.a[2][1]=0;
        if(n==1)
        {
            printf("%lld\n",n%m);
            continue;
        }
        ll k=n;
        n=(n-1)/2;
        base=mul(ji,ou);
        res=pow( base );

        if(!(k&1))
        {
            res=mul(ou,res);
        }
        ans=mul(res,ans);

        printf("%lld\n",ans.a[1][1]);

    }


}