generator 1（2019牛客暑期多校训练营（第五场）B，十进制矩阵倍增乘法）-CSDN博客

本文链接：https://blog.csdn.net/The___Flash/article/details/98178828

一.题目链接：

generator 1

二.题目大意：

给出 $x_{0},\; x_{1},\; a,\;b,\;n,\;mod .$

有递推式： $x_{i} = bx_{i - 2} + ax_{i - 1} (mod)$

求第 n 项. （ $1 \leq n \leq 10^{10^{6}}$ ）

三.分析：

这题如果直接用二进制倍增的话会 TLE. （要将 n 转换为 2 进制）

但如果十进制倍增的话，总共需要操作 $10^{6} \times 10$ 次. （貌似还不错的样子）

又学会了一个新知识点（从来没听说过。。。）

~~出题人：这题我小时候做过（orzzzzzzzz）~~

四.代码实现：

#include <set>
#include <map>
#include <ctime>
#include <queue>
#include <cmath>
#include <stack>
#include <bitset>
#include <vector>
#include <cstdio>
#include <sstream>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#define eps 1e-8
#define lc k * 2
#define rc k * 2 + 1
#define pi acos(-1.0)
#define ll long long int
using namespace std;

const int M = (int)1e6;
//const ll mod = (ll)1e9 + 7;
const ll inf = 0x3f3f3f3f3f;

struct node
{
    ll D[2][2];
}a, b1, b2;;

ll mod;
char s[M + 5];

node mul(struct node a, struct node b)
{
    struct node c;
    memset(c.D, 0, sizeof(c.D));
    for(int i = 0; i < 2; ++i)
    {
        for(int j = 0; j < 2; ++j)
        {
            for(int k = 0; k < 2; ++k)
            {
                c.D[i][j] = (c.D[i][j] + a.D[i][k] * b.D[k][j] % mod) % mod;
            }
        }
    }
    return c;
}

int main()
{
    scanf("%lld %lld %lld %lld", &a.D[0][0], &a.D[0][1], &b1.D[1][1], &b1.D[0][1]);
    b1.D[1][0] = 1;
    scanf("%s", s + 1);
    scanf("%lld", &mod);
    int len = strlen(s + 1);
    for(int i = len; i; --i)
    {
        b2 = b1;
        for(int j = 1; j <= 9; ++j)
        {
            if(s[i] - '0' == j)
                a = mul(a, b2);
            b2 = mul(b2, b1);
        }
        b1 = b2;
    }
    printf("%lld\n", a.D[0][0]);
    return 0;
}