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-10 of 19 results. Next

A130234 Minimal index k of a Fibonacci number such that Fibonacci(k) >= n (the 'upper' Fibonacci Inverse).

Original entry on oeis.org

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

Views

Author

Hieronymus Fischer, May 17 2007

Keywords

Comments

Inverse of the Fibonacci sequence (A000045), nearly, since a(Fibonacci(n)) = n except for n = 2 (see A130233 for another version). a(n+1) is equal to the partial sum of the Fibonacci indicator sequence (see A104162).

Examples

			a(10) = 7, since Fibonacci(7) = 13 >= 10 but Fibonacci(6) = 8 < 10.
		

Crossrefs

Partial sums: A130236.
Other related sequences: A000045, A130233, A130256, A130260, A104162, A108852.
Lucas inverse: A130241 - A130248.

Programs

  • Maple
    A130234 := proc(n)
        local i;
        for i from 0 do
            if A000045(i) >= n then
                return i;
            end if;
        end do:
    end proc: # R. J. Mathar, Jan 31 2015
  • Mathematica
    a[n_] := For[i = 0, True, i++, If[Fibonacci[i] >= n, Return[i]]];
    a /@ Range[0, 80] (* Jean-François Alcover, Apr 13 2020 *)
  • PARI
    a(n)=my(k);while(fibonacci(k)Charles R Greathouse IV, Feb 03 2014, corrected by M. F. Hasler, Apr 07 2021

Formula

a(n) = ceiling(log_phi((sqrt(5)*n + sqrt(5*n^2-4))/2)) = ceiling(arccosh(sqrt(5)*n/2)/log(phi)) where phi = (1+sqrt(5))/2, the golden ratio, for n > 0.
a(n) = A130233(n-1) + 1 for n > 0.
G.f.: x/(1-x) * Sum_{k >= 0} x^Fibonacci(k).
a(n) = ceiling(log_phi(sqrt(5)*n - 1)) for n > 0, where phi is the golden ratio. - Hieronymus Fischer, Jul 02 2007
a(n) = A108852(n-1). - R. J. Mathar, Jan 31 2015

A130235 Partial sums of the 'lower' Fibonacci Inverse A130233.

Original entry on oeis.org

0, 2, 5, 9, 13, 18, 23, 28, 34, 40, 46, 52, 58, 65, 72, 79, 86, 93, 100, 107, 114, 122, 130, 138, 146, 154, 162, 170, 178, 186, 194, 202, 210, 218, 227, 236, 245, 254, 263, 272, 281, 290, 299, 308, 317, 326, 335, 344, 353, 362, 371, 380, 389, 398, 407, 417, 427
Offset: 0

Views

Author

Hieronymus Fischer, May 17 2007

Keywords

Crossrefs

Programs

  • Magma
    m:=120;
    f:= func< x | (&+[x^Fibonacci(j): j in [1..Floor(3*Log(3*m+1))]])/(1-x)^2 >;
    R:=PowerSeriesRing(Rationals(), m+1);
    [0] cat Coefficients(R!( f(x) )); // G. C. Greubel, Mar 17 2023
    
  • Mathematica
    nmax = 90; CoefficientList[Series[Sum[x^Fibonacci[k], {k, 1, 1 + Log[3/2 + Sqrt[5]*nmax]/Log[GoldenRatio]}]/(1-x)^2, {x, 0, nmax}], x] (* Vaclav Kotesovec, Apr 14 2020 *)
  • SageMath
    m=120
    def f(x): return sum( x^fibonacci(j) for j in range(1, int(3*log(3*m+1))))/(1-x)^2
    def A130235_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( f(x) ).list()
    A130235_list(m) # G. C. Greubel, Mar 17 2023

Formula

a(n) = Sum_{k=0..n} A130233(k) = (n+1)*A130233(n) - Fib(A130233(n)+2) + 1.
G.f.: 1/(1-x)^2 * Sum_{k>=1} x^Fib(k). [corrected by Joerg Arndt, Apr 14 2020]

A130240 Partial sums of A130239.

Original entry on oeis.org

0, 2, 4, 6, 9, 12, 15, 18, 21, 25, 29, 33, 37, 41, 45, 49, 53, 57, 61, 65, 69, 73, 77, 81, 85, 90, 95, 100, 105, 110, 115, 120, 125, 130, 135, 140, 145, 150, 155, 160, 165, 170, 175, 180, 185, 190, 195, 200, 205, 210, 215, 220, 225, 230, 235, 240, 245, 250, 255, 260
Offset: 0

Views

Author

Hieronymus Fischer, May 17 2007

Keywords

Crossrefs

Programs

Formula

a(n) = Sum_{k=0..n} A130239(k).
a(n) = (n+1)*A130233(sqrt(n)) - Fib(A130233(sqrt(n)) + 1) * Fib(A130232(sqrt(n))).
G.f.: (1/(1-x)^2) * Sum_{k>=1} x^(Fib(k)^2).

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).

