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

A007051 a(n) = (3^n + 1)/2.

Original entry on oeis.org

1, 2, 5, 14, 41, 122, 365, 1094, 3281, 9842, 29525, 88574, 265721, 797162, 2391485, 7174454, 21523361, 64570082, 193710245, 581130734, 1743392201, 5230176602, 15690529805, 47071589414, 141214768241, 423644304722, 1270932914165, 3812798742494, 11438396227481
Offset: 0

Views

Author

Colin Mallows, N. J. A. Sloane, Simon Plouffe, Robert G. Wilson v

Keywords

Comments

Number of ordered trees with n edges and height at most 4.
Number of palindromic structures using a maximum of three different symbols. - Marks R. Nester
Number of compositions of all even natural numbers into n parts <= 2 (0 is counted as a part), see example. - Adi Dani, May 14 2011
Consider the mapping f(a/b) = (a + 2*b)/(2*a + b). Taking a = 1, b = 2 to start with, and carrying out this mapping repeatedly on each new (reduced) rational number gives the sequence 1/2, 4/5, 13/14, 40/41, ... converging to 1. The sequence contains the denominators = (3^n+1)/2. The same mapping for N, i.e., f(a/b) = (a + N*b)/(a+b) gives fractions converging to N^(1/2). - Amarnath Murthy, Mar 22 2003
Second binomial transform of the expansion of cosh(x). - Paul Barry, Apr 05 2003
The sequence (1, 1, 2, 5, ...) = 3^n/6 + 1/2 + 0^n/3 has binomial transform A007581. - Paul Barry, Jul 20 2003
Number of (s(0), s(1), ..., s(2n+2)) such that 0 < s(i) < 6 and |s(i) - s(i-1)| = 1 for i = 1, 2, ..., 2n+2, s(0) = 1, s(2n+2) = 1. - Herbert Kociemba, Jun 10 2004
Density of regular language L over {1,2,3}^* (i.e., number of strings of length n in L) described by regular expression 11*+11*2(1+2)*+11*2(1+2)*3(1+2+3)*. - Nelma Moreira, Oct 10 2004
Sums of rows of the triangle in A119258. - Reinhard Zumkeller, May 11 2006
Number of n-words from the alphabet A = {a,b,c} which contain an even number of a's. - Fung Cheok Yin (cheokyin_restart(AT)yahoo.com.hk), Aug 30 2006
Let P(A) be the power set of an n-element set A. Then a(n) = the number of pairs of elements {x,y} of P(A) for which either 0) x and y are disjoint and for which x is not a subset of y and y is not a subset of x, or 1) x = y. - Ross La Haye, Jan 10 2008
a(n+1) gives the number of primitive periodic multiplex juggling sequences of length n with base state <2>. - Steve Butler, Jan 21 2008
a(n) is also the number of idempotent order-preserving and order-decreasing partial transformations (of an n-chain). - Abdullahi Umar, Oct 02 2008
Equals row sums of triangle A147292. - Gary W. Adamson, Nov 05 2008
Equals leftmost column of A071919^3. - Gary W. Adamson, Apr 13 2009
A010888(a(n))=5 for n >= 2, that is, the digital root of the terms >= 5 equals 5. - Parthasarathy Nambi, Jun 03 2009
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*charpoly(A,2). - Milan Janjic, Jan 27 2010
Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=6, (i>1), A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n>=1, a(n)=(-1)^(n-1)*charpoly(A,3). - Milan Janjic, Feb 21 2010
It appears that if s(n) is a rational sequence of the form s(1)=2, s(n)= (2*s(n-1)+1)/(s(n-1)+2), n>1 then s(n)=a(n)/(a(n-1)-1).
Form an array with m(1,n)=1 and m(i,j) = Sum_{k=1..i-1} m(k,j) + Sum_{k=1..j-1} m(i,k), which is the sum of the terms to the left of m(i,j) plus the sum above m(i,j). The sum of the terms in antidiagonal(n-1) = a(n). - J. M. Bergot, Jul 16 2013
From Peter Bala, Oct 29 2013: (Start)
An Engel expansion of 3 to the base b := 3/2 as defined in A181565, with the associated series expansion 3 = b + b^2/2 + b^3/(2*5) + b^4/(2*5*14) + .... Cf. A034472.
More generally, for a positive integer n >= 3, the sequence [1, n - 1, n^2 - n - 1, ..., ( (n - 2)*n^k + 1 )/(n - 1), ...] is an Engel expansion of n/(n - 2) to the base n/(n - 1). Cases include A007583 (n = 4), A083065 (n = 5) and A083066 (n = 6). (End)
Diagonal elements (and one more than antidiagonal elements) of the matrix A^n where A=(2,1;1,2). - David Neil McGrath, Aug 17 2014
From M. Sinan Kul, Sep 07 2016: (Start)
a(n) is equal to the number of integer solutions to the following equation when x is equal to the product of n distinct primes: 1/x = 1/y + 1/z where 0 < x < y <= z.
If z = k*y where k is a fraction >= 1 then the solutions can be given as: y = ((k+1)/k)*x and z = (k+1)*x.
Here k can be equal to any divisor of x or to the ratio of two divisors.
For example for x = 2*3*5 = 30 (product of three distinct primes), k would have the following 14 values: 1, 6/5, 3/2, 5/3, 2, 5/2, 3, 10/3, 5, 6, 15/2, 10, 15, 30.
As an example for k = 10/3, we would have y=39, z=130 and 1/39 + 1/130 = 1/30.
Here finding the number of fractions would be equivalent to distributing n balls (distinct primes) to two bins (numerator and denominator) with no empty bins which can be found using Stirling numbers of the second kind. So another definition for a(n) is: a(n) = 2^n + Sum_{i=2..n} Stirling2(i,2)*binomial(n,i).
(End)
a(n+1) is the smallest i for which the Catalan number C(i) (see A000108) is divisible by 3^n for n > 0. This follows from the rule given by Franklin T. Adams-Watters for determining the multiplicity with which a prime divides C(n). We need to find the smallest number in base 3 to achieve a given count. Applied to prime 3, 1 is the smallest digit that counts but requires to be followed by 2 which cannot be at the end to count. Therefore the number in base 3 of the form 1{n-1 times}20 = (3^(n+1) + 1)/2 + 1 = a(n+1)+1 is the smallest number to achieve count n which implies the claim. - Peter Schorn, Mar 06 2020
Let A be a Toeplitz matrix of order n, defined by: A[i,j]=1, if ij; A[i,i]=2. Then, for n>=1, a(n) = det A. - Dmitry Efimov, Oct 28 2021
a(n) is the least number k such that A065363(k) = -(n-1), for n > 0. - Amiram Eldar, Sep 03 2022

