2013年省赛H题

2013年省赛H题
你不能每次都快速幂算A^x,优化就是预处理，把10^9预处理成10^5和10^4。
想法真的是非常巧妙啊
N=100000
构造两个数组，f1[N],间隔为A
f2[1e4]间隔为A^N,中间用f1来填补
f[x]=f1[x%N]*f2[x/N]%P;

 #include<iostream>
 #include<cstdio>
 #include<queue>
 #include<algorithm>
 #include<cmath>
 #include<ctime>
 #include<set>
 #include<map>
 #include<stack>
 #include<cstring>
 #define inf 2147483647
 #define ls rt<<1
 #define rs rt<<1|1
 #define lson ls,nl,mid,l,r
 #define rson rs,mid+1,nr,l,r
 #define N 100000
 #define For(i,a,b) for(long long i=a;i<=b;i++)
 #define p(a) putchar(a)
 #define g() getchar()

 using namespace std;
 long long T,ans,cnt;
 long long n,A,K,a,b,m,P;
 long long f1[],f2[],f[];

 void in(long long &x){
     long long y=;
     char c=g();x=;
     while(c<''||c>''){
         if(c=='-')y=-;
         c=g();
     }
     while(c<=''&&c>=''){
         x=(x<<)+(x<<)+c-'';c=g();
     }
     x*=y;
 }
 void o(long long x){
     if(x<){
         p('-');
         x=-x;
     }
     if(x>)o(x/);
     p(x%+'');
 }

 void predeal(){
     f1[]=f2[]=;
     For(i,,N)
         f1[i]=f1[i-]*A%P;
     For(i,,1e4)
         f2[i]=f2[i-]*f1[N]%P;
 }

 int main(){
     in(T);
     while(T--){
         in(n);in(A);in(K);in(a);in(b);in(m);in(P);
         predeal();
         f[]=K;
         For(i,,n)
             f[i]=(a*f[i-]+b)%m;
         ans=;
         For(i,,n){
             ans+=f1[f[i]%N]*f2[f[i]/N]%P;
             ans%=P;
         }
         cout<<"Case #"<<++cnt<<": "<<ans<<endl;
     }
     return ;
 }