JZOJ-senior-5941. 【NOIP2018模拟11.01】乘

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/HuangXinyue1017/article/details/83821665

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$ 的 $0$ ~ $10^6$ 次方和 $1*10^6$ ~ $10^6*10^6$
询问时将 $b$ 分成 $b$ $div$ $n$ 和 $b$ $mod$ $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);
	}
}