[Sdoi2013]随机数生成器（BSGS）

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<map>
using namespace std;
typedef long long ll;
struct Thash{
    static const int MOD=;
    static const int MAXN=1e6+;
    int tot,head[MOD+],next[MAXN],h[MAXN],val[MAXN];
    inline void clear(){tot=;memset(head,,sizeof head);}
    inline void insert(int H,int VAL){
        for(int i=head[H%MOD];i;i=next[i]) if(h[i]==H){val[i]=VAL;return ;}
        next[++tot]=head[H%MOD];head[H%MOD]=tot;h[tot]=H;val[tot]=VAL;
    }
    inline int get(int H){
        for(int i=head[H%MOD];i;i=next[i]) if(h[i]==H) return val[i];
        return ;
    }
}M;
int T;
ll p,a,b,x1,t;
inline int read(){
    int x=,f=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
    return x*f;
}
ll fpow(ll a,ll p,ll mod){
    ll res=;
    for(;p;p>>=,a=a*a%mod) if(p&) res=res*a%mod;
    return res;
}
void exgcd(ll a,ll b,ll &d,ll &x,ll &y){
    if(!b){d=a;x=;y=;return ;}
    exgcd(b,a%b,d,y,x);
    y-=x*(a/b);
}
ll BSGS(ll A,ll B,ll mod){
    A%=mod;
    if(!A){
        if(!B) return ;
        return -;
    }
    ll m=sqrt(mod)+,ni=fpow(A,mod-m-,mod);
    ll t=,y=;
    M.clear();
    M.insert(,m+);
    for(int i=;i<m;i++){
        t=t*A%mod;
        if(!M.get(t)) M.insert(t,i);
    }
    for(int i=;i<m;i++){
        int u=M.get(B*y%mod);
        if(u){
            if(u==m+) u=;
            return i*m+u;
        }
        y=y*ni%mod;
    }
    return -;
}
ll calc1(){
    ll z=(t-x1+p)%p,d,x,y;
    exgcd(b,p,d,x,y);
    if(z%d) return -;z/=d;
    x=x*z%p;
    if(x<) x+=p;
    return x+;
}
ll calc2(){
    ll c=fpow(a-,p-,p),A=(x1+b*c)%p,z=(b*c+t)%p,d,x,y;
    exgcd(A,p,d,x,y);
    if(z%d) return -;z/=d;
    if(x<p) x=x%p+p;
    d=BSGS(a,x*z%p,p);
    if(~d) return d+;
    return -;
}
ll solve(){
    if(x1==t) return ;
    if(!a){
        if(b==t) return ;
        else return -;
    }
    if(a==) return calc1();
    return calc2();
}
int main(){
    freopen("random.in","r",stdin);
    freopen("random.out","w",stdout);
    for(T=read();T--;){
        p=read();a=read();b=read();x1=read();t=read();
        printf("%lld\n",solve());
    }
    return ;
}