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.

Previous Showing 51-60 of 587 results. Next

A046717 a(n) = 2*a(n-1) + 3*a(n-2), a(0) = a(1) = 1.

Original entry on oeis.org

1, 1, 5, 13, 41, 121, 365, 1093, 3281, 9841, 29525, 88573, 265721, 797161, 2391485, 7174453, 21523361, 64570081, 193710245, 581130733, 1743392201, 5230176601, 15690529805, 47071589413, 141214768241, 423644304721, 1270932914165, 3812798742493, 11438396227481
Offset: 0

Views

Author

Gervais Deroo and M. Deroo

Keywords

Comments

Form the digraph with matrix A = [0,1,1,1; 1,0,1,1; 1,1,0,1; 1,0,1,1]. Then the sequence 0,1,1,5,... or (3^(n-1)-(-1)^n)/2+0^n/3 with g.f. x(1-x)/(1-2x-3x^2) corresponds to the (1,2) term of A^n. - Paul Barry, Oct 02 2004
3*a(n+1) + a(n) = 4*A060925(n); a(n+1) = A015518(n) + A060925(n); a(n+1) - 6*A015518(n) = (-1)^n. - Creighton Dement, Nov 15 2004
The sequence corresponds to the (1,1) term of the matrix [1,2;2,1]^n. - Simone Severini, Dec 04 2004
The same sequence may be obtained by the following process. Starting a priori with the fraction 1/1, the numerators of fractions built according to the rule: add top and bottom to get the new bottom, add top and 4 times the bottom to get the new top. The limit of the sequence of fractions is 2. - Cino Hilliard, Sep 25 2005
a(n)^2 + (2*A015518(n))^2 = a(2n). E.g., a(3) = 13, 2*A015518(3) = 14, A046717(6) = 365. 13^2 + 14^2 = 365. - Gary W. Adamson, Jun 17 2006
Equals INVERTi transform of A104934: (1, 2, 8, 28, 100, 356, 1268, ...). - Gary W. Adamson, Jul 21 2010
a(n) is the number of compositions of n when there are 1 type of 1 and 4 types of other natural numbers. - Milan Janjic, Aug 13 2010
An elephant sequence, see A175655. For the central square just one A[5] vector, with decimal value 341, leads to this sequence (without the first leading 1). For the corner squares this vector leads to the companion sequence A015518 (without the leading 0). - Johannes W. Meijer, Aug 15 2010
Pisano period lengths: 1, 1, 2, 1, 4, 2, 6, 4, 2, 4, 10, 2, 6, 6, 4, 8, 16, 2, 18, 4, ... - R. J. Mathar, Aug 10 2012
a(n) is the number of words of length n over a ternary alphabet whose position in the lexicographic order is a multiple of two. - Alois P. Heinz, Apr 13 2022
a(n) is the sum, for k=0..3, of the number of walks of length n between two vertices at distance k of the cube graph. - Miquel A. Fiol, Mar 09 2024

References

  • John Derbyshire, Prime Obsession, Joseph Henry Press, April 2004, see p. 16.

Crossrefs

The first difference sequence of A015518.
Row sums of triangle A080928.
The following sequences (and others) belong to the same family: A001333, A000129, A026150, A002605, A046717, A015518, A084057, A063727, A002533, A002532, A083098, A083099, A083100, A015519.
Cf. A015518.
Cf. A104934. - Gary W. Adamson, Jul 21 2010

Programs

  • Magma
    [n le 2 select 1 else 2*Self(n-1)+3*Self(n-2): n in [1..35]]; // Vincenzo Librandi, Jul 21 2013
    
  • Magma
    [(3^n + (-1)^n)/2: n in [0..30]]; // G. C. Greubel, Jan 07 2018
  • Maple
    a[0]:=1:a[1]:=1:for n from 2 to 50 do a[n]:=2*a[n-1]+3*a[n-2] od: seq(a[n], n=0..33); # Zerinvary Lajos, Dec 14 2008
    seq(denom(((-2)^(2*n)+6^(2*n))/((-2)^n+6^n)),n=0..26)
  • Mathematica
    Table[(3^n + (-1)^n)/2, {n, 0, 30}] (* Artur Jasinski, Dec 10 2006 *)
    CoefficientList[ Series[(1 - x)/(1 - 2x - 3x^2), {x, 0, 30}], x]  (* Robert G. Wilson v, Apr 04 2011 *)
    Table[ MatrixPower[{{1, 2}, {1, 1}}, n][[1, 1]], {n, 0, 30}] (* Robert G. Wilson v, Apr 04 2011 *)
  • PARI
    {a(n) = (3^n+(-1)^n)/2};
    for(n=0,30, print1(a(n), ", ")) /* modified by G. C. Greubel, Jan 07 2018 */
    
  • PARI
    x='x+O('x^30); Vec((1-x)/((1+x)*(1-3*x))) \\ G. C. Greubel, Jan 07 2018
    
  • Sage
    [lucas_number2(n,2,-3)/2 for n in range(0, 27)] # Zerinvary Lajos, Apr 30 2009
    

Formula

