A347818 Smallest n-digit brilliant number.
4, 10, 121, 1003, 10201, 100013, 1018081, 10000043, 100140049, 1000000081, 10000600009, 100000000147, 1000006000009, 10000000000073, 100000380000361, 1000000000000003, 10000001400000049, 100000000000000831, 1000000014000000049, 10000000000000000049, 100000000380000000361
Offset: 1
Examples
a(6) = 100013 = 103 * 971. a(7) = 1018081 = 1009 * 1009. a(8) = 10000043 = 2089 * 4787. a(9) = 100140049 = 10007 * 10007.
Links
- Dario Alejandro Alpern, Brilliant numbers
- World of Numbers, Smallest n-digit prp
Programs
-
Mathematica
Join[{4,10},Table[Module[{k=1},While[PrimeOmega[10^n+k]!=2||Length[ Union[ IntegerLength/@ FactorInteger[ 10^n+k][[;;,1]]]]!=1,k+=2];10^n+k],{n,2,20}]] (* Harvey P. Dale, Jan 09 2024 *)
-
PARI
isA078972(n)=my(f=factor(n)); (#f[, 1]==1 && f[1, 2]==2) || (#f[, 1]==2 && f[1, 2]==1 && f[2, 2]==1 && #Str(f[1, 1])==#Str(f[2, 1])) A084476(n)=for(k=0,10^n,if(isA078972(10^(2*n-1)+k),return(k))) a(n)=if(n%2,nextprime(10^((n-1)/2))^2,10^(n-1)+A084476(n/2)) \\ after Charles R Greathouse IV in A078972
Comments