[瞎搞]JZOJ 5941 乘

最新推荐文章于 2021-03-05 10:31:55 发布

weixin_33778778

最新推荐文章于 2021-03-05 10:31:55 发布

阅读量73

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/mastervan/p/9901264.html

本文介绍了一种处理大数幂取模运算的高效算法，通过将大指数分解为两部分，利用预计算的幂值数组，实现了快速计算。文章提供了完整的C++代码实现，并在pow.in和pow.out文件中进行了输入输出示例。

Description

Input

Output

Sample Input

Sample Input1：
4 3 9 6
5 8 7 7

Sample Input2：
见下发文件中的 ex_pow2.in/out.

Sample Output

Sample Output1：
0

Sample Output2：

Data Constraint

分析

我们发现b最大就10^12，那么我们把它拆成两个部分，一部分是10^6以下的幂，一部分是%10^6的幂，那这样两个数组就能存下来，计算的时候简单地乘起来就好了

#include <iostream>
#include <cstdio>
using namespace std;
typedef long long ll;
ll a,p,ans,b,m,c,l,pow[1000010],powup[1000010];
int q,k;

int main() {
    freopen("pow.in","r",stdin);
    freopen("pow.out","w",stdout);
    scanf("%lld%lld%d%d",&a,&p,&q,&k);
    scanf("%lld%lld%lld%lld",&b,&l,&m,&c);
    b%=l;m%=l;c%=l;
    pow[0]=powup[0]=1;
    for (int i=1;i<=1000000;i++) pow[i]=pow[i-1]*a%p;
    powup[1]=pow[1000000];
    for (int i=2;i<=l/1000000;i++) powup[i]=powup[i-1]*powup[1]%p;
    for (int i=1;i<=q;i++) {
        b=(b*m%l+c)%l;
        ans^=pow[b%1000000]*powup[b/1000000]%p;
        if (!(i%k)) printf("%lld\n",ans);
    }
}