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-8 of 8 results.

A130254 Duplicate of A130252.

Original entry on oeis.org

0, 1, 4, 7, 11, 15, 20, 25, 30, 35, 40, 45, 51, 57, 63, 69, 75, 81, 87, 93, 99, 105, 112, 119, 126, 133, 140, 147, 154, 161, 168, 175, 182, 189, 196, 203, 210, 217, 224, 231, 238, 245, 252, 259, 267, 275, 283, 291, 299, 307, 315, 323, 331, 339, 347, 355, 363, 371
Offset: 0

Views

Author

Keywords

A130248 Partial sums of the Lucas Inverse A130247.

Original entry on oeis.org

1, 1, 3, 6, 9, 12, 16, 20, 24, 28, 33, 38, 43, 48, 53, 58, 63, 69, 75, 81, 87, 93, 99, 105, 111, 117, 123, 129, 136, 143, 150, 157, 164, 171, 178, 185, 192, 199, 206, 213, 220, 227, 234, 241, 248, 255, 263, 271, 279, 287, 295, 303, 311, 319, 327, 335, 343, 351
Offset: 1

Views

Author

Hieronymus Fischer, May 19 2007

Keywords

Crossrefs

Other related sequences: A000032, A130241, A130242, A130243, A130244, A130245, A130246, A130251, A130252, A130257, A130261. Fibonacci inverse see A130233 - A130240, A104162.

Programs

  • Mathematica
    Join[{1, 1}, Table[Sum[Floor[Log[GoldenRatio, k + 1/2]], {k, 1, n}], {n, 3, 50}]] (* G. C. Greubel, Dec 24 2017 *)

Formula

a(n)=sum{1<=k<=n, A130247(k)}=2+(n+1)*A130247(n)-A000032(A130247(n)+2) for n>=3. G.f.: g(x)=1/(1-x)^2*(sum{k>=1, x^Lucas(k)}-x^2).

A130253 Number of Jacobsthal numbers (A001045) <=n.

Original entry on oeis.org

1, 3, 3, 4, 4, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9
Offset: 0

Views

Author

Hieronymus Fischer, May 20 2007

Keywords

Comments

Partial sums of the Jacobsthal indicator sequence (A105348).
For n<>1, we have a(A001045(n))=n+1.

Examples

			a(9)=5 because there are 5 Jacobsthal numbers <=9 (0,1,1,3 and 5).
		

Crossrefs

For partial sums see A130252. Other related sequences A001045, A130249, A130250, A130253, A105348. Also A130233, A130235, A130241, A108852, A130245.

Programs

Formula

a(n) = floor(log_2(3n+1)) + 1 = ceiling(log_2(3n+2)).
a(n) = A130249(n) + 1 = A130250(n+1).
G.f.: 1/(1-x)*(Sum_{k>=0} x^A001045(k)).

A130236 Partial sums of the 'upper' Fibonacci Inverse A130234.

Original entry on oeis.org

0, 1, 4, 8, 13, 18, 24, 30, 36, 43, 50, 57, 64, 71, 79, 87, 95, 103, 111, 119, 127, 135, 144, 153, 162, 171, 180, 189, 198, 207, 216, 225, 234, 243, 252, 262, 272, 282, 292, 302, 312, 322, 332, 342, 352, 362, 372, 382, 392, 402, 412, 422, 432, 442, 452, 462, 473
Offset: 0

Views

Author

Hieronymus Fischer, May 17 2007

Keywords

Crossrefs

Programs

  • Magma
    m:=120;
    f:= func< x | x*(&+[x^Fibonacci(j): j in [0..Floor(3*Log(3*m+1))]])/(1-x)^2 >;
    R:=PowerSeriesRing(Rationals(), m+1);
    [0] cat Coefficients(R!( f(x) )); // G. C. Greubel, Mar 18 2023
    
  • Mathematica
    b[n_]:= For[i=0, True, i++, If[Fibonacci[i] >= n, Return[i]]];
    b/@ Range[0, 56]//Accumulate (* Jean-François Alcover, Apr 13 2020 *)
  • SageMath
    m=120
    def f(x): return x*sum( x^fibonacci(j) for j in range(1+int(3*log(3*m+1))))/(1-x)^2
    def A130236_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( f(x) ).list()
    A130236_list(m) # G. C. Greubel, Mar 18 2023

