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 16 results. Next

A239564 a(n) = (round(c^prime(n)) - 1)/prime(n), where c is the pentanacci constant (A103814).

Original entry on oeis.org

154, 504, 5758, 19912, 245714, 11251030, 40679232, 1967728552, 26525975822, 97753187576, 1335948880418, 68398141417510, 3547322151373882, 13260715720748120, 697034813138756392, 9825603574709578482, 36935066391752894480, 1970457739485406707872
Offset: 5

Views

Author

Keywords

Comments

For n>=5, round(c^prime(n)) == 1 (mod 2*prime(n)). Proof in Shevelev link. In particular, all terms are even.

Crossrefs

A138372 Count of post-period decimal digits up to which the rounded n-th convergent to A103814 agrees with the exact value.

Original entry on oeis.org

1, 3, 3, 4, 3, 4, 5, 6, 6, 7, 8, 8, 9, 10, 11, 12, 13, 13, 14, 16, 16, 18, 21, 23, 23, 24, 26, 27, 27, 29, 29, 30, 30, 31, 32, 34, 37, 36, 39, 40, 41, 42, 43, 43, 46, 47, 48, 48, 49, 50, 50, 52, 53, 54, 55, 56, 56, 58, 58, 59, 62, 63, 65, 66, 67, 68, 69, 70, 70, 71, 73, 73, 75
Offset: 1

Views

Author

Artur Jasinski, Mar 17 2008

Keywords

Comments

This is a measure of the quality of the n-th convergent to the Pentanacci constant A103814 if the convergent and the exact value are compared rounded to an increasing number of digits. The sequence of rounded values of A103814 is 2, 2.0, 1.97, 1.966, 1.9659, 1.96595, 1.965948, 1.9659482 etc. The n-th convergents are 2 (n=1), 57/29 (n=2), 116/59 (n=3), 173/88 (n=4), 462/235 (n=5) etc, each with associated rounded decimal expansions.
a(n) is the maximum number of post-period digits of the two expansions if compared at the same level of rounding.

Examples

			For n=5, the 5th convergent is 462/235 = 1.96595744.., with a sequence of rounded representations 2, 2.0, 1.97, 1.966, 1.9660, 1.96596, etc.
Rounded to 1, 2, or 3 post-period decimal digits, this is the same as the rounded version of the exact value, but disagrees if both are rounded to 4 decimal digits, where 1.9659 <> 1.9660.
So a(5) = 3 (digits), the maximum rounding level with agreement.
		

Crossrefs

Extensions

Definition and values replaced as defined via continued fractions - R. J. Mathar, Oct 01 2009

A058265 Decimal expansion of the tribonacci constant t, the real root of x^3 - x^2 - x - 1.

Original entry on oeis.org

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

Views

Author

Robert G. Wilson v, Dec 07 2000

Keywords

Comments

