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

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, 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, 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

Offset

1,2

Links

Antti Karttunen, Table of n, a(n) for n = 1..10000

Mathematica

Table[GCD[n, Times @@ MapIndexed[Prime[First[#2]]^#1 &, Sort[FactorInteger[n][[All, -1]], Greater]]], {n, 98}] (* Michael De Vlieger, May 21 2017 *)

Prog

(PARI)
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
A083258(n) = gcd(n, A046523(n)); \\ Antti Karttunen, May 21 2017

Crossrefs

Cf. A046523, A071364, A083255-A083260, A083261.

Sequence in context: A340346 A193267 A327832 * A083259 A342242 A214060

Adjacent sequences: A083255 A083256 A083257 * A083259 A083260 A083261

Keyword

nonn

Author

Labos Elemer, May 09 2003

Status

approved