题意:给出一个数L,现在要求一个数A,使得L|A,并且A的每一位都是8,问A的最少位数是多少。
思路:
首先要把A转化一下,把A表示成,则(p是一个整数),令,那么就存在一个整数k使得,因此,就是求。同余方程有解的条件是(10,m)=1,简单证明:若(10,m)的素因子为2 or 5,但是和10互质,因此没有2、5的素因子,即方程无解。
通过欧拉公式可以知道,但此时并不一定是最小解。设,d是使同余方程成立的最小解,那么就有,k是正整数。也就是说d为的因子,因此枚举的因子,然后判断同余方程是否成立即可。
代码:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<stack>
#include<set>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
ll gcd(ll a,ll b) { return b==0?a:gcd(b,a%b);}
int factor[500],num[500],cnt;
ll res;
ll euler_phi(ll n)
{
int m=sqrt(n+0.5);
ll ans=n;
for(int i=2;i<=m;++i)
{
if(n%i==0)
{
ans=ans-ans/i;
while(n%i==0) n/=i;
}
}
if(n>1) ans=ans-ans/n;
return ans;
}
ll mul_mod(ll a,ll b,ll m)
{
ll res=0;
a%=m;b%=m;
while(b)
{
if(b&1) res=(res+a)%m;
a<<=1;
if(a>=m) a%=m;
b>>=1;
}
return res;
}
ll pow_mod(ll a,ll n,ll m)
{
ll res=1;
while(n)
{
if(n&1) res=mul_mod(res,a,m);
a=mul_mod(a,a,m);
n>>=1;
}
return res;
}
void dfs(int pos,ll now,ll m)
{
if(pos==cnt)
{
if(pow_mod(10,now,m)==1)
res=min(res,now);
return ;
}
ll tmp=1;
for(int i=0;i<=num[pos];++i)
{
dfs(pos+1,tmp*now,m);
tmp*=factor[pos];
}
}
void solve(ll n,ll m)
{
int M=sqrt(n+0.5);
cnt=0;
for(int i=2;i<=M;++i)
{
if(n%i==0)
{
factor[cnt]=i;
num[cnt]=0;
while(n%i==0) {n/=i;num[cnt]++;}
cnt++;
}
}
if(n>1) {factor[cnt]=n;num[cnt++]=1;}
res=m;
dfs(0,1,m);
}
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
int tcase=0;
ll L,m,d;
while(~scanf("%I64d",&L))
{
if(L==0) break;
tcase++;
printf("Case %d: ",tcase);
m=L/gcd(8,L)*9;
if(gcd(10,m)!=1)
{
printf("0\n");
continue;
}
d=euler_phi(m);
solve(d,m);
printf("%I64d\n",res);
}
return 0;
}