Examples

			From _Adi Dani_, May 14 2011: (Start)
a(3)=14 because all compositions of even natural numbers into 3 parts <=2 are
for 0: (0,0,0)
for 2: (0,1,1), (1,0,1), (1,1,0), (0,0,2), (0,2,0), (2,0,0)
for 4: (0,2,2), (2,0.2), (2,2,0), (1,1,2), (1,2,1), (2,1,1)
for 6: (2,2,2).
(End)

References

  • J. M. Borwein, D. H. Bailey and R. Girgensohn, Experimentation in Mathematics, A K Peters, Ltd., Natick, MA, 2004. x+357 pp. See p. 47.
  • Adi Dani, Quasicompositions of natural numbers, Proceedings of III congress of mathematicians of Macedonia, Struga Macedonia 29 IX -2 X 2005 pages 225-238.
  • R. K. Guy, Unsolved Problems in Number Theory, Springer, 1st edition, 1981. See section E11.
  • M. R. Nester (1999). Mathematical investigations of some plant interaction designs. PhD Thesis. University of Queensland, Brisbane, Australia. [See A056391 for pdf file of Chap. 2]
  • P. Ribenboim, The Book of Prime Number Records. Springer-Verlag, NY, 2nd ed., 1989, p. 60.
  • P. Ribenboim, The Little Book of Big Primes, Springer-Verlag, NY, 1991, p. 53.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Links

Crossrefs

Cf. A056449, A064881-A064886, A008277, A007581, A056272, A056273, A000392, A000079, A034472, A147292, A003462, A065363, A071919, A007583, A083065, A083066.
A row of the array in A278984.

Programs

Formula

