hdu 1005

题目:http://acm.hdu.edu.cn/showproblem.php?pid=1005
/*
    本题的n很大不用快速幂肯定超时,根据递推公式f(n) =A*f(n-1)+Bf(n-2)构造矩阵
 */
#include<iostream>
#include <stdio.h>
#include<algorithm>

using namespace std;
int a,b,n;
struct Mat{
    int num[3][3];
};
Mat mult(Mat a,Mat b){
    Mat c;
    for(int i=1;i<=2;i++)
        for(int j=1;j<=2;j++){
            c.num[i][j]=0;
            for(int k=1;k<=2;k++)
                c.num[i][j]=(c.num[i][j]+a.num[i][k]*b.num[k][j])%7;
        }
    return c;
}
Mat my_pow(Mat x,int y){
    Mat ans;
    for(int i=1;i<=2;i++)
        for(int j=1;j<=2;j++){
            if(i==j)
                ans.num[i][j]=1;
            else
                ans.num[i][j]=0;
        }
    for(;y;y>>=1){
        if(y&1)
            ans=mult(ans,x);
        x=mult(x,x);
    } return ans;
}
int main(){
    while(scanf("%d%d%d",&a,&b,&n),(a||b||n)){
        Mat A;
        A.num[1][1]=a;
        A.num[1][2]=b;
        A.num[2][1]=1;
        A.num[2][2]=0;
        if(n ==1){
            printf("%d\n",a);//不判会超时
            continue;
        }
        if(n ==2){
            printf("%d\n",b);
            continue;
        }
        A=my_pow(A,n-2);
        printf("%d\n",(A.num[1][1]+A.num[1][2])%7);//虽然在上面mod7了,结果竟然忘了,WA了几次
    }
    return 0;
}

/*
     还有一个更好的办法,代码简洁,但是很不容易想到,利用数论的知识,结果mod7之后是有周
     期性的(7*7=49),结果是在0--6之间的数,直接打表就可以了。
*/
#include<iostream>
#include <stdio.h>
#include<algorithm>

using namespace std;
int a,b,n,i;
int arry[60];
int main(){
    while(scanf("%d%d%d",&a,&b,&n),(a||b||n)){
        arry[1]= arry[2] = 1;
        for(i =3;i < 60;i++)
            arry[i] = (a*arry[i-1] +b*arry[i-2])%7;
        n%=49;
        printf("%d\n",arry[n]);
    }
    return 0;   
}



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值