A130255 Maximal index k of an odd Fibonacci number (A001519) such that A001519(k) = Fibonacci(2k-1) <= n (the 'lower' odd Fibonacci Inverse).

Original entry on oeis.org

1, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6
Offset: 1

Views

Author

Hieronymus Fischer, May 24 2007, Jul 02 2007

Keywords

Comments

Inverse of the odd Fibonacci sequence (A001519), nearly, since a(A001519(n))=n except for n=0 (see A130256 for another version). a(n)+1 is the number of odd Fibonacci numbers (A001519) <= n (for n >= 1).

Examples

			a(10)=3 because A001519(3) = 5 <= 10, but A001519(4) = 13 > 10.
		

Crossrefs

Cf. partial sums A130257. Other related sequences: A000045, A130233, A130237, A130239, A130256, A130259, A104160. Lucas inverse: A130241 - A130248.

Programs

  • Magma
    phi:=(1+Sqrt(5))/2; [Floor((1 +Argsinh(Sqrt(5)*n/2)/Log(phi))/2): n in [1..100]]; // G. C. Greubel, Sep 09 2018
  • Mathematica
    Table[Floor[(1 +ArcSinh[Sqrt[5]*n/2]/Log[GoldenRatio])/2], {n, 1, 100}] (* G. C. Greubel, Sep 09 2018 *)
  • PARI
    phi=(1+sqrt(5))/2; vector(100, n, floor((1 +asinh(sqrt(5)*n/2)/log(phi))/2)) \\ G. C. Greubel, Sep 09 2018
    

Formula

a(n) = floor((1 + arcsinh(sqrt(5)*n/2)/log(phi))/2).
a(n) = floor((1 + arccosh(sqrt(5)*n/2)/log(phi))/2).
a(n) = floor((1 + log_phi(sqrt(5)*n))/2) for n >= 1, where phi = (1 + sqrt(5))/2.
G.f.: g(x) = 1/(1-x)*Sum_{k>=1} x^Fibonacci(2k-1).
a(n) = floor((1/2)*(1 + log_phi(sqrt(5)*n + 1))) for n >= 1.

A130257 Partial sums of the 'lower' odd Fibonacci Inverse A130255.

Original entry on oeis.org

1, 3, 5, 7, 10, 13, 16, 19, 22, 25, 28, 31, 35, 39, 43, 47, 51, 55, 59, 63, 67, 71, 75, 79, 83, 87, 91, 95, 99, 103, 107, 111, 115, 120, 125, 130, 135, 140, 145, 150, 155, 160, 165, 170, 175, 180, 185, 190, 195, 200, 205, 210, 215, 220, 225, 230, 235, 240, 245, 250
Offset: 1

Views

Author

Hieronymus Fischer, May 24 2007

Keywords

Crossrefs

Programs

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

Formula