Formula

a(n) = Sum_{k=0..n} A130234(k).
a(n) = n*A130234(n) - Fibonacci(A130234(n)+1) + 1.
G.f.: (x/(1-x)^2) * Sum_{k>=0} x^Fibonacci(k).

A130246 Partial sums of A130245.

Original entry on oeis.org

0, 1, 3, 6, 10, 14, 18, 23, 28, 33, 38, 44, 50, 56, 62, 68, 74, 80, 87, 94, 101, 108, 115, 122, 129, 136, 143, 150, 157, 165, 173, 181, 189, 197, 205, 213, 221, 229, 237, 245, 253, 261, 269, 277, 285, 293, 301, 310, 319, 328, 337, 346, 355, 364, 373, 382, 391
Offset: 0

Views

Author

Hieronymus Fischer, May 19 2007

Keywords

Crossrefs

Other related sequences: A000032, A130241, A130243, A130244, A130248, A130251, A130252, A130255, A130257, A130261. Fibonacci inverse see A130233 - A130240, A104162.

Programs

  • Magma
    [0] cat [(&+[1+Floor(Log((2*k+1)/2)/Log((1+Sqrt(5))/2)): k in [1..n]]): n in [1..100]]; // G. C. Greubel, Sep 09 2018
  • Mathematica
    Table[Sum[1 + Floor[Log[GoldenRatio, (2*k + 1)/2]], {k, 1, n}], {n, 0, 100}] (* G. C. Greubel, Sep 09 2018 *)
  • PARI
    for(n=0, 100, print1(sum(k=1,n, 1 + floor(log((2*k+1)/2)/log((1+sqrt(5))/2))), ", ")) \\ G. C. Greubel, Sep 09 2018
    

Formula

a(n) = Sum_{k=1..n} A130245(k).
a(n) = 1 +(n+1)*A130245(n) - A000032(A130245(n)+1) for n=0 or n >= 2.
G.f.: 1/(1-x)^2*Sum_{k>=0} x^A000032(k).

A130251 Partial sums of A130249.

Original entry on oeis.org

0, 2, 4, 7, 10, 14, 18, 22, 26, 30, 34, 39, 44, 49, 54, 59, 64, 69, 74, 79, 84, 90, 96, 102, 108, 114, 120, 126, 132, 138, 144, 150, 156, 162, 168, 174, 180, 186, 192, 198, 204, 210, 216, 223, 230, 237, 244, 251, 258, 265, 272, 279, 286, 293, 300, 307, 314, 321
Offset: 0

Views

Author

Hieronymus Fischer, May 20 2007

Keywords

Examples

			G.f. = 2*x + 4*x^2 + 7*x^3 + 10*x^4 + 14*x^5 + 18*x^6 + 22*x^7 + ... - _Michael Somos_, Sep 17 2018
		

Crossrefs

