POJ 3696 The Luckiest number(08合肥数論)
2838 ワード
転載は出典を明記してください.ありがとうございます.http://blog.csdn.net/acm_cxlove/articale/detail/7854526 by---cxlove
題目:一つの数nを与えて、彼の倍数を探して、しかもこの数は一人当たり8で、最小のそれの桁数を見つけて、さもなくば0を出力します.
http://poj.org/problem?id=3696
実は面白い数論です.コードの難しさもあまりないです.
残念ながら、できません.
8888888は必ず8*(10^k-1)/9=m*nと書くことができます.nは元の数字で、kは桁数です.
8*(10^k-1)=9***mを少し整理して、令r=gcd(8,n)
8/r*(10^k-1)=9*n/r*mです.8/r*(10^k-1)を9*n/rの倍数にし、8/rと9*n/rを互質にするため、10^k-1==0(MOD 9*n/r)
つまり10^k==1(MOD 9*n/r)です.Ola定理によって知っていますが、もし10と9*n/rが互質であれば、10^phi(9*n/r)=1(mod(9*n/r)です.
したがって、解がない場合は非質です.
次にphi(9*n/r)を質因子に分解し、すべての約数を得て、順次小さいときから式を満足するかどうかを判断する.
注意9*n/rの最大範囲は1.8*10^10です.オーバーフローに注意して、やはり高精度のアナログ乗算をしました.
題目:一つの数nを与えて、彼の倍数を探して、しかもこの数は一人当たり8で、最小のそれの桁数を見つけて、さもなくば0を出力します.
http://poj.org/problem?id=3696
実は面白い数論です.コードの難しさもあまりないです.
残念ながら、できません.
8888888は必ず8*(10^k-1)/9=m*nと書くことができます.nは元の数字で、kは桁数です.
8*(10^k-1)=9***mを少し整理して、令r=gcd(8,n)
8/r*(10^k-1)=9*n/r*mです.8/r*(10^k-1)を9*n/rの倍数にし、8/rと9*n/rを互質にするため、10^k-1==0(MOD 9*n/r)
つまり10^k==1(MOD 9*n/r)です.Ola定理によって知っていますが、もし10と9*n/rが互質であれば、10^phi(9*n/r)=1(mod(9*n/r)です.
したがって、解がない場合は非質です.
次にphi(9*n/r)を質因子に分解し、すべての約数を得て、順次小さいときから式を満足するかどうかを判断する.
注意9*n/rの最大範囲は1.8*10^10です.オーバーフローに注意して、やはり高精度のアナログ乗算をしました.
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#define N 80005
#define maxn 150005
#define LL long long
#define pb(a) push_back(a)
using namespace std;
bool flag[N]={0};
int cnt=0,prime[N];
LL gcd(LL a,LL b){
return b==0?a:gcd(b,a%b);
}
void Prime(){
for(int i=2;i<N;i++){
if(flag[i]) continue;
prime[cnt++]=i;
for(int j=2;j*i<N;j++)
flag[i*j]=true;
}
}
LL get_eular(LL n){
LL ret=1;
for(int i=0;i<cnt&&prime[i]*prime[i]<=n;i++){
if(n%prime[i]==0){
ret*=prime[i]-1;
n/=prime[i];
while(n%prime[i]==0){
n/=prime[i];
ret*=prime[i];
}
}
}
if(n>1) ret*=n-1;
return ret;
}
LL fac[N][2],tot;
void Split(LL n){
tot=0;
for(int i=0;i<cnt&&prime[i]*prime[i]<=n;i++){
if(n%prime[i]==0){
fac[tot][0]=prime[i];
fac[tot][1]=0;
while(n%prime[i]==0){
n/=prime[i];
fac[tot][1]++;
}
tot++;
}
}
if(n>1){fac[tot][0]=n;fac[tot++][1]=1;}
}
vector<LL>fact;
void dfs(int idx,LL num){
if(idx>=tot){
fact.pb(num);
return;
}
LL tmp=1;
for(int i=0;i<=fac[idx][1];i++,tmp*=fac[idx][0])
dfs(idx+1,num*tmp);
}
LL MultMod(LL a,LL b,LL MOD){
a%=MOD;
b%=MOD;
LL ret=0;
while(b){
if(b&1){
ret+=a;
if(ret>=MOD) ret-=MOD;
}
a=a<<1;
if(a>=MOD) a-=MOD;
b=b>>1;
}
return ret;
}
LL PowMod(LL a,LL b,LL MOD){
LL ret=1;
while(b){
if(b&1) ret=MultMod(ret,a,MOD);
a=MultMod(a,a,MOD);
b>>=1;
}
return ret;
}
int main(){
LL n;
int cas=0;
Prime();
while(scanf("%I64d",&n)!=EOF&&n){
LL p=(LL)9*n/gcd(8,n);
printf("Case %d: ",++cas);
if(gcd(10,p)!=1){
printf("0
");
continue;
}
LL phi=get_eular(p);
Split(phi);
fact.clear();
dfs(0,1);
sort(fact.begin(),fact.end());
for(int i=0;i<fact.size();i++){
if(PowMod(10,fact[i],p)==1){
printf("%I64d
",fact[i]);
break;
}
}
}
return 0;
}