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 91-100 of 587 results. Next

A345918 Numbers k such that the k-th composition in standard order (row k of A066099) has reverse-alternating sum > 0.

Original entry on oeis.org

1, 2, 4, 6, 7, 8, 11, 12, 14, 16, 19, 20, 21, 22, 24, 26, 27, 28, 30, 31, 32, 35, 37, 38, 40, 42, 44, 47, 48, 51, 52, 54, 56, 59, 60, 62, 64, 67, 69, 70, 72, 73, 74, 76, 79, 80, 82, 83, 84, 86, 87, 88, 91, 92, 93, 94, 96, 99, 100, 101, 102, 104, 106, 107, 108
Offset: 1

Views

Author

Gus Wiseman, Jul 09 2021

Keywords

Comments

The reverse-alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(k-i) y_i.
The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

Examples

			The initial terms and the corresponding compositions:
     1: (1)        26: (1,2,2)        52: (1,2,3)
     2: (2)        27: (1,2,1,1)      54: (1,2,1,2)
     4: (3)        28: (1,1,3)        56: (1,1,4)
     6: (1,2)      30: (1,1,1,2)      59: (1,1,2,1,1)
     7: (1,1,1)    31: (1,1,1,1,1)    60: (1,1,1,3)
     8: (4)        32: (6)            62: (1,1,1,1,2)
    11: (2,1,1)    35: (4,1,1)        64: (7)
    12: (1,3)      37: (3,2,1)        67: (5,1,1)
    14: (1,1,2)    38: (3,1,2)        69: (4,2,1)
    16: (5)        40: (2,4)          70: (4,1,2)
    19: (3,1,1)    42: (2,2,2)        72: (3,4)
    20: (2,3)      44: (2,1,3)        73: (3,3,1)
    21: (2,2,1)    47: (2,1,1,1,1)    74: (3,2,2)
    22: (2,1,2)    48: (1,5)          76: (3,1,3)
    24: (1,4)      51: (1,3,1,1)      79: (3,1,1,1,1)
		

Crossrefs

