poj3696 The luckiest number

$问题转化:$
$$\begin{array}{c}x\\ \overbrace{888\dots 8}\end{array}=\frac{8\ast (10^x-1)}{9}$$
$$令\gcd (8,L)=d,$$
$$L|\frac{8(10^x-1)}{9}⟺\frac{9L}{d}|10^x-1⟺10^x\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}\frac{9L}{d})$$

$可以通过反证法证明:$
$$若a\perp n,则最小的x使a^x\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}n)的x满足x|\phi (n)$$
$所以求出了\phi (n)之后，就可以枚举它的因数快速幂判断满不满足$

$tip：此处的快速幂中的乘法注意防止溢出。$

#include<iostream>
#include<cstdio>

using namespace std;
long long l;

long long mul(long long x,long long y)
{
	long long ret=0;
	while(y)
	{
		if(y&1)(ret+=x)%=l;
		y>>=1;
		x=2*x%l;
	}
	return ret%l;
}

long long qpow(long long x,long long y)
{
	long long ret=1;
	while(y)
	{
		if(y&1)ret=mul(ret,x);
		y>>=1;
		x=mul(x,x);
	}
	return ret%l;
}
int main()
{
	int tot=0;
	while(scanf("%lld",&l)&&l)
	{
		for(int i=1;i<=3;++i)if(!(l&1))l>>=1;
		l*=9;
		
		long long tl=l;
		long long phi=l;
		for(long long i=2;i*i<=tl;++i)
		{
			if(tl%i==0)
			{
				while(tl%i==0)tl/=i;
				phi/=i;
				phi*=i-1;
			}
		}
		if(tl!=1)
		{
			phi/=tl;
			phi*=(tl-1);
		}
		long long ans=1e16;
		for(long long i=1;i*i<=phi;++i)
		{
			if(phi%i==0)
			{
				if(qpow(10,i)==1)
				{
					ans=min(ans,i);
				}
				if(qpow(10,phi/i)==1)
				{
					ans=min(ans,phi/i);
				}
			}
		}
		printf("Case %d: ",++tot);
		if(ans==1e16)printf("0\n");
		else printf("%lld\n",ans);
	}	
	return 0;
}