cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A083258 a(n) = gcd(A046523(n), n).

This page as a plain text file.
%I A083258 #15 May 21 2017 14:31:58
%S A083258 1,2,1,4,1,6,1,8,1,2,1,12,1,2,3,16,1,6,1,4,3,2,1,24,1,2,1,4,1,30,1,32,
%T A083258 3,2,1,36,1,2,3,8,1,6,1,4,3,2,1,48,1,2,3,4,1,6,1,8,3,2,1,60,1,2,3,64,
%U A083258 1,6,1,4,3,10,1,72,1,2,3,4,1,6,1,16,1,2,1,12,1,2,3,8,1,30,1,4,3,2,1,96,1,2
%N A083258 a(n) = gcd(A046523(n), n).
%H A083258 Antti Karttunen, <a href="/A083258/b083258.txt">Table of n, a(n) for n = 1..10000</a>
%t A083258 Table[GCD[n, Times @@ MapIndexed[Prime[First[#2]]^#1 &, Sort[FactorInteger[n][[All, -1]], Greater]]], {n, 98}] (* _Michael De Vlieger_, May 21 2017 *)
%o A083258 (PARI)
%o A083258 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 A083258 A083258(n) = gcd(n,A046523(n)); \\ _Antti Karttunen_, May 21 2017
%Y A083258 Cf. A046523, A071364, A083255-A083260, A083261.
%K A083258 nonn
%O A083258 1,2
%A A083258 _Labos Elemer_, May 09 2003