login
A377990
a(n) = sigma(sigma(A350388(n))) * sigma(sigma(A350389(n))), where A350388 and A350389 are the largest unitary divisor of n that is a square, and the largest unitary divisor of n that is an exponentially odd number, respectively.
2
1, 4, 7, 8, 12, 28, 15, 24, 14, 39, 28, 56, 24, 60, 60, 32, 39, 56, 42, 96, 63, 91, 60, 168, 32, 96, 90, 120, 72, 195, 63, 104, 124, 120, 124, 112, 60, 168, 120, 234, 96, 252, 84, 224, 168, 195, 124, 224, 80, 128, 195, 192, 120, 360, 195, 360, 186, 234, 168, 480, 96, 252, 210, 128, 224, 403, 126, 312, 252, 403, 195
OFFSET
1,2
COMMENTS
Differs from A051027 at 52, 98, 156, 164, 245, ..., = A377991.
FORMULA
a(n) = A051027(A350388(n)) * A051027(A350389(n)).
a(n) = sigma(A351568(n)) * sigma(A351569(n)).
PROG
(PARI)
A350388(n) = { my(m=1, f=factor(n)); for(k=1, #f~, if(0==(f[k, 2]%2), m *= (f[k, 1]^f[k, 2]))); (m); };
A350389(n) = { my(m=1, f=factor(n)); for(k=1, #f~, if(1==(f[k, 2]%2), m *= (f[k, 1]^f[k, 2]))); (m); };
A377990(n) = (sigma(sigma(A350388(n))) * sigma(sigma(A350389(n))));
KEYWORD
nonn
AUTHOR
Antti Karttunen, Nov 15 2024
STATUS
approved