Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A336697
a(n) = gcd((s+1)/2, 1+sigma(s)), where s is the n-th odd square, (2n-1)^2.
5
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 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, 1, 1, 1, 5, 1, 1, 13, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
OFFSET
1,17
COMMENTS
See comments in A336700, A337337, and A337339.
FORMULA
For all n >= 1, a(A048673(n)) = A337337(n).
MATHEMATICA
Array[GCD[(# + 1)/2, 1 + DivisorSigma[1, #]] &[(2 # - 1)^2] &, 120] (* Michael De Vlieger, Aug 24 2020 *)
PROG
(PARI) A336697(n) = { my(s=((n+n-1)^2)); gcd((s+1)/2, 1+sigma(s)); };
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Aug 18 2020
STATUS
approved