%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, !SquareFreeQ[ GCD[ #, Numerator[ BernoulliB[2# ]] ]] &] (* _Robert G. Wilson v_, Aug 28 2004 *)
%Y Cf. A097704.
%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 _Jean-François Alcover_, Dec 04 2013