【HDU 1005】 Number Sequence

【题目链接】

点击打开链接

【算法】

矩阵乘法快速幂，即可

【代码】

#include<bits/stdc++.h>
using namespace std;

int a,b,n;

struct Matrix
{
        int mat[][];
} ans;
inline void multipy(Matrix &a,Matrix b)
{
        int i,j,k;
        Matrix ans;
        ans.mat[][] = ans.mat[][] = ans.mat[][] = ans.mat[][] = ;
        for (i = ; i <= ; i++)
        {
                for (j = ; j <= ; j++)
                {
                        for (k = ; k <= ; k++)
                        {
                                ans.mat[i][j] = (ans.mat[i][j] + a.mat[i][k] * b.mat[k][j]) % ;
                        }
                }
        }
        a = ans;
}
inline Matrix power(Matrix a,int n)
{
        Matrix ans,b = a;
        ans.mat[][] = ; ans.mat[][] = ;
        ans.mat[][] = ; ans.mat[][] = ;
        while (n)
        {
                if (n & ) multipy(ans,b);
                n >>= ;
                multipy(b,b);
        }
        return ans;
}
int main()
{

        while (scanf("%d%d%d",&a,&b,&n) != EOF)
        {
                if (!a && !b && !n) break;
                if (n <= )
                {
                     printf("%d\n",);
                     continue;
                }
                ans.mat[][] = a; ans.mat[][] = b;
                ans.mat[][] = ; ans.mat[][] = ;
                ans = power(ans,n-);
                printf("%d\n",(ans.mat[][]+ans.mat[][])%);
        }

        return ;

}