HDU - 2855 Fibonacci Check-up 矩阵快速幂

题目大意:这里写图片描述
F函数是fibonacci函数,F(0) = 0,F(1) = 1

解题思路:这题我也不会推,打表可以得到结论。。。
G(n) = 3 * G(n-1) - G(n-2),G函数表示这里写图片描述
这题要注意最后求得的可能是负数,所以最后要处理一下

#include<cstdio>
typedef long long ll;
const int N = 2;
struct Matrix {
    int mat[N][N];
}A, B, tmp;
int n, m;
void init() {
    A.mat[0][0] = 3;
    A.mat[0][1] = 1;
    A.mat[1][0] = -1;
    A.mat[1][1] = 0;

    B.mat[0][0] = B.mat[1][1] = 1;
    B.mat[0][1] = B.mat[1][0] = 0;
}

Matrix matMul(Matrix x, Matrix y) {
    for(int i = 0; i < N; i++)
        for(int j = 0; j < N; j++) {
            tmp.mat[i][j] = 0;
            for(int k = 0; k < N; k++)
                tmp.mat[i][j] = (tmp.mat[i][j] + x.mat[i][k] * y.mat[k][j]) % m;
        }
    return tmp;
}

void solve() {
    while(n) {
        if(n & 1)
            B = matMul(B,A);
        A = matMul(A,A);
        n >>= 1;
    }
}

int main() {
    int test;
    scanf("%d", &test);
    while(test--) {
        scanf("%d%d", &n, &m);
        if(n == 0) {
            printf("0\n");
        }
        else {
            init();
            n--;
            solve();
            printf("%d\n", (B.mat[0][0] + m ) % m);
        }
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值