G.f.: (1-x)/((1+x)*(1-3*x)).
a(n) = (3^n + (-1)^n)/2.
a(n) = Sum_{k=0..n} binomial(n, 2k)2^(2k). - Paul Barry, Feb 26 2003
Binomial transform of A000302 (powers of 4) with interpolated zeros. Inverse binomial transform of A081294. - Paul Barry, Mar 17 2003
E.g.f.: exp(x)cosh(2x). - Paul Barry, Mar 17 2003
a(n) = ceiling(3^n/4) + floor(3^n/4) = ceiling(3^n/4)^2 - floor(3^n/4)^2. - Paul Barry, Jan 17 2005
a(n) = Sum_{k=0..n} Sum_{j=0..n} C(n,j)C(n-j,k)*(1+(-1)^(j-k))/2. - Paul Barry, May 21 2006
a(n) = Sum_{k=0..n} A098158(n,k)*4^(n-k). - Philippe Deléham, Dec 26 2007
a(n) = (3^n + (-1)^n)/2. - M. F. Hasler, Mar 20 2008
a(n) = 2 A015518(n) + (-1)^n; for n > 0, a(n) = A080925(n). - M. F. Hasler, Mar 20 2008
((1 + sqrt4)^n + (1 - sqrt4)^n)/2. The offset is 0. a(3)=13. - Al Hakanson (hawkuu(AT)gmail.com), Nov 22 2008
If p[1]=1 and p[i]=4 (i > 1), and if A is 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) = det A. - Milan Janjic, Apr 29 2010
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - x*(4*k-1)/(x*(4*k+3) - 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 26 2013
G.f.: G(0)/2, where G(k) = 1 + (-1)^k/(3^k - 3*9^k*x/(3*3^k*x + (-1)^k/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Oct 17 2013

Extensions

Description corrected by and more terms from Michael Somos

A000378 Sums of three squares: numbers of the form x^2 + y^2 + z^2.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 14, 16, 17, 18, 19, 20, 21, 22, 24, 25, 26, 27, 29, 30, 32, 33, 34, 35, 36, 37, 38, 40, 41, 42, 43, 44, 45, 46, 48, 49, 50, 51, 52, 53, 54, 56, 57, 58, 59, 61, 62, 64, 65, 66, 67, 68, 69, 70, 72, 73, 74, 75, 76, 77, 78, 80, 81, 82, 83
Offset: 1

Views

Author

Keywords

Comments

An equivalent definition: numbers of the form x^2 + y^2 + z^2 with x,y,z >= 0.
Bourgain studies "the spatial distribution of the representation of a large integer as a sum of three squares, on the small and critical scale as well as their electrostatic energy. The main results announced give strong evidence to the thesis that the solutions behave randomly. This is in sharp contrast to what happens with sums of two or four or more square." Sums of two nonzero squares are A000404. - Jonathan Vos Post, Apr 03 2012
The multiplicities for a(n) (if 0 <= x <= y <= z) are given as A000164(a(n)), n >= 1. Compare with A005875(a(n)) for integer x, y and z, and order taken into account. - Wolfdieter Lang, Apr 08 2013
a(n)^k is a member of this sequence for any k > 1. - Boris Putievskiy, May 05 2013
The selection rule for the planes with Miller indices (hkl) to undergo X-ray diffraction in a simple cubic lattice is h^2+k^2+l^2 = N where N is a term of this sequence. See A004014 for f.c.c. lattice. - Mohammed Yaseen, Nov 06 2022

Examples

			a(1) = 0 = 0^2 + 0^2 + 0^2. A005875(0) = 1 = A000164(0).
a(9) = 9 = 0^2 + 0^2 + 3^2 =  1^2 +  2^2 + 2^2. A000164(9) = 2. A000164(9) = 30 = 2*3 + 8*3 (counting signs and order). - _Wolfdieter Lang_, Apr 08 2013
		

References

  • J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 107.
  • E. Grosswald, Representations of Integers as Sums of Squares. Springer-Verlag, NY, 1985, p. 37.
  • R. K. Guy, Unsolved Problems in Number Theory, Springer, 1st edition, 1981. See section C20.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 311.

Crossrefs

Union of A000290, A000404 and A000408 (common elements).
Union of A000290, A000415 and A000419 (disjunct sets).
Complement of A004215.
Cf. A005875 (number of representations if x, y and z are integers).

Programs

  • Maple
    isA000378 := proc(n) # return true or false depending on n being in the list
        local x,y ;
        for x from 0 do
            if 3*x^2 > n then
                return false;
            end if;
            for y from x do
                if x^2+2*y^2 > n then
                    break;
                else
                    if issqr(n-x^2-y^2) then
                        return true;
                    end if;
                end if;
            end do:
        end do:
    end proc:
    A000378 := proc(n) # generate A000378(n)
        option remember;
        local a;
        if n = 1 then
            0;
        else
            for a from procname(n-1)+1 do
                if isA000378(a) then
                    return a;
                end if;
            end do:
        end if;
    end proc:
    seq(A000378(n),n=1..100) ; # R. J. Mathar, Sep 09 2015
  • Mathematica
    okQ[n_] := If[EvenQ[k = IntegerExponent[n, 2]], m = n/2^k; Mod[m, 8] != 7, True]; Select[Range[0, 100], okQ] (* Jean-François Alcover, Feb 08 2016, adapted from PARI *)
  • PARI
    isA000378(n)=my(k=valuation(n, 2)); if(k%2==0, n>>=k; n%8!=7, 1)
    
  • PARI
    list(lim)=my(v=List(),k,t); for(x=0,sqrtint(lim\=1), for(y=0, min(sqrtint(lim-x^2),x), k=x^2+y^2; for(z=0,min(sqrtint(lim-k), y), listput(v,k+z^2)))); Set(v) \\ Charles R Greathouse IV, Sep 14 2015
    
  • Python
    def valuation(n, b):
        v = 0
        while n > 1 and n%b == 0: n //= b; v += 1
        return v
    def ok(n): return n//4**valuation(n, 4)%8 != 7
    print(list(filter(ok, range(84)))) # Michael S. Branicky, Jul 15 2021
    
  • Python
    from itertools import count, islice
    def A000378_gen(): # generator of terms
        return filter(lambda n:n>>2*(bin(n)[:1:-1].index('1')//2) & 7 < 7, count(1))
    A000378_list = list(islice(A000378_gen(),30)) # Chai Wah Wu, Jun 27 2022
    
  • Python
    def A000378(n):
        def f(x): return n-1+sum(((x>>(i<<1))-7>>3)+1 for i in range(x.bit_length()>>1))
        m, k = n-1, f(n-1)
        while m != k: m, k = k, f(k)
        return m # Chai Wah Wu, Feb 14 2025

Formula

Legendre: a nonnegative integer is a sum of three squares iff it is not of the form 4^k m with m == 7 (mod 8).
n^(2k+1) is in the sequence iff n is in the sequence. - Ray Chandler, Feb 03 2009
Complement of A004215; complement of A000302(i)*A004771(j), i,j>=0. - Boris Putievskiy, May 05 2013
a(n) = 6n/5 + O(log n). - Charles R Greathouse IV, Mar 14 2014

Extensions

More terms from Ray Chandler, Sep 05 2004

A024036 a(n) = 4^n - 1.

Original entry on oeis.org

0, 3, 15, 63, 255, 1023, 4095, 16383, 65535, 262143, 1048575, 4194303, 16777215, 67108863, 268435455, 1073741823, 4294967295, 17179869183, 68719476735, 274877906943, 1099511627775, 4398046511103, 17592186044415, 70368744177663, 281474976710655
Offset: 0

Views

Author

Keywords

Comments

This sequence is the normalized length per iteration of the space-filling Peano-Hilbert curve. The curve remains in a square, but its length increases without bound. The length of the curve, after n iterations in a unit square, is a(n)*2^(-n) where a(n) = 4*a(n-1)+3. This is the sequence of a(n) values. a(n)*(2^(-n)*2^(-n)) tends to 1, the area of the square where the curve is generated, as n increases. The ratio between the number of segments of the curve at the n-th iteration (A015521) and a(n) tends to 4/5 as n increases. - Giorgio Balzarotti, Mar 16 2006
Numbers whose base-4 representation is 333....3. - Zerinvary Lajos, Feb 03 2007
From Eric Desbiaux, Jun 28 2009: (Start)
It appears that for a given area, a square n^2 can be divided into n^2+1 other squares.
It's a rotation and zoom out of a Cartesian plan, which creates squares with side
= sqrt( (n^2) / (n^2+1) ) --> A010503|A010532|A010541... --> limit 1,
and diagonal sqrt(2*sqrt((n^2)/(n^2+1))) --> A010767|... --> limit A002193.
(End)
Also the total number of line segments after the n-th stage in the H tree, if 4^(n-1) H's are added at the n-th stage to the structure in which every "H" is formed by 3 line segments. A164346 (the first differences of this sequence) gives the number of line segments added at the n-th stage. - Omar E. Pol, Feb 16 2013
a(n) is the cumulative number of segment deletions in a Koch snowflake after (n+1) iterations. - Ivan N. Ianakiev, Nov 22 2013
Inverse binomial transform of A005057. - Wesley Ivan Hurt, Apr 04 2014
For n > 0, a(n) is one-third the partial sums of A002063(n-1). - J. M. Bergot, May 23 2014
Also the cyclomatic number of the n-Sierpinski tetrahedron graph. - Eric W. Weisstein, Sep 18 2017

Examples

			G.f. = 3*x + 15*x^2 + 63*x^3 + 255*x^4 + 1023*x^5 + 4095*x^6 + ...
		

References

  • Graham Everest, Alf van der Poorten, Igor Shparlinski, and Thomas Ward, Recurrence Sequences, Amer. Math. Soc., 2003; see esp. p. 255.

Crossrefs

Programs

  • Haskell
    a024036 = (subtract 1) . a000302
    a024036_list = iterate ((+ 3) . (* 4)) 0
    -- Reinhard Zumkeller, Oct 03 2012
    
  • Maple
    A024036:=n->4^n-1; seq(A024036(n), n=0..30); # Wesley Ivan Hurt, Apr 04 2014
  • Mathematica
    Array[4^# - 1 &, 50, 0] (* Vladimir Joseph Stephan Orlovsky, Nov 03 2009 *)
    (* Start from Eric W. Weisstein, Sep 19 2017 *)
    Table[4^n - 1, {n, 0, 20}]
    4^Range[0, 20] - 1
    LinearRecurrence[{5, -4}, {0, 3}, 20]
    CoefficientList[Series[3 x/(1 - 5 x + 4 x^2), {x, 0, 20}], x]
    (* End *)
  • PARI
    for(n=0, 100, print1(4^n-1, ", ")) \\ Felix Fröhlich, Jul 04 2014
  • Sage
    [gaussian_binomial(2*n,1, 2) for n in range(21)] # Zerinvary Lajos, May 28 2009
    
  • Sage
    [stirling_number2(2*n+1, 2) for n in range(21)] # Zerinvary Lajos, Nov 26 2009
    

Formula

a(n) = 3*A002450(n). - N. J. A. Sloane, Feb 19 2004
G.f.: 3*x/((-1+x)*(-1+4*x)) = 1/(-1+x) - 1/(-1+4*x). - R. J. Mathar, Nov 23 2007
E.g.f.: exp(4*x) - exp(x). - Mohammad K. Azarian, Jan 14 2009
a(n) = A000051(n)*A000225(n). - Reinhard Zumkeller, Feb 14 2009
A079978(a(n)) = 1. - Reinhard Zumkeller, Nov 22 2009
a(n) = A179857(A000225(n)), for n > 0; a(n) > A179857(m), for m < A000225(n). - Reinhard Zumkeller, Jul 31 2010
a(n) = 4*a(n-1) + 3, with a(0) = 0. - Vincenzo Librandi, Aug 01 2010
A000120(a(n)) = 2*n. - Reinhard Zumkeller, Feb 07 2011
a(n) = (3/2)*A020988(n). - Omar E. Pol, Mar 15 2012
a(n) = (Sum_{i=0..n} A002001(i)) - 1 = A178789(n+1) - 3. - Ivan N. Ianakiev, Nov 22 2013
a(n) = n*E(2*n-1,1)/B(2*n,1), for n > 0, where E(n,x) denotes the Euler polynomials and B(n,x) the Bernoulli polynomials. - Peter Luschny, Apr 04 2014
a(n) = A000302(n) - 1. - Sean A. Irvine, Jun 18 2019
Sum_{n>=1} 1/a(n) = A248721. - Amiram Eldar, Nov 13 2020
a(n) = A080674(n) - A002450(n). - Elmo R. Oliveira, Dec 02 2023

Extensions

More terms Wesley Ivan Hurt, Apr 04 2014

A053738 If k is in sequence then 2*k and 2*k+1 are not (and 1 is in the sequence); numbers with an odd number of digits in binary.

Original entry on oeis.org

1, 4, 5, 6, 7, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109
Offset: 1

Views

Author

Henry Bottomley, Apr 06 2000

Keywords

Comments

Runs of successive numbers have lengths which are powers of 4.
Apparently, for any m>=1, 2^m is the largest power of 2 dividing sum(k=1,n,binomial(2k,k)^m) if and only if n is in the sequence. - Benoit Cloitre, Apr 27 2003
Numbers that begin with a 1 in base 4. - Michel Marcus, Dec 05 2013
The lower and upper asymptotic densities of this sequence are 1/3 and 2/3, respectively. - Amiram Eldar, Feb 01 2021

Crossrefs

Complement of A053754.

Programs

  • Maple
    seq(seq(i,i=4^k..2*4^k-1),k=0..5); # Robert Israel, Dec 28 2016
  • Mathematica
    Select[Range[110],OddQ[IntegerLength[#,2]]&] (* Harvey P. Dale, Sep 29 2012 *)
  • PARI
    isok(n, b=4) = digits(n, b)[1] == 1; \\ Michel Marcus, Dec 05 2013
    
  • PARI
    a(n) = n + 1<Kevin Ryde, Mar 27 2021

Formula

G.f.: x/(1-x)^2 + Sum_{k>=1} 2^(2k-1)*x^((4^k+2)/3)/(1-x). - Robert Israel, Dec 28 2016

A058622 a(n) = 2^(n-1) - ((1+(-1)^n)/4)*binomial(n, n/2).

Original entry on oeis.org

0, 1, 1, 4, 5, 16, 22, 64, 93, 256, 386, 1024, 1586, 4096, 6476, 16384, 26333, 65536, 106762, 262144, 431910, 1048576, 1744436, 4194304, 7036530, 16777216, 28354132, 67108864, 114159428, 268435456, 459312152, 1073741824, 1846943453
Offset: 0

Views

Author

Yong Kong (ykong(AT)curagen.com), Dec 29 2000

Keywords

Comments

a(n) is the number of n-digit binary sequences that have more 1's than 0's. - Geoffrey Critzer, Jul 16 2009
Maps to the number of walks that end above 0 on the number line with steps being 1 or -1. - Benjamin Phillabaum, Mar 06 2011
Chris Godsil observes that a(n) is the independence number of the (n+1)-folded cube graph; proof is by a Cvetkovic's eigenvalue bound to establish an upper bound and a direct construction of the independent set by looking at vertices at an odd (resp., even) distance from a fixed vertex when n is odd (resp., even). - Stan Wagon, Jan 29 2013
Also the number of subsets of {1,2,...,n} that contain more odd than even numbers. For example, for n=4, a(4)=5 and the 5 subsets are {1}, {3}, {1,3}, {1,2,3}, {1,3,4}. See A014495 when same number of even and odd numbers. - Enrique Navarrete, Feb 10 2018
Also half the number of length-n binary sequences with a different number of zeros than ones. This is also the number of integer compositions of n with nonzero alternating sum, where the alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i. Also the number of integer compositions of n+1 with alternating sum <= 0, ranked by A345915 (reverse: A345916). - Gus Wiseman, Jul 19 2021

Examples

			G.f. = x + x^2 + 4*x^3 + 5*x^4 + 16*x^5 + 22*x^6 + 64*x^7 + 93*x^8 + ...
From _Gus Wiseman_, Jul 19 2021: (Start)
The a(1) = 1 through a(5) = 16 compositions with nonzero alternating sum:
  (1)  (2)  (3)      (4)      (5)
            (1,2)    (1,3)    (1,4)
            (2,1)    (3,1)    (2,3)
            (1,1,1)  (1,1,2)  (3,2)
                     (2,1,1)  (4,1)
                              (1,1,3)
                              (1,2,2)
                              (1,3,1)
                              (2,1,2)
                              (2,2,1)
                              (3,1,1)
                              (1,1,1,2)
                              (1,1,2,1)
                              (1,2,1,1)
                              (2,1,1,1)
                              (1,1,1,1,1)
(End)
		

References

  • A. P. Prudnikov, Yu. A. Brychkov and O.I. Marichev, "Integrals and Series", Volume 1: "Elementary Functions", Chapter 4: "Finite Sums", New York, Gordon and Breach Science Publishers, 1986-1992, Eq. (4.2.1.7)

Crossrefs

The odd bisection is A000302.
The even bisection is A000346.
The following relate to compositions with nonzero alternating sum:
- The complement is counted by A001700 or A138364.
- The version for alternating sum > 0 is A027306.
- The unordered version is A086543 (even bisection: A182616).
- The version for alternating sum < 0 is A294175.
- These compositions are ranked by A345921.
A011782 counts compositions.
A097805 counts compositions by alternating (or reverse-alternating) sum.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A345197 counts compositions by length and alternating sum.
Compositions of n, 2n, or 2n+1 with alternating/reverse-alternating sum k:
- k = 0: counted by A088218, ranked by A344619/A344619.
- k = 1: counted by A000984, ranked by A345909/A345911.
- k = -1: counted by A001791, ranked by A345910/A345912.
- k = 2: counted by A088218, ranked by A345925/A345922.
- k = -2: counted by A002054, ranked by A345924/A345923.
- k >= 0: counted by A116406, ranked by A345913/A345914.
- k > 0: counted by A027306, ranked by A345917/A345918.
- k < 0: counted by A294175, ranked by A345919/A345920.
- k even: counted by A081294, ranked by A053754/A053754.
- k odd: counted by A000302, ranked by A053738/A053738.

Programs

  • Magma
    [(2^n -(1+(-1)^n)*Binomial(n, Floor(n/2))/2)/2: n in [0..40]]; // G. C. Greubel, Aug 08 2022
    
  • Mathematica
    Table[Sum[Binomial[n, Floor[n/2 + i]], {i, 1, n}], {n, 0, 32}] (* Geoffrey Critzer, Jul 16 2009 *)
    a[n_] := If[n < 0, 0, (2^n - Boole[EvenQ @ n] Binomial[n, Quotient[n, 2]])/2]; (* Michael Somos, Nov 22 2014 *)
    a[n_] := If[n < 0, 0, n! SeriesCoefficient[(Exp[2 x] - BesselI[0, 2 x])/2, {x, 0, n}]]; (* Michael Somos, Nov 22 2014 *)
    Table[2^(n - 1) - (1 + (-1)^n) Binomial[n, n/2]/4, {n, 0, 40}] (* Eric W. Weisstein, Mar 21 2018 *)
    CoefficientList[Series[2 x/((1-2x)(1 + 2x + Sqrt[(1+2x)(1-2x)])), {x, 0, 40}], x] (* Eric W. Weisstein, Mar 21 2018 *)
    ats[y_]:=Sum[(-1)^(i-1)*y[[i]],{i,Length[y]}];Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],ats[#]!=0&]],{n,0,15}] (* Gus Wiseman, Jul 19 2021 *)
  • PARI
    a(n) = 2^(n-1) - ((1+(-1)^n)/4)*binomial(n, n\2); \\ Michel Marcus, Dec 30 2015
    
  • PARI
    my(x='x+O('x^100)); concat(0, Vec(2*x/((1-2*x)*(1+2*x+((1+2*x)*(1-2*x))^(1/2))))) \\ Altug Alkan, Dec 30 2015
    
  • Python
    from math import comb
    def A058622(n): return (1<>1)>>1) if n else 0 # Chai Wah Wu, Aug 25 2025
  • SageMath
    [(2^n - binomial(n, n//2)*((n+1)%2))/2 for n in (0..40)] # G. C. Greubel, Aug 08 2022
    

Formula

a(n) = 2^(n-1) - ((1+(-1)^n)/4)*binomial(n, n/2).
a(n) = Sum_{i=0..floor((n-1)/2)} binomial(n, i).
G.f.: 2*x/((1-2*x)*(1+2*x+((1+2*x)*(1-2*x))^(1/2))). - Vladeta Jovovic, Apr 27 2003
E.g.f: (e^(2x)-I_0(2x))/2 where I_n is the Modified Bessel Function. - Benjamin Phillabaum, Mar 06 2011
Logarithmic derivative of the g.f. of A210736 is a(n+1). - Michael Somos, Nov 22 2014
Even index: a(2n) = 2^(n-1) - A088218(n). Odd index: a(2n+1) = 2^(2n). - Gus Wiseman, Jul 19 2021
D-finite with recurrence n*a(n) +2*(-n+1)*a(n-1) +4*(-n+1)*a(n-2) +8*(n-2)*a(n-3)=0. - R. J. Mathar, Sep 23 2021
a(n) = 2^n-A027306(n). - R. J. Mathar, Sep 23 2021
A027306(n) - a(n) = A126869(n). - R. J. Mathar, Sep 23 2021

A002697 a(n) = n*4^(n-1).

Original entry on oeis.org

0, 1, 8, 48, 256, 1280, 6144, 28672, 131072, 589824, 2621440, 11534336, 50331648, 218103808, 939524096, 4026531840, 17179869184, 73014444032, 309237645312, 1305670057984, 5497558138880, 23089744183296
Offset: 0

Views

Author

Keywords

Comments

Coefficient of x^(2n-2) in Chebyshev polynomial T(2n) is -a(n).
Let M_n be the n X n matrix m_(i,j) = 1 + 2*abs(i-j); then det(M_n) = (-1)^(n-1)*a(n-1). - Benoit Cloitre, May 28 2002
Number of subsequences 00 in all words of length n+1 on the alphabet {0,1,2,3}. Example: a(2)=8 because we have 000,001,002,003,100,200,300 (the other 57=A125145(3) words of length 3 have no subsequences 00). a(n) = Sum_{k=0..n} k*A128235(n+1, k). - Emeric Deutsch, Feb 27 2007
Let P(A) be the power set of an n-element set A. Then a(n) = the sum of the size of the symmetric difference of x and y for every subset {x,y} of P(A). - Ross La Haye, Dec 30 2007 (See the comment from Bernard Schott below.)
Let P(A) be the power set of an n-element set A and B be the Cartesian product of P(A) with itself. Then remove (y,x) from B when (x,y) is in B and x != y and call this R35. Then a(n) = the sum of the size of the symmetric difference of x and y for every (x,y) of R35. [proposed edit of comment just above; by Ross La Haye]
The numbers in this sequence are the Wiener indices of the graphs of n-cubes (Boolean hypercubes). For example, the 3-cube is the graph of the standard cube whose Wiener index is 48. - K.V.Iyer, Feb 26 2009
From Gary W. Adamson, Sep 06 2009: (Start)
Starting (1, 8, 48, ...) = 4th binomial transform of [1, 4, 0, 0, 0, ...].
Equals the sum of terms in 2^n X 2^n semi-magic square arrays in which each row and column is composed of a binomial frequency of terms in the set (1, 3, 5, 7, ...).
The first few such arrays = [1] [1,3; 3,1]; /Q.
[1, 3, 5, 3;
3, 1, 3, 5;
5, 3, 1, 3;
3, 5, 3, 1]
(sum of terms = 48, with a binomial frequency of (1, 2, 1) as to (1, 3, 5) in each row and column)
[1, 3, 5, 3, 5, 7, 5, 3;
3, 1, 3, 5, 7, 5, 3, 5;
5, 3, 1, 3, 5, 3, 5, 7;
3, 5, 3, 1, 3, 5, 7, 5;
5, 7, 5, 3, 1, 3, 5, 3;
7, 5, 3, 5, 3, 1, 3, 5;
5, 3, 5, 7, 5, 3, 1, 3;
3, 5, 7, 5, 3, 5, 3, 1]
(sum of terms = 256, with each row and column composed of one 1, three 3's, three 5's, and one 7)
... (End)
Let P(A) be the power set of an n-element set A and B be the Cartesian product of P(A) with itself. Then a(n) = the sum of the size of the intersection of x and y for every (x,y) of B. - Ross La Haye, Jan 05 2013
Following the last comment of Ross, A002699 is the similar sequence when "intersection" is replaced by "symmetric difference" and A212698 is the similar sequence when "intersection" is replaced by "union". - Bernard Schott, Jan 04 2013
Also, following the first comment of Ross, A082134 is the similar sequence when "symmetric difference" is replaced by "intersection" and A133224 is the similar sequence when "symmetric difference" is replaced by "union". - Bernard Schott, Jan 15 2013
Let [n] denote the set {1,2,3,...,n} and denote an n-permutation of the elements of [n] by p = p(1)p(2)p(3)...p(n), where p(i) is the i-th entry in the linear order given by p. Then (p(i),p(j)) is an inversion of p if i < j but p(i) > p(j). Denote the number of inversions of p by inv(p) and call a 2n-permutation p = p(1)p(2)...p(2n) 2-ordered if p(1) < p(3) < ... < p(2n-1) and p(2) < p(4) < ... < p(2n). Then Sum(inv(p)) = n*4^(n-1), where the sum is taken over all 2-ordered 2n-permutations of p. See Bona reference below. - Ross La Haye, Jan 21 2014
Sum over all peaks of Dyck paths of semilength n of the product of the x and y coordinates. - Alois P. Heinz, May 29 2015
Sum of the number of all edges over all j-dimensional subcubes of the boolean hypercube graph of dimension n, Q_n, for all j, so a(n) = Sum_{j=1..n} binomial(n,j)*2^(n-j) * j*2^(j-1). - Constantinos Kourouzides, Mar 24 2024

Examples

			From _Bernard Schott_, Jan 04 2013: (Start)
See the comment about intersection of X and Y.
If A={b,c}, then in P(A) we have:
{b}Inter{b}={b},
{b}Inter{b,c}={b},
{c}Inter{c}={c},
{c}Inter{b,c}={c},
{b,c}Inter{b}={b},
{b,c}Inter{c}={c},
{b,c}Inter{b,c}={b,c}
and : #{b}+ #{b}+ #{c}+ #{c}+ #{b}+ #{c}+ #{b,c} = 8 = 2*4^(2-1) = a(2).
The other intersections are empty.
(End)
		

References

  • Miklos Bona, Combinatorics of Permutations, Chapman and Hall/CRC, 2004, pp. 1, 43, 64.
  • C. Lanczos, Applied Analysis. Prentice-Hall, Englewood Cliffs, NJ, 1956, p. 516.
  • 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

Programs

Formula

a(n) = n*4^(n-1).
G.f.: x/(1-4x)^2. a(n+1) is the convolution of powers of 4 (A000302). - Wolfdieter Lang, May 16 2003
Third binomial transform of n. E.g.f.: x*exp(4x). - Paul Barry, Jul 22 2003
a(n) = Sum_{k=0..n} k*binomial(2*n, 2*k). - Benoit Cloitre, Jul 30 2003
For n>=0, a(n+1) = Sum_{i+j+k+l=n} binomial(2i, i)*binomial(2j, j)*binomial(2k, k)*binomial(2l, l). - Philippe Deléham, Jan 22 2004
a(n) = Sum_{k=0..n} 4^(n-k)*binomial(n-k+1, k)*binomial(1, (k+1)/2)*(1-(-1)^k)/2. - Paul Barry, Oct 15 2004
Sum_{n>0} 1/a(n) = 8*log(2) - 4*log(3). - Jaume Oliver Lafont, Sep 11 2009
a(0) = 0, a(n) = 4*a(n-1) + 4^(n-1). - Vincenzo Librandi, Dec 31 2010
a(n+1) is the convolution of A000984 with A002457. - Rui Duarte, Oct 08 2011
a(0) = 0, a(1) = 1, a(n) = 8*a(n-1) - 16*a(n-2). - Harvey P. Dale, Jan 18 2012
a(n) = A002699(n)/2 = A212698(n)/3. - Bernard Schott, Jan 04 2013
G.f.: W(0)*x/2 , where W(k) = 1 + 1/( 1 - 4*x*(k+2)/( 4*x*(k+2) + (k+1)/W(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 19 2013
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*log(5/4). - Amiram Eldar, Oct 28 2020
a(n) = (1/2)*Sum_{k=0..n} k*binomial(2*n, k). Compare this with the formula of Benoit Cloitre above. - Wolfdieter Lang, Nov 12 2021
a(n) = (-1)^(n-1)*det(M(n)) for n > 0, where M(n) is the n X n symmetric Toeplitz matrix whose first row consists of 1, 3, ..., 2*n-1. - Stefano Spezia, Aug 04 2022

A236913 Number of partitions of 2n of type EE (see Comments).

Original entry on oeis.org

1, 1, 3, 6, 12, 22, 40, 69, 118, 195, 317, 505, 793, 1224, 1867, 2811, 4186, 6168, 9005, 13026, 18692, 26613, 37619, 52815, 73680, 102162, 140853, 193144, 263490, 357699, 483338, 650196, 870953, 1161916, 1544048, 2044188, 2696627, 3545015, 4644850, 6066425
Offset: 0

Views

Author

Clark Kimberling, Feb 01 2014

Keywords

Comments

The partitions of n are partitioned into four types:
EO, even # of odd parts and odd # of even parts, A236559;
OE, odd # of odd parts and even # of even parts, A160786;
EE, even # of odd parts and even # of even parts, A236913;
OO, odd # of odd parts and odd # of even parts, A236914.
A236559 and A160786 are the bisections of A027193;
A236913 and A236914 are the bisections of A027187.

Examples

			The partitions of 4 of type EE are [3,1], [2,2], [1,1,1,1], so that a(2) = 3.
type/k . 1 .. 2 .. 3 .. 4 .. 5 .. 6 .. 7 .. 8 ... 9 ... 10 .. 11
EO ..... 0 .. 1 .. 0 .. 2 .. 0 .. 5 .. 0 .. 10 .. 0 ... 20 .. 0
OE ..... 1 .. 0 .. 2 .. 0 .. 4 .. 0 .. 8 .. 0 ... 16 .. 0 ... 29
EE ..... 0 .. 1 .. 0 .. 3 .. 0 .. 6 .. 0 .. 12 .. 0 ... 22 .. 0
OO ..... 0 .. 0 .. 1 .. 0 .. 3 .. 0 .. 7 .. 0 ... 14 .. 0 ... 27
From _Gus Wiseman_, Feb 09 2021: (Start)
This sequence counts even-length partitions of even numbers, which have Heinz numbers given by A340784. For example, the a(0) = 1 through a(4) = 12 partitions are:
  ()  (11)  (22)    (33)      (44)
            (31)    (42)      (53)
            (1111)  (51)      (62)
                    (2211)    (71)
                    (3111)    (2222)
                    (111111)  (3221)
                              (3311)
                              (4211)
                              (5111)
                              (221111)
                              (311111)
                              (11111111)
(End)
		

Crossrefs

Note: A-numbers of ranking sequences are in parentheses below.
The ordered version is A000302.
The case of odd-length partitions of odd numbers is A160786 (A340931).
The Heinz numbers of these partitions are (A340784).
A027187 counts partitions of even length/maximum (A028260/A244990).
A034008 counts compositions of even length.
A035363 counts partitions into even parts (A066207).
A047993 counts balanced partitions (A106529).
A058695 counts partitions of odd numbers (A300063).
A058696 counts partitions of even numbers (A300061).
A067661 counts strict partitions of even length (A030229).
A072233 counts partitions by sum and length.
A339846 counts factorizations of even length.
A340601 counts partitions of even rank (A340602).
A340785 counts factorizations into even factors.
A340786 counts even-length factorizations into even factors.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, [1, 0$3],
          `if`(i<1, [0$4], b(n, i-1)+`if`(i>n, [0$4], (p->
          `if`(irem(i, 2)=0, [p[3], p[4], p[1], p[2]],
              [p[2], p[1], p[4], p[3]]))(b(n-i, i)))))
        end:
    a:= n-> b(2*n$2)[1]:
    seq(a(n), n=0..40);  # Alois P. Heinz, Feb 16 2014
  • Mathematica
    z = 25; m1 = Map[Length[Select[Map[{Count[#, True], Count[#, False]} &,  OddQ[IntegerPartitions[2 #]]], EvenQ[(*Odd*)First[#]] && OddQ[(*Even*)Last[#]] &]] &, Range[z]]; m2 = Map[Length[Select[Map[{Count[#, True], Count[#, False]} &,      OddQ[IntegerPartitions[2 # - 1]]], OddQ[(*Odd*)First[#]] && EvenQ[(*Even*)Last[#]] &]] &, Range[z]]; m3 = Map[Length[Select[Map[{Count[#, True], Count[#, False]} &,
    OddQ[IntegerPartitions[2 #]]], EvenQ[(*Odd*)First[#]] && EvenQ[(*Even*)Last[#]] &]] &, Range[z]] ; m4 = Map[Length[Select[Map[{Count[#, True], Count[#, False]} &,
    OddQ[IntegerPartitions[2 # - 1]]], OddQ[(*Odd*)First[#]] && OddQ[(*Even*)Last[#]] &]] &, Range[z]];
    m1 (* A236559, type EO*)
    m2 (* A160786, type OE*)
    m3 (* A236913, type EE*)
    m4 (* A236914, type OO*)
    (* Peter J. C. Moses, Feb 03 2014 *)
    b[n_, i_] := b[n, i] = If[n == 0, {1, 0, 0, 0}, If[i < 1, {0, 0, 0, 0}, b[n, i - 1] + If[i > n, {0, 0, 0, 0}, Function[p, If[Mod[i, 2] == 0, p[[{3, 4, 1, 2}]], p[[{2, 1, 4, 3}]]]][b[n - i, i]]]]]; a[n_] := b[2*n, 2*n][[1]]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Oct 27 2015, after Alois P. Heinz *)
    Table[Length[Select[IntegerPartitions[2n],EvenQ[Length[#]]&]],{n,0,15}] (* Gus Wiseman, Feb 09 2021 *)

Extensions

More terms from Alois P. Heinz, Feb 16 2014

A345197 Concatenation of square matrices A(n), each read by rows, where A(n)(k,i) is the number of compositions of n of length k with alternating sum i, where 1 <= k <= n, and i ranges from -n + 2 to n in steps of 2.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 03 2021

Keywords

Comments

The alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i.

Examples

			The matrices for n = 1..7:
  1   0 1   0 0 1   0 0 0 1   0 0 0 0 1   0 0 0 0 0 1   0 0 0 0 0 0 1
      1 0   1 1 0   1 1 1 0   1 1 1 1 0   1 1 1 1 1 0   1 1 1 1 1 1 0
            0 1 0   0 1 2 0   0 1 2 3 0   0 1 2 3 4 0   0 1 2 3 4 5 0
                    0 1 0 0   0 2 2 0 0   0 3 4 3 0 0   0 4 6 6 4 0 0
                              0 0 1 0 0   0 0 2 3 0 0   0 0 3 6 6 0 0
                                          0 0 1 0 0 0   0 0 3 3 0 0 0
                                                        0 0 0 1 0 0 0
Matrix n = 5 counts the following compositions:
           i=-3:        i=-1:          i=1:            i=3:        i=5:
        -----------------------------------------------------------------
   k=1: |    0            0             0               0          (5)
   k=2: |   (14)         (23)          (32)            (41)         0
   k=3: |    0          (131)       (221)(122)   (311)(113)(212)    0
   k=4: |    0       (1211)(1112)  (2111)(1121)         0           0
   k=5: |    0            0          (11111)            0           0
		

Crossrefs

The number of nonzero terms in each matrix appears to be A000096.
The number of zeros in each matrix appears to be A000124.
Row sums and column sums both appear to be A007318 (Pascal's triangle).
The matrix sums are A131577.
Antidiagonal sums appear to be A163493.
The reverse-alternating version is also A345197 (this sequence).
Antidiagonals are A345907.
Traces are A345908.
A000041 counts partitions of 2n with alternating sum 0, ranked by A000290.
A011782 counts compositions.
A097805 counts compositions by alternating (or reverse-alternating) sum.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A316524 gives the alternating sum of prime indices (reverse: A344616).
A344610 counts partitions by sum and positive reverse-alternating sum.
A344611 counts partitions of 2n with reverse-alternating sum >= 0.
Other tetrangles: A318393, A318816, A320808, A321912.
Compositions of n, 2n, or 2n+1 with alternating/reverse-alternating sum k:
- k = 0: counted by A088218, ranked by A344619/A344619.
- k = 1: counted by A000984, ranked by A345909/A345911.
- k = -1: counted by A001791, ranked by A345910/A345912.
- k = 2: counted by A088218, ranked by A345925/A345922.
- k = -2: counted by A002054, ranked by A345924/A345923.
- k >= 0: counted by A116406, ranked by A345913/A345914.
- k <= 0: counted by A058622(n-1), ranked by A345915/A345916.
- k > 0: counted by A027306, ranked by A345917/A345918.
- k < 0: counted by A294175, ranked by A345919/A345920.
- k != 0: counted by A058622, ranked by A345921/A345921.
- k even: counted by A081294, ranked by A053754/A053754.
- k odd: counted by A000302, ranked by A053738/A053738.

Programs

  • Mathematica
    ats[y_]:=Sum[(-1)^(i-1)*y[[i]],{i,Length[y]}];
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],Length[#]==k&&ats[#]==i&]],{n,0,6},{k,1,n},{i,-n+2,n,2}]

A053754 If k is in the sequence then 2*k and 2*k+1 are not (and 0 is in the sequence); when written in binary k has an even number of bits (0 has 0 digits).

Original entry on oeis.org

0, 2, 3, 8, 9, 10, 11, 12, 13, 14, 15, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148
Offset: 1

Views

Author

Henry Bottomley, Apr 06 2000

Keywords

Comments

Runs of successive terms with same number of bits have length twice powers of 4 (A081294). [Clarified by Michel Marcus, Oct 21 2020]
The sequence A081294 counts compositions of even numbers - Gus Wiseman, Aug 12 2021
A031443 is a subsequence; A179888 is the intersection of this sequence and A032925. - Reinhard Zumkeller, Jul 31 2010
The lower and upper asymptotic densities of this sequence are 1/3 and 2/3, respectively. - Amiram Eldar, Feb 01 2021
From Gus Wiseman, Aug 10 2021: (Start)
Also numbers k such that the k-th composition in standard order (row k of A066099) has even sum. The terms and corresponding compositions begin:
0: () 2: (2) 8: (4)
3: (1,1) 9: (3,1)
10: (2,2)
11: (2,1,1)
12: (1,3)
13: (1,2,1)
14: (1,1,2)
15: (1,1,1,1)
The following pertain to compositions in standard order: A000120, A029837, A070939, A066099, A124767.
(End)

Crossrefs

Positions of even terms in A029837 with offset 0.
The complement (the odd version) is A053738, counted by A000302.
The version for Heinz numbers of partitions is A300061, counted by A058696.

Programs

  • Haskell
    a053754 n = a053754_list !! (n-1)
    a053754_list = 0 : filter (even . a070939) [1..]
    -- Reinhard Zumkeller, Apr 18 2015
    
  • Mathematica
    Select[Range[0, 150], EvenQ @ IntegerLength[#, 2] &] (* Amiram Eldar, Feb 01 2021 *)
  • PARI
    lista(nn) = {my(va = vector(nn)); for (n=2, nn, my(k=va[n-1]+1); while (#select(x->(x==k\2), va), k++); va[n] = k;); va;} \\ Michel Marcus, Oct 20 2020
    
  • PARI
    a(n) = n-1 + (1<Kevin Ryde, Apr 30 2021

Extensions

Offset corrected by Reinhard Zumkeller, Jul 30 2010

A002001 a(n) = 3*4^(n-1), n>0; a(0)=1.

Original entry on oeis.org

1, 3, 12, 48, 192, 768, 3072, 12288, 49152, 196608, 786432, 3145728, 12582912, 50331648, 201326592, 805306368, 3221225472, 12884901888, 51539607552, 206158430208, 824633720832, 3298534883328, 13194139533312, 52776558133248, 211106232532992, 844424930131968
Offset: 0

Views

Author

N. J. A. Sloane, Dec 11 1996

Keywords

Comments

Second binomial transform of (1,1,4,4,16,16,...) = (3*2^n+(-2)^n)/4. - Paul Barry, Jul 16 2003
Number of vertices (or sides) formed after the (n-1)-th iterate towards building a Koch's snowflake. - Lekraj Beedassy, Jan 24 2005
For n >= 1, a(n) is the number of functions f:{1,2,...,n}->{1,2,3,4} such that for a fixed x in {1,2,...,n} and a fixed y in {1,2,3,4} we have f(x) <> y. - Aleksandar M. Janjic and Milan Janjic, Mar 27 2007
a(n) = (n+1) terms in the sequence (1, 2, 3, 3, 3, ...) dot (n+1) terms in the sequence (1, 1, 3, 12, 48, ...). Example: a(4) = 192 = (1, 2, 3, 3, 3) dot (1, 1, 3, 12, 48) = (1 + 2 + 9 + 36 + 144). - Gary W. Adamson, Aug 03 2010
a(n) is the number of compositions of n when there are 3 types of each natural number. - Milan Janjic, Aug 13 2010
See A178789 for the number of acute (= exterior) angles of the Koch snowflake referred to in the above comment by L. Beedassy. - M. F. Hasler, Dec 17 2013
After 1, subsequence of A033428. - Vincenzo Librandi, May 26 2014
a(n) counts walks (closed) on the graph G(1-vertex; 1-loop x3, 2-loop x3, 3-loop x3, 4-loop x3, ...). - David Neil McGrath, Jan 01 2015
For n > 1, a(n) are numbers k such that (2^(k-1) mod k)/(2^k mod k) = 2; 2^(a(n)-1) == 2^(2n-1) (mod a(n)) and 2^a(n) == 2^(2n-2) (mod a(n)). - Thomas Ordowski, Apr 22 2020
For n > 1, a(n) is the number of 4-colorings of the Hex graph of size 2 X (n-1). More generally, for q > 2, the number of q-colorings of the Hex graph of size 2 X n is given by q*(q - 1)*(q - 2)^(2*n - 2). - Sela Fried, Sep 25 2023
For n > 1, a(n) is the number of pixels in the HEALPix discretization of the sphere of order n-2; HEALPix is a common sphere pixellization scheme in astronomy, cosmology, and nuclear engineering. - Jayson R. Vavrek, Aug 08 2024

Crossrefs

First difference of 4^n (A000302).

Programs

Formula

From Paul Barry, Apr 20 2003: (Start)
a(n) = (3*4^n + 0^n)/4 (with 0^0=1).
E.g.f.: (3*exp(4*x) + 1)/4. (End)
With interpolated zeros, this has e.g.f. (3*cosh(2*x) + 1)/4 and binomial transform A006342. - Paul Barry, Sep 03 2003
a(n) = Sum_{j=0..1} Sum_{k=0..n} C(2n+j, 2k). - Paul Barry, Nov 29 2003
G.f.: (1-x)/(1-4*x). The sequence 1, 3, -12, 48, -192, ... has g.f. (1+7*x)/(1+4*x). - Paul Barry, Feb 12 2004
a(n) = 3*Sum_{k=0..n-1} a(k). - Adi Dani, Jun 24 2011
G.f.: 1/(1-3*Sum_{k>=1} x^k). - Joerg Arndt, Jun 24 2011
Row sums of triangle A134316. - Gary W. Adamson, Oct 19 2007
a(n) = A011782(n) * A003945(n). - R. J. Mathar, Jul 08 2009
If p(1)=3 and p(i)=3 for i > 1, and if A is the Hessenberg matrix of order n defined by A(i,j) = p(j-i+1) when i <= j, A(i,j)=-1 when i=j+1, and A(i,j) = 0 otherwise, then, for n >= 1, a(n) = det A. - Milan Janjic, Apr 29 2010
a(n) = 4*a(n-1), a(0)=1, a(1)=3. - Vincenzo Librandi, Dec 31 2010
G.f.: 1 - G(0) where G(k) = 1 - 1/(1-3*x)/(1-x/(x-1/G(k+1))); (recursively defined continued fraction). - Sergei N. Gladkovskii, Jan 25 2013
G.f.: x+2*x/(G(0)-2), where G(k) = 1 + 1/(1 - x*(3*k+1)/(x*(3*k+4) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 26 2013
a(n) = ceiling(3*4^(n-1)). - Wesley Ivan Hurt, Dec 17 2013
Construct the power matrix T(n,j) = [A(n)^*j]*[S(n)^*(j-1)] where A(n)=(3,3,3,...) and S(n)=(0,1,0,0,...). (* is convolution operation.) Then T(n,j) counts n-walks containing j loops on the single vertex graph above and a(n) = Sum_{j=1..n} T(n,j). (S(n)^*0=I.) - David Neil McGrath, Jan 01 2015
Previous Showing 51-60 of 587 results. Next