[瞎搞]JZOJ 5941 乘

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);
    }
}