cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A248633 Least k such that 20/27- sum{(h^2)/4^h, h = 1..k} < 1/8^n.

Original entry on oeis.org

3, 5, 7, 9, 11, 12, 14, 16, 17, 19, 21, 22, 24, 25, 27, 29, 30, 32, 33, 35, 36, 38, 40, 41, 43, 44, 46, 47, 49, 50, 52, 53, 55, 57, 58, 60, 61, 63, 64, 66, 67, 69, 70, 72, 73, 75, 76, 78, 80, 81, 83, 84, 86, 87, 89, 90, 92, 93, 95, 96, 98, 99, 101, 102, 104
Offset: 1

Views

Author

Clark Kimberling, Oct 11 2014

Keywords

Comments

This sequence provides insight into the manner of convergence of sum{(h^2)/4^h, h = 1..k} to 20/27. The difference sequence of A248633 entirely of 1s and 2s, so that A248634 and A248635 partition the positive integers.

Examples

			Let s(n) = 20/27 - sum{(h^2)/4^h, h = 1..n}.  Approximations follow:
n ... s(n) ........ 1/8^n
1 ... 0.49074 ... 0.125000
2 ... 0.24074 ... 0.015625
3 ... 0.10011 ... 0.001953
4 ... 0.03761 ... 0.000244
5 ... 0.01320 ... 0.000030
a(2) = 5 because s(5) < 1/8^2 < s(2).
		

Crossrefs

Programs

  • Mathematica
    z= 300; p[k_] := p[k] = Sum[(h^2/4^h), {h, 1, k}];
    d = N[Table[20/27 - p[k], {k, 1, z/5}], 12];
    f[n_] := f[n] = Select[Range[z], 20/27 - p[#] < 1/8^n &, 1];
    u = Flatten[Table[f[n], {n, 1, z}]]  (* A248633 *)
    d = Differences[u]
    Flatten[Position[d, 1]]  (* A248634 *)
    Flatten[Position[d, 2]]  (* A248635 *)

A248634 Numbers k such that A248633(k+1) = A248633(k) + 1.

Original entry on oeis.org

5, 8, 11, 13, 16, 18, 20, 23, 25, 27, 29, 31, 34, 36, 38, 40, 42, 44, 46, 49, 51, 53, 55, 57, 59, 61, 63, 65, 67, 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 90, 92, 94, 96, 98, 101, 103, 105, 107, 109, 111, 113, 115, 117, 119, 121, 123, 125, 127, 129, 131, 133
Offset: 1

Views

Author

Clark Kimberling, Oct 11 2014

Keywords

Examples

			(A248633(k+1) = A248633(k)) = (2, 2, 2, 2, 1, 2, 2, 1, 2, 2, 1, 2, 1, 2, 2, 1,... ), so that A248634 = (5, 8, 11, 13, ... ) and A248635 = (1, 2, 3, 4, 6, 7, 9, 10, 12, ...).
		

Crossrefs

Programs

  • Mathematica
    z= 300; p[k_] := p[k] = Sum[(h^2/4^h), {h, 1, k}];
    d = N[Table[20/27 - p[k], {k, 1, z/5}], 12];
    f[n_] := f[n] = Select[Range[z], 20/27 - p[#] < 1/8^n &, 1];
    u = Flatten[Table[f[n], {n, 1, z}]]  (* A248633 *)
    d = Differences[u]
    Flatten[Position[d, 1]]  (* A248634 *)
    Flatten[Position[d, 2]]  (* A248635 *)

A248637 Numbers k such that A248636(k+1) = A248636(k) + 1.

Original entry on oeis.org

2, 4, 7, 8, 10, 12, 14, 15, 17, 19, 20, 22, 23, 25, 26, 28, 29, 31, 32, 34, 35, 37, 38, 40, 41, 43, 44, 46, 47, 48, 50, 51, 53, 54, 56, 57, 59, 60, 61, 63, 64, 66, 67, 69, 70, 71, 73, 74, 76, 77, 78, 80, 81, 83, 84, 86, 87, 88, 90, 91, 93, 94, 95, 97, 98
Offset: 1

Views

Author

Clark Kimberling, Oct 11 2014

Keywords

Examples

			(A248636(k+1) = A248636(k)) = (2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 1, 1, 2,... ), so that A248637 = (2, 4, 7, 8, 10, 12, 14, ... ) and A248638 = (1, 3, 5, 6, 9, 11, 13, ...).
		

Crossrefs

Programs

  • Mathematica
    z = 300; p[k_] := p[k] = Sum[(h^3/3^h), {h, 1, k}];
    d = N[Table[33/8 - p[k], {k, 1, z/5}], 12]
    f[n_] := f[n] = Select[Range[z], 33/8 - p[#] < 1/4^n &, 1];
    u = Flatten[Table[f[n], {n, 1, z}]]  (* A248636 *)
    d = Differences[u]
    v = Flatten[Position[d, 1]] (* A248637 *)
    w = Flatten[Position[d, 2]] (* A248638 *)
Showing 1-3 of 3 results.