a(n) = 3*a(n-1) - 1.
Binomial transform of Chebyshev coefficients A011782. - Paul Barry, Mar 16 2003
From Paul Barry, Mar 16 2003: (Start)
a(n) = 4*a(n-1) - 3*a(n-2) for n > 1, a(0)=1, a(1)=2.
G.f.: (1 - 2*x)/((1 - x)*(1 - 3*x)). (End)
E.g.f.: exp(2*x)*cosh(x). - Paul Barry, Apr 05 2003
a(n) = Sum_{k=0..floor(n/2)} binomial(n, 2*k)*2^(n-2*k). - Paul Barry, May 08 2003
This sequence is also the partial sums of the first 3 Stirling numbers of second kind: a(n) = S(n+1, 1) + S(n+1, 2) + S(n+1, 3) for n >= 0; alternatively it is the number of partitions of [n+1] into 3 or fewer parts. - Mike Zabrocki, Jun 21 2004
For c=3, a(n) = (c^n)/c! + Sum_{k=1..c-2}((k^n)/k!*(Sum_{j=2..c-k}(((-1)^j)/j!))) or = Sum_{k=1..c} g(k, c)*k^n where g(1, 1) = 1, g(1, c) = g(1, c-1) + ((-1)^(c-1))/(c-1)! for c > 1, and g(k, c) = g(k-1, c-1)/k for c > 1 and 2 <= k <= c. - Nelma Moreira, Oct 10 2004
The i-th term of the sequence is the entry (1, 1) in the i-th power of the 2 X 2 matrix M = ((2, 1), (1, 2)). - Simone Severini, Oct 15 2005
If p[i]=fibonacci(2i-3) 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
INVERT transform of A001519: [1, 1, 2, 5, 13, 34, ...]. - Gary W. Adamson, Jun 13 2011
a(n) = M^n*[1,1,1,0,0,0,...], leftmost column term; where M = an infinite bidiagonal matrix with all 1's in the superdiagonal and (1,2,3,...) in the main diagonal and the rest zeros. - Gary W. Adamson, Jun 23 2011
a(n) = M^n*[1,1,1,0,0,0,...], top entry term; where M is an infinite bidiagonal matrix with all 1's in the superdiagonal, (1,2,3,...) as the main diagonal, and the rest zeros. - Gary W. Adamson, Jun 24 2011
a(n) = A201730(n,0). - Philippe Deléham, Dec 05 2011
a(n) = A006342(n) + A006342(n-1). - Yuchun Ji, Sep 19 2018
From Dmitry Efimov, Oct 29 2021: (Start)
a(2*m+1) = Product_{k=-m..m} (2+i*tan(Pi*k/(2*m+1))),
a(2*m) = Product_{k=-m..m-1} (2+i*tan(Pi*(2*k+1)/(4*m))),
where i is the imaginary unit. (End)

A007583 a(n) = (2^(2*n + 1) + 1)/3.

Original entry on oeis.org

1, 3, 11, 43, 171, 683, 2731, 10923, 43691, 174763, 699051, 2796203, 11184811, 44739243, 178956971, 715827883, 2863311531, 11453246123, 45812984491, 183251937963, 733007751851, 2932031007403, 11728124029611, 46912496118443, 187649984473771, 750599937895083
Offset: 0

Views

Author

Simon Plouffe

Keywords

Comments

Let u(k), v(k), w(k) be the 3 sequences defined by u(1)=1, v(1)=0, w(1)=0 and u(k+1)=u(k)+v(k)-w(k), v(k+1)=u(k)-v(k)+w(k), w(k+1)=-u(k)+v(k)+w(k); let M(k)=Max(u(k),v(k),w(k)); then a(n)=M(2n)=M(2n-1). - Benoit Cloitre, Mar 25 2002
Also the number of words of length 2n generated by the two letters s and t that reduce to the identity 1 by using the relations ssssss=1, tt=1 and stst=1. The generators s and t along with the three relations generate the dihedral group D6=C2xD3. - Jamaine Paddyfoot (jay_paddyfoot(AT)hotmail.com) and John W. Layman, Jul 08 2002
Binomial transform of A025192. - Paul Barry, Apr 11 2003
Number of walks of length 2n+1 between two adjacent vertices in the cycle graph C_6. Example: a(1)=3 because in the cycle ABCDEF we have three walks of length 3 between A and B: ABAB, ABCB and AFAB. - Emeric Deutsch, Apr 01 2004
Numbers of the form 1 + Sum_{i=1..m} 2^(2*i-1). - Artur Jasinski, Feb 09 2007
Prime numbers of the form 1+Sum[2^(2n-1)] are in A000979. Numbers x such that 1+Sum[2^(2n-1)] is prime for n=1,2,...,x is A127936. - Artur Jasinski, Feb 09 2007
Related to A024493(6n+1), A131708(6n+3), A024495(6n+5). - Paul Curtz, Mar 27 2008
Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=-6, (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,2). - Milan Janjic, Feb 21 2010
Number of toothpicks in the toothpick structure of A139250 after 2^n stages. - Omar E. Pol, Feb 28 2011
Numbers whose binary representation is "10" repeated (n-1) times with "11" appended on the end, n >= 1. For example 171 = 10101011 (2). - Omar E. Pol, Nov 22 2012
a(n) is the smallest number for which A072219(a(n)) = 2*n+1. - Ramasamy Chandramouli, Dec 22 2012
An Engel expansion of 2 to the base b := 4/3 as defined in A181565, with the associated series expansion 2 = b + b^2/3 + b^3/(3*11) + b^4/(3*11*43) + .... Cf. A007051. - Peter Bala, Oct 29 2013
The positive integer solution (x,y) of 3*x - 2^n*y = 1, n>=0, with smallest x is (a(n/2), 2) if n is even and (a((n-1)/2), 1) if n is odd. - Wolfdieter Lang, Feb 15 2014
The smallest positive number that requires at least n additions and subtractions of powers of 2 to be formed. See Puzzling StackExchange link. - Alexander Cooke Jul 16 2023

