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[100010],f2[100010],f[1000010];

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

void predeal(){
    f1[0]=f2[0]=1;
    For(i,1,N)
        f1[i]=f1[i-1]*A%P;
    For(i,1,1e4)
        f2[i]=f2[i-1]*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[1]=K;
        For(i,2,n)
            f[i]=(a*f[i-1]+b)%m;
        ans=0;
        For(i,1,n){
            ans+=f1[f[i]%N]*f2[f[i]/N]%P;
            ans%=P;
        }
        cout<<"Case #"<<++cnt<<": "<<ans<<endl;
    }
    return 0;
}