login
a(n) = [nr+kr]-[nr]-[kr], where r=sqrt(2), k=3, [ ]=floor.
4

%I #11 Sep 08 2022 08:45:56

%S 0,1,0,0,0,0,1,0,0,0,0,1,0,1,0,0,0,0,1,0,0,0,0,1,0,1,0,0,0,0,1,0,0,0,

%T 0,1,0,0,0,0,1,0,1,0,0,0,0,1,0,0,0,0,1,0,1,0,0,0,0,1,0,0,0,0,1,0,0,0,

%U 0,1,0,1,0,0,0,0,1,0,0,0,0,1,0,1,0,0,0,0,1,0,0,0,0,1,0,1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1,0,1,0,0,0,0,1,0,0,0,0,1,0

%N a(n) = [nr+kr]-[nr]-[kr], where r=sqrt(2), k=3, [ ]=floor.

%C See A187950.

%H G. C. Greubel, <a href="/A187969/b187969.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = [(n+3)*r] - [n*r] - [3*r], where r=sqrt(2).

%t r=2^(1/2);

%t seqA=Table[Floor[(n+3)r]-Floor[n*r]-Floor[3r], {n,1,220}] (* A187969 *)

%t Flatten[Position[seqA,0] ] (* A187970 *)

%t Flatten[Position[seqA,1] ] (* A187971 *)

%o (PARI) for(n=1,30, print1(floor((n+3)*sqrt(2)) - floor(n*sqrt(2)) - floor(3*sqrt(2)), ", ")) \\ _G. C. Greubel_, Jan 31 2018

%o (Magma) [Floor((n+3)*Sqrt(2)) - Floor(n*Sqrt(2)) - Floor(3*Sqrt(2)): n in [1..30]]; // _G. C. Greubel_, Jan 31 2018

%Y Cf. A187950, A187970, A187971.

%K nonn

%O 1

%A _Clark Kimberling_, Mar 17 2011