OFFSET
1,1
COMMENTS
REFERENCES
Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004, Section B46, p. 154.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
41 is prime, 40/5 = 8 = 2^3, hence 41 is not a term.
101 is prime, 100/5 = 20 = 2^2 * 5, hence 101 is a term.
491 is prime, 490/7 = 70 = 2 * 5 * 7, hence 491 is not a term.
521 is prime, 520/13 = 40 = 2^3 * 5, hence 521 is a term.
MAPLE
alias(pf = NumberTheory:-PrimeFactors): gpf := n -> max(pf(n)):
is_a := n -> isprime(n) and pf((n-1)/gpf(n-1)) = {2, 5}:
select(is_a, [$5..3371]); # Peter Luschny, Dec 13 2020
MATHEMATICA
q[n_] := Divisible[n, 10] && ((PrimeQ[(r = n/2^IntegerExponent[n, 2]/5^(e = IntegerExponent[n, 5]))] && r > 5) || (r == 1 && e > 1)); Select[Range[3500], PrimeQ[#] && q[# - 1] &] (* Amiram Eldar, Dec 13 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Bernard Schott, Dec 13 2020
STATUS
approved