References

  • H. W. Gould, Combinatorial Identities, Morgantown, 1972, (1.77), page 10.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Links

Crossrefs

Cf. also A006054, A006356, A005578, A075680.
Partial sums of A081294.
Cf. A002450, A004171, A007051, A083065, A083066, A083884.
Cf. A000978, A000979, A124400, A124401, A127936, A127955, A127956, A127957, A127958.
Cf. location of records in A007302.

Programs

  • GAP
    List([0..25], n-> (2^(2*n+1) + 1)/3); # G. C. Greubel, Dec 25 2019
  • Haskell
    a007583 = (`div` 3) . (+ 1) . a004171
-- Reinhard Zumkeller, Jan 09 2013
  • Magma
    [(2^(2*n+1) + 1)/3: n in [0..30] ]; // Vincenzo Librandi, Apr 28 2011
  • Maple
    a[0]:=1:for n from 1 to 50 do a[n]:=4*a[n-1]-1 od: seq(a[n], n=0..23); # Zerinvary Lajos, Feb 22 2008, with correction by K. Spage, Aug 20 2014
A007583 := proc(n)
    (2^(2*n+1)+1)/3 ;
end proc: # R. J. Mathar, Feb 19 2015
  • Mathematica
    (* From Michael De Vlieger, Aug 22 2016 *)
Table[(2^(2n+1) + 1)/3, {n, 0, 23}]
Table[1 + 2Sum[4^k, {k, 0, n-1}], {n, 0, 23}]
NestList[4# -1 &, 1, 23]
Table[Sum[Binomial[n+k, 2k]/2^(k-n), {k, 0, n}], {n, 0, 23}]
CoefficientList[Series[(1-2x)/(1-5x+4x^2), {x, 0, 23}], x] (* End *)
  • PARI
    a(n)=sum(k=-n\3,n\3,binomial(2*n+1,n+1+3*k))
  • PARI
    a=1; for(n=1,23, print1(a,", "); a=bitor(a,3*a)) \\ K. Spage, Aug 20 2014
  • PARI
    Vec((1-2*x)/(1-5*x+4*x^2) + O(x^30)) \\ Altug Alkan, Dec 08 2015
  • PARI
    apply( {A007583(n)=2<<(2*n)\/3}, [0..25]) \\ M. F. Hasler, Nov 30 2021
  • Sage
    [(2^(2*n+1) + 1)/3 for n in (0..25)] # G. C. Greubel, Dec 25 2019

Formula

a(n) = 2*A002450(n) + 1.
From Wolfdieter Lang, Apr 24 2001: (Start)
a(n) = Sum_{m = 0..n} A060920(n, m) = A002450(n+1) - 2*A002450(n).
G.f.: (1-2*x)/(1-5*x+4*x^2). (End)
a(n) = Sum_{k = 0..n} binomial(n+k, 2*k)/2^(k - n).
a(n) = 4*a(n-1) - 1, n > 0.
From Paul Barry, Mar 17 2003: (Start)
a(n) = 1 + 2*Sum_{k = 0..n-1} 4^k;
a(n) = A001045(2n+1). (End)
a(n) = A020988(n-1) + 1 = A039301(n+1) - 1 = A083584(n-1) + 2. - Ralf Stephan, Jun 14 2003
a(0) = 1; a(n+1) = a(n) * 4 - 1. - Regis Decamps (decamps(AT)users.sf.net), Feb 04 2004 (correction to lead index by K. Spage, Aug 20 2014)
a(n) = Sum_{i + j + k = n; 0 <= i, j, k <= n} (n+k)!/i!/j!/(2*k)!. - Benoit Cloitre, Mar 25 2004
a(n) = 5*a(n-1) - 4*a(n-2). - Emeric Deutsch, Apr 01 2004
a(n) = 4^n - A001045(2*n). - Paul Barry, Apr 17 2004
a(n) = 2*(A001045(n))^2 + (A001045(n+1))^2. - Paul Barry, Jul 15 2004
a(n) = left and right terms in M^n * [1 1 1] where M = the 3X3 matrix [1 1 1 / 1 3 1 / 1 1 1]. M^n * [1 1 1] = [a(n) A002450(n+1) a(n)] E.g. a(3) = 43 since M^n * [1 1 1] = [43 85 43] = [a(3) A002450(4) a(3)]. - Gary W. Adamson, Dec 18 2004
a(n) = A072197(n) - A020988(n). - Creighton Dement, Dec 31 2004
a(n) = A139250(2^n). - Omar E. Pol, Feb 28 2011
a(n) = A193652(2*n+1). - Reinhard Zumkeller, Aug 08 2011
a(n) = Sum_{k = -floor(n/3)..floor(n/3)} binomial(2*n, n+3*k)/2. - Mircea Merca, Jan 28 2012
a(n) = 2^(2*(n+1)) - A072197(n). - Vladimir Pletser, Apr 12 2014
a(n) == 2*n + 1 (mod 3). Indeed, from Regis Decamps' formula (Feb 04 2004) we have a(i+1) - a(i) == -1 (mod 3), i= 0, 1, ..., n - 1. Summing, we have a(n) - 1 == -n (mod 3), and the formula follows. - Vladimir Shevelev, May 20 2015
For n > 0 a(n) = A133494(0) + 2 * (A133494(n) + Sum_{x = 1..n - 1}Sum_{k = 0..x - 1}(binomial(x - 1, k)*(A133494(k+1) + A133494(n-x+k)))). - J. Conrad, Dec 06 2015
a(n) = Sum_{k = 0..2n} (-2)^k == 1 + Sum_{k = 1..n} 2^(2k-1). - Bob Selcoe, Aug 21 2016
E.g.f.: (1 + 2*exp(3*x))*exp(x)/3. - Ilya Gutkovskiy, Aug 21 2016
A075680(a(n)) = 1, for n > 0. - Ralf Stephan, Jun 17 2025

A002253 Numbers k such that 3*2^k + 1 is prime.

Original entry on oeis.org

1, 2, 5, 6, 8, 12, 18, 30, 36, 41, 66, 189, 201, 209, 276, 353, 408, 438, 534, 2208, 2816, 3168, 3189, 3912, 20909, 34350, 42294, 42665, 44685, 48150, 54792, 55182, 59973, 80190, 157169, 213321, 303093, 362765, 382449, 709968, 801978, 916773, 1832496, 2145353
Offset: 1

Views

Author

N. J. A. Sloane

Keywords

Comments

From Zak Seidov, Mar 08 2009: (Start)
List is complete up to 3941000 according to the list of RB & WK.
So far there are only 4 primes: 2, 5, 41, 353. (End)

References

  • D. E. Knuth, The Art of Computer Programming. Addison-Wesley, Reading, MA, Vol. 2, p. 614.
  • H. Riesel, "Prime numbers and computer methods for factorization", Progress in Mathematics, Vol. 57, Birkhauser, Boston, 1985, Chap. 4, see pp. 381-384.
  • 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).

