A189458 a(n) = n+[nr/s]+[nt/s]; r=2, s=sqrt(2), t=1+sqrt(2).
3, 7, 12, 15, 20, 24, 27, 32, 36, 41, 44, 48, 53, 56, 61, 65, 70, 73, 77, 82, 85, 90, 94, 97, 102, 106, 111, 114, 119, 123, 126, 131, 135, 140, 143, 147, 152, 155, 160, 164, 167, 172, 176, 181, 184, 189, 193, 196, 201, 205, 210, 213, 217, 222, 225, 230, 234, 239, 242, 246, 251, 254, 259, 263, 266, 271, 275, 280, 283, 287, 292, 295, 300, 304, 309, 312, 316, 321, 324, 329, 333, 336, 341
Offset: 1
Keywords
Links
- G. C. Greubel, Table of n, a(n) for n = 1..10000
Programs
-
Magma
[n+Floor((2*n)/Sqrt(2))+Floor((n*(1+Sqrt(2)))/Sqrt(2)): n in [1..120]]; // G. C. Greubel, Aug 19 2018
-
Mathematica
r=2; s=2^(1/2); t=1+2^(1/2); a[n_] := n + Floor[n*s/r] + Floor[n*t/r]; b[n_] := n + Floor[n*r/s] + Floor[n*t/s]; c[n_] := n + Floor[n*r/t] + Floor[n*s/t]; Table[a[n], {n, 1, 120}] (*A189457*) Table[b[n], {n, 1, 120}] (*A189458*) Table[c[n], {n, 1, 120}] (*A186222, conjectured*) Table[n+Floor[(2n)/Sqrt[2]]+Floor[(n(1+Sqrt[2]))/Sqrt[2]],{n,90}] (* Harvey P. Dale, Feb 04 2015 *) -
PARI
vector(120, n, n+floor((2*n)/sqrt(2))+floor((n*(1+sqrt(2)))/sqrt(2))) \\ G. C. Greubel, Aug 19 2018
Comments