|
%I #15 Sep 08 2022 08:45:57
%S 1,1,0,1,1,0,1,0,0,1,0,0,0,1,1,0,1,1,0,1,0,0,1,0,0,0,1,1,0,1,1,0,1,0,
%T 0,1,0,0,0,1,1,0,1,1,0,1,1,0,1,0,0,0,1,1,0,1,1,0,1,1,0,1,0,0,0,1,1,0,
%U 1,1,0,1,1,0,0,0,0,0,1,1,1,1,1,0,1,1,0,0,0,0,0,1,1,1,1,1,0,1,1,0,0,0,0,0,0,1,1,1,1,0,1,1,0,0,0,0,0,0,1,1,1,1,0,1,1,0,0,1,0,0,0
%N a(n) = [n*u+n*v] - [n*u] - [n*v], where u=sin(Pi/8), v=cos(Pi/8), and []=floor.
%H G. C. Greubel, <a href="/A190233/b190233.txt">Table of n, a(n) for n = 1..10000</a>
%t u = Sin[Pi/8]; v = Cos[Pi/8];
%t f[n_] := Floor[n*u + n*v] - Floor[n*u] - Floor[n*v]
%t t = Table[f[n], {n, 1, 120}] (* A190233 *)
%t Flatten[Position[t, 0]] (* A190234 *)
%t Flatten[Position[t, 1]] (* A190235 *)
%o (PARI) for(n=1,30, print1(floor(n*(sin(Pi/8) + cos(Pi/8))) - floor(n*cos(Pi/8)) - floor(n*sin(Pi/8)), ", ")) \\ _G. C. Greubel_, Dec 27 2017
%o (Magma) C<i> := ComplexField(); [Floor(n*(Sin(Pi(C)/8) + Cos(Pi(C)/8))) - Floor(n*Sin(Pi(C)/8)) - Floor(n*Cos(Pi(C)/8)): n in [1..30]]; // _G. C. Greubel_, Dec 27 2017
%Y Cf. A190234, A190235.
%K nonn
%O 1
%A _Clark Kimberling_, May 06 2011
|
