题意:给你a、b、p(素数,<=1e6+3)、x(<=1e12),
求满足n (1 ≤ n ≤ x) 且
的n个数。
题解:等式左边分开考虑,n%p显然最小循环周期是p,而a^n%p最小循环周期是p-1,因为n=0时,a^n%p=1,由费马小定理,a^(p-1)%p=1,得证。
还有一个事情就是n*a^n%p的最小循环周期是p*(p-1),这个事情大概可以这样考虑:显然LCM(p,(p-1))=p*(p-1)是一个循环周期,然后也找不到更小的了,得证…………感觉有点草率OTZ。
对于每一个n,因为a^n的最小循环周期是p-1,因此我们可以在1..p-1范围枚举n,对于每一个a^n独立求解,那么我们先计算出a^n,用逆元除到右边去,使得这个等式变成 n和b/a^n同余,因为我们固定了a^n,所以我们的n每次只能加上p-1,这样给n带来的影响就是变成了(n+p-1)%p=(n-1)%p,那么我们可以很轻松的得到最小的满足条件的正整数n,再结合x的范围,以及周期=p*(p-1)可以得到这一个a^n对答案的总贡献……就做完了。。。
脑子是个好东西……可惜我没有……
Code:
include<bits/stdc++.h>
using namespace std;
typedef long long LL;
LL a,b,p,x;
LL quick(LL x,LL y,LL p){
LL res=1;
while (y){
if (y&1)res = res*x%p;
x = x*x%p;
y/=2;
}
return res;
}
int main(){
cin>>a>>b>>p>>x;
LL ans =0;
LL temp = a;
for (LL i=1;i<p;i++,temp=temp*a%p){
LL c = b*quick(temp,p-2,p)%p;
LL k = (i-c+p)%p;
LL minn = i+k*(p-1);
if (minn>x)continue;
ans += max(0LL,(x-minn)/(p*(p-1))+1);
}
cout<<ans<<endl;
return 0;
}