A306507 a(n) = gcd(n!^2+1, sigma(n!)), where sigma() denotes the sum of the divisors.
1, 1, 1, 1, 1, 13, 1, 17, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 61, 1, 1, 1, 1, 1, 1, 1, 1, 61, 1, 1, 1, 193, 1, 1, 1, 757, 61, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 109, 1, 1, 1, 181, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 113
Offset: 1
Keywords
Programs
-
Mathematica
Table[GCD[(n!)^2+1,DivisorSigma[1,n!]],{n,90}] (* Harvey P. Dale, Jun 03 2021 *)
-
PARI
a(n) = gcd(n!^2+1, sigma(n!)); \\ Michel Marcus, Feb 20 2019
Extensions
More terms from Michel Marcus, Feb 20 2019
Comments