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

A062308 Largest prime factor of 11^n+1 (A034524).

Original entry on oeis.org

2, 3, 61, 37, 7321, 13421, 1117, 1623931, 6304673, 590077, 212601841, 58367, 20113, 59583967, 55527473, 7537711, 447600088289, 2218331, 3138426605161, 1390636259, 1120648576818041, 3421169496361, 2649263870814793
Offset: 0

Views

Author

Jason Earls, Jul 12 2001

Keywords

Crossrefs

Cf. similar sequences listed in A274903.

Programs

  • Magma
    [Maximum(PrimeDivisors(11^n+1)): n in [0..40]]; // Vincenzo Librandi, Jul 12 2016
  • Mathematica
    Table[FactorInteger[11^n + 1][[-1, 1]], {n, 0, 20}] (* Vincenzo Librandi, Jul 12 2016 *)
  • PARI
    for(n=0,22,print(factor(11^n+1)))
    
  • PARI
    { for (n=0, 80, f=factor(11^n + 1)~; write("b062308.txt", n, " ", f[1, length(f)]) ) } \\ Harry J. Smith, Aug 04 2009
    

Formula

a(n) = A006530(A034524(n)). - Vincenzo Librandi, Jul 12 2016

Extensions

Terms to a(80) in b-file from Harry J. Smith, Jun 01 2010
a(81)-a(301) in b-file from Amiram Eldar, Feb 08 2020
a(302)-a(325) in b-file from Max Alekseyev, Apr 25 2022, Oct 11 2023

A000051 a(n) = 2^n + 1.

Original entry on oeis.org

2, 3, 5, 9, 17, 33, 65, 129, 257, 513, 1025, 2049, 4097, 8193, 16385, 32769, 65537, 131073, 262145, 524289, 1048577, 2097153, 4194305, 8388609, 16777217, 33554433, 67108865, 134217729, 268435457, 536870913, 1073741825, 2147483649, 4294967297, 8589934593
Offset: 0

Views

Author

Keywords

Comments

Same as Pisot sequence L(2,3).
Length of the continued fraction for Sum_{k=0..n} 1/3^(2^k). - Benoit Cloitre, Nov 12 2003
See also A004119 for a(n) = 2a(n-1)-1 with first term = 1. - Philippe Deléham, Feb 20 2004
From the second term on (n>=1), in base 2, these numbers present the pattern 1000...0001 (with n-1 zeros), which is the "opposite" of the binary 2^n-2: (0)111...1110 (cf. A000918). - Alexandre Wajnberg, May 31 2005
Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=5, (i>1), A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n>=1, a(n-1)=(-1)^(n-1)* charpoly(A,3). - Milan Janjic, Jan 27 2010
First differences of A006127. - Reinhard Zumkeller, Apr 14 2011
The odd prime numbers in this sequence form A019434, the Fermat primes. - David W. Wilson, Nov 16 2011
Pisano period lengths: 1, 1, 2, 1, 4, 2, 3, 1, 6, 4, 10, 2, 12, 3, 4, 1, 8, 6, 18, 4, ... . - R. J. Mathar, Aug 10 2012
Is the mentioned Pisano period lengths (see above) the same as A007733? - Omar E. Pol, Aug 10 2012
Only positive integers that are not 1 mod (2k+1) for any k>1. - Jon Perry, Oct 16 2012
For n >= 1, a(n) is the total length of the segments of the Hilbert curve after n iterations. - Kival Ngaokrajang, Mar 30 2014
Frénicle de Bessy (1657) proved that a(3) = 9 is the only square in this sequence. - Charles R Greathouse IV, May 13 2014
a(n) is the number of distinct possible sums made with at most two elements in {1,...,a(n-1)} for n > 0. - Derek Orr, Dec 13 2014
For n > 0, given any set of a(n) lattice points in R^n, there exist 2 distinct members in this set whose midpoint is also a lattice point. - Melvin Peralta, Jan 28 2017
Also the number of independent vertex sets, irredundant sets, and vertex covers in the (n+1)-star graph. - Eric W. Weisstein, Aug 04 and Sep 21 2017
Also the number of maximum matchings in the 2(n-1)-crossed prism graph. - Eric W. Weisstein, Dec 31 2017
Conjecture: For any integer n >= 0, a(n) is the permanent of the (n+1) X (n+1) matrix with M(j, k) = -floor((j - k - 1)/(n + 1)). This conjecture is inspired by the conjecture of Zhi-Wei Sun in A036968. - Peter Luschny, Sep 07 2021

