OFFSET
1,2
COMMENTS
If k is a term then every number with the same prime signature as k is a term of A374904.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
lps = Cases[Import["https://oeis.org/A025487/b025487.txt", "Table"], {_, _}][[;; , 2]]; f[p_, e_] := e/2 + 1; d[1] = 1; d[n_] := Denominator[Plus @@ f @@@ FactorInteger[n]]; Select[lps, d[#] == 1 &]
PROG
(PARI) is(n) = {my(f = factor(n), p = f[, 1], e = f[, 2]); n == 1 || (prime(#p) == p[#p] && e == vecsort(e, , 4) && denominator(vecprod(apply(x -> x/2 +1, e))) == 1); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jul 23 2024
STATUS
approved
