Codeforces919E

前言：

惭愧，昨天竟然没有做出来……

题解：

其实挺简单。首先像这种题目肯定是有个周期的，对于这题，一个周期的长度就是$p\times(p-1)$。因为$n\%p$的周期是$p$，$a^n\%p$的周期是$p-1$（根据费马小定理或欧拉定理得出），而且，对于一个周期内的所有$n\%p=x$(意思是$n=x+k\times p$)，只会有一个解，因为它对应的$a^n\%p$都是不一样的，所以我们只要求出一个第一个周期内的解就可以了。对于一个枚举的$n$，我们先把$a^n$除掉，也就是$b$乘上它的逆元，然后不断地调整这个$n$，因为$n$每加上一个$p-1$，$a^n\%p$是不变的，而$n\%p$却减了$1$，发现这点后就很简单了。

可能是昨天没想出来的原因：

1、昨天想到把两边除掉$n$，没想到除掉$a^n$。
2、想着hack人。

代码：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
LL read()
{
    LL x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x*f;
}
LL a,b,p,x;
LL Pow(LL x,LL y,LL mod)
{
    if(!y)return 1;if(y==1)return x%mod;
    LL t=Pow(x,y>>1,mod),ans=t*t%mod;
    if(y&1)ans=ans*x%mod;
    return ans;
}
int main()
{
    a=read(),b=read(),p=read(),x=read();
    LL P=p*(p-1),ans=0;
    for(LL i=1;i<p;i++)
    {
        LL r=b*Pow(Pow(a,i,p),p-2,p)%p,c;
        if(i>=r)c=(i+(i-r)*(p-1))%P;
        else c=(i+(i+p-r)*(p-1))%P;
        if(!c)c=P;
        ans+=((x/P)+((x%P>=c)?1:0));
    }
    printf("%lld",ans);
}