%I #26 Jan 01 2024 14:01:33
%S 25,49,75,98,121,125,169,175,196,225,242,245,275,289,325,338,343,361,
%T 363,375,392,425,475,484,490,507,525,529,539,575,578,625,637,675,676,
%U 686,722,725,726,775,784,825,833,841,845,847,867,875,925,931,961,968
%N Numbers n such that gcd(n, numerator(B(2n))) is not squarefree where B(2n) are the Bernoulli numbers.
%C All members so far are themselves not squarefree.  _Robert G. Wilson v_, Aug 28 2004
%t Select[ Range[1000], !SquareFreeQ[ GCD[ #, Numerator[ BernoulliB[2# ]] ]] &] (* _Robert G. Wilson v_, Aug 28 2004 *)
%Y Cf. A097704.
%K nonn
%O 1,1
%A _Benoit Cloitre_, Feb 06 2002
%E More terms from _Robert G. Wilson v_, Aug 28 2004
%E 'Needs["NumberTheory`"]' (no longer needed) removed from Mathematica code by _JeanFrançois Alcover_, Dec 04 2013
