login
A326062
a(1) = gcd((sigma(n)-A032742(n))-n, n-A032742(n)), where A032742 gives the largest proper divisor of n.
15
1, 1, 2, 1, 4, 3, 6, 1, 1, 1, 10, 2, 12, 1, 2, 1, 16, 3, 18, 2, 2, 1, 22, 12, 1, 1, 2, 14, 28, 3, 30, 1, 2, 1, 2, 1, 36, 1, 2, 10, 40, 3, 42, 2, 6, 1, 46, 4, 1, 1, 2, 2, 52, 3, 2, 4, 2, 1, 58, 6, 60, 1, 2, 1, 2, 3, 66, 2, 2, 1, 70, 3, 72, 1, 2, 2, 2, 3, 78, 2, 1, 1, 82, 14, 2, 1, 2, 4, 88, 9, 2, 2, 2, 1, 2, 12, 96, 1, 6, 1, 100, 3, 102, 2, 2
OFFSET
1,3
COMMENTS
See comments in A326063 and A326064.
LINKS
FORMULA
a(1) = 1; for n > 1, a(n) = gcd(A060681(n), A318505(n)).
a(n) = gcd((A000203(n)-A032742(n))-n, n-A032742(n)).
PROG
(PARI)
A032742(n) = if(1==n, n, n/vecmin(factor(n)[, 1]));
A326062(n) = gcd((sigma(n)-A032742(n))-n, n-A032742(n));
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jun 06 2019
STATUS
approved