"The tribonacci constant, the only real solution to the equation x^3 - x^2 - x - 1 = 0, which is related to tribonacci sequences (in which U_n = U_n-1 + U_n-2 + U_n-3) as the Golden Ratio is related to the Fibonacci sequence and its generalizations. This ratio also appears when a snub cube is inscribed in an octahedron or a cube, by analogy once again with the appearance of the Golden Ratio when an icosahedron is inscribed in an octahedron. [John Sharp, 1997]"
The tribonacci constant corresponds to the Golden Section in a tripartite division 1 = u_1 + u_2 + u_3 of a unit line segment; i.e., if 1/u_1 = u_1/u_2 = u_2/u_3 = c, c is the tribonacci constant. - Seppo Mustonen, Apr 19 2005
The other two polynomial roots are the complex-conjugated pair -0.4196433776070805662759262... +- i* 0.60629072920719936925934... - R. J. Mathar, Oct 25 2008
For n >= 3, round(q^prime(n)) == 1 (mod 2*prime(n)). Proof in Shevelev link. - Vladimir Shevelev, Mar 21 2014
Concerning orthogonal projections, the tribonacci constant is the ratio of the diagonal of a square to the width of a rhombus projected by rotating a square along its diagonal in 3D until the angle of rotation equals the apparent apex angle at approximately 57.065 degrees (also the corresponding angle in the formula generating A256099). See illustration in the links. - Peter M. Chema, Jan 02 2017
From Wolfdieter Lang, Aug 10 2018: (Start)
Real eigenvalue t of the tribonacci Q-matrix <<1, 1, 1>,<1, 0, 0>,<0, 1, 0>>.
Limit_{n -> oo} T(n+1)/T(n) = t (from the T recurrence), where T = {A000073(n+2)}_{n >= 0}. (End)
The nonnegative powers of t are t^n = T(n)*t^2 + (T(n-1) + T(n-2))*t + T(n-1)*1, for n >= 0, with T(n) = A000073(n), with T(-1) = 1 and T(-2) = -1, This follows from the recurrences derived from t^3 = t^2 + t + 1. See the examples below. For the negative powers see A319200. - Wolfdieter Lang, Oct 23 2018
Note that we have: t + t^(-3) = 2, and the k-nacci constant approaches 2 when k approaches infinity (Martin Gardner). - Bernard Schott, May 16 2022
The roots of this cubic are found from those of y^3 - (4/3)*y - 38/27, after adding 1/3. - Wolfdieter Lang, Aug 24 2022
The algebraic number t - 1 has minimal polynomial x^3 + 2*x^2 - 2 over Q. The roots coincide with those of y^3 - (4/3)*y - 38/27, after subtracting 2/3. - Wolfdieter Lang, Sep 20 2022
The value of the ratio R/r of the radius R of a uniform ball to the radius r of a spherical hole in it with a common point of contact, such that the center of gravity of the object lies on the surface of the spherical hole (Schmidt, 2002). - Amiram Eldar, May 20 2023

Examples

			1.8392867552141611325518525646532866004241787460975922467787586394042032220\
    81966425738435419428307014141979826859240974164178450746507436943831545\
    820499513796249655539644613666121540277972678118941041...
From _Wolfdieter Lang_, Oct 23 2018: (Start)
The coefficients of t^2, t, 1 for t^n begin, for n >= 0:
    n     t^2    t    1
    -------------------
    0      0     0    1
    1      0     1    0
    2      1     0    0
    1      1     1    1
    4      2     2    1
    5      4     3    2
    6      7     6    4
    7     13    11    7
    8     24    20   13
    9     44    37   24
   10     81    68   44
...  (End)
		

References

  • Steven R. Finch, Mathematical Constants, Cambridge, 2003, Section 1.2.2.
  • Martin Gardner, The Second Scientific American Book of Mathematical Puzzles and Diversions, "Phi: The Golden Ratio", Chapter 8, p. 101, Simon & Schuster, NY, 1961.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Revised Edition, Penguin Books, London, England, 1997, page 23.

Crossrefs

Cf. A000073, A019712 (continued fraction), A133400, A254231, A158919 (spectrum = floor(n*t)), A357101 (x^3-2*x^2-2).
Cf. A192918 (reciprocal), A276800 (square), A276801 (cube), A319200.
k-nacci constants: A001622 (Fibonacci), this sequence (tribonacci), A086088 (tetranacci), A103814 (pentanacci), A118427 (hexanacci), A118428 (heptanacci).

Programs

  • Maple
    Digits:=200; fsolve(x^3=x^2+x+1); # N. J. A. Sloane, Mar 16 2019
  • Mathematica
    RealDigits[ N[ 1/3 + 1/3*(19 - 3*Sqrt[33])^(1/3) + 1/3*(19 + 3*Sqrt[33])^(1/3), 100]] [[1]]
    RealDigits[Root[x^3-x^2-x-1,1],10,120][[1]] (* Harvey P. Dale, Mar 23 2019 *)
  • Maxima
    set_display(none)$ fpprec:100$ bfloat(rhs(solve(t^3-t^2-t-1,t)[3])); /* Dimitri Papadopoulos, Nov 09 2023 */
  • PARI
    default(realprecision, 20080); x=solve(x=1, 2, x^3 - x^2 - x - 1); for (n=1, 20000, d=floor(x); x=(x-d)*10; write("b058265.txt", n, " ", d));  \\ Harry J. Smith, May 30 2009
    
  • PARI
    q=(1+sqrtn(19+3*sqrt(33),3)+sqrtn(19-3*sqrt(33),3))/3 \\ Use \p# to set 'realprecision'. - M. F. Hasler, Mar 23 2014
    

