A322362 - OEIS

A322362

a(n) = gcd(n, A166590(n)), where A166590 is completely multiplicative with a(p) = p+2 for prime p.

1, 2, 1, 4, 1, 2, 1, 8, 1, 2, 1, 4, 1, 2, 5, 16, 1, 2, 1, 4, 3, 2, 1, 8, 1, 2, 1, 4, 1, 10, 1, 32, 1, 2, 7, 4, 1, 2, 3, 8, 1, 6, 1, 4, 5, 2, 1, 16, 1, 2, 1, 4, 1, 2, 1, 8, 3, 2, 1, 20, 1, 2, 9, 64, 5, 2, 1, 4, 1, 14, 1, 8, 1, 2, 5, 4, 1, 6, 1, 16, 1, 2, 1, 12, 1, 2, 1, 8, 1, 10, 1, 4, 3, 2, 1, 32, 1, 2, 1, 4, 1, 2, 1, 8, 105

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,2

LINKS

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

Antti Karttunen, Data supplement: n, a(n) computed for n = 1..100000

FORMULA

a(n) = gcd(n, A166590(n)).

a(A037074(n)) = A006512(n).

MATHEMATICA

a[n_] := If[n == 1, 1, GCD[n, Times@@ ((First[#]+2)^Last[#] &/@FactorInteger[n])]]; Array[a, 120] (* Amiram Eldar, Dec 05 2018~ *)

PROG