login
A346099
a(n) = gcd(n, A346098(n)).
4
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.
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. A346090 (positions of ones).
Sequence in context: A278057 A366876 A330750 * A345941 A088492 A345931
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jul 07 2021
STATUS
approved