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 11-20 of 196 results. Next

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

A120452 Number of partitions of n-1 boys and one girl with no couple.

Original entry on oeis.org

1, 1, 3, 5, 9, 14, 23, 34, 52, 75, 109, 153, 216, 296, 407, 549, 739, 981, 1300, 1702, 2224, 2879, 3716, 4761, 6083, 7721, 9774, 12306, 15450, 19307, 24064, 29867, 36978, 45614, 56130, 68846, 84250, 102793, 125148, 151955, 184123, 222553, 268482
Offset: 1

Views

Author

Yasutoshi Kohmoto, Jul 20 2006

Keywords

Comments

From Gus Wiseman, Jun 08 2021: (Start)
Also the number of:
- integer partitions of 2n with reverse-alternating sum 2;
- reversed integer partitions of 2n with alternating sum 2;
- integer partitions of 2n with exactly two odd parts, one of which is the greatest;
- odd-length integer partitions of 2n whose conjugate partition has exactly two odd parts.
Note that integer partitions of 2n with alternating or reverse-alternating sum 0 are counted by A000041, ranked by A000290.
(End)

Examples

			n=5:
If partitions have no pair "o*", then a(5)=9 ("o" means a boy, "*" means a girl): {o, o, o, o, *}, {o, o, *, oo}, {*, oo, oo}, {o, *, ooo}, {o, o, oo*}, {oo, oo*}, {*, oooo}, {o, ooo*}, {oooo*}.
From _Gus Wiseman_, Jun 08 2021: (Start)
The a(1) = 1 through a(6) = 14 partitions of 2n with reverse-alternating sum 2:
  (2)  (211)  (222)    (332)      (442)        (552)
              (321)    (431)      (541)        (651)
              (21111)  (22211)    (22222)      (33222)
                       (32111)    (32221)      (33321)
                       (2111111)  (33211)      (43221)
                                  (43111)      (44211)
                                  (2221111)    (54111)
                                  (3211111)    (2222211)
                                  (211111111)  (3222111)
                                               (3321111)
                                               (4311111)
                                               (222111111)
                                               (321111111)
                                               (21111111111)
For example, the partition (43221) has reverse-alternating sum 1 - 2 + 2 - 3 + 4 = 2, so is counted under a(6).
The a(1) = 1 through a(6) = 14 partitions of 2n with exactly two odd parts, one of which is the greatest:
  (11)  (31)  (33)   (53)    (55)     (75)
              (51)   (71)    (73)     (93)
              (321)  (332)   (91)     (111)
                     (521)   (532)    (543)
                     (3221)  (541)    (552)
                             (721)    (732)
                             (3322)   (741)
                             (5221)   (921)
                             (32221)  (5322)
                                      (5421)
                                      (7221)
                                      (33222)
                                      (52221)
                                      (322221)
(End)
		

Crossrefs

A diagonal of A103919.
A diagonal of A344612.
A000097 counts partitions of 2n with alternating sum 2.
A001700/A088218 appear to count compositions with reverse-alternating sum 2.
A058696 counts partitions of 2n, ranked by A300061.
A344610 counts partitions of 2n by sum and positive reverse-alternating sum.
A344611 counts partitions of 2n with reverse-alternating sum >= 0.
A344741 counts partitions of 2n with reverse-alternating sum -2.

Programs

  • Mathematica
    a[n_] := Total[PartitionsP[Range[0, n-3]]] + PartitionsP[n-1];
    Array[a, 50] (* Jean-François Alcover, Jun 05 2021 *)

Formula

a(n) = A000070(n-2) + A002865(n-1). - Fung Cheok Yin (cheokyin_restart(AT)yahoo.com.hk), Aug 15 2006
a(n) = A000070(n-1) - A000041(n-2) = A000070(n-3) + A000041(n-1). - Max Alekseyev, Aug 23 2006
a(n) ~ exp(Pi*sqrt(2*n/3)) / (2^(3/2)*Pi*sqrt(n)) * (1 - 37*Pi/(24*sqrt(6*n))). - Vaclav Kotesovec, Oct 25 2016

Extensions

More terms from Fung Cheok Yin (cheokyin_restart(AT)yahoo.com.hk), Aug 15 2006
More terms from Max Alekseyev, Aug 23 2006

A238349 Triangle T(n,k) read by rows: T(n,k) is the number of compositions of n with k parts p at position p (fixed points), n>=0, 0<=k<=n.

Original entry on oeis.org

1, 0, 1, 1, 1, 0, 2, 1, 1, 0, 3, 4, 1, 0, 0, 6, 7, 3, 0, 0, 0, 11, 16, 4, 1, 0, 0, 0, 22, 29, 12, 1, 0, 0, 0, 0, 42, 60, 23, 3, 0, 0, 0, 0, 0, 82, 120, 47, 7, 0, 0, 0, 0, 0, 0, 161, 238, 100, 12, 1, 0, 0, 0, 0, 0, 0, 316, 479, 198, 30, 1, 0, 0, 0, 0, 0, 0, 0, 624, 956, 404, 61, 3, 0, 0, 0, 0, 0, 0, 0, 0, 1235, 1910, 818, 126, 7, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Joerg Arndt and Alois P. Heinz, Feb 25 2014

