CF#460 div2 E - 费马小+逆元

解题思路:

根据费马小定理可以将a^n变为a^(n%(p-1)),那么n丶范围就变成了0-(p-2),那么原式就可以化为(i+k(p-1))a^i，k为非负整数，i就是取余后的数，那么问题就变成当k取什么值时，这个式子对p取模会为b。我们不妨首先去找最小的k是多少，设a^i的逆元c，t = b*c%p，那么i+k(p-1) == t + jp,k和j都是非负整数，当i==t时，k=0；当i>t时，k=i-t；当i<t时，k=p-t+i(根据k<=p可以全部推出)。有了最小值k，就有了最小式子的值y，设y1>y,且也满足y1*a^i%p == b,又y1 = y + m(p-1)，m为非负整数，那么就又要确定m的最小值，m最小值其实就是(p-1)*a^i和p的最小公倍数除以(p-1)*a^i,有了m就知道其他取值是2m,3m,4m.....问题就解决了。

#include<bits/stdc++.h>
using namespace std;
const int mx = 1e6+10;
typedef long long ll;
ll x,a,b,p;
ll qpow(ll x,ll y)
{
    ll ans = 1;
    while(y){
        if(y&1) ans = ans*x%p;
        y >>= 1;
        x = x*x%p;
    }
    return ans;
}
int main()
{
    scanf("%I64d%I64d%I64d%I64d",&a,&b,&p,&x);
    ll p1 = b % p,k,ans = 0;
    for(int i=1;i<p;i++){
        ll po = qpow(a,i%(p-1)),ret = 0;
        ll t = p1*qpow(po,p-2)%p;
        if(i>t) k = i - t;
        else if(i==t) k = 0;
        else k = p - t + i;
        ll y = i + k*(p-1);
        if(y<=x){
            t = (p-1)*po%p;
            k = t*p/__gcd(t,p)/t;
            ret = (x-y)/(k*(p-1)) + 1;
        }
        ans += ret;
    }
    printf("%I64d\n",ans);
    return 0,0;
}