Links

Crossrefs

Cf. A002254, A004119, A181565.
See A039687 for the actual primes.

Programs

  • PARI
    is(n)=isprime(3*2^n+1) \\ Charles R Greathouse IV, Feb 17 2017
  • PARI
    A2253=[1]; A002253(n)=for(k=#A2253, n-1, my(m=A2253[k]); until(ispseudoprime(3<M. F. Hasler, Mar 03 2023

Formula

a(n) = log_2((A039687(n)-1)/3) = floor(log_2(A039687(n)/3)). - M. F. Hasler, Mar 03 2023

Extensions

Corrected and extended according to the list of Ray Ballinger and Wilfrid Keller by Zak Seidov, Mar 08 2009
Edited by N. J. A. Sloane, Mar 13 2009
a(47) and a(48) from the Ballinger & Keller web page, Joerg Arndt, Apr 07 2013
a(49) from https://t5k.org/primes/page.php?id=116922, Fabrice Le Foulher, Mar 09 2014
Terms moved from Data to b-file (Links), and additional term appended to b-file, by Jeppe Stig Nielsen, Oct 30 2020

A087289 a(n) = 2^(2*n+1) + 1.

Original entry on oeis.org

3, 9, 33, 129, 513, 2049, 8193, 32769, 131073, 524289, 2097153, 8388609, 33554433, 134217729, 536870913, 2147483649, 8589934593, 34359738369, 137438953473, 549755813889, 2199023255553, 8796093022209, 35184372088833, 140737488355329, 562949953421313, 2251799813685249
Offset: 0

Views

Author

W. Edwin Clark, Aug 29 2003

Keywords

Comments

Number of pairs of polynomials (f,g) in GF(2)[x] satisfying deg(f) <= n, deg(g) <= n and gcd(f,g) = 1.
An unpublished result due to Stephen Suen, David desJardins, and W. Edwin Clark. This is the case k = 2, q = 2 of their formula q^((n+1)*k) * (1 - 1/q^(k-1) + (q-1)/q^((n+1)*k)) for the number of ordered k-tuples (f_1, ..., f_k) of polynomials in GF(q)[x] such that deg(f_i) <= n for all i and gcd(f_1, ..., f_k) = 1.
Apparently the same as A084508 shifted left.
Terms in binary are palindromes of the form 1x1 where x is a string of 2*n zeros (A152577). - Brad Clardy, Sep 01 2011
For n > 0, a(n) is the number k such that the number of iterations of the map k -> (3k +1)/8 == 4 (mod 8) until reaching (3k +1)/8 <> 4 (mod 8) equals n. (see the Collatz problem: the start of the parity trajectory of a(n) is n times {100} = 100100100100...100abcd...). - Michel Lagneau, Jan 23 2012
An Engel expansion of 2 to the base 4 as defined in A181565, with the associated series expansion 2 = 4/3 + 4^2/(3*9) + 4^3/(3*9*33) + 4^4/(3*9*33*129) + .... Cf. A199561 and A207262. - Peter Bala, Oct 29 2013
For x = A083420(n), y = A000079(n+1), z = a(n) then x^2 + 2*y^2 = z^2. - Vincenzo Librandi, Jun 09 2014
A254046(n+1) is the 3-adic valuation of a(n). - Fred Daniel Kline, Jan 11 2017

Examples

			a(0) = 3 since there are three pairs, (0,1), (1,0) and (1,1) of polynomials (f,g) in GF(2)[x] of degree at most 0 such that gcd(f,g) = 1.

Links

Crossrefs

Cf. A087290, A087291, A087292, A099393.
Equals A004171 + 1.
Cf. also A007583, A181565, A199561, A207262, A254046, A283070.
Cf. A000079, A083420, A084508, A152577.

Programs

  • Magma
    [2^(2*n+1) + 1: n in [0..30]]; // Vincenzo Librandi, May 16 2011
  • Mathematica
    Table[2^(2 n + 1) + 1, {n, 0, 20}] (* or *) 3 NestList[4 # - 1 &, 1, 20]
(* or *) CoefficientList[Series[(3 - 6 x)/((1 - x) (1 - 4 x)), {x, 0, 20}], x] (* Michael De Vlieger, Mar 03 2017 *)
  • PARI
    a(n)=2^(2*n+1)+1 \\ Charles R Greathouse IV, Sep 24 2015

Formula

G.f.: (3-6*x)/((1-x)*(1-4*x)).
a(n) = 3*A007583(n).
a(n) = 4*a(n-1) - 3. - Lekraj Beedassy, Apr 29 2005
a(n) = A099393(n+1) - 2*A099393(n). - Brad Clardy, Sep 01 2011
a(n) = 2^(2*n + 1)*a(-1-n) for all n in Z. - Michael Somos, Jan 11 2017
a(n) = A283070(n) - 1. - Peter M. Chema, Mar 02 2017
From Elmo R. Oliveira, Feb 22 2025: (Start)
E.g.f.: exp(x)*(2*exp(3*x) + 1).
a(n) = 5*a(n-1) - 4*a(n-2). (End)

A004119 a(0) = 1; thereafter a(n) = 3*2^(n-1) + 1.

Original entry on oeis.org

1, 4, 7, 13, 25, 49, 97, 193, 385, 769, 1537, 3073, 6145, 12289, 24577, 49153, 98305, 196609, 393217, 786433, 1572865, 3145729, 6291457, 12582913, 25165825, 50331649, 100663297, 201326593, 402653185, 805306369, 1610612737, 3221225473, 6442450945, 12884901889
Offset: 0

Views

Author

N. J. A. Sloane, R. K. Guy

Keywords

Comments

Also Pisot sequence L(4,7) (cf. A008776).
Alternatively, define the sequence S(a(1),a(2)) by: a(n+2) is the least integer such that a(n+2)/a(n+1) > a(n+1)/a(n) for n > 0. This is S(4,7).
a(n) = number of terms of the arithmetic progression with first term 2^(2n-1) and last term 2^(2n+1). So common difference is 2^n. E.g., a(2)=7 corresponds to (8,12,16,20,24,28,32). - Augustine O. Munagi, Feb 21 2007
Equals binomial transform of [1, 3, 0, 3, 0, 3, 0, 3, ...]. - Gary W. Adamson, Aug 27 2010

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Links

Crossrefs

Cf. A049988, A049775, A000051, A008776.
A181565 is an essentially identical sequence.
For primes see A002253 and A039687.

Programs

  • Magma
    [1] cat [n le 1 select 4 else 2*Self(n-1)-1: n in [1..40]]; // Vincenzo Librandi, Dec 16 2015
  • Maple
    A004119:=-(-1-z+3*z**2)/(2*z-1)/(z-1); # Simon Plouffe in his 1992 dissertation
  • Mathematica
    s=4;lst={1,s};Do[s=s+(s-1);AppendTo[lst,s],{n,5!}];lst (* Vladimir Joseph Stephan Orlovsky, Nov 30 2009 *)
Prepend[Table[3*2^n + 1, {n, 0, 32}], 1] (* or *)
{1}~Join~LinearRecurrence[{3, -2}, {4, 7}, 33] (* Michael De Vlieger, Dec 16 2015 *)
  • PARI
    a(n)=3<Charles R Greathouse IV, Sep 28 2015

Formula

a(n) = 3a(n-1) - 2a(n-2).
For n>3, a(3)=13, a(n)= a(n-1)+2*floor(a(n-1)/2). - Benoit Cloitre, Aug 14 2002
For n>=1, a(n) = A049775(n+1)/2^(n-2). For n>=2, a(n) = 2a(n-1)-1; see also A000051. - Philippe Deléham, Feb 20 2004
O.g.f.: -(-1-x+3*x^2)/((2*x-1)*(x-1)). - R. J. Mathar, Nov 23 2007
For n>0, a(n) = 2*a(n-1)-1. - Vincenzo Librandi, Dec 16 2015
E.g.f.: exp(x)*(1 + 3*sinh(x)). - Stefano Spezia, May 06 2023

Extensions

Edited by N. J. A. Sloane, Dec 16 2015 at the suggestion of Bruno Berselli

A039687 Primes of the form 3*2^k + 1.

Original entry on oeis.org

7, 13, 97, 193, 769, 12289, 786433, 3221225473, 206158430209, 6597069766657, 221360928884514619393, 2353913150770005286438421033702874906038383291674012942337
Offset: 1

Views

Author

Felice Russo

Keywords

Comments

Primes of the form 6m+1 (A002476) of A074781. - Bernard Schott, Dec 14 2020

Links

Crossrefs

For more terms see A002253. These are the primes in A004119 (or A181565).
Cf. A002476, A074781.
Subsequence of A081091.

Programs

Formula

a(n) = 3*2^A002253(n) + 1. - M. F. Hasler, Mar 03 2023

A252753 Tree of Eratosthenes: a(0) = 1, a(1) = 2; after which, a(2n) = A250469(a(n)), a(2n+1) = 2 * a(n).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 9, 8, 7, 10, 15, 12, 25, 18, 21, 16, 11, 14, 27, 20, 35, 30, 33, 24, 49, 50, 51, 36, 55, 42, 45, 32, 13, 22, 39, 28, 65, 54, 57, 40, 77, 70, 87, 60, 85, 66, 69, 48, 121, 98, 147, 100, 125, 102, 105, 72, 91, 110, 123, 84, 115, 90, 93, 64, 17, 26, 63, 44, 95, 78, 81, 56, 119, 130, 159, 108, 145, 114, 117, 80
Offset: 0

Views

Author

Antti Karttunen, Jan 02 2015

Keywords

Comments

This sequence can be represented as a binary tree. Each child to the left is obtained by applying A250469 to the parent, and each child to the right is obtained by doubling the parent:
1
|
...................2...................
3 4
5......../ \........6 9......../ \........8
/ \ / \ / \ / \
/ \ / \ / \ / \
/ \ / \ / \ / \
7 10 15 12 25 18 21 16
11 14 27 20 35 30 33 24 49 50 51 36 55 42 45 32
etc.
Sequence A252755 is the mirror image of the same tree. A253555(n) gives the distance of n from 1 in both trees.

Links

Crossrefs

Inverse: A252754.
Row sums: A253787, products: A253788.
Fixed points of a(n-1): A253789.
Similar permutations: A005940, A252755, A054429, A250245.
Cf. also A001248, A001511, A055396, A083221, A181565, A250469, A253555.

Programs

  • Mathematica
    (* b = A250469 *)
b[1] = 1; b[n_] := If[PrimeQ[n], NextPrime[n], m1 = p1 = FactorInteger[n][[ 1, 1]]; For[k1 = 1, m1 <= n, m1 += p1; If[m1 == n, Break[]]; If[ FactorInteger[m1][[1, 1]] == p1, k1++]]; m2 = p2 = NextPrime[p1]; For[k2 = 1, True, m2 += p2, If[FactorInteger[m2][[1, 1]] == p2, k2++]; If[k1+2 == k2, Return[m2]]]];
a[0] = 1; a[1] = 2; a[n_] := a[n] = If[EvenQ[n], b[a[n/2]], 2 a[(n-1)/2]];
Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Mar 08 2016 *)

Formula

a(0) = 1, a(1) = 2; after which, a(2n) = A250469(a(n)), a(2n+1) = 2 * a(n).
As a composition of related permutations:
a(n) = A252755(A054429(n)).
a(n) = A250245(A005940(1+n)).
Other identities. For all n >= 1:
A055396(a(n)) = A001511(n). [A005940 has the same property.]
a(A003945(n)) = A001248(n) for n>=1. - Peter Luschny, Jan 13 2015

A083575 a(0) = 6; for n>0, a(n) = 2*a(n-1) - 1.

Original entry on oeis.org

6, 11, 21, 41, 81, 161, 321, 641, 1281, 2561, 5121, 10241, 20481, 40961, 81921, 163841, 327681, 655361, 1310721, 2621441, 5242881, 10485761, 20971521, 41943041, 83886081, 167772161, 335544321, 671088641, 1342177281, 2684354561, 5368709121, 10737418241
Offset: 0

Views

Author

N. J. A. Sloane, Jun 15 2003

Keywords

Comments

The primes in this sequence are listed in A050526. - M. F. Hasler, Oct 30 2010
An Engel expansion of 2/5 to the base 2 as defined in A181565, with the associated series expansion 2/5 = 2/6 + 2^2/(6*11) + 2^3/(6*11*21) + 2^4/(6*11*21*41) + ... . - Peter Bala, Oct 29 2013

Links

Crossrefs

Cf. A020737, A050526, A083683, A083686, A083705, A168596, A181565, A195744.

Programs

  • Magma
    [5*2^n+1 : n in [0..30]]; // Vincenzo Librandi, Nov 03 2011
  • Mathematica
    NestList[2#-1&,6,40] (* Harvey P. Dale, Jun 23 2017 *)
  • PARI
    a(n)=5<M. F. Hasler, Oct 30 2010
  • PARI
    Vec((6-7*x)/((1-x)*(1-2*x)) + O(x^40)) \\ Colin Barker, Sep 20 2016

Formula

a(n) = 5*2^n + 1. - M. F. Hasler, Oct 30 2010
a(n) = 3*a(n-1) - 2*a(n-2), n>1. - Vincenzo Librandi, Nov 03 2011
G.f.: (6-7*x) / ((1-x)*(1-2*x)). - Colin Barker, Sep 20 2016
E.g.f.: exp(x)*(1 + 5*exp(x)). - Stefano Spezia, Oct 08 2022
Product_{n>=0} (1 + 1/a(n)) = 7/5. - Amiram Eldar, Aug 04 2024

A083066 5th row of number array A083064.

Original entry on oeis.org

1, 5, 29, 173, 1037, 6221, 37325, 223949, 1343693, 8062157, 48372941, 290237645, 1741425869, 10448555213, 62691331277, 376147987661, 2256887925965, 13541327555789, 81247965334733, 487487792008397, 2924926752050381
Offset: 0

Views

Author

Paul Barry, Apr 21 2003

Keywords

Comments

Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=8, (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,2). - Milan Janjic, Feb 21 2010
An Engel expansion of 3/2 to the base b := 6/5 as defined in A181565, with the associated series expansion 3/2 = b + b^2/5 + b^3/(5*29) + b^4/(5*29*173) + .... Cf. A007051. - Peter Bala, Oct 29 2013

Links

Crossrefs

Cf. A083064, A083065, A083067, A007051, A007583.

Programs

Formula

a(n) = (4*6^n+1)/5.
G.f.: (1-2*x)/((1-6*x)*(1-x)).
E.g.f.: (4*exp(6*x)+exp(x))/5.
a(n) = 6*a(n-1)-1 with n>0, a(0)=1. - Vincenzo Librandi, Aug 08 2010
a(n) = 7*a(n-1)-6*a(n-2). - Vincenzo Librandi, Nov 04 2011
a(n) = 6^n - Sum_{i=0..n-1} 6^i for n>0. - Bruno Berselli, Jun 20 2013

A136412 a(n) = (5*4^n + 1)/3.

Original entry on oeis.org

2, 7, 27, 107, 427, 1707, 6827, 27307, 109227, 436907, 1747627, 6990507, 27962027, 111848107, 447392427, 1789569707, 7158278827, 28633115307, 114532461227, 458129844907, 1832519379627, 7330077518507, 29320310074027
Offset: 0

Views

Author

Paul Curtz, Mar 31 2008

Keywords

Comments

An Engel expansion of 4/5 to the base b := 4/3 as defined in A181565, with the associated series expansion 4/5 = b/2 + b^2/(2*7) + b^3/(2*7*27) + b^4/(2*7*27*107) + .... Cf. A199115 and A140660. - Peter Bala, Oct 29 2013

Links

Crossrefs

Sequences of the form (m*4^n + 1)/3: A007583 (m=2), this sequence (m=5), A199210 (m=11), A199210 (m=11), A206373 (m=14).
Cf. A007302, A140660, A181565, A199115.

Programs

Formula

a(n) = 4*a(n-1) - 1.
a(n) = A199115(n)/3.
O.g.f.: (2-3*x)/((1-x)*(1-4*x)). - R. J. Mathar, Apr 04 2008
a(n) = 5*a(n-1) - 4*a(n-2). - Vincenzo Librandi, Nov 04 2011
E.g.f.: (1/3)*(5*exp(4*x) + exp(x)). - G. C. Greubel, Jan 19 2023

Extensions

Formula in definition and more terms from R. J. Mathar, Apr 04 2008
Showing 1-10 of 39 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.