Keywords

Comments

T(n*(n+3)/2,n) = A227682(n).
From Vaclav Kotesovec, Sep 07 2014: (Start)
In general, column k is asymptotic to c(k) * 2^n. The constants c(k) numerically:
c(0) = 0.144394047543301210639449860964615390044455952420342... = A048651/2
c(1) = 0.231997216225445223894202367545783700531838988546098... = c(0)*A065442
c(2) = 0.104261929557371534733906196116707679501974368826074...
c(3) = 0.017956317806894073430249112172514186063327165575720...
c(4) = 0.001343254222922697613125145839110293324517874530073...
c(5) = 0.000046459767012163920051487037952792359225887287888...
c(6) = 0.000000768651747857094917953943327540619110335556499...
c(7) = 0.000000006200599904985793344094393321042983316604040...
c(8) = 0.000000000024656652167851516076173236693314090168122...
c(9) = 0.000000000000048633746319332356416193899916110113745...
c(10)= 0.000000000000000047750743608910618576944191079881479...
c(20)= 1.05217230403079700467566...*10^(-63)
For big k is c(k) ~ m * 2^(-k*(k+1)/2), where m = 1/(4*c(0)) = 1/(2*A048651) = 1.7313733097275318...
(End)

Examples

			Triangle starts:
00:  1,
01:  0, 1,
02:  1, 1, 0,
03:  2, 1, 1, 0,
04:  3, 4, 1, 0, 0,
05:  6, 7, 3, 0, 0, 0,
06:  11, 16, 4, 1, 0, 0, 0,
07:  22, 29, 12, 1, 0, 0, 0, 0,
08:  42, 60, 23, 3, 0, 0, 0, 0, 0,
09:  82, 120, 47, 7, 0, 0, 0, 0, 0, 0,
10:  161, 238, 100, 12, 1, 0, 0, 0, 0, 0, 0,
11:  316, 479, 198, 30, 1, 0, 0, 0, 0, 0, 0, 0,
12:  624, 956, 404, 61, 3, 0, 0, 0, 0, 0, 0, 0, 0,
13:  1235, 1910, 818, 126, 7, 0, 0, 0, 0, 0, 0, 0, 0, 0,
14:  2449, 3817, 1652, 258, 16, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
15:  4864, 7633, 3319, 537, 30, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
...
From _Gus Wiseman_, Apr 03 2022: (Start)
Row n = 5 counts the following compositions (empty columns indicated by dots):
  (5)     (14)     (113)   .  .  .
  (23)    (32)     (122)
  (41)    (131)    (1211)
  (212)   (221)
  (311)   (1112)
  (2111)  (1121)
          (11111)
(End)
		

References

  • M. Archibald, A. Blecher and A. Knopfmacher, Fixed points in compositions and words, accepted by the Journal of Integer Sequences.

Crossrefs

