A351084 - OEIS

A351084

a(n) = gcd(n, A328572(n)), where A328572 converts the primorial base expansion of n into its prime product form, but with 1 subtracted from all nonzero digits.

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 5, 1, 1, 1, 1, 25, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 25, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 5, 1, 7, 1, 1, 5, 1, 1, 1, 7, 5, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 1, 49, 1, 1, 1, 1, 1, 1, 35

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

OFFSET

0,16

LINKS

Antti Karttunen, Table of n, a(n) for n = 0..65537

Index entries for sequences related to primorial base

FORMULA

a(n) = gcd(n, A328572(n)) = gcd(A324198(n), A351083(n)).

a(n) = gcd(n, A085731(A276086(n))) = gcd(n, A276086(n), A327860(n)).

PROG

(PARI)

A328572(n) = { my(m=1, p=2); while(n, if(n%p, m *= p^((n%p)-1)); n = n\p; p = nextprime(1+p)); (m); };