【POJ】3243 Clever Y

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <iostream>
typedef long long ll;
using namespace std;
int gcd(int a, int b) { return b?gcd(b, a%b):a; }
void exgcd(ll a, ll b, ll &d, ll &x, ll &y) { if(!b) { d=a; x=1; y=0; return; } exgcd(b, a%b, d, y, x); y-=a/b*x; }
int ni(int a, int b) {
	static ll x, y, d;
	exgcd(a, b, d, x, y);
	return (x+b)%b;
}
int ipow(int a, int b, int c) { int x=1; for(; b; b>>=1, a=(ll)a*a%c) if(b&1) x=(ll)x*a%c; return x; }
struct H {
	static const int md=3999997;
	bool vis[md]; int dt[md], foo[md], s[md], top;
	void clr() { while(top) { int x=s[top--]; vis[x]=0, dt[x]=-1; } }
	void add(int a, int b) {
		int x=a%md;
		while(1) { if(!vis[x] || dt[x]==a) { if(!vis[x]) s[++top]=x; vis[x]=1; dt[x]=a; foo[x]=b; break; } ++x; if(x==md) x=0; }
	}
	int find(int a) {
		int x=a%md;
		while(1) { if(!vis[x]) return -1; if(dt[x]==a) return foo[x]; ++x; if(x==md) x=0; }
	}
}h;
int a, b, mo;
bool spj() {
	if(b==1) puts("0");
	else if(!a && !b) puts("1");
	else if(!b) puts("No Solution");
	else return 0;
	return 1;
}
void work() {
	b%=mo;
	if(spj()) return;
	for(int i=0, t=1; i<30; ++i, t=(ll)t*a%mo) if(t==b) { printf("%d\n", i); return; }
	int a1=0, d=1, g;
	while((g=gcd(a, mo))!=1) { if(b%g) { puts("No Solution"); return; } ++a1, mo/=g, b/=g, d=(ll)a/g*d%mo; }
	d=ni(d, mo);
	int m=sqrt(0.5+mo);
	h.clr();
	for(int i=0, t=(ll)b*d%mo; i<m; ++i, t=(ll)t*a%mo) h.add(t, i);
	for(int i=0, wn=ipow(a, m, mo), t=1; i<=m; ++i, t=(ll)t*wn%mo) { int x=h.find(t); if(x!=-1) { printf("%d\n", a1+i*m-x); return; } }
	puts("No Solution");
}
int main() {
	while(scanf("%d%d%d", &a, &mo, &b) && (a|b|mo)) work();
	return 0;
}