Row sums are A011782.
The version for permutations is A008290.
The version with all zeros removed is A238350.
The version for reversed partitions is A238352.
The corresponding rank statistic is A352512, nonfixed A352513.
The version for nonfixed points is A352523, A352520 (k=1).
Below: comps = compositions, first = column k=0, stat = rank statistic.
- A352521 counts comps by strong nonexcedances, first A219282, stat A352514.
- A352522 counts comps by weak nonexcedances, first A238874, stat A352515.
- A352524 counts comps by strong excedances, first A008930, stat A352516.
- A352525 counts comps by weak excedances, A177510 (k=1), stat A352517.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, expand(
           add(b(n-j, i+1)*`if`(i=j, x, 1), j=1..n)))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(n, 1)):
    seq(T(n), n=0..15);
  • Mathematica
    b[n_, i_] := b[n, i] = If[n == 0, 1, Expand[Sum[b[n-j, i+1]*If[i == j, x, 1], {j, 1, n}]]]; T[n_] := Function[{p}, Table[Coefficient[p, x, i], {i, 0, n}]][b[n, 1]]; Table[T[n], {n, 0, 15}] // Flatten (* Jean-François Alcover, Jan 06 2015, after Alois P. Heinz *)
    pq[y_]:=Length[Select[Range[Length[y]],#==y[[#]]&]];
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],pq[#]==k&]],{n,0,9},{k,0,n}] (* Gus Wiseman, Apr 03 2022 *)

A114121 Expansion of (sqrt(1 - 4*x) + (1 - 2*x))/(2*(1 - 4*x)).

Original entry on oeis.org

1, 2, 7, 26, 99, 382, 1486, 5812, 22819, 89846, 354522, 1401292, 5546382, 21977516, 87167164, 345994216, 1374282019, 5461770406, 21717436834, 86392108636, 343801171354, 1368640564996, 5450095992964, 21708901408216, 86492546019214
Offset: 0

Views

Author

Paul Barry, Feb 13 2006

Keywords

Comments

Second binomial transform of A032443 with interpolated zeros.
a(n) is the total number of lattice points, taken over all Dyck n-paths (A000108), that (i) lie on or above ground level and (ii) lie on or directly below a peak. For example with n = 2, UUDD has 1 peak contributing 3 lattice points--(2, 0), (2, 1) and (2, 2) when the path starts at the origin--and UDUD has 2 peaks, each contributing 2 lattice points and so a(2) = 3 + 4 = 7. - David Callan, Jul 14 2006
Hankel transform is binomial(n + 2, 2). - Paul Barry, Dec 04 2007
Image of (-1)^n under the Riordan array ((1/2)(1/(1 - 4x) + 1/sqrt(1 - 4x)), c(x) - 1), c(x) the g.f. of A000108. - Paul Barry, Jun 15 2008
From Gus Wiseman, Jun 21 2021: (Start)
Also the even bisection of A116406 = number compositions of n with alternating sum >= 0, where the alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i. The a(3) = 26 compositions are:
(6) (33) (114) (1122) (11112) (111111)
(42) (123) (1131) (11121)
(51) (132) (1221) (11211)
(213) (2112) (12111)
(222) (2121) (21111)
(231) (2211)
(312) (3111)
(321)
(411)
(End)

Examples

			G.f. = 1 + 2*x + 7*x^2 + 26*x^3 + 99*x^4 + 382*x^5 + 1486*x^6 + 5812*x^7 + ...
		

Crossrefs

The case of alternating sum = 0 is A001700.
The case of alternating sum < 0 is A008549.
This is the even bisection of A116406.
The restriction to reversed partitions is A344611.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A124754 gives the alternating sum of standard compositions.
A316524 is the alternating sum of the prime indices of n.
A344611 counts partitions of 2n with reverse-alternating sum >= 0.

Programs

  • Maple
    seq(sum(binomial(2*n,2*k+irem(n,2)),k=0..floor((1/2)*n)),n=0..20)
    seq(binomial(2*n-1,n)+4^(n-1)-(1/4)*0^n,n=0..20)
  • Mathematica
    a[ n_] := SeriesCoefficient[((1 + 1/Sqrt[1 - 4 x])/2)^2, {x, 0, n}] (* Michael Somos, Dec 31 2013 *)
    ats[y_]:=Sum[(-1)^(i-1)*y[[i]],{i,Length[y]}];Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],ats[#]>=0&]],{n,0,15,2}] (* Gus Wiseman, Jun 21 2021 *)

Formula

a(n) = Sum_{k=0..n} C(n, k)*2^(n-k-2)*(2^k + C(k, k/2))*(1 + (-1)^k).
a(n) = (A000984(n) + A081294(n))/2.
From Paul Barry, Jun 15 2008: (Start)
G.f.: (1 - 4*x + (1 - 2*x)*sqrt(1 - 4*x))/(2*(1 - 4*x)^(3/2)).
a(n) = Sum_{k=0..n} ( Sum_{j=0..n} C(2*n, n-k-j)*(-1)^j ). (End)
a(n) = Sum_{k=0..n} C(2*n, n-k)*(1 + (-1)^k)/2. - Paul Barry, Aug 06 2009
From Paul Barry, Sep 07 2009: (Start)
a(n) = C(2*n-1, n-1) + (4^n + 3*0^n)/4.
Integral representation a(n) = (1/(2*pi))*(Integral_{x=0..4} x^n/sqrt(x(4 - x))) + (4^n + 0^n)/4. (End)
a(n) = Sum_{k=0..floor(n/2)} C(2*n, 2*k + (n mod 2)). - Mircea Merca, Jun 21 2011
Conjecture: n*a(n) + 2*(3 - 4*n)*a(n-1) + 8*(2*n-3)*a(n-2) = 0. - R. J. Mathar, Nov 07 2012
Conjecture verified using the differential equation (16*x^2-8*x+1)*g'(x) + (8*x-2)*g(x)-2*x=0 satisfied by the G.f. - Robert Israel, Jul 27 2020
a(n) = Sum_{i=0..n} (sum_{j=0..n} binomial(n, i+j)*binomial(n, j-i)). - Yalcin Aktar, Jan 07 2013.
G.f.: (1 + (1 - 4*x)^(-1/2))^2 / 4. Convolution square of A088218. - Michael Somos, Dec 31 2013
0 = (1 + 2*n)*b(n) - (5 + 4*n)*b(n+1) + (4 + 2*n)*b(n+2) if n > 0 where b(n) = a(n) / 4^n. - Michael Somos, Dec 31 2013
0 = b(n+3) * (2*b(n+2) - 7*b(n+1) + 5*b(n)) + b(n+2) * (-b(n+2) + 7*b(n+1) - 7*b(n)) + b(n+1) * (-b(n+1) + 2*b(n)) if n > 0 where b(n) = a(n) / 4^n. - Michael Somos, Dec 31 2013
For n > 0, a(n) = 2^(2n-1) - A008549(n). - Gus Wiseman, Jun 27 2021
a(n) = [x^n] 1/((1-2*x) * (1-x)^(n-1)). - Seiichi Manyama, Apr 10 2024

A026810 Number of partitions of n in which the greatest part is 4.

Original entry on oeis.org

0, 0, 0, 0, 1, 1, 2, 3, 5, 6, 9, 11, 15, 18, 23, 27, 34, 39, 47, 54, 64, 72, 84, 94, 108, 120, 136, 150, 169, 185, 206, 225, 249, 270, 297, 321, 351, 378, 411, 441, 478, 511, 551, 588, 632, 672, 720, 764, 816, 864, 920, 972, 1033, 1089, 1154, 1215, 1285, 1350
Offset: 0

Views

Author

Keywords

Comments

Also number of partitions of n into exactly 4 parts.
Also the number of weighted cubic graphs on 4 nodes (=the tetrahedron) with weight n. - R. J. Mathar, Nov 03 2018
From Gus Wiseman, Jun 27 2021: (Start)
Also the number of strict integer partitions of 2n with alternating sum 4, or (by conjugation) partitions of 2n covering an initial interval of positive integers with exactly 4 odd parts. The strict partitions with alternating sum 4 are:
(4) (5,1) (6,2) (7,3) (8,4) (9,5) (10,6)
(5,2,1) (5,3,2) (5,4,3) (6,5,3) (7,6,3)
(6,3,1) (6,4,2) (7,5,2) (8,6,2)
(7,4,1) (8,5,1) (9,6,1)
(6,3,2,1) (6,4,3,1) (6,5,4,1)
(7,4,2,1) (7,4,3,2)
(7,5,3,1)
(8,5,2,1)
(6,4,3,2,1)
(End)

Examples

			From _Gus Wiseman_, Jun 27 2021: (Start)
The a(4) = 1 through a(10) = 9 partitions of length 4:
  (1111)  (2111)  (2211)  (2221)  (2222)  (3222)  (3322)
                  (3111)  (3211)  (3221)  (3321)  (3331)
                          (4111)  (3311)  (4221)  (4222)
                                  (4211)  (4311)  (4321)
                                  (5111)  (5211)  (4411)
                                          (6111)  (5221)
                                                  (5311)
                                                  (6211)
                                                  (7111)
(End)
		

References

  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 275.
  • D. E. Knuth, The Art of Computer Programming, vol. 4,fascicle 3, Generating All Combinations and Partitions, Section 7.2.1.4., p. 56, exercise 31.

Crossrefs

Cf. A001400, A026811, A026812, A026813, A026814, A026815, A026816, A069905 (3 positive parts), A002621 (partial sums), A005044 (first differences).
A non-strict version is A000710 or A088218.
This is column k = 2 of A152146.
A reverse version is A343941.

Programs

  • Magma
    [Round((n^3+3*n^2-9*n*(n mod 2))/144): n in [0..60]]; // Vincenzo Librandi, Oct 14 2015
  • Maple
    A049347 := proc(n)
            op(1+(n mod 3),[1,-1,0]) ;
    end proc:
    A056594 := proc(n)
            op(1+(n mod 4),[1,0,-1,0]) ;
    end proc:
    A026810 := proc(n)
            1/288*(n+1)*(2*n^2+4*n-13+9*(-1)^n) ;
            %-A049347(n)/9 ;
            %+A056594(n)/8 ;
    end proc: # R. J. Mathar, Jul 03 2012
  • Mathematica
    Table[Count[IntegerPartitions[n], {4, _}], {n, 0, 60}]
    LinearRecurrence[{1, 1, 0, 0, -2, 0, 0, 1, 1, -1}, {0, 0, 0, 0, 1, 1, 2, 3, 5, 6}, 60] (* Vincenzo Librandi, Oct 14 2015 *)
    Table[Length[IntegerPartitions[n, {4}]], {n, 0, 60}] (* Eric Rowland, Mar 02 2017 *)
    CoefficientList[Series[x^4/Product[1 - x^k, {k, 1, 4}], {x, 0, 60}], x] (* Robert A. Russell, May 13 2018 *)
  • PARI
    for(n=0, 60, print(n, " ", round((n^3 + 3*n^2 -9*n*(n % 2))/144))); \\ Washington Bomfim, Jul 03 2012
    
  • PARI
    x='x+O('x^60); concat([0, 0, 0, 0], Vec(x^4/((1-x)*(1-x^2)*(1-x^3)*(1-x^4)))) \\ Altug Alkan, Oct 14 2015
    
  • PARI
    vector(60, n, n--; (n+1)*(2*n^2+4*n-13+9*(-1)^n)/288 + real(I^n)/8 - ((n+2)%3-1)/9) \\ Altug Alkan, Oct 26 2015
    
  • PARI
    print1(0,", "); for(n=1,60,j=0;forpart(v=n,j++,,[4,4]); print1(j,", ")) \\ Hugo Pfoertner, Oct 01 2018
    

Formula

G.f.: x^4/((1-x)*(1-x^2)*(1-x^3)*(1-x^4)) = x^4/((1-x)^4*(1+x)^2*(1+x+x^2)*(1+x^2)).
a(n+4) = A001400(n). - Michael Somos, Apr 07 2012
a(n) = round( (n^3 + 3*n^2 -9*n*(n mod 2))/144 ). - Washington Bomfim, Jan 06 2021 and Jul 03 2012
a(n) = (n+1)*(2*n^2+4*n-13+9*(-1)^n)/288 -A049347(n)/9 +A056594(n)/8. - R. J. Mathar, Jul 03 2012
From Gregory L. Simay, Oct 13 2015: (Start)
a(n) = (n^3 + 3*n^2 - 9*n)/144 + a(m) - (m^3 + 3*m^2 - 9*m)/144 if n = 12k + m and m is odd. For example, a(23) = a(12*1 + 11) = (23^3 + 3*23^2 - 9*23)/144 + a(11) - (11^3 + 3*11^2 - 9*11)/144 = 94.
a(n) = (n^3 + 3*n^2)/144 + a(m) - (m^3 + 3*m^2)/144 if n = 12k + m and m is even. For example, a(22) = a(12*1 + 10) = (22^3 + 3*22^2)/144 + a(10) - (10^3 + 3*10^2)/144 = 84. (End)
a(n) = A008284(n,4). - Robert A. Russell, May 13 2018
From Gregory L. Simay, Jul 28 2019: (Start)
a(2n+1) = a(2n) + a(n+1) - a(n-3) and
a(2n) = a(2n-1) + a(n+2) - a(n-2). (End)

A347437 Number of factorizations of n with integer alternating product.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Sep 06 2021

Keywords

Comments

A factorization of n is a weakly increasing sequence of positive integers > 1 with product n.
We define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)).

