A354946 a(n) = gcd(A001923(n),n^n).

1, 1, 1, 32, 1, 1, 1, 4, 1, 1, 1, 8, 1, 1, 9, 4, 17, 3, 19, 16, 7, 1, 1, 4, 1, 1, 1, 8, 1, 1, 1, 4, 3, 17, 1, 768, 1, 1, 1, 4, 1, 1, 1, 8, 1, 1, 1, 4, 1, 125, 3, 16, 1, 9, 1, 4, 1, 1, 1, 8, 1, 1, 1, 4, 5, 1, 1, 32, 9, 1, 1, 12, 1, 1, 5, 8, 1, 1, 1, 4, 1, 1

Offset: 1

Author

José María Grau Ribas, Jun 12 2022

Keywords

• nonn

Links

• Michel Marcus, Table of n, a(n) for n = 1..5000

Crossrefs

Cf. A000312, A001923.

Programs

Mathematica

A001923[n_]:=Sum[k^k,{k,1,n}]; A[n_]:=GCD[A001923[n],n^n]; Table[A[n],{n,1,33}]

PARI

a(n) = gcd(sum(k=1, n, k^k), n^n); \\ Michel Marcus, Jul 11 2022