Programs

  • Magma
    [0] cat [(&+[Floor(Log(3*k+1)/Log(2)) : k in [1..n]]): n in [1..100]]; // G. C. Greubel, Sep 09 2018
    
  • Mathematica
    Join[{0}, Table[Sum[Floor[Log[2, 3*k + 1]], {k, 1, n}], {n, 1, 2500}]] (* G. C. Greubel, Sep 09 2018 *)
  • PARI
    for(n=0,100, print1(sum(k=1,n, floor(log(3*k+1)/log(2))), ", ")) \\ G. C. Greubel, Sep 09 2018
    
  • Python
    def A130251(n): return (n+1)*((m:=(3*n+1).bit_length())-1)-(((1<>1) # Chai Wah Wu, Apr 17 2025

Formula

a(n) = Sum_{k=0..n} A130249(k).
a(n) = (n+1)*floor(log_2(3*n+1)) - (1/2)*(A001045(floor(log_2(3*n+1))+2)-1).
G.f.: (1/(1-x)^2)*Sum_{k>=1} x^A001045(k).

A130244 Partial sums of the 'upper' Lucas Inverse A130242.

Original entry on oeis.org

0, 0, 0, 2, 5, 9, 13, 17, 22, 27, 32, 37, 43, 49, 55, 61, 67, 73, 79, 86, 93, 100, 107, 114, 121, 128, 135, 142, 149, 156, 164, 172, 180, 188, 196, 204, 212, 220, 228, 236, 244, 252, 260, 268, 276, 284, 292, 300, 309, 318, 327, 336, 345, 354, 363, 372, 381, 390
Offset: 0

Views

Author

Hieronymus Fischer, May 19 2007

Keywords

Crossrefs

Other related sequences: A000032, A130241, A130243, A130245, A130246, A130248, A130252, A130258, A130262. Fibonacci inverse see A130233 - A130240, A104162.

Programs

  • Magma
    [0,0] cat [(&+[Ceiling(Log(k + 1/2)/Log((1+Sqrt(5))/2)) : k in [0..n]]): n in [1..50]]; // G. C. Greubel, Sep 12 2018
  • Mathematica
    Join[{0, 0}, Table[Sum[Ceiling[Log[GoldenRatio, k + 1/2]], {k, 0, n}], {n, 1, 50}]] (* G. C. Greubel, Sep 12 2018 *)
  • PARI
    for(n=-1,50, print1(if(n==-1, 0, if(n==0, 0, sum(k=0, n, ceil(log(k + 1/2)/log((1+sqrt(5))/2))))), ", ")) \\ G. C. Greubel, Sep 12 2018
    

Formula

a(n) = Sum_{k=0..n} A130242(k).
a(n) = n*A130242(n) - A000032(A130242(n) +1) for n>=3.
G.f.: x/(1-x)^2*(2*x^2 + Sum{k>=2, x^Lucas(k)}).

A130250 Minimal index k of a Jacobsthal number such that A001045(k) >= n (the 'upper' Jacobsthal inverse).

Original entry on oeis.org

0, 1, 3, 3, 4, 4, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9
Offset: 0

Views

Author

Hieronymus Fischer, May 20 2007

Keywords

Comments

Inverse of the Jacobsthal sequence (A001045), nearly, since a(A001045(n))=n except for n=2 (see A130249 for another version). a(n+1) is equal to the partial sum of the Jacobsthal indicator sequence (see A105348).

Examples

			a(10)=5 because A001045(5) = 11 >= 10, but A001045(4) = 5 < 10.
		

Crossrefs

For partial sums see A130252.

Programs

  • Magma
    [0] cat [Ceiling(Log(2,3*n-1)): n in [1..120]]; // G. C. Greubel, Mar 18 2023
    
  • Mathematica
    Table[If[n==0, 0, Ceiling[Log[2, 3*n-1]]], {n,0,120}] (* G. C. Greubel, Mar 18 2023 *)
  • Python
    def A130250(n): return (3*n-2).bit_length() if n else 0 # Chai Wah Wu, Apr 17 2025
  • SageMath
    def A130250(n): return 0 if (n==0) else ceil(log(3*n-1, 2))
    [A130250(n) for n in range(121)] # G. C. Greubel, Mar 18 2023
    

Formula

a(n) = ceiling(log_2(3n-1)) = 1 + floor(log_2(3n-2)) for n >= 1.
a(n) = A130249(n-1) + 1 = A130253(n-1) for n >= 1.
G.f.: (x/(1-x))*Sum_{k>=0} x^A001045(k).
Showing 1-8 of 8 results.