Examples

			The factorizations for n = 4, 16, 36, 48, 54, 64, 108:
  (4)   (16)      (36)      (48)        (54)    (64)          (108)
  (2*2) (4*4)     (6*6)     (2*4*6)     (2*3*9) (8*8)         (2*6*9)
        (2*2*4)   (2*2*9)   (3*4*4)     (3*3*6) (2*4*8)       (3*6*6)
        (2*2*2*2) (2*3*6)   (2*2*12)            (4*4*4)       (2*2*27)
                  (3*3*4)   (2*2*2*2*3)         (2*2*16)      (2*3*18)
                  (2*2*3*3)                     (2*2*4*4)     (3*3*12)
                                                (2*2*2*2*4)   (2*2*3*3*3)
                                                (2*2*2*2*2*2)
		

Crossrefs

Positions of 1's are A005117, complement A013929.
Allowing any alternating product <= 1 gives A339846.
Allowing any alternating product > 1 gives A339890.
The restriction to powers of 2 is A344607.
The even-length case is A347438, also the case of alternating product 1.
The reciprocal version is A347439.
Allowing any alternating product < 1 gives A347440.
The odd-length case is A347441.
The reverse version is A347442.
The additive version is A347446, ranked by A347457.
Allowing any alternating product >= 1 gives A347456.
The restriction to perfect squares is A347458, reciprocal A347459.
The ordered version is A347463.
A001055 counts factorizations.
A046099 counts factorizations with no alternating permutations.
A071321 gives the alternating sum of prime factors of n (reverse: A071322).
A273013 counts ordered factorizations of n^2 with alternating product 1.
A347460 counts possible alternating products of factorizations.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    altprod[q_]:=Product[q[[i]]^(-1)^(i-1),{i,Length[q]}];
    Table[Length[Select[facs[n],IntegerQ@*altprod]],{n,100}]
  • PARI
    A347437(n, m=n, ap=1, e=0) = if(1==n, if(e%2, 1==denominator(ap), 1==numerator(ap)), sumdiv(n, d, if((d>1)&&(d<=m), A347437(n/d, d, ap * d^((-1)^e), 1-e)))); \\ Antti Karttunen, Oct 22 2023

