题解 | generator 1-2019牛客暑期多校训练营第五场B题

最新推荐文章于 2019-09-09 16:32:07 发布

牛客竞赛

最新推荐文章于 2019-09-09 16:32:07 发布

阅读量160

点赞数

分类专栏：题解 19牛客暑期多校训练营第五场文章标签：题解 19牛客暑期多校训练营

本文链接：https://blog.csdn.net/lastbye/article/details/100521686

版权

题解同时被 2 个专栏收录

97 篇文章 0 订阅

订阅专栏

19牛客暑期多校训练营第五场

10 篇文章 0 订阅

订阅专栏

该博客主要解析了2019年牛客暑期多校训练营第五场的一道题目，提供了详细的输入输出描述及示例，并给出了题解和代码实现。推荐访问牛客竞赛区获取更多题目和详细题解。

摘要由CSDN通过智能技术生成

题目来源于牛客竞赛：https://ac.nowcoder.com/acm/contest/discuss
题目描述：
在这里插入图片描述
输入描述：

输出描述：

示例1：

示例2：

题解：

代码：

#include<cstdio>
#include<cstring>
#include<tuple>
#include <utility>
typedef unsigned long long ULL;
const int SIZE = 3000010;
ULL MOD;
char s[SIZE];
void mul(ULL* c1, ULL* c2, ULL *res){
    res[0] = (c1[0] * c2[0] + c1[1] * c2[2]) % MOD;
    res[1] = (c1[0] * c2[1] + c1[1] * c2[3]) % MOD;
    res[2] = (c1[2] * c2[0] + c1[3] * c2[2]) % MOD;
    res[3] = (c1[3] * c2[3] + c1[2] * c2[1]) % MOD;
}
int main() {
    int a,b;
    int x1,x2;
    scanf("%d%d%d%d", &x1, &x2, &a, &b);
    scanf("%s%llu",s, &MOD);
    int len = 0;
    for(; s[len]; len++);
    s[len-1]--;
    for(int i = len - 1; i >= 0 && s[i] < '0'; i--){
        s[i] = '9';
        s[i-1]--;
    }
    ULL now0 = x1, now1 = x2;
    ULL d[4][4];
    d[0][0] = 0;
    d[0][1] = 1;
    d[0][2] = b;
    d[0][3] = a;
    for(int it = len - 1; it >= 0; it--){
        memset(d[1], 0, sizeof(ULL) * 12);
        for(int p = 1; p < 4; p++){
            mul(d[p-1], d[p-1], d[p]);
        }
        s[it] -= '0';
        for(int p = 0; p < 4; p++){
            if((s[it] >> p) & 1){
                ULL* ml = d[p];
                std::tie(now0, now1) = std::make_pair((ml[0] * now0 + ml[1] * now1) % MOD,(ml[2] * now0 + ml[3] * now1) % MOD);
            }
        }
        mul(d[1], d[3], d[0]);
    }
    printf("%llu\n", now1);
    return 0;
}

更多问题，更详细题解可关注牛客竞赛区，一个刷题、比赛、分享的社区。
传送门：https://ac.nowcoder.com/acm/contest/discuss