OFFSET
1,1
COMMENTS
Larger values not certified.
Some of the larger entries may only correspond to probable primes.
a(25) > 16000. - Michael S. Branicky, Apr 27 2025
EXAMPLE
5 is a term because (10-5)*(10-4)*(10-3)*(10-2)*(10-1) + 1 = 15121 is prime.
10 is a term because (100-10)*(10-9)*..*(10-1) + 1 = 32659201 is prime.
MAPLE
b:= proc(n) option remember; `if`(n=0, 1,
(10^length(n)-n)*b(n-1))
end:
q:= n-> isprime(1+b(n)):
select(q, [$1..300])[]; # Alois P. Heinz, Aug 13 2025
MATHEMATICA
f[n_] := 10^Length[IntegerDigits[n]] - n; p = 1; Do[p *= f[n]; If[PrimeQ[p + 1], Print[n]], {n, 4635}] (* Ryan Propper, May 20 2006 *)
CROSSREFS
KEYWORD
nonn,base,hard,more
AUTHOR
Jason Earls, Aug 01 2005
EXTENSIONS
a(14)-a(20) from Ryan Propper, May 20 2006
a(21)-a(24) from Michael S. Branicky, Apr 25 2025
STATUS
approved