Formula

t = (1/3)*(1+(19+3*sqrt(33))^(1/3)+(19-3*sqrt(33))^(1/3)). - Zak Seidov, Jun 08 2005
t = 1 - Sum_{k>=1} A057597(k+2)/(T_k*T_(k+1)), where T_n = A000073(n+1). - Vladimir Shevelev, Mar 02 2013
1/t + 1/t^2 + 1/t^3 = 1/A058265 + 1/A276800 + 1/A276801 = 1. - N. J. A. Sloane, Oct 28 2016
t = (4/3)*cosh((1/3)*arccosh(19/8)) + 1/3. - Wolfdieter Lang, Aug 24 2022
t = 2 - Sum_{k>=0} binomial(4*k + 2, k)/((k + 1)*2^(4*k + 3)). - Antonio Graciá Llorente, Oct 28 2024

A001591 Pentanacci numbers: a(n) = a(n-1) + a(n-2) + a(n-3) + a(n-4) + a(n-5), a(0)=a(1)=a(2)=a(3)=0, a(4)=1.

Original entry on oeis.org

0, 0, 0, 0, 1, 1, 2, 4, 8, 16, 31, 61, 120, 236, 464, 912, 1793, 3525, 6930, 13624, 26784, 52656, 103519, 203513, 400096, 786568, 1546352, 3040048, 5976577, 11749641, 23099186, 45411804, 89277256, 175514464, 345052351, 678355061, 1333610936, 2621810068
Offset: 0

Views

Author

Keywords

Comments

Number of permutations satisfying -k <= p(i) - i <= r, i=1..n-4, with k=1, r=4. - Vladimir Baltic, Jan 17 2005
a(n) is the number of compositions of n-4 with no part greater than 5. - Vladimir Baltic, Jan 17 2005
The pentanomial (A035343(n)) transform of a(n) is a(5n+4), n >= 0. - Bob Selcoe, Jun 10 2014
a(n) is the number of ways to tile a strip of length n-4 with squares, dominoes, trominoes (of length 3), and rectangles with length 4 (tetraminoes) and length 5 (pentaminoes). - Wajdi Maaloul, Jun 21 2022

Examples

			n=2: a(14) = (1*1 + 2*1 + 3*2 + 4*4 + 5*8 + 4*16 + 3*31 + 2*61 + 1*120) = 464. - _Bob Selcoe_, Jun 10 2014
G.f. = x^4 + x^5 + 2*x^6 + 4*x^7 + 8*x^8 + 16*x^9 + 31*x^10 + 120*x^11 + ...
		

References

  • Silvia Heubach and Toufik Mansour, Combinatorics of Compositions and Words, CRC Press, 2010.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Row 5 of arrays A048887 and A092921 (k-generalized Fibonacci numbers).
Cf. A106303 (Pisano period lengths).
Cf. A035343 (pentanomial coefficients).

