Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #15 May 29 2021 20:01:10
%S 1,1,1,7,1,1,1,5,13,1,1,7,1,1,1,31,1,13,1,7,1,1,1,5,31,1,10,7,1,1,1,
%T 21,1,1,1,91,1,1,1,5,1,1,1,7,13,1,1,31,57,31,1,7,1,10,1,5,1,1,1,7,1,1,
%U 13,127,1,1,1,7,1,1,1,65,1,1,31,7,1,1,1,31,121,1,1,7,1,1,1,5,1,13,1,7,1,1,1,21,1,57,13
%N a(n) = sigma(n) / gcd(sigma(n), A001615(n)).
%C This is not multiplicative. The first point where a(m*n) = a(m)*a(n) does not hold for coprime m and n is 108 = 4*27, where a(108) = 35, although a(4) = 7 and a(27) = 10. See A344702.
%H Antti Karttunen, <a href="/A344696/b344696.txt">Table of n, a(n) for n = 1..10404</a>
%H Antti Karttunen, <a href="/A344696/a344696.txt">Data supplement: n, a(n) computed for n = 1..65537</a>
%H <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a>
%F a(n) = A000203(n) / A344695(n).
%o (PARI)
%o A001615(n) = if(1==n,n, my(f=factor(n)); prod(i=1, #f~, f[i, 1]^f[i, 2] + f[i, 1]^(f[i, 2]-1))); \\ After code in A001615
%o A344696(n) = { my(u=sigma(n)); (u/gcd(u,A001615(n))); };
%Y Cf. A000203, A001615, A005117 (positions of ones), A344695, A344697, A344698, A344702.
%Y Cf. also A344756.
%K nonn
%O 1,4
%A _Antti Karttunen_, May 26 2021