amax = 25; km0 = 10; Clear[cc]; L[a_, s_, km_] := Sum[ JacobiSymbol[ -a, 2 k + 1]/(2 k + 1)^s, {k, 0, km}]; c[1, n_, km_] := 2 (2 n)! L[1, 2 n + 1, km] (2/Pi)^(2 n + 1) // Round; c[a_ /; a > 1, n_, km_] := (2 n)! L[a, 2 n + 1, km] (2 a/Pi)^(2 n + 1)/Sqrt[a] // Round; cc[km_] := cc[km] = Table[ c[a, 3, km], {a, 1, amax} ]; cc[km0]; cc[km = 2 km0]; While[cc[km] != cc[km/2, km = 2 km]];
A000508 = cc[km] (*
Jean-François Alcover, Feb 09 2016 *)