OFFSET
1,1
COMMENTS
The asymptotic density of this sequence is A336065 = 0.848957... (Schinzel and Šalát, 1994).
REFERENCES
József Sándor and Borislav Crstici, Handbook of Number theory II, Kluwer Academic Publishers, 2004, chapter 3, p. 331.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Andrzej Schinzel and Tibor Šalát, Remarks on maximum and minimum exponents in factoring, Mathematica Slovaca, Vol. 44, No. 5 (1994), pp. 505-514.
EXAMPLE
4 = 2^2 is a term since A051903(4) = 2 is a divisor of 4.
MATHEMATICA
H[1] = 0; H[n_] := Max[FactorInteger[n][[;; , 2]]]; Select[Range[2, 100], Divisible[#, H[#]] &]
PROG
(PARI) isok(m) = if (m>1, (m % vecmax(factor(m)[, 2])) == 0); \\ Michel Marcus, Jul 08 2020
CROSSREFS
A005117 (except for 1) is subsequence.
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jul 07 2020
STATUS
approved