The Luckiest Number[POJ3696]

欢迎大家访问我的老师的OJ———caioj.cn

题面描述

传送门

思路

一句话题意：求至少多少个$8$连在一次组成的正整数是$L$的倍数。
$\because x个8$实际上就是$(10^x-1)/9\ast 8$，现在我们可以转化一下题意:
求最小的$x$，满足
$$L\mid (10^x-1)/9\ast 8$$

$$L\mid (10^x-1)/9\ast 8⟷L\ast 9\mid (10^x-1)\ast 8⟷\frac{9L}{\gcd (8,L)}\mid 10^x-1⟷10^x\equiv 1(\mathrm{mod}\frac{9L}{\gcd (8,L)})$$

引理:

若正整数a,n互质，且满足$a^x\equiv 1(\mathrm{mod}n)$的最小正整数$x_0$是$\phi (n)$的约数。

证明:

假设满足$a^x\equiv 1(\mathrm{mod}n)$的最小正整数$x_0$不能整除$\phi (n)$。

设$\phi (n)=q{x}_0+r(0<r<x_0)$.

因为$a^{x_0}\equiv 1(\mathrm{mod}n)$，

则$a^{q{x}_0}\equiv 1(\mathrm{mod}n)$，

又$a^{\phi (n)}\equiv 1(\mathrm{mod}n)$，

故$a^r\equiv 1(\mathrm{mod}n)$，

与假设$x_0$最小相矛盾。故假设不成立。原命题成立。

证毕。

根据以上引理，我们只要求出$\phi (9L/gcd(L,8))$,枚举它的所有约数，用快速幂逐一检查是否满足条件即可。

AC code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#define ll long long
#define gc getchar()
using namespace std;
const int N=15;
ll p[N];int c[N];int t=0;
inline void qr(ll &x)
{
	x=0;int f=1;char c=gc;
	while(c<'0'||c>'9'){if(c=='-')f=-1;c=gc;}
	while(c>='0'&&c<='9'){x=x*10+(c^48);c=gc;}
	x*=f;
}
ll L;
inline void mul(ll &a,ll b,ll p)
{
	ll c=(long double)a*b/p;
	a=a*b-c*p;
}
inline ll pow_mod(ll a,ll b,ll p)
{
	ll ans=1%p;a%=p;
	while(b)
	{
		if(b&1)mul(ans,a,p);
		b>>=1;mul(a,a,p);
	}
	return ans;
}
inline ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
inline void div(ll x)
{
	p[0]=0;
	for(ll i=2;i*i<=x;i++)
		if(x%i==0)
		{
			p[++p[0]]=i;
			c[p[0]]=0;
			while(x%i==0)x/=i,++c[p[0]];
		}
	if(x>1)p[++p[0]]=x,c[p[0]]=1;
}
bool dfs(int x,ll ans)
{
	if(x>p[0])
	{
		if(pow_mod(10ll,ans,L)==1)return printf("Case %d: %lld\n",t,ans),true;
		return false;
	}
	ll y=1;
	for(int i=0;i<=c[x];i++)
	{
		if(dfs(x+1,ans*y))return true;
		y*=p[x];
	}
	return false;
}
inline ll phi(ll x)
{
	ll ans=x;
	for(ll i=2;i*i<=x;i++)
	{
		if(x%i==0)
		{
			ans=ans/i*(i-1);
			while(x%i==0)x/=i;
		}
	}
	if(x>1)ans=ans/x*(x-1);
	return ans;
}
int main()
{
	while(qr(L),L)
	{
		L=L*9/gcd(8,L);++t;
		if(gcd(10ll,L)>1){printf("Case %d: 0\n",t);continue;}
		ll phx=phi(L);
		div(phx);
		dfs(1,1);
	}
	return 0;
}