|
%I #19 Jul 08 2023 18:09:19
%S 3,0,21,3,201,13,18081,43,140049,81,600009,147,6000009,73,380000361,3,
%T 1400000049,831,14000000049,49,380000000361,987,600000000009,691,
%U 78000000001521,183,740000000001369,4153,6200000000000961,279
%N Least k such that 10^n+k is a brilliant number (cf. A078972).
%C If n is an even positive exponent, then a(n) is the first prime greater than 10^(n/2) squared less 10^n.
%H Chai Wah Wu, <a href="/A083289/b083289.txt">Table of n, a(n) for n = 0..42</a>
%H Dario Alejandro Alpern, <a href="https://www.alpertron.com.ar/BRILLIANT.HTM">Brilliant numbers</a>
%t NextPrim[n_] := Block[{k = n + 1}, While[ !PrimeQ[k], k++ ]; k]; LengthBase10[n_] := Floor[ Log[10, n] + 1]; f[n_] := Block[{k = 0}, If[ EvenQ[n] && n > 1, NextPrim[ 10^(n/2)]^2 - 10^n, While[fi = FactorInteger[10^n + k]; Plus @@ Flatten[ Table[ # [[2]], {1}] & /@ fi] != 2 || Length[ Union[ LengthBase10 /@ Flatten[ Table[ # [[1]], {1}] & /@ fi]]] != 1, k++ ]; k]]; Table[ f[n], {n, 0, 30}]
%o (Python)
%o from sympy import nextprime, factorint
%o def A083289(n):
%o a, b = divmod(n,2)
%o c, d = 10**n, 10**a
%o if b == 0: return nextprime(d)**2-c
%o k = 0
%o while True:
%o fs = factorint(c+k,multiple=True)
%o if len(fs) == 2 and min(fs) >= d:
%o return k
%o k += 1 # _Chai Wah Wu_, Sep 28 2021
%Y Cf. A078972, A084475, A084476.
%K base,nonn
%O 0,1
%A _Jason Earls_, Jun 03 2003
%E Edited and extended by _Robert G. Wilson v_, Jun 27 2003
|
