OFFSET
1,1
COMMENTS
Presumably, all odd terms are in A000466.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(6) = 12 is a term because 12 = 1*12 with 12 - 1 = 11 prime.
MAPLE
filter:= k -> ormap(d -> d^2 > k and isprime(d - k/d), numtheory:-divisors(k)):
select(filter, [$1..200]); # Robert Israel, Jun 30 2025
MATHEMATICA
A384488Q[k_] := AnyTrue[Divisors[k], PrimeQ[# - k/#] &];
Select[Range[200], A384488Q] (* Paolo Xausa, Jun 30 2025 *)
PROG
(Magma) [k: k in [1..180] | not #[d: d in Divisors (k) | IsPrime(d-(k div d))] eq 0];
(PARI) isok(k) = fordiv(k, d, if (isprime(d - k/d), return(1))); \\ Michel Marcus, Jun 01 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, May 30 2025
STATUS
approved
