login
A331291
a(n) = A001222(n) * A329348(n).
3
1, 1, 1, 4, 1, 4, 1, 3, 2, 4, 1, 12, 1, 4, 12, 8, 1, 6, 1, 12, 12, 4, 1, 12, 4, 4, 3, 12, 1, 15, 1, 5, 12, 4, 16, 16, 1, 4, 12, 4, 1, 3, 1, 12, 3, 4, 1, 5, 2, 12, 12, 12, 1, 8, 8, 32, 12, 4, 1, 12, 1, 4, 9, 12, 26, 36, 1, 12, 12, 15, 1, 15, 1, 4, 15, 12, 4, 36, 1, 10, 4, 4, 1, 8, 22, 4, 12, 32, 1, 8, 12, 12, 12, 4
OFFSET
1,4
COMMENTS
Records occur at n = 1, 4, 12, 30, 35, 56, 66, 132, 145, 208, 276, 279, 469, ...
FORMULA
a(n) = A001222(n) * A329348(n) = A001222(n) * (A331188(n) mod A117366(n)).
a(p) = 1 for all primes p.
PROG
(PARI)
A002110(n) = prod(i=1, n, prime(i));
A331291(n) = if(1==n, 1, my(f=factor(n), p=nextprime(1+vecmax(f[, 1]))); ((prod(i=1, #f~, A002110(primepi(f[i, 1]))^(f[i, 2]-(#f~==i)))%p)*bigomega(n)));
CROSSREFS
KEYWORD
nonn,less
AUTHOR
Antti Karttunen, Jan 15 2020
STATUS
approved