A194383 a(n) = (1/2) * A194382(n).
2, 3, 5, 6, 8, 9, 11, 12, 14, 15, 17, 20, 23, 26, 29, 32, 52, 55, 58, 61, 64, 67, 69, 70, 72, 73, 75, 76, 78, 79, 81, 82, 84, 85, 86, 88, 89, 91, 92, 94, 95, 97, 98, 100, 101, 102, 104, 105, 107, 108, 110, 111, 113, 114, 116, 117, 119, 122, 125, 128, 131, 134
Offset: 1
Keywords
Programs
-
Mathematica
r = Sqrt[8]; c = 1/2; x[n_] := Sum[FractionalPart[k*r], {k, 1, n}] y[n_] := Sum[FractionalPart[c + k*r], {k, 1, n}] t1 = Table[If[y[n] < x[n], 1, 0], {n, 1, 100}]; Flatten[Position[t1, 1]] (* A194381 *) t2 = Table[If[y[n] == x[n], 1, 0], {n, 1, 300}]; Flatten[Position[t2, 1]] (* A194382 *) %/2 (* A194383 *) t3 = Table[If[y[n] > x[n], 1, 0], {n, 1, 600}]; Flatten[Position[t3, 1]] (* A194384 *)
Comments