Time Limits: 2000 ms Memory Limits: 262144 KB
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
Solution
多次给出同一个数的指数,求这个数的幂
因为是同一个数嘛,肯定不可以每次都用快速幂,这太浪费时间了
于是我们可以先预处理
a
a
a 的
0
0
0 ~
1
0
6
10^6
106 次方和
1
∗
1
0
6
1*10^6
1∗106 ~
1
0
6
∗
1
0
6
10^6*10^6
106∗106
询问时将
b
b
b 分成
b
b
b
d
i
v
div
div
n
n
n 和
b
b
b
m
o
d
mod
mod
n
n
n 两部分求解即可
Code
#include<algorithm>
#include<cstdio>
#define fo(i,a,b) for(int i=a;i<=b;++i)
#define fd(i,a,b) for(int i=a;i>=b;--i)
#define ll long long
using namespace std;
const int N=1e6;
ll a,p,q,k,b,l,m,c,d[N+5],e[N+5];
int main()
{
freopen("pow.in","r",stdin);
freopen("pow.out","w",stdout);
scanf("%lld%lld%lld%lld",&a,&p,&q,&k);
scanf("%lld%lld%lld%lld",&b,&l,&m,&c);
d[0]=1; fo(i,1,N) d[i]=d[i-1]*a%p;
e[0]=1; fo(i,1,N) e[i]=e[i-1]*d[N]%p;
ll ans=0;
fo(i,1,q/k)
{
ll sum=0;
fo(j,(i-1)*k+1,i*k)
{
b=(m*b%l+c)%l;
sum^=d[b%N]*e[b/N]%p;
}
ans^=sum;
printf("%lld\n",ans);
}
}