Codeforces 919E 小学数学入门

题意：给你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;
}