A328323 a(n) = gcd(A328316(n),A328316(1+n)).

1, 1, 1, 3, 1, 1, 1, 1, 2401, 861018571670257

Offset: 0

Author

Antti Karttunen, Oct 14 2019

Keywords

• nonn

Comments

All terms are odd, because the parity of A328316 alternates.

a(10) is 208 decimal digits long, and can be found in the b-file.

Links

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

Crossrefs

Cf. A276086, A324198, A328316, A328317.

Programs

PARI

A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A328316(n) = if(!n,0,A276086(A328316(n-1)));
A328323(n) = gcd(A328316(n),A328316(1+n));
\\ Or alternatively, as:
A324198(n) = gcd(n, A276086(n));
A328323(n) = A324198(A328316(n));

Formula

a(n) = A324198(A328316(n)) = gcd(A328316(n),A328316(1+n)).