OFFSET
1,2
COMMENTS
MATHEMATICA
r = Sqrt[6]; 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, 500}];
Flatten[Position[t1, 1]] (* empty *)
t2 = Table[If[y[n] == x[n], 1, 0], {n, 1, 400}];
Flatten[Position[t2, 1]] (* A194376 *)
t3 = Table[If[y[n] > x[n], 1, 0], {n, 1, 100}];
Flatten[Position[t3, 1]] (* A194377 *)
PROG
(PARI) is(n)=my(r=sqrt(6), f=x->x-x\1); sum(k=1, n, f(1/2+k*r)-f(k*r))>0 \\ Charles R Greathouse IV, Jul 25 2012
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Aug 23 2011
STATUS
approved