【第22期】观点:IT 行业加班,到底有没有价值?

递推+矩阵快速幂

原创 2016年08月28日 20:54:53

题目描述
由于长度为一的方块只有一种方案,长度为二的有四种方案(不包含长度为一中的情况),长度为三的有两种方案(不包含长度为二中的情况),得递推式:
f[i] = f[i-1] + f[i-2] * 4 + f[i-3] * 2;
由于n <= 10 ^ 18, 考虑使用矩阵快速幂,用如下矩阵存状态:
f[i], 0, 0
f[i+1], 0, 0
f[i+2], 0, 0
用如下矩阵转移状态:
0, 1, 0
0, 0, 1
2, 4, 1
代码:

#include <cstdio>
#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstring>

#define For(i,j,k) for(LL i = j;i <= k;i ++)
const int Mod = 1000000007;
typedef long long LL;

struct Matrix{
    const static int N = 3;
    LL M[N][N];

    Matrix(int c = 0){
        memset(M, 0, sizeof(M));
        if(c == 1)
            For(i,0,N-1)
                M[i][i] = 1;
        if(c == 2){
            M[0][1] = M[1][2] = M[2][2] = 1;
            M[2][0] = 2, M[2][1] = 4;
        }
        if(c == 3){
            M[0][0] = M[1][0] = 1;
            M[2][0] = 5;
        }
    }

    Matrix operator * (const Matrix &B) const{
        Matrix C;
        For(i,0,N-1)
            For(j,0,N-1)
                For(k,0,N-1)
                    C.M[i][j] = (C.M[i][j] + M[i][k] * B.M[k][j]) % Mod;
        return C;
    }

    Matrix operator ^ (LL exp) const{
        Matrix Ans(1), T = *this;
        while(exp){
            if(exp & 1) Ans = Ans * T;
            T = T * T;
            exp >>= 1;
        }
        return Ans;
    }

    void print() const{
        For(i,0,N-1)
            For(j,0,N-1)
                printf("%lld%c", M[i][j], j == N - 1 ? '\n' : ' ');
    }

};

LL n;
int main(){
    scanf("%lld", &n);
    Matrix Ans = (Matrix(2) ^ n) * 3;
    printf("%lld\n", Ans.M[0][0]);
    return 0;
}
版权声明:本文为博主原创文章,转载需注明原文地址。 举报

相关文章推荐

HDU 5950 Recursive sequence 【递推+矩阵快速幂】 (2016ACM/ICPC亚洲区沈阳站)

Recursive sequence Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Ot...

【递推+矩阵快速幂】【HDU2604】【Queuing】

Queuing Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total S...

矩阵快速幂在常系数线性递推关系中的应用

先引入一下,知乎上有一个问题 关于斐波拉契数列的一个低级问题 。题主询问了关于求解斐波拉契数列第n项对10007取模的结果。而这个n,可以达到10610^6甚至10910^9 。 解法已经在排名第一...
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)