HDU 5451 Best Solver （2015年沈阳赛区网络赛B题）

1.题目描述：点击打开链接

2.解题思路：第一次做这类通过矩阵求解整数部分的题目，学习了。首先是如何求(5+2sqrt(6))^n的整数部分，这里可以参考如下的博客链接：

点击打开链接

接下来，本题还有另外一个难点就是如何降低指数的次数。由于我们已经把问题变换为了通过矩阵快速幂来求解，而且根据线性递推数列取模后均为周期序列的结论，本题的矩阵也存在周期性，看到网上有的是根据群论的知识直接算出了周期，但是我比较倾向于用暴力的方式来找循环节，毕竟比赛的时候还是要坚持KISS原则的。

3.代码：

#include<iostream>
#include<algorithm>
#include<cassert>
#include<string>
#include<sstream>
#include<set>
#include<bitset>
#include<vector>
#include<stack>
#include<map>
#include<queue>
#include<deque>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<cctype>
#include<complex>
#include<functional>
#pragma comment(linker, "/STACK:1024000000,1024000000")

using namespace std;

#define rep(i,n) for(int i=0;i<(n);i++)
#define me(s) memset(s,0,sizeof(s))
#define pb push_back
#define lid (id<<1)
#define rid (id<<1|1)
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> P;

const int N=46340;

typedef int Matrix[2][2];

struct State
{
    int d,p;
    State(int d=0,int p=0):d(d),p(p){}
    bool operator<(const State&u)const
    {
        return d<u.d;
    }
};

int n,m;
Matrix A,B;
vector<State>g[N];

 void mul(Matrix A,Matrix B,Matrix&s)
{
    Matrix c;
    memset(c,0,sizeof(c));
    rep(i,2)
        rep(j,2)
            rep(k,2)
    c[i][j]=(c[i][j]+A[i][k]*B[k][j]%m)%m;
    memcpy(s,c,sizeof(c));
}

int solve()
{
    int start=0,cnt=1;
    while(!start)
    {
        int a=A[0][0],b=A[1][0];
        State u=State(b,cnt);
        for(int i=0;i<g[a].size();i++)
            if(g[a][i].d==b)
            start=g[a][i].p;
        mul(A,A,A);
        g[a].push_back(u);
        cnt++;
    }
    cnt-=start+1;
    memcpy(A,B,sizeof(B));
    for(int i=0;i<min(start,n);i++)
        mul(A,A,A);
    if(n>start)
    {
        n-=start;
        n=n%cnt;
        for(int i=0;i<n;i++)
            mul(A,A,A);
    }
    mul(A,B,A);
    return (2*A[0][0]-1+m)%m;
}

int main()
{
    int T;
    int rnd=0;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        A[0][0]=A[1][1]=5%m;
        A[1][0]=2%m,A[0][1]=12%m;
        memcpy(B,A,sizeof(A));
        for(int i=0;i<m;i++)g[i].clear();
        printf("Case #%d: %d\n",++rnd,solve());
    }
}