login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A351076
Numbers k such that the maximal exponent in the prime factorization of A327860(k) is greater than or equal to A051903(k), the maximal exponent in the prime factorization of k.
4
1, 3, 5, 7, 8, 10, 11, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 31, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 61, 62, 65, 66, 67, 69, 70, 71, 72, 73, 74, 77, 78, 79, 80, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 97, 98
OFFSET
1,2
COMMENTS
Numbers k for which A328391(k) >= A051903(k).
PROG
(PARI)
A051903(n) = if((1==n), 0, vecmax(factor(n)[, 2]));
A327860(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= (p^e); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); };
isA351076(n) = (A051903(A327860(n)) >= A051903(n));
CROSSREFS
Cf. A003415, A276086, A051903, A327860, A328391, A351075 (complement), A351077 (the characteristic function).
Positions of nonnegative terms in A351074.
Cf. also A350076.
Sequence in context: A360392 A189829 A175146 * A171508 A049068 A185602
KEYWORD
nonn
AUTHOR
Antti Karttunen, Feb 01 2022
STATUS
approved