cp's OEIS Frontend

A342459 a(n) = gcd(A048250(n), A342001(n)).

%I A342459 #10 Mar 29 2021 15:05:01
%S A342459 1,1,1,1,1,1,1,3,2,1,1,4,1,3,8,1,1,1,1,6,2,1,1,1,2,3,1,8,1,1,1,1,2,1,
%T A342459 12,2,1,3,8,1,1,1,1,12,1,1,1,2,2,9,4,14,1,3,8,1,2,1,1,2,1,3,1,3,6,1,1,
%U A342459 18,2,1,1,1,1,3,1,20,6,1,1,2,4,1,1,2,2,3,8,1,1,1,4,24,2,1,24,1,1,1,1,2,1,1,1,1,1
%N A342459 a(n) = gcd(A048250(n), A342001(n)).
%H A342459 Antti Karttunen, <a href="/A342459/b342459.txt">Table of n, a(n) for n = 1..65537</a>
%F A342459 a(n) = A342458(n) / A003557(n) = gcd(A048250(n), A342001(n)).
%F A342459 a(n) = A342001(n) / A342919(n).
%o A342459 (PARI)
%o A342459 A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
%o A342459 A003557(n) = (n/factorback(factorint(n)[, 1]));
%o A342459 A342001(n) = (A003415(n) / A003557(n));
%o A342459 A048250(n) = factorback(apply(p -> p+1,factor(n)[,1]));
%o A342459 A342459(n) = gcd(A048250(n), A342001(n));
%Y A342459 Cf. A001615, A003415, A003557, A048250, A342001, A342458, A342919.
%Y A342459 Cf. also A066086, A342416.
%K A342459 nonn
%O A342459 1,8
%A A342459 _Antti Karttunen_, Mar 28 2021