Sample Input1： 4 3 9 6 5 8 7 7
Sample Output
Sample Output1： 0
做法（转自JZOJ）：考虑 a 是定值, 而 b ≤ 1012 , 我们可以预处理 a 的 0...106 次方与 1 ? 106 , 2 ? 106 ...106 ? 106 次 方, 询问时把 b 切成两半, 拼起来就是答案. 这样询问就是 O(1) 的. 复杂度 O(q + 106 ).
标签：预处理 freopen 次方 def hide define out stdin input
Sample Input1： 4 3 9 6 5 8 7 7
Sample Output1： 0
做法（转自JZOJ）：考虑 a 是定值, 而 b ≤ 1012 , 我们可以预处理 a 的 0...106 次方与 1 ? 106 , 2 ? 106 ...106 ? 106 次 方, 询问时把 b 切成两半, 拼起来就是答案. 这样询问就是 O(1) 的. 复杂度 O(q + 106 ).
1 #include <cstdio> 2 #include <cstring> 3 #include <iostream> 4 #include <algorithm> 5 #include <string> 6 #define LL long long 7 #define N 1000007 8 using namespace std; 9 LL a,p,q,k; 10 LL b,m,l,C,b1; 11 LL ans,mia[N],la[N]; 12 13 LL ksm(LL a,LL b){ 14 LL root=1,base=a; 15 while(b){ 16 if (b&1) root*=base,root%=p; 17 base*=base; 18 base%=p; 19 b>>=1; 20 } 21 return root; 22 } 23 24 25 int main(){ 26 freopen("pow.in","r",stdin); 27 freopen("pow.out","w",stdout); 28 scanf("%lld%lld%lld%lld",&a,&p,&q,&k); 29 scanf("%lld%lld%lld%lld",&b,&l,&m,&C); 30 mia[0]=1; 31 la[0]=1; 32 LL g=1000000; 33 for(int i=1;i<=g;++i) 34 mia[i]=(mia[i-1]*a)%p,la[i]=ksm(mia[i],g)%p; 35 for(int i=1;i<=q;++i){ 36 b1=(b*m+C)%l; 37 ans=ans^((mia[b1%g]*la[b1/g])%p); 38 if (i%k==0) printf("%lld\n",ans); 39 b=b1; 40 } 41 }
标签：预处理 freopen 次方 def hide define out stdin input
原文地址：https://www.cnblogs.com/traveller-ly/p/9892586.html
