|
%I #10 Feb 02 2018 16:20:17
%S 0,2,1,3,0,2,6,2,7,1,6,12,5,12,4,11,1,8,17,6,15,2,12,23,9,21,5,17,0,
%T 12,26,8,22,2,16,32,11,28,5,22,40,16,35,10,29,2,21,42,14,36,6,28,51,
%U 20,44,11,35,1,25,51,16,42,5,32,60,22,51,11,40,70,29,60,18,49,5,36
%N Number of k in [1,n] for which <n*r>+<k*r> > 1, where < > = fractional part and r = sqrt(2); row sums of A164664.
%H G. C. Greubel, <a href="/A194665/b194665.txt">Table of n, a(n) for n = 1..10000</a>
%t r = Sqrt[2]; z = 13;
%t p[x_] := FractionalPart[x]; f[x_] := Floor[x];
%t w[n_, k_] := p[r^n] + p[r^k] - p[r^n + r^k]
%t Flatten[Table[w[n, k], {n, 1, z}, {k, 1, n}]]
%t TableForm[Table[w[n, k], {n, 1, z}, {k, 1, n}]]
%t s[n_] := Sum[w[n, k], {k, 1, n}] (* A194663 *)
%t Table[s[n], {n, 1, 100}]
%t h[n_, k_] := f[p[n*r] + p[k*r]]
%t Flatten[Table[h[n, k], {n, 1, z}, {k, 1, n}]]
%t (* A194664 *)
%t TableForm[Table[h[n, k], {n, 1, z}, {k, 1, n}]]
%t t[n_] := Sum[h[n, k], {k, 1, n}]
%t Table[t[n], {n, 1, 100}] (* A194665 *)
%Y Cf. A194664.
%K nonn
%O 1,2
%A _Clark Kimberling_, Sep 01 2011
|
