bzoj2875 [Noi2012]随机数生成器 [矩阵+快乘]

http://www.lydsy.com/JudgeOnline/problem.php?id=2875

Description

给出 a,c,g,mod,x0,n,xn=(axn1+c) ,求 xn%g

Input

包含6个用空格分割的m,a,c,X0,n和g,其中a,c,X0是非负整数,m,n,g是正整数。

Output

输出一个数,即Xn mod g

Sample Input

11 8 7 1 5 3

Sample Output

2

Solution

比较裸。

#include <bits/stdc++.h>

const int MAXN = 5;

using namespace std;

long long m, a, c, x0, n, g;
struct matrix {
    long long x, y;
    matrix() {}
    matrix(long long a, long long b) : x(a), y(b) {}
}origin;

long long mult(long long a, long long b) {
    long long ans = 0;
    while (b) {
        if (b & 1) ans = (ans + a) % m;
        a = (a + a) % m;
        b >>= 1;
    }
    return ans;
}

inline matrix operator * (matrix a, matrix b) {
    return matrix(mult(a.x, b.x), (mult(a.x, b.y) + a.y) % m);
}

inline matrix pow(matrix a, long long x) {
    matrix res = a;
    x--;
    while (x) {
        if (x & 1) {
            res = res * a;
        }
        x >>= 1;
        a = a * a;
    }
    return res;
}

int main() {
    scanf("%lld %lld %lld %lld %lld %lld", &m, &a, &c, &x0, &n, &g);
    origin = matrix(a, c);
    origin = pow(origin, n);
    printf("%lld\n", (mult(origin.x, x0) + origin.y) % m % g);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值