The version for prime indices is A000037.
The version for Heinz numbers of partitions is A026424, counted by A027193.
These compositions are counted by A027306.
These are the positions of terms > 0 in A344618.
The weak (k >= 0) version is A345914.
The version for unreversed alternating sum is A345917.
The opposite (k < 0) version is A345920.
A011782 counts compositions.
A097805 counts compositions by alternating (or reverse-alternating) sum.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A236913 counts partitions of 2n with reverse-alternating sum <= 0.
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.
A345197 counts compositions by sum, length, and alternating sum.
Standard compositions: A000120, A066099, A070939, A228351, A124754, A344618.
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
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    sats[y_]:=Sum[(-1)^(i-Length[y])*y[[i]],{i,Length[y]}];
    Select[Range[0,100],sats[stc[#]]>0&]

A159991 Powers of 60: a(n) = 60^n.

Original entry on oeis.org

1, 60, 3600, 216000, 12960000, 777600000, 46656000000, 2799360000000, 167961600000000, 10077696000000000, 604661760000000000, 36279705600000000000, 2176782336000000000000, 130606940160000000000000, 7836416409600000000000000, 470184984576000000000000000
Offset: 0

Views

Author

Reinhard Zumkeller, May 01 2009

Keywords

Examples

			G.f. = 1 + 60*x + 3600*x^2 + 216000*x^3 + 12960000*x^4 + 77600000*x^5 + ... - _Michael Somos_, Jan 01 2019
		

Crossrefs

Programs

Formula

a(n) = A000400(n)*A011557(n) = A000351(n)*A001021(n) = A000302(n)*A001024(n) = A000244(n)*A009964(n). (Corrected by Robert B Fowler, Jan 25 2023)
From Muniru A Asiru, Nov 21 2018: (Start)
a(n) = 60^n.
a(n) = 60*a(n-1) for n > 0, a(0) = 1.
G.f.: 1/(1-60*x).
E.g.f: exp(60*x). (End)
a(n) = 1/a(-n) for all n in Z. - Michael Somos, Jan 01 2019

A345920 Numbers k such that the k-th composition in standard order (row k of A066099) has reverse-alternating sum < 0.

Original entry on oeis.org

5, 9, 17, 18, 23, 25, 29, 33, 34, 39, 45, 49, 57, 65, 66, 68, 71, 75, 77, 78, 81, 85, 89, 90, 95, 97, 98, 103, 105, 109, 113, 114, 119, 121, 125, 129, 130, 132, 135, 139, 141, 142, 149, 153, 154, 159, 161, 169, 177, 178, 183, 189, 193, 194, 199, 205, 209, 217
Offset: 1

Views

Author

Gus Wiseman, Jul 09 2021

Keywords

Comments

The reverse-alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(k-i) y_i.
The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

Examples

			The initial terms and the corresponding compositions:
      5: (2,1)         68: (4,3)
      9: (3,1)         71: (4,1,1,1)
     17: (4,1)         75: (3,2,1,1)
     18: (3,2)         77: (3,1,2,1)
     23: (2,1,1,1)     78: (3,1,1,2)
     25: (1,3,1)       81: (2,4,1)
     29: (1,1,2,1)     85: (2,2,2,1)
     33: (5,1)         89: (2,1,3,1)
     34: (4,2)         90: (2,1,2,2)
     39: (3,1,1,1)     95: (2,1,1,1,1,1)
     45: (2,1,2,1)     97: (1,5,1)
     49: (1,4,1)       98: (1,4,2)
     57: (1,1,3,1)    103: (1,3,1,1,1)
     65: (6,1)        105: (1,2,3,1)
     66: (5,2)        109: (1,2,1,2,1)
		

Crossrefs

The version for prime indices is {}.
The version for Heinz numbers of partitions is A119899.
These compositions are counted by A294175 (even bisection: A008549).
These are the positions of terms < 0 in A344618.
The complement is A345914.
The weak (k <= 0) version is A345916.
The opposite (k > 0) version is A345918.
The version for unreversed alternating sum is A345919.
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).
A236913 counts partitions of 2n with reverse-alternating sum <= 0.
A316524 gives the alternating sum of prime indices (reverse: A344616).
A345197 counts compositions by sum, length, and alternating sum.
Standard compositions: A000120, A066099, A070939, A228351, A124754, A344618.
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
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    sats[y_]:=Sum[(-1)^(i-Length[y])*y[[i]],{i,Length[y]}];
    Select[Range[0,100],sats[stc[#]]<0&]

A345921 Numbers k such that the k-th composition in standard order (row k of A066099) has alternating sum != 0.

Original entry on oeis.org

1, 2, 4, 5, 6, 7, 8, 9, 11, 12, 14, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 37, 38, 39, 40, 42, 44, 45, 47, 48, 49, 51, 52, 54, 56, 57, 59, 60, 62, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81
Offset: 1

Views

Author

Gus Wiseman, Jul 10 2021

Keywords

Comments

The alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i.
The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.
Also numbers k such that the k-th composition in standard order has reverse-alternating sum != 0.

Examples

			The initial terms and the corresponding compositions:
     1: (1)        20: (2,3)          35: (4,1,1)
     2: (2)        21: (2,2,1)        37: (3,2,1)
     4: (3)        22: (2,1,2)        38: (3,1,2)
     5: (2,1)      23: (2,1,1,1)      39: (3,1,1,1)
     6: (1,2)      24: (1,4)          40: (2,4)
     7: (1,1,1)    25: (1,3,1)        42: (2,2,2)
     8: (4)        26: (1,2,2)        44: (2,1,3)
     9: (3,1)      27: (1,2,1,1)      45: (2,1,2,1)
    11: (2,1,1)    28: (1,1,3)        47: (2,1,1,1,1)
    12: (1,3)      29: (1,1,2,1)      48: (1,5)
    14: (1,1,2)    30: (1,1,1,2)      49: (1,4,1)
    16: (5)        31: (1,1,1,1,1)    51: (1,3,1,1)
    17: (4,1)      32: (6)            52: (1,2,3)
    18: (3,2)      33: (5,1)          54: (1,2,1,2)
    19: (3,1,1)    34: (4,2)          56: (1,1,4)
		

Crossrefs

The version for Heinz numbers of partitions is A000037.
These compositions are counted by A058622.
These are the positions of terms != 0 in A124754.
The complement (k = 0) is A344619.
The positive (k > 0) version is A345917 (reverse: A345918).
The negative (k < 0) version is A345919 (reverse: A345920).
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).
A345197 counts compositions by sum, length, and alternating sum.
Standard compositions: A000120, A066099, A070939, A228351, A124754, A344618.
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
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    ats[y_]:=Sum[(-1)^(i-1)*y[[i]],{i,Length[y]}];
    Select[Range[0,100],ats[stc[#]]!=0&]

A002063 a(n) = 9*4^n.

Original entry on oeis.org

9, 36, 144, 576, 2304, 9216, 36864, 147456, 589824, 2359296, 9437184, 37748736, 150994944, 603979776, 2415919104, 9663676416, 38654705664, 154618822656, 618475290624, 2473901162496, 9895604649984, 39582418599936, 158329674399744, 633318697598976
Offset: 0

Views

Author

Keywords

Comments

a(n) is twice the area of the trapezoid created by the four points (2^n,2^(n+1)), (2^(n+1), 2^n), (2^(n+1), 2^(n+2)), (2^(n+2), 2^(n+1)). - J. M. Bergot, May 23 2014
These are squares that can be expressed as sum of exactly two distinct powers of two. For instance, a(4) = 9*4^4 = 2304 = 2^11 + 2^8 . It is conjectured that these are the only squares with this characteristic (tested on squares up to (10^7)^2). - Andres Cicuttin, Apr 23 2016
Conjecture is true. It is equivalent to prove that the Diophantine equation m^2 = 2^k*(1+2^h), where h>0, has solutions only when h=3. Dividing by 2^k we must obtain an odd square on the left, since 1+2^h is odd, so we can write (2*r+1)^2 = 1+2^h. Expanding, we have 4*r*(r+1) = 2^h, from which it follows that r must be equal to 1 and thus h=3, since r and r+1 must be powers of 2. - Giovanni Resta, Jul 27 2017

Crossrefs

Essentially the same as A055841. First differences of A002001.
Cf. A000302.

Programs

Formula

From Philippe Deléham, Nov 23 2008: (Start)
a(n) = 4*a(n-1), n > 0; a(0)=9.
G.f.: 9/(1-4*x). (End)
a(n) = 9*A000302(n). - Michel Marcus, Apr 23 2016
E.g.f.: 9*exp(4*x). - Ilya Gutkovskiy, Apr 23 2016
a(n) = 2^(2*n+3) + 2^(2*n). - Andres Cicuttin, Apr 26 2016
a(n) = A004171(n+1) + A000302(n). - Zhandos Mambetaliyev, Nov 19 2016

A047653 Constant term in expansion of (1/2) * Product_{k=-n..n} (1 + x^k).

Original entry on oeis.org

1, 2, 4, 10, 26, 76, 236, 760, 2522, 8556, 29504, 103130, 364548, 1300820, 4679472, 16952162, 61790442, 226451036, 833918840, 3084255128, 11451630044, 42669225172, 159497648600, 597950875256, 2247724108772, 8470205600640, 31991616634296, 121086752349064
Offset: 0

Views

Author

Keywords

Comments

Or, constant term in expansion of Product_{k=1..n} (x^k + 1/x^k)^2. - N. J. A. Sloane, Jul 09 2008
Or, maximal coefficient of the polynomial (1+x)^2 * (1+x^2)^2 *...* (1+x^n)^2.
a(n) = A000302(n) - A181765(n).
From Gus Wiseman, Apr 18 2023: (Start)
Also the number of subsets of {1..2n} that are empty or have mean n. The a(0) = 1 through a(3) = 10 subsets are:
{} {} {} {}
{1} {2} {3}
{1,3} {1,5}
{1,2,3} {2,4}
{1,2,6}
{1,3,5}
{2,3,4}
{1,2,3,6}
{1,2,4,5}
{1,2,3,4,5}
Also the number of subsets of {-n..n} with no 0's but with sum 0. The a(0) = 1 through a(3) = 10 subsets are:
{} {} {} {}
{-1,1} {-1,1} {-1,1}
{-2,2} {-2,2}
{-2,-1,1,2} {-3,3}
{-3,1,2}
{-2,-1,3}
{-2,-1,1,2}
{-3,-1,1,3}
{-3,-2,2,3}
{-3,-2,-1,1,2,3}
(End)

Crossrefs

Cf. A025591.
Cf. A053632; variant: A127728.
For median instead of mean we have A079309(n) + 1.
Odd bisection of A133406.
A000980 counts nonempty subsets of {1..2n-1} with mean n.
A007318 counts subsets by length, A327481 by mean.

Programs

  • Maple
    f:=n->coeff( expand( mul((x^k+1/x^k)^2,k=1..n) ),x,0);
    # second Maple program:
    b:= proc(n, i) option remember; `if`(n>i*(i+1)/2, 0,
          `if`(i=0, 1, 2*b(n, i-1)+b(n+i, i-1)+b(abs(n-i), i-1)))
        end:
    a:=n-> b(0, n):
    seq(a(n), n=0..40);  # Alois P. Heinz, Mar 10 2014
  • Mathematica
    b[n_, i_] := b[n, i] = If[n>i*(i+1)/2, 0, If[i == 0, 1, 2*b[n, i-1]+b[n+i, i-1]+b[Abs[n-i], i-1]]]; a[n_] := b[0, n]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Mar 10 2014, after Alois P. Heinz *)
    nmax = 26; d = {1}; a1 = {};
    Do[
      i = Ceiling[Length[d]/2];
      AppendTo[a1, If[i > Length[d], 0, d[[i]]]];
      d = PadLeft[d, Length[d] + 2 n] + PadRight[d, Length[d] + 2 n] +
        2 PadLeft[PadRight[d, Length[d] + n], Length[d] + 2 n];
    , {n, nmax}];
    a1 (* Ray Chandler, Mar 15 2014 *)
    Table[Length[Select[Subsets[Range[2n]],Length[#]==0||Mean[#]==n&]],{n,0,6}] (* Gus Wiseman, Apr 18 2023 *)
  • PARI
    a(n)=polcoeff(prod(k=-n,n,1+x^k),0)/2
    
  • PARI
    {a(n)=sum(k=0,n*(n+1)/2,polcoeff(prod(m=1,n,1+x^m+x*O(x^k)),k)^2)} \\ Paul D. Hanna, Nov 30 2010

Formula

Sum of squares of coefficients in Product_{k=1..n} (1+x^k):
a(n) = Sum_{k=0..n(n+1)/2} A053632(n,k)^2. - Paul D. Hanna, Nov 30 2010
a(n) = A000980(n)/2.
a(n) ~ sqrt(3) * 4^n / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Sep 11 2014
From Gus Wiseman, Apr 18 2023 (Start)
a(n) = A133406(2n+1).
a(n) = A212352(n) + 1.
a(n) = A362046(2n) + 1.
(End)

Extensions

More terms from Michael Somos, Jun 10 2000

A056450 a(n) = (3*2^n - (-2)^n)/2.

Original entry on oeis.org

1, 4, 4, 16, 16, 64, 64, 256, 256, 1024, 1024, 4096, 4096, 16384, 16384, 65536, 65536, 262144, 262144, 1048576, 1048576, 4194304, 4194304, 16777216, 16777216, 67108864, 67108864, 268435456, 268435456, 1073741824, 1073741824, 4294967296
Offset: 0

Views

Author

Keywords

Comments

Number of palindromes of length n using a maximum of four different symbols.
Number of achiral rows of n colors using up to four colors. - Robert A. Russell, Nov 09 2018
Interleaving of A000302 and 4*A000302.
Unsigned version of A141125.
Binomial transform is A164907. Second binomial transform is A164908. Third binomial transform is A057651. Fourth binomial transform is A016129.

Examples

			At length n=1 there are a(1)=4 palindromes, A, B, C, D.
At length n=2, there are a(2)=4 palindromes, AA, BB, CC, DD.
At length n=3, there are a(3)=16 palindromes, AAA, BBB, CCC, DDD, ABA, BAB, ... , CDC, DCD.
		

References

  • 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]

Crossrefs

Column k=4 of A321391.
Cf. A016116.
Essentially the same as A213173.
Cf. A000302 (oriented), A032121 (unoriented), A032087(n>1) (chiral).

Programs

  • Magma
    [ (3*2^n-(-2)^n)/2: n in [0..31] ];
    
  • Magma
    [4^Floor((n+1)/2): n in [0..40]]; // Vincenzo Librandi, Aug 16 2011
    
  • Mathematica
    Table[4^Ceiling[n/2], {n,0,40}] (* or *)
    CoefficientList[Series[(1 + 4 x)/((1 + 2 x) (1 - 2 x)), {x, 0, 31}], x] (* or *)
    LinearRecurrence[{0, 4}, {1, 4}, 40] (* Robert A. Russell, Nov 07 2018 *)
  • PARI
    a(n)=4^((n+1)\2) \\ Charles R Greathouse IV, Apr 08 2012
    
  • PARI
    a(n)=(3*2^n-(-2)^n)/2 \\ Charles R Greathouse IV, Oct 03 2016

Formula

a(n) = 4^floor((n+1)/2).
a(n) = 4*a(n-2) for n > 1; a(0) = 1, a(1) = 4.
G.f.: (1+4*x) / (1-4*x^2). - R. J. Mathar, Jan 19 2011 [Adapted to offset 0 by Robert A. Russell, Nov 07 2018]
a(n+3) = a(n+2)*a(n+1)/a(n). - Reinhard Zumkeller, Mar 04 2011
a(n) = 4*abs(A164111(n-1)). - R. J. Mathar, Jan 19 2011
a(n) = C(4,0)*A000007(n) + C(4,1)*A057427(n) + C(4,2)*A056453(n) + C(4,3)*A056454(n) + C(4,4)*A056455(n). - Robert A. Russell, Nov 08 2018

Extensions

a(0)=1 prepended by Robert A. Russell, Nov 07 2018
Edited by N. J. A. Sloane, Sep 29 2019

A091042 Triangle of even numbered entries of odd numbered rows of Pascal's triangle A007318.

Original entry on oeis.org

1, 1, 3, 1, 10, 5, 1, 21, 35, 7, 1, 36, 126, 84, 9, 1, 55, 330, 462, 165, 11, 1, 78, 715, 1716, 1287, 286, 13, 1, 105, 1365, 5005, 6435, 3003, 455, 15, 1, 136, 2380, 12376, 24310, 19448, 6188, 680, 17, 1, 171, 3876, 27132, 75582, 92378, 50388, 11628, 969, 19, 1, 210, 5985, 54264, 203490, 352716, 293930, 116280, 20349, 1330, 21
Offset: 0

Views

Author

Wolfdieter Lang, Jan 23 2004

Keywords

Comments

The row polynomials Pe(n, x) := Sum_{m=0..n} a(n, m)*x^m appear as numerators of the generating functions for the even numbered column sequences of array A034870.
Elements have the same parity as those of Pascal's triangle.
All zeros of polynomial Pe(n, x) are negative. They are -tan^2(Pi/2*n+1), -tan^2(2*Pi/2*n+1), ..., -tan^2(n*Pi/2*n+1). Moreover, for m >= 1, Pe(m, -x^2) is the characteristic polynomial of the linear difference equation with constant coefficients for differences between multiples of 2*m+1 with even and odd digit sum in base 2*m in the interval [0,(2*m)^n). - Vladimir Shevelev and Peter J. C. Moses, May 22 2012
Row reverse of A103327. - Peter Bala, Jul 29 2013
The row polynomial Pe(d, x), multiplied by (2*d)!/d! = A001813(d), gives the numerator polynomial of the o.g.f. of the sequence of the diagonal d, for d >= 0, of the Sheffer triangle Lah[4,1] given in A048854. - Wolfdieter Lang, Oct 12 2017

Examples

			Triangle a(n, m) begins:
n\m  0   1    2     3      4      5      6      7     8    9  10 ...
0:   1
1:   1   3
2:   1  10    5
3:   1  21   35     7
4:   1  36  126    84      9
5:   1  55  330   462    165     11
6:   1  78  715  1716   1287    286     13
7:   1 105 1365  5005   6435   3003    455     15
8:   1 136 2380 12376  24310  19448   6188    680    17
9:   1 171 3876 27132  75582  92378  50388  11628   969   19
10:  1 210 5985 54264 203490 352716 293930 116280 20349 1330  21
... reformatted and extended. - _Wolfdieter Lang_, Oct 12 2017
From _Peter Bala_, Jan 30 2022: (Start)
(1/2)*(N^2 + N) = Sum_{j = 1..N} j.
(1/2)*(N^2 + N)^3 = Sum_{j = 1..N} j^3 + 3*Sum_{j = 1..N} j^5.
(1/2)*(N^2 + N)^5 = Sum_{j = 1..N} j^5 + 10*Sum_{j = 1..N} j^7 + 5*Sum_{j = 1..N} j^9.
(1/2)*(N^2 + N)^7 = Sum_{j = 1..N} j^7 + 21*Sum_{j = 1..N} j^9 + 35*Sum_{j = 1..N} j^11 + 7*Sum_{j = 1..N} j^13. (End)
		

References

  • A. M. Yaglom and I. M. Yaglom, An elementary proof of the Wallis, Leibniz and Euler formulas for pi. Uspekhi Matem. Nauk, VIII (1953), 181-187(in Russian).

Crossrefs

Cf. A212500, A038754. A000302 (row sums), A085478, A103327 (row reverse), A048854, A103328.

Programs

  • GAP
    Flat(List([0..12], n-> List([0..n], k-> Binomial(2*n+1, 2*k) ))); # G. C. Greubel, Aug 01 2019
    
  • Magma
    [[Binomial(2*n+1,2*k): k in [0..n]]: n in [0..12]]; // G. C. Greubel, Aug 01 2019
    
  • Maple
    f := (x, t) -> cosh(sqrt(x)*t)*sinh(t); seq(seq(coeff(((2*n+1)!*coeff(series(f(x,t),t,2*n+2),t,2*n+1)),x,k),k=0..n),n=0..9); # Peter Luschny, Jul 29 2013
  • Mathematica
    T[n_, k_] /; 0 <= k <= n := T[n, k] = 2T[n-1, k] + 2T[n-1, k-1] + 2T[n-2, k-1] - T[n-2, k] - T[n-2, k-2]; T[0, 0] = T[1, 0] = 1; T[1, 1] = 3; T[, ] = 0;
    Table[T[n, k], {n, 0, 10}, {k, 0, n}]//Flatten (* Jean-François Alcover, Jul 29 2018, after Philippe Deléham *)
    Table[Binomial[2*n+1, 2*k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Aug 01 2019 *)
  • PARI
    T(n,k) = binomial(2*n+1, 2*k); \\ G. C. Greubel, Aug 01 2019
    
  • Python
    from math import comb, isqrt
    def A091042(n): return comb((r:=(m:=isqrt(k:=n+1<<1))-(k<=m*(m+1)))<<1|1,n-comb(r+1,2)<<1) # Chai Wah Wu, Apr 30 2025
  • Sage
    [[binomial(2*n+1, 2*k) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Aug 01 2019
    

Formula

T(n, m) = binomial(2*n+1, 2*m) = A007318(2*n+1, 2*m), n >= m >= 0, otherwise 0.
From Peter Bala, Jul 29 2013: (Start)
E.g.f.: sinh(t)*cosh(sqrt(x)*t) = t + (1 + 3*x)*t^3/3! + (1 + 10*x + 5*x^2)*t^5/5! + (1 + 21*x + 35*x^2 + 7*x^3)*t^7/7! + ....
O.g.f.: A(x,t) = (1 + (x - 1)*t)/( (1 + (x - 1)*t)^2 - 4*t*x ) = 1 + (1 + 3*x)*t + (1 + 10*x + 5*x^2)*t^2 + ...
The function A( x/(x + 4), t*(x + 4)/4 ) = 1 + (1 + x)*t + (1 + 3*x + x^2)*t^2 + ... is the o.g.f. for A085478.
O.g.f. for n-th diagonal: ( Sum_{k = 0..n} binomial(2*n,2*k)*x^k )/(1 - x)^(2*n).
n-th row polynomial R(n,x) = (1/2)*( (1 + sqrt(x))^(2*n+1) - (sqrt(x) - 1)^(2*n+1) ).
Row sums A000302. (End)
T(n, k) = 2*T(n-1,k) + 2*T(n-1,k-1) + 2*T(n-2,k-1) - T(n-2,k) - T(n-2,k-2) with T(0,0)=T(1,0)=1, T(1,1)=3, T(n,k)=0 if k < 0 or if k > n. - Philippe Deléham, Nov 26 2013
From Peter Bala, Jan 31 2022: (Start)
Define S(r,N) = Sum_{j = 1..N} j^r. Then the following identity holds for n >= 0: (1/2)*(N^2 + N)^(2*n+1) = T(n,0)*S(2*n+1,N) + T(n,1)*S(2*n+3,N) + ... + T(n,n)* S(4*n+1,N). Some examples are given below. (End)

A097073 Expansion of (1-x+2*x^2)/((1+x)*(1-2*x)).

Original entry on oeis.org

1, 0, 4, 4, 12, 20, 44, 84, 172, 340, 684, 1364, 2732, 5460, 10924, 21844, 43692, 87380, 174764, 349524, 699052, 1398100, 2796204, 5592404, 11184812, 22369620, 44739244, 89478484, 178956972, 357913940, 715827884, 1431655764, 2863311532
Offset: 0

Views

Author

Paul Barry, Jul 22 2004

Keywords

Comments

Partial sums are A097074.
Pairwise sums are {1, 1, 4, 16, 32, ...} or 2^n -Sum_{k=0..n} binomial(n,k)*(-1)^(n+k)*k.

Crossrefs

Cf. A001045, A078008 (form a(n)=2^n-a(n-1)).

Programs

Formula

a(n) = (2*2^n + 4*(-1)^n)/3 - 0^n.
a(n) = A001045(n+1) + (-1)^n - 0^n.
a(n) = 2*A078008(n) - 0^n.
a(2*n+1) + a(2*n+2) = A000302(n+1). - Paul Curtz, Jun 30 2008
G.f.: 1 - x + x*Q(0), where Q(k) = 1 + 2*x^2 + (4*k+5)*x - x*(4*k+1 + 2*x)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Oct 07 2013
E.g.f.: (1/3)*( 2*exp(2*x) + 4*exp(-x) - 3 ). - G. C. Greubel, Aug 19 2022

Extensions

Obscure variable k in Orlovsky comment replaced with a(n) by R. J. Mathar, Apr 23 2009

A256890 Triangle T(n,k) = t(n-k, k); t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 2.

Original entry on oeis.org

1, 2, 2, 4, 12, 4, 8, 52, 52, 8, 16, 196, 416, 196, 16, 32, 684, 2644, 2644, 684, 32, 64, 2276, 14680, 26440, 14680, 2276, 64, 128, 7340, 74652, 220280, 220280, 74652, 7340, 128, 256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172, 256, 512, 72076, 1637860, 10978444, 27227908, 27227908, 10978444, 1637860, 72076, 512
Offset: 0

Views

Author

Dale Gerdemann, Apr 12 2015

Keywords

Comments

Related triangles may be found by varying the function f(x). If f(x) is a linear function, it can be parameterized as f(x) = a*x + b. With different values for a and b, the following triangles are obtained:
a\b 1.......2.......3.......4.......5.......6
The row sums of these, and similarly constructed number triangles, are shown in the following table:
a\b 1.......2.......3.......4.......5.......6.......7.......8.......9
The formula can be further generalized to: t(n,m) = f(m+s)*t(n-1,m) + f(n-s)*t(n,m-1), where f(x) = a*x + b. The following table specifies triangles with nonzero values for s (given after the slash).
a\b 0 1 2 3
-2 A130595/1
-1
0
With the absolute value, f(x) = |x|, one obtains A038221/3, A038234/4, A038247/5, A038260/6, A038273/7, A038286/8, A038299/9 (with value for s after the slash).
If f(x) = A000045(x) (Fibonacci) and s = 1, the result is A010048 (Fibonomial).
In the notation of Carlitz and Scoville, this is the triangle of generalized Eulerian numbers A(r, s | alpha, beta) with alpha = beta = 2. Also the array A(2,1,4) in the notation of Hwang et al. (see page 31). - Peter Bala, Dec 27 2019

Examples

			Array, t(n, k), begins as:
   1,    2,      4,        8,        16,         32,          64, ...;
   2,   12,     52,      196,       684,       2276,        7340, ...;
   4,   52,    416,     2644,     14680,      74652,      357328, ...;
   8,  196,   2644,    26440,    220280,    1623964,    10978444, ...;
  16,  684,  14680,   220280,   2643360,   27227908,   251195000, ...;
  32, 2276,  74652,  1623964,  27227908,  381190712,  4677894984, ...;
  64, 7340, 357328, 10978444, 251195000, 4677894984, 74846319744, ...;
Triangle, T(n, k), begins as:
    1;
    2,     2;
    4,    12,      4;
    8,    52,     52,       8;
   16,   196,    416,     196,      16;
   32,   684,   2644,    2644,     684,      32;
   64,  2276,  14680,   26440,   14680,    2276,     64;
  128,  7340,  74652,  220280,  220280,   74652,   7340,   128;
  256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172,   256;
		

Crossrefs

Programs

  • Magma
    A256890:= func< n,k | (&+[(-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n: j in [0..k]]) >;
    [A256890(n,k): k in [0..n], n in [0..10]]; // G. C. Greubel, Oct 18 2022
    
  • Mathematica
    Table[Sum[(-1)^(k-j)*Binomial[j+3, j] Binomial[n+4, k-j] (j+2)^n, {j,0,k}], {n,0, 9}, {k,0,n}]//Flatten (* Michael De Vlieger, Dec 27 2019 *)
  • PARI
    t(n,m) = if ((n<0) || (m<0), 0, if ((n==0) && (m==0), 1, (m+2)*t(n-1, m) + (n+2)*t(n, m-1)));
    tabl(nn) = {for (n=0, nn, for (k=0, n, print1(t(n-k, k), ", ");); print(););} \\ Michel Marcus, Apr 14 2015
    
  • SageMath
    def A256890(n,k): return sum((-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n for j in range(k+1))
    flatten([[A256890(n,k) for k in range(n+1)] for n in range(11)]) # G. C. Greubel, Oct 18 2022

Formula

T(n,k) = t(n-k, k); t(0,0) = 1, t(n,m) = 0 if n < 0 or m < 0 else t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 2.
Sum_{k=0..n} T(n, k) = A001715(n).
T(n,k) = Sum_{j = 0..k} (-1)^(k-j)*binomial(j+3,j)*binomial(n+4,k-j)*(j+2)^n. - Peter Bala, Dec 27 2019
Modified rule of Pascal: T(0,0) = 1, T(n,k) = 0 if k < 0 or k > n else T(n,k) = f(n-k) * T(n-1,k-1) + f(k) * T(n-1,k), where f(x) = x + 2. - Georg Fischer, Nov 11 2021
From G. C. Greubel, Oct 18 2022: (Start)
T(n, n-k) = T(n, k).
T(n, 0) = A000079(n). (End)
Previous Showing 91-100 of 587 results. Next