hdu 1005 Number Sequence（矩阵快速幂,找规律,模版更通用）

题目

第一次做是看了大牛的找规律结果，如下：

//显然我看了答案，循环节点是48，但是为什么是48，据说是高手打表出来的
#include<stdio.h>
int main()
{
    int f[],a,b,i,n;
    f[]=;f[]=;
    while(scanf("%d%d%d",&a,&b,&n)!=EOF)
    {
        if(a==&&b==&&n==)break;

        for(i=;i<;i++)
        {
            f[i]=(a*f[i-])%+(b*f[i-])%;
        }
        printf("%d\n",f[n%]%);
    }
    return ;
}

第二次做是学了矩阵快速幂，这是经典的矩阵快速幂简单题

//简单的矩阵快速幂
//f(1) = 1, f(2) = 1, f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7.
//注意origin矩阵中的第二行的a和b的位置摆放
#include<stdio.h>
#include<string.h>
int num=,mod=;
struct matrix
{
    int a[][];
};
matrix multiply(matrix x,matrix y)//矩阵乘法
{
    matrix temp;
    memset(temp.a,,sizeof(temp.a));
    for(int i=;i<num;i++)
    {
        for(int k=;k<num;k++)
        {
            for(int j=;j<num;j++)
            {
                temp.a[i][j]=(temp.a[i][j]+x.a[i][k]*y.a[k][j])%mod;
            }
        }
    }
    return temp;
}
matrix calc(matrix a,int n)//矩阵快速幂——a^n
{
    if(n==)return a;
    matrix e;
    for(int i=;i<num;i++)
        for(int j=;j<num;j++)
            e.a[i][j]=(i==j);

    while(n)
    {
        if(n&)
            e=multiply(e,a);
        n>>=;
        a=multiply(a,a);
    }
    return e;
}
int main()
{
    int n,a,b;
    while(scanf("%d%d%d",&a,&b,&n)!=EOF)
    {
        if(a==&&b==&&n==)break;
        matrix origin= {,};
        origin.a[][]=b;origin.a[][]=a;
        matrix answ={,,
                     ,};
        if(n>)
            answ=multiply(calc(origin,n-),answ);
        printf("%d\n",answ.a[][]);
    }
    return ;
}