OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..100
MATHEMATICA
usigma[n_] := Block[{d = Divisors[n]}, Plus @@ Select[d, GCD[ #, n/# ] == 1 &]]; lmt = 1296000; t = (Select[ Range[ lmt], DivisorSigma[1, # ] == 2usigma[ # ] &] - 108)/216; u = (Select[ Range[ Floor[(lmt - 108)/432]], !SquareFreeQ[ GCD[ #, Numerator[ BernoulliB[ 2# ]] ]] &] -1)/2; v = Table[ 3k - 2, {k, Floor[(lmt - 108)/216]}]; Complement[ Range[ Floor[ (lmt - 108)/216]], t, u, v]
q[n_] := Mod[n, 3] != 1 && (Divisible[2*n + 1, 3] || (! Divisible[2*n + 1, 3] && ! SquareFreeQ[2*n + 1])) && SquareFreeQ[GCD[2*n + 1, BernoulliB[4*n + 2]]]; Select[Range[10^4, q] (* Amiram Eldar, Aug 31 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Ralf Stephan and Robert G. Wilson v, Sep 15 2004
EXTENSIONS
More terms from Amiram Eldar, Aug 31 2024
STATUS
approved