【打表计划】【小澳的坐标系】 矩阵快速幂

【打表计划】

题目描述

CYR 想打 n个表，他一次可以打一个、两个、三个或四个表，求他打 n 个表的方案总数。

输入格式

一行一个整数 n，表示 CYR 想打的表的数量。

输出格式

一行一个整数，表示 CYR 打表的方案数，答案对2769433 取模。

测试样例

样例输入
4

样例输出
8

数据范围与提示

对于 30% 的数据，n≤5；
对于 50% 的数据，n≤10 ^ 5​​；
对于 80% 的数据，n≤10 ^ 9；
对于 100% 的数据，n≤10 ^ 16​​。

我本来以为这道题是跟跳台阶一样的水题【naive】，结果是矩阵快速幂，可是我矩乘和快速幂都不会啊【滑稽*n】

//代码主要是抄Menci和Rachel的qnq #include <cstdio>#include <cstring>#include <iostream>using namespace std;const int MOD = 2769433;long long n;struct Matrix {    long long a[4][4];    Matrix ()    {        memset(a, 0, sizeof(a));    }};Matrix operator* (const Matrix &a, const Matrix &b){    Matrix res;    for (int i = 0; i < 4; i++)    {        for (int j = 0; j < 4; j++)        {            for (int k = 0; k < 4; k++)            {                res.a[i][j] = (res.a[i][j] % MOD + (a.a[i][k] % MOD * b.a[k][j] % MOD) % MOD) % MOD;            }        }    }    return res;}Matrix pow(Matrix a, long long n){    Matrix res;    if (n == 0)    {        res.a[0][0] = res.a[1][1] = res.a[2][2] = res.a[3][3] = 1; //单位矩阵 相当于1，也就是每个数乘它都得到自己本身         return res;    }    else if (n % 2 == 1)    {        return pow(a * a, n / 2) * a;    }    else    {        return pow(a * a, n / 2);    }}int main(){    scanf("%lld", &n);    if (n == 0)    {        printf("1\n");    }    else if (n == 1)    {        printf("1\n");    }    else if (n == 2)    {        printf("2\n");    }    else if (n == 3)    {        printf("4\n");    }    else if (n == 4)    {        printf("8\n");    }    else    {        Matrix t, s;        s.a[0][0] = 1, s.a[1][0] = 2, s.a[2][0] = 4, s.a[3][0] = 8;        t.a[0][0] = 0, t.a[0][1] = 1, t.a[0][2] = 0, t.a[0][3] = 0;        t.a[1][0] = 0, t.a[1][1] = 0, t.a[1][2] = 1, t.a[1][3] = 0;        t.a[2][0] = 0, t.a[2][1] = 0, t.a[2][2] = 0, t.a[2][3] = 1;        t.a[3][0] = 1, t.a[3][1] = 1, t.a[3][2] = 1, t.a[3][3] = 1;        s = pow(t, n - 4) * s;        long long ans = s.a[3][0];        printf("%lld\n", ans);    }}

我一开始想的是开两个结构体，因为有两个不同形状的矩阵嘛【滑稽】，但其实开一个结构体就行了，就是大的那个矩阵的形状，小的也用那个，不足的就直接补零了，方便了好多。

【小澳的坐标系】

题目描述

小澳的梦境中出现了一个平面直角坐标系，自原点，向四方无限延伸。

小澳在坐标系的原点，他可以向上、向左或者向右走。他可以走 n 步，但不能经过相同的点。

小澳想知道他有多少种走法。

输入格式

输入文件仅第一行一个正整数 n ，表示小澳可以走的步数。

输出格式

输出文件共一行，输出一个正整数，表示答案（对10 ^ 9 + 7取模）。

测试样例

样例输入 1
2

样例输出 1
7

样例输入 2
3

样例输出 2
17

数据范围与提示

测试点编号
n
1−2 n≤10
3−4 n≤100
5−6 n≤1000
7−8 n≤10 ^ 6​​
9−10 n≤10 ^ 9

//代码还是主要抄Menci和Rachel的qnq 我已经是个废鱼了qnq #include <cstdio>#include <iostream>#include <algorithm>#include <cstring>using namespace std;const int MOD = 1e9 + 7;struct Matrix {    long long a[2][2];    Matrix()    {        memset(a, 0, sizeof(a));    }};Matrix operator*(const Matrix &a, const Matrix &b){    Matrix res;    for (int i = 0; i < 2; i++)    {        for (int j = 0; j < 2; j++)        {            for (int k = 0; k < 2; k++)            {                res.a[i][j] = (res.a[i][j] % MOD + (a.a[i][k] % MOD * b.a[k][j] % MOD) % MOD) % MOD;            }        }    }    return res;}Matrix pow(Matrix a, long long n){    Matrix res;    if (n == 0)    {        res.a[0][0] = res.a[1][1] = 1;        return res;    }    else    {        if (n % 2 == 1)        {            return pow(a * a, n / 2) * a;        }        else        {            return pow(a * a, n / 2);        }    }}int main(){    //freopen("coordinate.in", "r", stdin);    //freopen("coordinate.out", "w", stdout);     long long n;    scanf("%lld", &n);    if (n == 1) printf("3\n");    else if (n == 2) printf("7\n");    else    {        Matrix t, s;        t.a[0][0] = 3; t.a[1][0] = 7;        s.a[0][0] = 0; s.a[0][1] = 1;        s.a[1][0] = 1; s.a[1][1] = 2;        t = pow(s, n - 1) * t; //不懂         long long ans = t.a[0][0] % MOD; //不懂         printf("%lld\n", ans);    }    return 0;}

两道题我都是照着题解写的【辣鸡】，转移矩阵自己没推，不过自己凑一个只要满足要求就行。就像小澳的坐标系，题解的转移矩阵是【1,1】【2,1】，i在上，i - 1在下，而Rachel的是【0,1】【1,2】，是i - 1再上，i在下，都能A。
当n=0时的单位矩阵是我没有想到的，就是找一个矩阵让它乘自己自己不变。
还有转移矩阵一般都是方的

pow还可以写成

orz pow(orz a, int u){    orz res;    res.a[0][0] = res.a[1][1] = 1;    for (; u; u >>= 1, a = a * a) if (u & 1) res = res * a;    return res;}

//代码大都不是我自己写的，侵删