References

  • Paul Bachmann, Niedere Zahlentheorie (1902, 1910), reprinted Chelsea, NY, 1968, vol. 2, p. 75.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See pp. 46, 60, 244.
  • 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).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 141.

Crossrefs

Apart from the initial 1, identical to A094373.
See A008776 for definitions of Pisot sequences.
Column 2 of array A103438.
Cf. A007583 (a((n-1)/2)/3 for odd n).

Programs

  • Haskell
    a000051 = (+ 1) . a000079
    a000051_list = iterate ((subtract 1) . (* 2)) 2
    -- Reinhard Zumkeller, May 03 2012
    
  • Magma
    [2^n+1: n in [0..40]]; // G. C. Greubel, Jan 18 2025
  • Maple
    A000051:=-(-2+3*z)/(2*z-1)/(z-1); # Simon Plouffe in his 1992 dissertation
    a := n -> add(binomial(n,k)*bernoulli(n-k,1)*2^(k+1)/(k+1),k=0..n); # Peter Luschny, Apr 20 2009
  • Mathematica
    Table[2^n + 1, {n,0,40}]
    2^Range[0,40] + 1 (* Eric W. Weisstein, Jul 17 2017 *)
    LinearRecurrence[{3, -2}, {2, 3}, 40] (* Eric W. Weisstein, Sep 21 2017 *)
  • PARI
    a(n)=2^n+1
    
  • PARI
    first(n) = Vec((2 - 3*x)/((1 - x)*(1 - 2*x)) + O(x^n)) \\ Iain Fox, Dec 31 2017
    
  • Python
    def A000051(n): return (1<Chai Wah Wu, Dec 21 2022
    

Formula

a(n) = 2*a(n-1) - 1 = 3*a(n-1) - 2*a(n-2).
G.f.: (2-3*x)/((1-x)*(1-2*x)).
First differences of A052944. - Emeric Deutsch, Mar 04 2004
a(0) = 1, then a(n) = (Sum_{i=0..n-1} a(i)) - (n-2). - Gerald McGarvey, Jul 10 2004
Inverse binomial transform of A007689. Also, V sequence in Lucas sequence L(3, 2). - Ross La Haye, Feb 07 2005
a(n) = A127904(n+1) for n>0. - Reinhard Zumkeller, Feb 05 2007
Equals binomial transform of [2, 1, 1, 1, ...]. - Gary W. Adamson, Apr 23 2008
a(n) = A000079(n)+1. - Omar E. Pol, May 18 2008
E.g.f.: exp(x) + exp(2*x). - Mohammad K. Azarian, Jan 02 2009
a(n) = A024036(n)/A000225(n). - Reinhard Zumkeller, Feb 14 2009
From Peter Luschny, Apr 20 2009: (Start)
A weighted binomial sum of the Bernoulli numbers A027641/A027642 with A027641(1)=1 (which amounts to the definition B_{n} = B_{n}(1)).
a(n) = Sum_{k=0..n} C(n,k)*B_{n-k}*2^(k+1)/(k+1). (See also A052584.) (End)
a(n) is the a(n-1)-th odd number for n >= 1. - Jaroslav Krizek, Apr 25 2009
From Reinhard Zumkeller, Feb 28 2010: (Start)
a(n)*A000225(n) = A000225(2*n).
a(n) = A173786(n,0). (End)
If p[i]=Fibonacci(i-4) and if A is the Hessenberg matrix of order n defined by: A[i,j]=p[j-i+1], (i<=j), A[i,j]=-1, (i=j+1), and A[i,j]=0 otherwise, then, for n>=1, a(n-1)= det A. - Milan Janjic, May 08 2010
a(n+2) = a(n) + a(n+1) + A000225(n). - Ivan N. Ianakiev, Jun 24 2012
a(A006521(n)) mod A006521(n) = 0. - Reinhard Zumkeller, Jul 17 2014
a(n) = 3*A007583((n-1)/2) for n odd. - Eric W. Weisstein, Jul 17 2017
Sum_{n>=0} 1/a(n) = A323482. - Amiram Eldar, Nov 11 2020

A034472 a(n) = 3^n + 1.

Original entry on oeis.org

2, 4, 10, 28, 82, 244, 730, 2188, 6562, 19684, 59050, 177148, 531442, 1594324, 4782970, 14348908, 43046722, 129140164, 387420490, 1162261468, 3486784402, 10460353204, 31381059610, 94143178828, 282429536482, 847288609444, 2541865828330, 7625597484988
Offset: 0

Views

Author

Keywords

Comments

Companion numbers to A003462.
a(n) = A024101(n)/A024023(n). - Reinhard Zumkeller, Feb 14 2009
Mahler exhibits this sequence with n>=2 as a proof that there exists an infinite number of x coprime to 3, such that x belongs to A005836 and x^2 belong to A125293. - Michel Marcus, Nov 12 2012
a(n-1) is the number of n-digit base 3 numbers that have an even number of digits 0. - Yifan Xie, Jul 13 2024

Examples

			a(3)=28 because 4*a(2)-3*a(1)=4*10-3*4=28 (28 is also 3^3 + 1).
G.f. = 2 + 4*x + 10*x^2 + 28*x^3 + 82*x^4 + 244*x^5 + 730*x^5 + ...
		

References

  • Knuth, Donald E., Satisfiability, Fascicle 6, volume 4 of The Art of Computer Programming. Addison-Wesley, 2015, pages 148 and 220, Problem 191.
  • P. Ribenboim, The Little Book of Big Primes, Springer-Verlag, NY, 1991, pp. 35-36, 53.

Crossrefs

Programs

  • Magma
    [3^n+1: n in [0..30]]; // Vincenzo Librandi, Jan 11 2017
  • Maple
    ZL:= [S, {S=Union(Sequence(Z), Sequence(Union(Z, Z, Z)))}, unlabeled]: seq(combstruct[count](ZL, size=n), n=0..25); # Zerinvary Lajos, Jun 19 2008
    g:=1/(1-3*z): gser:=series(g, z=0, 43): seq(coeff(gser, z, n)+1, n=0..31); # Zerinvary Lajos, Jan 09 2009
  • Mathematica
    Table[3^n + 1, {n, 0, 24}]
  • PARI
    a(n) = 3^n + 1
    
  • PARI
    Vec(2*(1-2*x)/((1-x)*(1-3*x)) + O(x^50)) \\ Altug Alkan, Nov 15 2015
    
  • Sage
    [lucas_number2(n,4,3) for n in range(27)] # Zerinvary Lajos, Jul 08 2008
    
  • Sage
    [sigma(3,n) for n in range(27)] # Zerinvary Lajos, Jun 04 2009
    
  • Sage
    [3^n+1 for n in range(30)] # Bruno Berselli, Jan 11 2017
    

Formula

a(n) = 3*a(n-1) - 2 = 4*a(n-1) - 3*a(n-2). (Lucas sequence, with A003462, associated to the pair (4, 3).)
G.f.: 2*(1-2*x)/((1-x)*(1-3*x)). Inverse binomial transforms yields 2,2,4,8,16,... i.e., A000079 with the first entry changed to 2. Binomial transform yields A063376 without A063376(-1). - R. J. Mathar, Sep 05 2008
E.g.f.: exp(x) + exp(3*x). - Mohammad K. Azarian, Jan 02 2009
a(n) = A279396(n+3,3). - Wolfdieter Lang, Jan 10 2017
a(n) = 2*A007051(n). - R. J. Mathar, Apr 07 2022

Extensions

Additional comments from Rick L. Shepherd, Feb 13 2002

A052539 a(n) = 4^n + 1.

Original entry on oeis.org

2, 5, 17, 65, 257, 1025, 4097, 16385, 65537, 262145, 1048577, 4194305, 16777217, 67108865, 268435457, 1073741825, 4294967297, 17179869185, 68719476737, 274877906945, 1099511627777, 4398046511105, 17592186044417
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Comments

The sequence is a Lucas sequence V(P,Q) with P = 5 and Q = 4, so if n is a prime number, then V_n(5,4) - 5 is divisible by n. The smallest pseudoprime q which divides V_q(5,4) - 5 is 15.
Also the edge cover number of the (n+1)-Sierpinski tetrahedron graph. - Eric W. Weisstein, Sep 20 2017
First bisection of A000051, A049332, A052531 and A014551. - Klaus Purath, Sep 23 2020

Crossrefs

Programs

  • GAP
    List([0..30], n-> 4^n+1); # G. C. Greubel, May 09 2019
  • Magma
    [4^n+1: n in [0..30] ]; // Vincenzo Librandi, Apr 30 2011
    
  • Maple
    spec := [S,{S=Union(Sequence(Union(Z,Z,Z,Z)),Sequence(Z))},unlabeled]: seq(combstruct[count](spec,size=n), n=0..30);
    A052539:=n->4^n + 1; seq(A052539(n), n=0..30); # Wesley Ivan Hurt, Jun 12 2014
  • Mathematica
    Table[4^n + 1, {n, 0, 30}]
    (* From Eric W. Weisstein, Sep 20 2017 *)
    4^Range[0, 30] + 1
    LinearRecurrence[{5, -4}, {2, 5}, 30]
    CoefficientList[Series[(2-5x)/(1-5x+4x^2), {x, 0, 30}], x] (* End *)
  • PARI
    a(n)=4^n+1 \\ Charles R Greathouse IV, Nov 20 2011
    
  • Sage
    [4^n+1 for n in (0..30)] # G. C. Greubel, May 09 2019
    

Formula

a(n) = 4^n + 1.
a(n) = 4*a(n-1) - 3 = 5*a(n-1) - 4*a(n-2).
G.f.: (2 - 5*x)/((1 - 4*x)*(1 - x)).
E.g.f.: exp(x) + exp(4*x). - Mohammad K. Azarian, Jan 02 2009
From Klaus Purath, Sep 23 2020: (Start)
a(n) = 3*4^(n-1) + a(n-1).
a(n) = (a(n-1)^2 + 9*4^(n-2))/a(n-2).
a(n) = A178675(n) - 3. (End)

A034474 a(n) = 5^n + 1.

Original entry on oeis.org

2, 6, 26, 126, 626, 3126, 15626, 78126, 390626, 1953126, 9765626, 48828126, 244140626, 1220703126, 6103515626, 30517578126, 152587890626, 762939453126, 3814697265626, 19073486328126, 95367431640626, 476837158203126
Offset: 0

Views

Author

Keywords

Comments

a(n) is the deficiency of 3*5^n (see A033879). - Patrick J. McNab, May 28 2017

Examples

			G.f. = 2 + 6*x + 26*x^2 + 126*x^3 + 626*x^4 + 3126*x^5 + 15626*x^6 + ...
		

Crossrefs

Programs

Formula

a(n) = 5*a(n-1) - 4 with a(0) = 2.
a(n) = 6*a(n-1) - 5*a(n-2) for n > 1.
From Mohammad K. Azarian, Jan 02 2009: (Start)
G.f.: 1/(1-x) + 1/(1-5*x) = (2-6*x)/((1-x)*(1-5*x)).
E.g.f.: exp(x) + exp(5*x). (End)
a(n) = A279396(n+5,5). - Wolfdieter Lang, Jan 10 2017
From Elmo R. Oliveira, Dec 06 2023: (Start)
a(n) = A000351(n) + 1.
a(n) = 2*A034478(n). (End)

A062394 a(n) = 6^n + 1.

Original entry on oeis.org

2, 7, 37, 217, 1297, 7777, 46657, 279937, 1679617, 10077697, 60466177, 362797057, 2176782337, 13060694017, 78364164097, 470184984577, 2821109907457, 16926659444737, 101559956668417, 609359740010497, 3656158440062977
Offset: 0

Views

Author

Henry Bottomley, Jun 22 2001

Keywords

Crossrefs

Sequences of the form m^n + 1: A000012 (m=0), A007395 (m=1), A000051 (m=2), A034472 (m=3), A052539 (m=4), A034474 (m=5), this sequence (m=6), A034491 (m=7), A062395 (m=8), A062396 (m=9), A062397 (m=10), A034524 (m=11), A178248 (m=12), A141012 (m=13), A228081 (m=64).
Cf. A000400.

Programs

  • Magma
    [6^n + 1: n in [0..30] ]; // Vincenzo Librandi, Apr 30 2011
    
  • Mathematica
    6^Range[0,30] +1
    LinearRecurrence[{7,-6},{2,7},30] (* Harvey P. Dale, Aug 11 2015 *)
  • PARI
    vector(20, n, n--; 6^n + 1) \\ Michel Marcus, Jun 11 2015
    
  • SageMath
    [6^n+1 for n in range(31)] # G. C. Greubel, Mar 11 2023

Formula

a(n) = 6*a(n-1) - 5.
a(n) = A000400(n) + 1.
a(n) = 7*a(n-1) - 6*a(n-2).
From Mohammad K. Azarian, Jan 02 2009: (Start)
G.f.: 1/(1-x) + 1/(1-6*x).
E.g.f.: exp(x) + exp(6*x). (End)

A034491 a(n) = 7^n + 1.

Original entry on oeis.org

2, 8, 50, 344, 2402, 16808, 117650, 823544, 5764802, 40353608, 282475250, 1977326744, 13841287202, 96889010408, 678223072850, 4747561509944, 33232930569602, 232630513987208, 1628413597910450, 11398895185373144
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

Formula

a(n) = 7*a(n-1) - 6.
a(n) = 8*a(n-1) - 7*a(n-2).
From Mohammad K. Azarian, Jan 02 2009: (Start)
G.f.: 1/(1-x) + 1/(1-7*x).
E.g.f.: exp(x) + exp(7*x). (End)
a(n) = 2*A034494(n). - G. C. Greubel, Mar 11 2023

A062395 a(n) = 8^n + 1.

Original entry on oeis.org

2, 9, 65, 513, 4097, 32769, 262145, 2097153, 16777217, 134217729, 1073741825, 8589934593, 68719476737, 549755813889, 4398046511105, 35184372088833, 281474976710657, 2251799813685249, 18014398509481985, 144115188075855873
Offset: 0

Views

Author

Henry Bottomley, Jun 22 2001

Keywords

Comments

Any number of the form b^k+1 is composite for b>2 and k odd since b+1 algebraically divides b^k+1. - Robert G. Wilson v, Aug 25 2002

References

  • D. M. Burton, Elementary Number Theory, Allyn and Bacon, Boston, MA, 1976, pp. 51.
  • G. Everest, A. van der Poorten, I. Shparlinski and T. Ward, Recurrence Sequences, Amer. Math. Soc., 2003; see esp. p. 255.

Crossrefs

Programs

  • Magma
    [8^n + 1: n in [0..40] ]; // Vincenzo Librandi, Apr 30 2011
  • Mathematica
    Table[8^n + 1, {n, 0, 20}]
    LinearRecurrence[{9,-8},{2,9},20] (* Harvey P. Dale, Jan 24 2019 *)
  • PARI
    for(n=0,22,print(8^n+1))
    

Formula

a(n) = 8a(n-1)-7 = A001018(n)+1 = 9a(n-1) - 8a(n-2).
G.f.: -(-2+9*x)/(-1+x)/(-1+8*x). - R. J. Mathar, Nov 16 2007
E.g.f.: e^x+e^(8*x). - Mohammad K. Azarian, Jan 02 2009

A062396 a(n) = 9^n + 1.

Original entry on oeis.org

2, 10, 82, 730, 6562, 59050, 531442, 4782970, 43046722, 387420490, 3486784402, 31381059610, 282429536482, 2541865828330, 22876792454962, 205891132094650, 1853020188851842, 16677181699666570, 150094635296999122
Offset: 0

Views

Author

Henry Bottomley, Jun 22 2001

Keywords

Crossrefs

Programs

Formula

a(n) = 9*a(n-1) - 8 = A001019(n) + 1 = 10*a(n-1) - 9*a(n-2).
From Mohammad K. Azarian, Jan 02 2009: (Start)
G.f.: 1/(1-x) + 1/(1-9*x).
E.g.f.: e^x + e^(9*x). (End)

A062397 a(n) = 10^n + 1.

Original entry on oeis.org

2, 11, 101, 1001, 10001, 100001, 1000001, 10000001, 100000001, 1000000001, 10000000001, 100000000001, 1000000000001, 10000000000001, 100000000000001, 1000000000000001, 10000000000000001, 100000000000000001
Offset: 0

Views

Author

Henry Bottomley, Jun 22 2001

Keywords

Comments

The first three terms (indices 0, 1 and 2) are the only known primes. Moreover, the terms not of the form a(2^k) are all composite, except for a(0). Indeed, for all n >= 0, a(2n+1) is divisible by 11, a(4n+2) is divisible by 101, a(8n+4) is divisible by 73, a(16n+8) is divisible by 17, a(32n+16) is divisible by 353, a(64n+32) is divisible by 19841, etc. - M. F. Hasler, Nov 03 2018 [Edited based on the comment by Jeppe Stig Nielsen, Oct 17 2019]
This sequence also results when each term is generated by converting the previous term into a Roman numeral, then replacing each letter with its corresponding decimal value, provided that the vinculum is used and numerals are written in a specific way for integers greater than 3999, e.g., IV with a vinculum over the I and V for 4000. - Jamie Robert Creasey, Apr 14 2021
By Mihăilescu's theorem, a(n) can never be a perfect power (see "Catalan's conjecture" in Links). - Marco Ripà, Mar 10 2025

Crossrefs

Except for the initial term, essentially the same as A000533. Cf. A054977, A007395, A000051, A034472, A052539, A034474, A062394, A034491, A062395, A062396, A007689, A063376, A063481, A074600-A074624, A034524, A178248, A228081 for numbers one more than powers, i.e., this sequence translated from base n (> 2) to base 10.
Cf. A038371 (smallest prime factor), A185121.

Programs

Formula

a(n) = 10*a(n-1) - 9 = A011557(n) + 1 = A002283(n) + 2.
From Mohammad K. Azarian, Jan 02 2009: (Start)
G.f.: 1/(1-x) + 1/(1-10*x).
E.g.f.: exp(x) + exp(10*x). (End)
Showing 1-10 of 22 results. Next