OFFSET
1,2
COMMENTS
For p prime, numbers of the form p^p satisfy the condition, hence A051674 is a subsequence. - Michel Marcus, May 19 2014
Also, numbers of the form p^q * q^p, with distinct primes p and q, satisfy the condition, hence A082949 is a subsequence. - Bernard Schott, Apr 11 2020
LINKS
Rémy Sigrist, Table of n, a(n) for n = 1..10000
Rémy Sigrist, PARI program for A054412
EXAMPLE
192 is included because 192 =2^6 *3^1 and 2*3 = 6*1.
MATHEMATICA
peppfQ[n_]:=Module[{f=Transpose[FactorInteger[n]]}, Times@@First[f] == Times@@Last[f]]; Select[Range[1.5*10^6], peppfQ] (* Harvey P. Dale, Oct 14 2015 *)
PROG
(PARI) isok(n) = my(f = factor(n)); prod(i=1, #f~, f[i, 2]) == prod(i=1, #f~, f[i, 1]); \\ Michel Marcus, May 19 2014
(PARI) See Links section.
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, May 09 2000
EXTENSIONS
More terms from James A. Sellers, May 23 2000
New name and three more terms from Michel Marcus, May 19 2014
STATUS
approved