Programs

  • Magma
    a:=[0,0,0,0,1]; [n le 5 select a[n] else Self(n-1) + Self(n-2) + Self(n-3) + Self(n-4) + Self(n-5): n in [1..40]]; // Marius A. Burtea, Oct 03 2019
    
  • Maple
    g:=1/(1-z-z^2-z^3-z^4-z^5): gser:=series(g, z=0, 49): seq((coeff(gser, z, n)), n=-4..32); # Zerinvary Lajos, Apr 17 2009
    # second Maple program:
    a:= n-> (<<0|1|0|0|0>, <0|0|1|0|0>, <0|0|0|1|0>, <0|0|0|0|1>, <1|1|1|1|1>>^n)[1, 5]:
    seq(a(n), n=0..44);  # Alois P. Heinz, Apr 09 2021
  • Mathematica
    CoefficientList[Series[x^4/(1 - x - x^2 - x^3 - x^4 - x^5), {x, 0, 50}], x]
    a[0] = a[1] = a[2] = a[3] = 0; a[4] = a[5] = 1; a[n_] := a[n] = 2 a[n - 1] - a[n - 6]; Array[a, 37, 0]
    LinearRecurrence[{1, 1, 1, 1, 1}, {0, 0, 0, 0, 1}, 50] (* Vladimir Joseph Stephan Orlovsky, May 25 2011 *)
  • Maxima
    a(n):=mod(floor(10^((n-4)*(n+1))*10^(5*(n+1))*(10^(n+1)-1)/(10^(6*(n+1))-2*10^(5*(n+1))+1)),10^n); /* Tani Akinari, Apr 10 2014 */
    
  • PARI
    a=vector(100);a[4]=a[5]=1;for(n=6,#a,a[n]=a[n-1]+a[n-2]+a[n-3]+a[n-4]+a[n-5]);concat(0, a) \\ Charles R Greathouse IV, Jul 15 2011
    
  • PARI
    A001591(n,m=5)=(matrix(m,m,i,j,i==j-1||i==m)^n)[1,m] \\ M. F. Hasler, Apr 20 2018
    
  • PARI
    a(n)= {my(x='x, p=polrecip(1 - x - x^2 - x^3 - x^4 - x^5)); polcoef(lift(Mod(x, p)^n), 4); }
    vector(41, n, a(n-1)) \\ Joerg Arndt, May 16 2021
    
  • Python
    def pentanacci():
        a, b, c, d, e = 0, 0, 0, 0, 1
        while True:
            yield a
            a, b, c, d, e = b, c, d, e, a + b + c + d + e
    f = pentanacci()
    print([next(f) for  in range(100)]) # _Reza K Ghazi Apr 09 2021

Formula

G.f.: x^4/(1 - x - x^2 - x^3 - x^4 - x^5). - Simon Plouffe in his 1992 dissertation.
G.f.: Sum_{n >= 0} x^(n+4) * (Product_{k = 1..n} (k + k*x + k*x^2 + k*x^3 + x^4)/(1 + k*x + k*x^2 + k*x^3 + k*x^4)). - Peter Bala, Jan 04 2015
Another form of the g.f.: f(z) = (z^4-z^5)/(1-2*z+z^6); then a(n) = Sum_{i=0..floor((n-4)/6)} ((-1)^i*binomial(n-4-5*i,i)*2^(n-4-6*i)) - Sum_{i=0..floor((n-5)/6)} ((-1)^i*binomial(n-5-5*i,i)*2^(n-5-6*i)) with convention Sum_{i=m..n} alpha(i) = 0 for m > n. - Richard Choulet, Feb 22 2010
a(n) = Sum_{k=1..n} (Sum_{r=0..k} (binomial(k,r) * Sum_{m=0..r} (binomial(r,m) * Sum_{j=0..m} (binomial(m,j)*binomial(j,n-m-k-j-r))))), n > 0. - Vladimir Kruchinin, Aug 30 2010
Sum_{k=0..4*n} a(k+b)*A035343(n,k) = a(5*n+b), b >= 0.
a(n) = 2*a(n-1) - a(n-6). - Vincenzo Librandi, Dec 19 2010
a(n) = (Sum_{i=0..n-1} a(i)*A074048(n-i))/(n-4) for n > 4. - Greg Dresden and Advika Srivastava, Oct 01 2019
For k>0 and n>0, a(n+5*k) = A074048(k)*a(n+4*k) - A123127(k-1)*a(n+3*k) + A123126(k-1)*a(n+2*k) - A074062(k)*a(n+k) + a(n). - Kai Wang, Sep 06 2020
lim n->oo a(n)/a(n-1) = A103814. - R. J. Mathar, Mar 11 2024

A086088 Decimal expansion of the limit of the ratio of consecutive terms in the tetranacci sequence A000078.

Original entry on oeis.org

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

Views

Author

Eric W. Weisstein, Jul 08 2003

Keywords

Comments

The tetranacci constant corresponds to the Golden Section in a quadripartite division 1 = u_1 + u_2 + u_3 + u_4 of a unit line segment, i.e., if 1/u_1 = u_1/u_2 = u_2/u_3 = u_3/u_4 = c, c is the tetranacci constant. - Seppo Mustonen, Apr 19 2005
The other 3 polynomial roots of 1+x+x^2+x^3-x^4 are -0.77480411321543385... and the complex-conjugated pair -0.07637893113374572508475 +- i * 0.814703647170386526841... - R. J. Mathar, Oct 25 2008
The continued fraction expansion starts 1, 1, 12, 1, 4, 7, 1, 21, 1, 2, 1, 4, 6, 1, 10, 1, 2, 2, 1, 7, 1, 1,... - R. J. Mathar, Mar 09 2012
For n>=4, round(c^prime(n)) == 1 (mod 2*prime(n)). Proof in Shevelev link. - Vladimir Shevelev, Mar 21 2014
Note that we have: c + c^(-4) = 2, and the k-nacci constant approaches 2 when k approaches infinity (Martin Gardner). - Bernard Schott, May 09 2022

Examples

			1.927561975...
		

References

  • Martin Gardner, The Second Scientific American Book Of Mathematical Puzzles and Diversions, "Phi: The Golden Ratio", Chapter 8, p. 101, Simon & Schuster, NY, 1961.

Crossrefs

Cf. A000078.
k-nacci constants: A001622 (Fibonacci), A058265 (tribonacci), this sequence (tetranacci), A103814 (pentanacci), A118427 (hexanacci), A118428 (heptanacci).

Programs

Formula

Equals 1/4 + sqrt(11/48 - s/72 + 7/s) + sqrt(11/24 + s/72 - 7/s + 1 / sqrt(704/507 - 128 * s/1521 + 7168 / (169 * s))) where s = (sqrt(177304464) + 7020)^(1/3). - Michal Paulovic, Oct 08 2022

A125899 Floor((Pentanacci ratio)^n).

Original entry on oeis.org

1, 3, 7, 14, 29, 57, 113, 223, 438, 862, 1695, 3333, 6552, 12882, 25326, 49791, 97887, 192440, 378328, 743774, 1462223, 2874655, 5651422, 11110404, 21842481, 42941187, 84420150, 165965646, 326279870, 641449337, 1261056193, 2479171198
Offset: 1

Views

Author

Artur Jasinski, Dec 13 2006

Keywords

Examples

			The Pentanacci ratio is A103814, a root of the polynomial x^5-x^4-x^3-x^2-x-1.
		

Crossrefs

Programs

  • Mathematica
    Table[Floor[(Root[ -1 - #1 - #1^2 - #1^3 - #1^4 + #1^5 &, 1])^n], {n, 1, 100}] (* Artur Jasinski *)

Extensions

Reference to A103814 added - R. J. Mathar, Sep 26 2010

A118427 Decimal expansion of hexanacci constant.

Original entry on oeis.org

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

Views

Author

Eric W. Weisstein, Apr 27 2006

Keywords

Comments

The continued fraction expansion starts 1, 1, 59, 1, 10, 2, 1, 6, 2, 1, 6, 1, 1, 7, 1, 71, 7, 1, 6, 8, ... - R. J. Mathar, Mar 09 2012
For n>=7, round(c^prime(n)) == 1 (mod 2*prime(n)). Proof in Shevelev link. - Vladimir Shevelev, Mar 21 2014
Note that we have: c + c^(-6) = 2, and the k-nacci constant approaches 2 when k approaches infinity (Martin Gardner). - Bernard Schott, May 06 2022

Examples

			1.9835828434243263303...
		

References

  • Martin Gardner, The Second Scientific American Book Of Mathematical Puzzles and Diversions, "Phi: The Golden Ratio", Chapter 8, Simon & Schuster, NY, 1961.

Crossrefs

Cf. A001592.
k-nacci constants: A001622 (Fibonacci), A058265 (tribonacci), A086088 (tetranacci), A103814 (pentanacci), this sequence (hexanacci), A118428 (heptanacci).

Programs

  • Mathematica
    RealDigits[ Root[ x^6 - x^5 - x^4 - x^3 - x^2 - x - 1, 2] , 10, 105] // First (* Jean-François Alcover, Feb 07 2013 *)

A184835 a(n) = n + floor(n/t) + floor(n/t^2) + floor(n/t^3) + floor(n/t^4), where t is the pentanacci constant.

Original entry on oeis.org

1, 3, 4, 7, 8, 10, 11, 15, 16, 18, 19, 22, 23, 25, 27, 31, 32, 34, 35, 38, 39, 41, 43, 46, 47, 49, 50, 53, 54, 57, 60, 62, 63, 65, 67, 69, 70, 73, 75, 77, 78, 80, 82, 84, 86, 89, 91, 93, 94, 96, 98, 100, 101, 104, 106, 108, 109, 112, 114, 116, 119, 121, 123, 124, 126, 128, 130, 131, 134, 136, 138, 139, 141, 143, 146, 148, 150, 152, 154, 155, 157, 159, 161, 163, 165, 167, 169, 170, 173, 175, 177, 179, 182, 183, 185, 186, 189, 190, 193, 194, 197, 198, 200, 201, 205, 206, 209, 210, 213, 214, 216, 217, 220, 222, 224, 227, 228
Offset: 1

Views

Author

Paul D. Hanna, Jan 23 2011

Keywords

Comments

This is one of five sequences that partition the positive integers.
Given t is the pentanacci constant, then the following sequences are disjoint:
. A184835(n) = n + [n/t] + [n/t^2] + [n/t^3] + [n/t^4],
. A184836(n) = n + [n*t] + [n/t] + [n/t^2] + [n/t^3],
. A184837(n) = n + [n*t] + [n*t^2] + [n/t] + [n/t^2],
. A184838(n) = n + [n*t] + [n*t^2] + [n*t^3] + [n/t],
. A184839(n) = n + [n*t] + [n*t^2] + [n*t^3] + [n*t^4], where []=floor.
This is a special case of Clark Kimberling's results given in A184812.

Examples

			Given t = pentanacci constant, then t = 1 + 1/t + 1/t^2 + 1/t^3 + 1/t^4,
t = 1.965948236645..., t^2 = 3.864952469169..., t^3 = 7.598296491482..., t^4 = 14.93785758893..., t^5 = 29.36705478623...
		

Crossrefs

Programs

  • Mathematica
    With[{pc=x/.FindRoot[x^5-x^4-x^3-x^2-x-1==0,{x,1.96},WorkingPrecision-> 100]}, Table[n+Total[Table[Floor[n/pc^i],{i,4}]],{n,150}]] (* Harvey P. Dale, Jun 21 2011 *)
  • PARI
    {a(n)=local(t=real(polroots(1+x+x^2+x^3+x^4-x^5)[1])); n+floor(n/t)+floor(n/t^2)+floor(n/t^3)+floor(n/t^4)}

Formula

Limit a(n)/n = t = 1.9659482366454853371899373...
a(n) = n + floor(n*p/u) + floor(n*q/u) + floor(n*r/u) + floor(n*s/u), where p=t, q=t^2, r=t^3, s=t^4, u=t^5, and t is the pentanacci constant.

A118428 Decimal expansion of heptanacci constant.

Original entry on oeis.org

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

Views

Author

Eric W. Weisstein, Apr 27 2006

Keywords

Comments

Other roots of the equation x^7 - x^6 - ... - x - 1 see in A239566. For n>=7, round(c^prime(n)) == 1 (mod 2*prime(n)). Proof in Shevelev link. - Vladimir Shevelev, Mar 21 2014
Note that we have: c + c^(-7) = 2, and the k-nacci constant approaches 2 when k approaches infinity (Martin Gardner). - Bernard Schott, May 07 2022

Examples

			1.9919641966050350210...
		

References

  • Martin Gardner, The Second Scientific American Book Of Mathematical Puzzles and Diversions, "Phi: The Golden Ratio", Chapter 8, p. 101, Simon & Schuster, NY, 1961.

Crossrefs

k-nacci constants: A001622 (Fibonacci), A058265 (tribonacci), A086088 (tetranacci), A103814 (pentanacci), A118427 (hexanacci), this sequence (heptanacci).

Programs

  • Mathematica
    RealDigits[x/.FindRoot[x^7+Total[-x^Range[0,6]]==0,{x,2}, WorkingPrecision-> 110]][[1]] (* Harvey P. Dale, Dec 13 2011 *)
  • PARI
    polrootsreal(x^7 - x^6 - x^5 - x^4 - x^3 - x^2 - x - 1)[1] \\ Charles R Greathouse IV, Feb 11 2025

A184836 a(n) = n + floor(n*t) + floor(n/t) + floor(n/t^2) + floor(n/t^3), where t is the pentanacci constant.

Original entry on oeis.org

2, 6, 9, 14, 17, 21, 24, 30, 33, 37, 40, 45, 48, 52, 55, 61, 64, 68, 71, 76, 79, 83, 87, 92, 95, 99, 102, 107, 110, 113, 118, 122, 125, 129, 133, 137, 140, 145, 149, 153, 156, 160, 164, 168, 171, 176, 180, 184, 187, 191, 195, 199, 202, 207, 211, 215, 218, 223, 226, 229, 234, 238, 242, 245, 249, 253, 257, 260, 265, 269, 273, 276, 280, 284, 288, 292, 296, 300, 304, 307, 311, 315, 319, 323, 327, 331, 335, 338, 342, 345, 349, 353, 358, 361, 365, 368, 373, 376, 381, 384, 389, 392, 396, 399, 404, 407, 412, 415, 420, 423
Offset: 1

Views

Author

Paul D. Hanna, Jan 23 2011

Keywords

Comments

This is one of five sequences that partition the positive integers.
Given t is the pentanacci constant, then the following sequences are disjoint:
. A184835(n) = n + [n/t] + [n/t^2] + [n/t^3] + [n/t^4],
. A184836(n) = n + [n*t] + [n/t] + [n/t^2] + [n/t^3],
. A184837(n) = n + [n*t] + [n*t^2] + [n/t] + [n/t^2],
. A184838(n) = n + [n*t] + [n*t^2] + [n*t^3] + [n/t],
. A184839(n) = n + [n*t] + [n*t^2] + [n*t^3] + [n*t^4], where []=floor.
This is a special case of Clark Kimberling's results given in A184812.

Examples

			Given t = pentanacci constant, then t^2 = 1 + t + 1/t + 1/t^2 + 1/t^3,
t = 1.965948236645..., t^2 = 3.864952469169..., t^3 = 7.598296491482..., t^4 = 14.93785758893..., t^5 = 29.36705478623...
		

Crossrefs

Programs

  • PARI
    {a(n)=local(t=real(polroots(1+x+x^2+x^3+x^4-x^5)[1])); n+floor(n*t)+floor(n/t)+floor(n/t^2)+floor(n/t^3)}

Formula

Limit a(n)/n = t^2 = 3.8649524691694932164414964...
a(n) = n + floor(n*p/s) + floor(n*q/s) + floor(n*r/s) + floor(n*u/s), where p=t, q=t^2, r=t^3, s=t^4, u=t^5, and t is the pentanacci constant.
Showing 1-10 of 16 results. Next