login
a(n) = gcd(A046523(n), A071364(n)).
4

%I #9 Aug 24 2017 13:38:32

%S 1,2,2,4,2,6,2,8,4,6,2,12,2,6,6,16,2,6,2,12,6,6,2,24,4,6,8,12,2,30,2,

%T 32,6,6,6,36,2,6,6,24,2,30,2,12,12,6,2,48,4,6,6,12,2,6,6,24,6,6,2,60,

%U 2,6,12,64,6,30,2,12,6,30,2,72,2,6,6,12,6,30,2,48,16,6,2,60,6,6,6,24,2,30,6

%N a(n) = gcd(A046523(n), A071364(n)).

%H Antti Karttunen, <a href="/A083260/b083260.txt">Table of n, a(n) for n = 1..10000</a>

%o (PARI)

%o A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); }; \\ This function from _Charles R Greathouse IV_, Aug 17 2011

%o A071364(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = prime(i)); factorback(f); }; \\ This function from _Michel Marcus_, Jun 13 2014

%o A083260(n) = gcd(A046523(n), A071364(n)); \\ _Antti Karttunen_, Aug 24 2017

%Y Cf. A046523, A071364, A083255, A083256, A083257, A083258, A083259, A083261.

%K nonn

%O 1,2

%A _Labos Elemer_, May 09 2003