Formula

a(2^n) = A344607(n).
a(n^2) = A347458(n).

Extensions

Data section extended up to a(108) by Antti Karttunen, Oct 22 2023

A347438 Number of unordered factorizations of n with alternating product 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Sep 06 2021

Keywords

Comments

Also the number of unordered factorizations of n with alternating sum 0.
Also the number of unordered factorizations of n with all even multiplicities.
This is the even-length case of A347437, the odd-length case being A347441.
An unordered factorization of n is a weakly increasing sequence of positive integers > 1 with product n.
We define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)).

Examples

			The a(n) factorizations for n = 16, 64, 144, 256, 576:
  4*4      8*8          12*12        16*16            24*24
  2*2*2*2  2*2*4*4      2*2*6*6      2*2*8*8          3*3*8*8
           2*2*2*2*2*2  3*3*4*4      4*4*4*4          4*4*6*6
                        2*2*2*2*3*3  2*2*2*2*4*4      2*2*12*12
                                     2*2*2*2*2*2*2*2  2*2*2*2*6*6
                                                      2*2*3*3*4*4
                                                      2*2*2*2*2*2*3*3
		

Crossrefs

Positions of zeros are A000037.
Positions of nonzero terms are A000290.
The restriction to perfect squares is A001055 (ordered: A273013).
The restriction to powers of 2 is A035363.
The additive version is A119620, ranked by A028982.
Positions of non-1's are A213367 \ {1}.
Positions of 1's are A280076 = {1} \/ A001248.
Sorted first positions are 1, 2, and all terms of A330972 squared.
Allowing any alternating product <= 1 gives A339846.
Allowing any alternating product > 1 gives A339890.
Allowing any integer alternating product gives A347437.
Allowing any integer reciprocal alternating product gives A347439.
Allowing any alternating product < 1 gives A347440.
Allowing any alternating product >= 1 gives A347456.
A046099 counts factorizations with no alternating permutations.
A071321 gives the alternating sum of prime factors (reverse: A071322).
A316524 gives the alternating sum of prime indices (reverse: A344616).
A344606 counts alternating permutations of prime factors.
A347441 counts odd-length factorizations with integer alternating product.
A347460 counts possible alternating products of factorizations.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    altprod[q_]:=Product[q[[i]]^(-1)^(i-1),{i,Length[q]}];
    Table[Length[Select[facs[n],altprod[#]==1&]],{n,100}]
  • PARI
    A347438(n, m=n, k=0, t=1) = if(1==n, (1==t), my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A347438(n/d, d, 1-k, t*(d^((-1)^k))))); (s)); \\ Antti Karttunen, Oct 30 2021

Formula

a(2^n) = A035363(n).
a(n^2) = A001055(n).

Extensions

Name and comments clarified (with unordered) by Jacob Sprittulla, Oct 05 2021

A024718 a(n) = (1/2)*(1 + Sum_{k=0..n} binomial(2*k, k)).

Original entry on oeis.org

1, 2, 5, 15, 50, 176, 638, 2354, 8789, 33099, 125477, 478193, 1830271, 7030571, 27088871, 104647631, 405187826, 1571990936, 6109558586, 23782190486, 92705454896, 361834392116, 1413883873976, 5530599237776, 21654401079326, 84859704298202, 332818970772254
Offset: 0

Views

Author

Keywords

Comments

Total number of leaves in all rooted ordered trees with at most n edges. - Michael Somos, Feb 14 2006
Also: Number of UH-free Schroeder paths of semilength n with horizontal steps only at level less than two [see Yan]. - R. J. Mathar, May 24 2008
Hankel transform is A010892. - Paul Barry, Apr 28 2009
Binomial transform of A005773. - Philippe Deléham, Dec 13 2009
Number of vertices all of whose children are leaves in all ordered trees with n+1 edges. Example: a(3) = 15; for an explanation see David Callan's comment in A001519. - Emeric Deutsch, Feb 12 2015

Examples

			G.f. = 1 + 2*x + 5*x^2 + 15*x^3 + 50*x^4 + 176*x^5 + 638*x^6 + ...
		

Crossrefs

Partial sums of A088218.
Bisection of A086905.
Second column of triangle A102541.

Programs

  • Maple
    A024718 := n -> (binomial(2*n, n)*hypergeom([1, -n], [1/2 - n], 1/4) + 1)/2:
    seq(simplify(A024718(n)), n = 0..26); # Peter Luschny, Dec 15 2024
  • Mathematica
    Table[Sum[Binomial[2k-1,k-1],{k,0,n}],{n,0,100}] (* Emanuele Munarini, May 18 2018 *)
  • PARI
    a(n) = (1 + sum(k=0, n, binomial(2*k, k)))/2; \\ Michel Marcus, May 18 2018

Formula

a(n) = A079309(n) + 1.
G.f.: 1/((1 - x)*(2 - C)), where C = g.f. for the Catalan numbers A000108. - N. J. A. Sloane, Aug 30 2002
Given g.f. A(x), then x * A(x - x^2) is the g.f. of A024494. - Michael Somos, Feb 14 2006
G.f.: (1 + 1 / sqrt(1 - 4*x)) / (2 - 2*x). - Michael Somos, Feb 14 2006
D-finite with recurrence: n*a(n) - (5*n-2)*a(n-1) + 2*(2*n-1)*a(n-2) = 0. - R. J. Mathar, Dec 02 2012
Remark: The above recurrence is true (it can be easily proved by differentiating the generating function). Notice that it is the same recurrence satisfied by the partial sums of the central binomial coefficients (A006134). - Emanuele Munarini, May 18 2018
0 = a(n)*(16*a(n+1) - 22*a(n+2) + 6*a(n+3)) + a(n+1)*(-18*a(n+1) + 27*a(n+2) - 7*a(n+3)) + a(n+2)*(-3*a(n+2) + a(n+3)) for all n in Z if a(n) = 1/2 for n < 0. - Michael Somos, Apr 23 2014
a(n) = ((1 - I/sqrt(3))/2 - binomial(2*n+1, n)*hypergeom([n+3/2, 1], [n+2], 4)). - Peter Luschny, May 18 2018
a(n) = [x^n] 1/((1-x+x^2) * (1-x)^n). - Seiichi Manyama, Apr 06 2024

Extensions

Name edited by Petros Hadjicostas, Aug 04 2020

A163493 Number of binary strings of length n which have the same number of 00 and 01 substrings.

Original entry on oeis.org

1, 2, 2, 3, 6, 9, 15, 30, 54, 97, 189, 360, 675, 1304, 2522, 4835, 9358, 18193, 35269, 68568, 133737, 260802, 509132, 995801, 1948931, 3816904, 7483636, 14683721, 28827798, 56637969, 111347879, 219019294, 431043814, 848764585, 1672056525, 3295390800, 6497536449
Offset: 0

Views

Author

Keywords

Comments

A variation of problem 11424 in the American Mathematical Monthly. Terms were brute-force calculated using Maple 10.
Proposed Problem 11610 in the Dec 2011 A.M.M.
From Gus Wiseman, Jul 27 2021: (Start)
Also the antidiagonal sums of the matrices counting integer compositions by length and alternating sum (A345197). So a(n) is the number of integer compositions of n + 1 of length (n - s + 3)/2, where s is the alternating sum of the composition. For example, the a(0) = 1 through a(6) = 7 compositions are:
(1) (2) (3) (4) (5) (6) (7)
(11) (21) (31) (41) (51) (61)
(121) (122) (123) (124)
(221) (222) (223)
(1112) (321) (322)
(1211) (1122) (421)
(1221) (1132)
(2112) (1231)
(2211) (2122)
(2221)
(3112)
(3211)
(11131)
(12121)
(13111)
For a bijection with the main (binary string) interpretation, take the run-lengths of each binary string of length n + 1 that satisfies the condition and starts with 1.
(End)

Examples

			1 + 2*x + 2*x^2 + 3*x^3 + 6*x^4 + 9*x^5 + 15*x^6 + 30*x^7 + 54*x^8 + 97*x^9 + ...
From _Gus Wiseman_, Jul 27 2021: (Start)
The a(0) = 1 though a(6) = 15 binary strings:
  ()  (0)  (1,0)  (0,0,1)  (0,0,1,0)  (0,0,1,1,0)  (0,0,0,1,0,1)
      (1)  (1,1)  (1,1,0)  (0,0,1,1)  (0,0,1,1,1)  (0,0,1,0,0,1)
                  (1,1,1)  (0,1,0,0)  (0,1,1,0,0)  (0,0,1,1,1,0)
                           (1,0,0,1)  (1,0,0,1,0)  (0,0,1,1,1,1)
                           (1,1,1,0)  (1,0,0,1,1)  (0,1,0,0,0,1)
                           (1,1,1,1)  (1,0,1,0,0)  (0,1,1,1,0,0)
                                      (1,1,0,0,1)  (1,0,0,1,1,0)
                                      (1,1,1,1,0)  (1,0,0,1,1,1)
                                      (1,1,1,1,1)  (1,0,1,1,0,0)
                                                   (1,1,0,0,1,0)
                                                   (1,1,0,0,1,1)
                                                   (1,1,0,1,0,0)
                                                   (1,1,1,0,0,1)
                                                   (1,1,1,1,1,0)
                                                   (1,1,1,1,1,1)
(End)
		

Crossrefs

Antidiagonal sums of the matrices A345197.
Row sums of A345907.
Taking diagonal instead of antidiagonal sums gives A345908.
A011782 counts compositions (or binary strings).
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).
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

  • Maple
    with(combinat): count := proc(n) local S, matches, A, k, i; S := subsets(\{seq(i, i=1..n)\}): matches := 0: while not S[finished] do A := S[nextvalue](): k := 0: for i from 1 to n-1 do: if not (i in A) and not (i+1 in A) then k := k + 1: fi: if not (i in A) and (i+1 in A) then k := k - 1: fi: od: if (k = 0) then matches := matches + 1: fi: end do; return(matches); end proc:
    # second Maple program:
    b:= proc(n, l, t) option remember; `if`(n-abs(t)<0, 0, `if`(n=0, 1,
          add(b(n-1, i, t+`if`(l=0, (-1)^i, 0)), i=0..1)))
        end:
    a:= n-> b(n, 1, 0):
    seq(a(n), n=0..36);  # Alois P. Heinz, Mar 20 2024
  • Mathematica
    a[0] = 1; a[n_] := Sum[Binomial[2*k - 1, k]*Binomial[n - 2*k, k] + Binomial[2*k, k]*Binomial[n - 2*k - 1, k], {k, 0, n/3}];
    Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Nov 28 2017, after Joel B. Lewis *)
    Table[Length[Select[Tuples[{0,1},n],Count[Partition[#,2,1],{0,0}]==Count[Partition[#,2,1],{0,1}]&]],{n,0,10}] (* Gus Wiseman, Jul 27 2021 *)
    a[0]:=1; a[n_]:=(1 + 3*HypergeometricPFQ[{1/2, 1-3*n/8, (1-n)/3, (2-n)/3, -n/3},{1, (1-n)/2, 1-n/2, -3*n/8}, -27])/2; Array[a,37,0] (* Stefano Spezia, Apr 26 2024 *)
  • Python
    from math import comb
    def A163493(n): return 2+sum((x:=comb((k:=m<<1)-1,m)*comb(n-k,m))+(x*(n-3*m)<<1)//(n-k) for m in range(1,n//3+1)) if n else 1 # Chai Wah Wu, May 01 2024

Formula

G.f.: 1/2/(1-x) + (1+2*x)/2/sqrt((1-x)*(1-2*x)*(1+x+2*x^2)). - Richard Stanley, corrected Apr 29 2011
G.f.: (1 + sqrt( 1 + 4*x / ((1 - x) * (1 - 2*x) * (1 + x + 2*x^2)))) / (2*(1 - x)). - Michael Somos, Jan 30 2012
a(n) = sum( binomial(2*k-1, k)*binomial(n-2*k,k) + binomial(2*k, k)*binomial(n-2*k-1, k), k=0..floor(n/3)). - Joel B. Lewis, May 21 2011
Conjecture: -n*a(n) +(2+n)*a(n-1) +(3n-12)*a(n-2) +(12-n)*a(n-3) +(2n-18)*a(n-4)+(56-12n)*a(n-5) +(8n-40)*a(n-6)=0. - R. J. Mathar, Nov 28 2011
G.f. y = A(x) satisfies x = (1 - x) * (1 - 2*x) * (1 + x + 2*x^2) * y * (y * (1 - x) - 1). - Michael Somos, Jan 30 2012
Sequence a(n) satisfies 0 = a(n) * (n^2-2*n) + a(n-1) * (-3*n^2+8*n-2) + a(n-2) * (3*n^2-10*n+2) + a(n-3) * (-5*n^2+18*n-6) + a(n-4) * (8*n^2-34*n+22) + a(n-5) * (-4*n^2+20*n-16) except if n=1 or n=2. - Michael Somos, Jan 30 2012
a(n) = (1 + 3*hypergeom([1/2, 1-3*n/8, (1-n)/3, (2-n)/3, -n/3],[1, (1-n)/2, 1-n/2, -3*n/8],-27))/2 for n > 0. - Stefano Spezia, Apr 26 2024
a(n) ~ 2^n / sqrt(Pi*n). - Vaclav Kotesovec, Apr 26 2024

A351005 Number of integer partitions of n into parts that are alternately equal and unequal.

Original entry on oeis.org

1, 1, 2, 1, 2, 2, 3, 2, 4, 2, 5, 4, 6, 4, 8, 5, 10, 6, 12, 8, 16, 9, 18, 12, 22, 14, 28, 16, 33, 20, 40, 24, 48, 28, 56, 34, 67, 40, 80, 46, 94, 56, 110, 64, 130, 75, 152, 88, 176, 102, 206, 118, 238, 138, 276, 159, 320, 182, 368, 210, 424, 242, 488, 276, 558
Offset: 0

Views

Author

Gus Wiseman, Jan 31 2022

Keywords

Comments

Also partitions whose multiplicities are all 2's, except possibly for the last, which may be 1.

Examples

			The a(1) = 1 through a(12) = 6 partitions (A..C = 10..12):
  1  2   3  4   5    6     7    8     9    A     B      C
     11     22  221  33    331  44    441  55    443    66
                     2211       332        442   551    552
                                3311       3322  33221  4422
                                           4411         5511
                                                        332211
		

Crossrefs

The even-length ordered version is A003242, ranked by A351010.
The even-length case is A035457.
Without equalities we have A122135, opposite A122129, even-length A122134.
The non-strict version is A351004, opposite A351003, even-length A035363.
The opposite version is A351006, even-length A351007.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],And@@Table[#[[i]]==#[[i+1]],{i,1,Length[#]-1,2}]&&And@@Table[#[[i]]!=#[[i+1]],{i,2,Length[#]-1,2}]&]],{n,0,30}]
Previous Showing 11-20 of 196 results. Next