A346099 a(n) = gcd(n, A346098(n)).

1, 2, 3, 4, 5, 3, 7, 1, 1, 5, 11, 3, 13, 7, 5, 1, 17, 1, 19, 5, 7, 11, 23, 3, 1, 13, 3, 7, 29, 5, 31, 2, 11, 17, 7, 1, 37, 19, 13, 1, 41, 1, 43, 11, 1, 23, 47, 1, 1, 1, 17, 13, 53, 9, 11, 7, 19, 29, 59, 1, 61, 31, 1, 4, 13, 11, 67, 17, 23, 1, 71, 3, 73, 37, 25, 19, 1, 13, 79, 1, 1, 41, 83, 1, 17, 43, 29, 11, 89

Offset

1,2

Comments

Only powers of primes (A000961) occur as terms. A346100 lists the exponents.

Links

Table of n, a(n) for n=1..89.

Index entries for sequences computed from indices in prime factorization

Index entries for sequences related to primorial base

Index entries for sequences related to primorial numbers

Formula

a(n) = gcd(n, A346098(n)) = gcd(n, A064989(A319626(A324886(n)))).

Prog

(PARI)
A064989(n) = { my(f = factor(n)); if((n>1 && f[1, 1]==2), f[1, 2] = 0); for (i=1, #f~, f[i, 1] = precprime(f[i, 1]-1)); factorback(f); };
A319626(n) = (n / gcd(n, A064989(n)));
A346099(n) = gcd(n, A064989(A319626(A324886(n)))); \\ Rest of program in A324886.

Crossrefs

Cf. A000961, A064989, A319626, A324886, A346095, A346096, A346097, A346098, A346100.

Cf. A346090 (positions of ones).

Sequence in context: A278057 A366876 A330750 * A345941 A088492 A345931

Adjacent sequences: A346096 A346097 A346098 * A346100 A346101 A346102

Keyword

nonn

Author

Antti Karttunen, Jul 07 2021

Status

approved