|
|
A373485
|
|
a(n) = gcd(A083345(n), A276085(n)), where A276085 is fully additive with a(p) = p#/p, and A083345 is the numerator of the fully additive function with a(p) = 1/p.
|
|
5
|
|
|
0, 1, 1, 1, 1, 1, 1, 3, 2, 7, 1, 4, 1, 1, 8, 2, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 8, 1, 1, 1, 5, 2, 1, 12, 1, 1, 1, 8, 1, 1, 1, 1, 4, 1, 1, 1, 1, 2, 1, 4, 2, 1, 1, 8, 1, 2, 1, 1, 1, 1, 1, 17, 3, 6, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 4, 6, 1, 1, 1, 4, 1, 1, 1, 2, 1, 8, 1, 1, 1, 20, 4, 2, 1, 12, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,8
|
|
COMMENTS
|
For all n >= 1, A373145(n) is a multiple of a(n).
|
|
LINKS
|
|
|
PROG
|
(PARI)
A083345(n) = { my(f=factor(n)); numerator(vecsum(vector(#f~, i, f[i, 2]/f[i, 1]))); };
A276085(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*prod(i=1, primepi(f[k, 1]-1), prime(i))); };
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|