题目链接:here
思路:一开始看到n是1e5,于是就直接用快速幂了,然后就TLE了,手动
既然nlogn不行,那么只能o(n)了; =-=
然后就只能打表了;
由于m是1e9,所以数组空间存不下,只能改变进制了;
将其改为33000进制,即f=k*a+b,k就是33000.
所以可以用2个k大小的数组储存A^i的值(i>0&&i<b),A^(i*a)的值(i>0&&i<k)
然后就直接分开打表了。
代码:
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <string>
#include <queue>
#include <map>
#include <vector>
using namespace std;
const int N=100005;
const int maxn=33000;
//const int mod=1e9+7;
int t,n,A,K,a,b,m,p;
long long num[40000];
long long mod[40000];
int dabiao()
{
mod[1]=A;
mod[0]=num[0]=1;
for(int i=2;i<=maxn;i++)
{
mod[i]=(mod[i-1]*A)%p;
}
long long c=mod[maxn];
for(int i=1;i<=maxn;i++)
{
num[i]=(c*num[i-1])%p;
}
return 0;
}
int main()
{
int Case=0;
scanf("%d",&t);
while(t--)
{
Case++;
scanf("%d %d %d %d %d %d %d",&n,&A,&K,&a,&b,&m,&p);
dabiao();
long long ans=0,fun=K;
for(int i=1;i<=n;i++)
{
ans+=num[fun/maxn]*mod[fun%maxn]%p;
ans%=p;
fun=(fun*a+b)%m;
}
printf("Case #%d: %lld\n", Case, ans);
}
return 0;
}