|
|
A126196
|
|
Numbers k such that gcd(A001008(k), A001008(floor(k/2))) > 1.
|
|
12
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
|
|
LINKS
|
|
|
MATHEMATICA
|
Select[Range[5000], GCD @@ Numerator@ HarmonicNumber@{#, Floor[#/2]} > 1 &] (* Giovanni Resta, May 13 2016 *)
|
|
PROG
|
(PARI) a001008(n)=numerator(sum(i=1, n, 1/i))
for(n=1, 1e6, if(gcd(a001008(n), a001008(n/2)) > 1, print1(n, ", "))) \\ Felix Fröhlich, Aug 08 2014
|
|
CROSSREFS
|
The corresponding GCDs are given by A126197.
|
|
KEYWORD
|
nonn,more
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|