a(n) = (n+1)*A130255(n) - A001906(A130255(n)).
a(n) = (n+1)*A130255(n) - Fib(2*A130255(n)).
G.f.: g(x)=1/(1-x)^2*sum(k>=1, x^Fib(2k-1)).

A130243 Partial sums of the 'lower' Lucas Inverse A130241.

Original entry on oeis.org

1, 2, 4, 7, 10, 13, 17, 21, 25, 29, 34, 39, 44, 49, 54, 59, 64, 70, 76, 82, 88, 94, 100, 106, 112, 118, 124, 130, 137, 144, 151, 158, 165, 172, 179, 186, 193, 200, 207, 214, 221, 228, 235, 242, 249, 256, 264, 272, 280, 288, 296, 304, 312, 320, 328, 336, 344, 352
Offset: 1

Views

Author

Hieronymus Fischer, May 19 2007

Keywords

Crossrefs

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

Programs

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

Formula

a(n) = Sum_{k=1..n} A130241(k).
a(n) = (n+1)*A130241(n) - A000032(A130241(n)+2) + 3.
G.f.: g(x) = 1/(1-x)^2*Sum_{k>=1} x^Lucas(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).

A130259 Maximal index k of an even Fibonacci number (A001906) such that A001906(k) = Fib(2k) <= n (the 'lower' even Fibonacci Inverse).

Original entry on oeis.org

0, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5
Offset: 0

Views

Author

Hieronymus Fischer, May 25 2007, Jul 02 2007

Keywords

Comments

Inverse of the even Fibonacci sequence (A001906), since a(A001906(n))=n (see A130260 for another version).
a(n)+1 is the number of even Fibonacci numbers (A001906) <=n.

Examples

			a(10)=3 because A001906(3)=8<=10, but A001906(4)=21>10.
		

Crossrefs

Cf. partial sums A130261. Other related sequences: A000045, A001519, A130233, A130237, A130239, A130255, A130260, A104160. Lucas inverse: A130241 - A130248.

Programs

  • Magma
    [Floor(Log((Sqrt(5)*n+1))/(2*Log((1+Sqrt(5))/2))): n in [0..100]]; // G. C. Greubel, Sep 12 2018
  • Mathematica
    Table[Floor[1/2*Log[GoldenRatio, (Sqrt[5]*n + 1)]], {n, 0, 100}] (* G. C. Greubel, Sep 12 2018 *)
  • PARI
    vector(100, n, n--; floor(log((sqrt(5)*n+1))/(2*log((1+sqrt(5))/2) ))) \\ G. C. Greubel, Sep 12 2018
    

Formula

a(n) = floor(arcsinh(sqrt(5)*n/2)/(2*log(phi))), where phi=(1+sqrt(5))/2.
a(n) = A130260(n+1) - 1.
G.f.: g(x) = 1/(1-x)*Sum_{k>=1} x^Fibonacci(2*k).
a(n) = floor(1/2*log_phi(sqrt(5)*n+1)) for n>=0.

A130261 Partial sums of the 'lower' even Fibonacci Inverse A130259.

Original entry on oeis.org

0, 1, 2, 4, 6, 8, 10, 12, 15, 18, 21, 24, 27, 30, 33, 36, 39, 42, 45, 48, 51, 55, 59, 63, 67, 71, 75, 79, 83, 87, 91, 95, 99, 103, 107, 111, 115, 119, 123, 127, 131, 135, 139, 143, 147, 151, 155, 159, 163, 167, 171, 175, 179, 183, 187, 192, 197, 202, 207, 212, 217
Offset: 0

Views

Author

Hieronymus Fischer, May 25 2007

Keywords

Crossrefs

Programs

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

Formula

a(n) = (n+1)*A130259(n) - A001519(A130259(n)+1) + 1.
a(n) = (n+1)*A130259(n) - Fib(2*A130259(n)+1) + 1.
G.f.: g(x) = 1/(1-x)^2*Sum_{k>=1} x^Fib(2*k).
Showing 1-10 of 19 results. Next