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

A005581 a(n) = (n-1)*n*(n+4)/6.

Original entry on oeis.org

0, 0, 2, 7, 16, 30, 50, 77, 112, 156, 210, 275, 352, 442, 546, 665, 800, 952, 1122, 1311, 1520, 1750, 2002, 2277, 2576, 2900, 3250, 3627, 4032, 4466, 4930, 5425, 5952, 6512, 7106, 7735, 8400, 9102, 9842, 10621, 11440, 12300, 13202, 14147, 15136, 16170
Offset: 0

Views

Author

Keywords

Comments

A class of Boolean functions of n variables and rank 2.
Also, number of inscribable triangles within a (n+4)-gon sharing with them its vertices but not its sides. - Lekraj Beedassy, Nov 14 2003
a(n) = A111808(n,3) for n > 2. - Reinhard Zumkeller, Aug 17 2005
If X is an n-set and Y a fixed 2-subset of X then a(n-2) is equal to the number of (n-3)-subsets of X intersecting Y. - Milan Janjic, Jul 30 2007
The sequence starting with offset 2 = binomial transform of [2, 5, 4, 1, 0, 0, 0, ...]. - Gary W. Adamson, Mar 20 2009
Let I=I_n be the n X n identity matrix and P=P_n be the incidence matrix of the cycle (1,2,3,...,n). Then, for n >= 4, a(n-4) is the number of (0,1) n X n matrices A <= P^(-1) + I + P having exactly two 1's in every row and column with perA=8. - Vladimir Shevelev, Apr 12 2010
Also arises as the number of triples of edges which can be chosen as the cut-points in the "three-opt" heuristic for a traveling salesman problem on (n+4) nodes. - James McDermott, Jul 10 2015
a(n) = risefac(n, 3)/3! - n is for n >= 1 also the number of independent components of a symmetric traceless tensor of rank 3 and dimension n. Here risefac is the rising factorial. - Wolfdieter Lang, Dec 10 2015
For n >= 2, a(n) is the number of characters in a word Q formed by concatenating all 'directed' ( left to right or vice versa), unrearranged subwords, from length 1 to (n-1), of a length (n-1) word q- allowing for the appearance of repeated subwords- and simply inserting an extra character for all subwords thus concatenated. - Christopher Hohl, May 30 2019

Examples

			In hexagon ABCDEF, the "interior" triangles are ACE and BDF, and a(6-4)=a(2)=2. - _Toby Gottfried_, Nov 12 2011
G.f. = 2*x^2 + 7*x^3 + 16*x^4 + 30*x^5 + 50*x^6 + 77*x^7 + 112*x^8 + ...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), Table 22.7, p. 797.
  • Joseph D. Konhauser, Dan Velleman and Stan Wagon,, Which Way Did the Bicycle Go?, MAA, 1996, p. 177.
  • V. S. Shevelyov (Shevelev), Extension of the Moser class of four-line Latin rectangles, DAN Ukrainy, Vol. 3 (1992), pp. 15-19. - Vladimir Shevelev, Apr 12 2010
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • A. M. Yaglom and I. M. Yaglom, Challenging Mathematical Problems with Elementary Solutions. Vol. I. Combinatorial Analysis and Probability Theory. New York: Dover Publications, Inc., 1987, p. 13, #51 (the case k=3) (First published: San Francisco: Holden-Day, Inc., 1964).

Crossrefs

Programs

Formula

G.f.: (x^2)*(2-x)/(1-x)^4.
a(n) = binomial(n+1, n-2) + binomial(n, n-2).
a(n) = A027907(n, 3), n >= 0 (fourth column of trinomial coefficients). - N. J. A. Sloane, May 16 2003
Convolution of {1, 2, 3, ...} with {2, 3, 4, ...}. - Jon Perry, Jun 25 2003
a(n+2) = 2*te(n) - te(n-1), e.g., a(5) = 2*te(3) - te(2) = 2*20 - 10 = 30, where te(n) are the tetrahedral numbers A000292. - Jon Perry, Jul 23 2003
a(n) is the coefficient of x^3 in the expansion of (1+x+x^2)^n. For example, a(1)=0 since (1+x+x^2)^1=1+x+x^2. - Peter C. Heinig (algorithms(AT)gmx.de), Apr 09 2007
E.g.f.: (x^2 + x^3/6) * exp(x). - Michael Somos, Apr 13 2007
a(n) = - A005586(-4-n) for all n in Z. - Michael Somos, Apr 13 2007
a(n) = C(4+n,3)-(n+4)*(n+1), since C(4+n,3) = number of all triangles in (n+4)-gon, and (n+4)*(n+1)=number of triangles with at least one of the edges included. Example: n=0,in a square, all 4 possible triangles include some of the square's edges and C(4+n,3)-(n+4)*(n+1)=4-4*1=0 = number of other triangles = a(0). - Toby Gottfried, Nov 12 2011
a(n) = 2*binomial(n,2) + binomial(n,3). - Vladimir Shevelev and Peter J. C. Moses, Jun 22 2012
a(0)=0, a(1)=0, a(2)=2, a(3)=7, a(n)=4*a(n-1)-6*a(n-2)+4*a(n-3)-a(n-4). - Harvey P. Dale, Sep 22 2012
a(n) = A000292(n-1) + A000217(n-1) for all n in Z. - Michael Somos, Jul 29 2015
a(n+2) = -A127672(6+n, n), n >= 0, with A127672 giving the coefficients of Chebyshev's C polynomials. See the Abramowitz-Stegun reference. - Wolfdieter Lang, Dec 10 2015
a(n) = GegenbauerC(N, -n, -1/2) where N = 3 if 3Peter Luschny, May 10 2016
From Amiram Eldar, Jan 09 2022: (Start)
Sum_{n>=2} 1/a(n) = 163/200.
Sum_{n>=2} (-1)^n/a(n) = 12*log(2)/5 - 253/200. (End)

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Jun 01 2000

A226290 Irregular triangle read by rows: T(n,k) is the number of binary pattern classes in the (3,n)-rectangular grid with k '1's and (3n-k) '0's: two patterns are in same class if one of them can be obtained by a reflection or 180-degree rotation of the other.

Original entry on oeis.org

1, 1, 2, 2, 1, 1, 2, 6, 6, 6, 2, 1, 1, 4, 13, 27, 39, 39, 27, 13, 4, 1, 1, 4, 22, 60, 139, 208, 252, 208, 139, 60, 22, 4, 1, 1, 6, 34, 129, 371, 794, 1310, 1675, 1675, 1310, 794, 371, 129, 34, 6, 1, 1, 6, 48, 218, 813, 2196, 4767, 8070, 11139, 12300, 11139, 8070, 4767, 2196, 813, 218, 48, 6, 1
Offset: 0

Views

Author

Yosu Yurramendi, Jun 02 2013

Keywords

Comments

Sum of rows (see example) gives A225827.
This triangle is to A225827 as Losanitsch's triangle A034851 is to A005418, and triangle A226048 to A225826.
By columns:
T(n,1) is A052928.
T(n,2) is A226292.
Also the number of equivalence classes of ways of placing k 1 X 1 tiles in an n X 3 rectangle under all symmetry operations of the rectangle. - Christopher Hunt Gribble, Apr 24 2015

Examples

			n\k 0 1  2   3   4    5    6    7     8     9    10   11   12
0   1
1   1 2  2   1
2   1 2  6   6   6    2    1
3   1 4 13  27  39   39   27   13     4     1
4   1 4 22  60 139  208  252  208   139    60    22    4    1
5   1 6 34 129 371  794 1310 1675  1675  1310   794  371  129    34   6   1
6   1 6 48 218 813 2196 4767 8070 11139 12300 11139 8070 4767  2196 813 218 48 6 1
...
The length of row n is 3*n+1.
		

Crossrefs

Programs

  • Mathematica
    T[n_, k_] := (Binomial[3n, k] + If[OddQ[n] || EvenQ[k], Binomial[Quotient[3 n, 2], Quotient[k, 2]], 0] + Sum[Binomial[n, k - 2i] Binomial[n, i] + Binomial[3 Mod[n, 2], k - 2i] Binomial[3 Quotient[n, 2], i], {i, 0, Quotient[k, 2]}])/4; Table[T[n, k], {n, 0, 6}, {k, 0, 3n}] // Flatten (* Jean-François Alcover, Oct 06 2017, after Andrew Howroyd *)
  • PARI
    T(n,k)={(binomial(3*n,k) + if(n%2==1||k%2==0,binomial(3*n\2,k\2),0) + sum(i=0,k\2, binomial(n,k-2*i) * binomial(1*n,i) + binomial(3*(n%2),k-2*i) * binomial(3*(n\2),i)))/4}
    for(n=0,6,for(k=0,3*n, print1(T(n,k), ", "));print) \\ Andrew Howroyd, May 30 2017

Extensions

Definition corrected by María Merino, May 19 2017

A004524 Three even followed by one odd.

Original entry on oeis.org

0, 0, 0, 1, 2, 2, 2, 3, 4, 4, 4, 5, 6, 6, 6, 7, 8, 8, 8, 9, 10, 10, 10, 11, 12, 12, 12, 13, 14, 14, 14, 15, 16, 16, 16, 17, 18, 18, 18, 19, 20, 20, 20, 21, 22, 22, 22, 23, 24, 24, 24, 25, 26, 26, 26, 27, 28, 28, 28, 29, 30, 30, 30, 31, 32, 32, 32, 33, 34, 34, 34, 35, 36, 36, 36, 37
Offset: 0

Views

Author

Keywords

Comments

Ignoring the first term, for n >= 0, n/2 rounded by the method called "banker's rounding", "statistician's rounding", or "round-to-even" gives 0, 0, 1, 2, 2, 2, 3, ..., where this method rounds k + 0.5 to k if positive integer k is even but rounds k + 0.5 to k + 1 when k + 1 is even. (If the method is indeed defined such that the above statement is also true with the word "positive" removed, then the first 0 term need not be ignored and this sequence can be further extended symmetrically with a(m) = -a(-m) for all integers m, an advantage over usual rounding.) The corresponding sequence for n/2 rounded by the common method is A004526 (considered as beginning with n = -1). - Rick L. Shepherd, Nov 16 2006
From Anthony Hernandez, Aug 08 2016: (Start)
Arrange the positive integers starting at 1 into a triangular array
1
2 3
4 5 6
7 8 9 10
11 12 13 14 15
16 17 18 19 20 21
22 23 24 25 26 27 28
29 30 31 32 33 34 35 36
and let e(n) count the even numbers in the n-th row of the array. Then e(n) = a(n+1). For example, e(6) = a(7) = 3 and there are three even numbers in the 6th row of the array. For the count of odd numbers, f(n), look at the sequence A004525. (End)
Also the domination number of the (n-1) X (n-1) white bishop graph. - Eric W. Weisstein, Jun 26 2017
Let (b(n)) be the p-INVERT of A010892 using p(S) = 1 - S^2; then b(n) = a(n+1) for n >= 0. See A292301. - Clark Kimberling, Sep 30 2017
Also the total domination number of the (n-2)-complete graph (for n>3), (n-2)-cycle graph (for n>4), and (n-2)-pan graph (for n>4). - Eric W. Weisstein, Apr 07 2018
The sequence is the interleaving of the duplicated even integers (A052928) with the nonnegative integers (A001477). - Guenther Schrack, Mar 05 2019

Examples

			G.f. = x^3 + 2*x^4 + 2*x^5 + 2*x^6 + 3*x^7 + 4*x^8 + 4*x^9 + 4*x^10 + ...
		

Crossrefs

Zero followed by partial sums of A021913.
First differences of A011848.

Programs

  • GAP
    List([0..79],n->Int(n/4)+Int((n+1)/4)); # Muniru A Asiru, Mar 06 2019
    
  • Haskell
    a004524 n = n `div` 4 + (n + 1) `div` 4
    a004524_list = 0 : 0 : 0 : 1 : map (+ 2) a004524_list
    -- Reinhard Zumkeller, Feb 22 2013, Jul 14 2012
    
  • Magma
    [Floor(n/4)+Floor((n+1)/4) : n in [0..80]]; // Wesley Ivan Hurt, Jul 21 2014
    
  • Maple
    A004524:=n->floor(n/4)+floor((n+1)/4): seq(A004524(n), n=0..50); # Wesley Ivan Hurt, Jul 21 2014
  • Mathematica
    Table[Floor[n/4] + Floor[(n + 1)/4], {n, 0, 80}] (* Wesley Ivan Hurt, Jul 21 2014 *)
    Flatten[Table[{n, n, n, n + 1}, {n, 0, 38, 2}]] (* Alonso del Arte, Aug 10 2016 *)
    Table[(n + Cos[n Pi/2] - 1)/2, {n, 0, 80}] (* Eric W. Weisstein, Apr 07 2018 *)
    Table[Floor[n/2 - 1] + Ceiling[n/4 - 1/2] - Floor[n/4 - 1/2], {n, 0, 80}] (* Eric W. Weisstein, Apr 07 2018 *)
    LinearRecurrence[{2, -2, 2, -1}, {0, 0, 1, 2}, {0, 80}] (* Eric W. Weisstein, Apr 07 2018 *)
    CoefficientList[Series[x^3/((1 - x)^2 (1 + x^2)), {x, 0, 80}], x] (* Eric W. Weisstein, Apr 07 2018 *)
    Table[Round[(n - 1)/2], {n, 0, 20}] (* Eric W. Weisstein, Jun 19 2024 *)
    Round[(Range[0, 20] - 1)/2] (* Eric W. Weisstein, Jun 19 2024 *)
    Table[PadRight[{},If[EvenQ[n],3,1],n],{n,0,40}]//Flatten (* Harvey P. Dale, Dec 11 2024 *)
  • PARI
    {a(n) = n\4 + (n+1)\4}; /* Michael Somos, Jul 19 2003 */
    
  • PARI
    concat([0,0,0], Vec(x^3/((1-x)^2*(1+x^2)) + O(x^80))) \\ Altug Alkan, Oct 31 2015
    
  • Python
    def A004524(n): return (n>>2)+(n+1>>2) # Chai Wah Wu, Jul 29 2022
  • Sage
    [floor(n/4)+floor((n+1)/4) for n in (0..80)] # G. C. Greubel, Mar 08 2019
    

Formula

a(n) = a(n-1) - a(n-2) + a(n-3) + 1 = (n-1) - A004525(n-1). - Henry Bottomley, Mar 08 2000
G.f.: x^3/((1 - x)^2*(1 + x^2)) = x^3*(1 - x^2)/((1 - x)^2*(1 - x^4)). - Michael Somos, Jul 19 2003
If the sequence is extended to negative arguments in the natural way, it satisfies a(n) = -a(2-n) for all n in Z. - Michael Somos, Jul 19 2003
a(n) = A092038(n-3) for n > 4. - Reinhard Zumkeller, Mar 28 2004
From Paul Barry, Oct 27 2004: (Start)
E.g.f.: (exp(x)*(x-1) + cos(x))/2.
a(n) = (n - 1 - cos(Pi*(n-2)/2))/2. (End)
a(n+3) = Sum_{k = 0..n} (1 + (-1)^C(n,2))/2. - Paul Barry, Mar 31 2008
a(n) = floor(n/4) + floor((n+1)/4). - Arkadiusz Wesolowski, Sep 19 2012
From Wesley Ivan Hurt, Jul 21 2014, Oct 31 2015: (Start)
a(n) = Sum_{i = 1..n-1} (floor(i/2) mod 2).
a(n) = n/2 - sqrt(n^2 mod 8)/2. (End)
Euler transform of length 4 sequence [2, -1, 0, 1]. - Michael Somos, Apr 03 2017
a(n) = (2*n - 2 + (1 + (-1)^n)*(-1)^(n*(n-1)/2))/4. - Guenther Schrack, Mar 04 2019
Sum_{n>=3} (-1)^(n+1)/a(n) = log(2) (A002162). - Amiram Eldar, Sep 29 2022

A047838 a(n) = floor(n^2/2) - 1.

Original entry on oeis.org

1, 3, 7, 11, 17, 23, 31, 39, 49, 59, 71, 83, 97, 111, 127, 143, 161, 179, 199, 219, 241, 263, 287, 311, 337, 363, 391, 419, 449, 479, 511, 543, 577, 611, 647, 683, 721, 759, 799, 839, 881, 923, 967, 1011, 1057, 1103, 1151, 1199, 1249, 1299, 1351, 1403
Offset: 2

Views

Author

Michael Somos, May 07 1999

Keywords

Comments

Define the organization number of a permutation pi_1, pi_2, ..., pi_n to be the following. Start at 1, count the steps to reach 2, then the steps to reach 3, etc. Add them up. Then the maximal value of the organization number of any permutation of [1..n] for n = 0, 1, 2, 3, ... is given by 0, 1, 3, 7, 11, 17, 23, ... (this sequence). This was established by Graham Cormode (graham(AT)research.att.com), Aug 17 2006, see link below, answering a question raised by Tom Young (mcgreg265(AT)msn.com) and Barry Cipra, Aug 15 2006
From Dmitry Kamenetsky, Nov 29 2006: (Start)
This is the length of the longest non-self-intersecting spiral drawn on an n X n grid. E.g., for n=5 the spiral has length 17:
1 0 1 1 1
1 0 1 0 1
1 0 1 0 1
1 0 0 0 1
1 1 1 1 1 (End)
It appears that a(n+1) is the maximum number of consecutive integers (beginning with 1) that can be placed, one after another, on an n-peg Towers of Hanoi, such that the sum of any two consecutive integers on any peg is a square. See the problem: http://online-judge.uva.es/p/v102/10276.html. - Ashutosh Mehra, Dec 06 2008
a(n) = number of (w,x,y) with all terms in {0,...,n} and w = |x+y-w|. - Clark Kimberling, Jun 11 2012
The same sequence also represents the solution to the "pigeons problem": maximal value of the sum of the lengths of n-1 line segments (connected at their end-points) required to pass through n trail dots, with unit distance between adjacent points, visiting all of them without overlaping two or more segments. In this case, a(0)=0, a(1)=1, a(2)=3, and so on. - Marco Ripà, Jan 28 2014
Also the longest path length in the n X n white bishop graph. - Eric W. Weisstein, Mar 27 2018
a(n) is the number of right triangles with sides n*(h-floor(h)), floor(h) and h, where h is the hypotenuse. - Andrzej Kukla, Apr 14 2021

Examples

			x^2 + 3*x^3 + 7*x^4 + 11*x^5 + 17*x^6 + 23*x^7 + 31*x^8 + 39*x^9 + 49*x^10 + ...
		

Crossrefs

Complement of A047839. First difference is A052928.
Partial sums: A213759(n-1) for n > 1. - Guenther Schrack, May 12 2018

Programs

  • Magma
    [Floor(n^2/2)-1 : n in [2..100]]; // Wesley Ivan Hurt, Aug 06 2015
  • Maple
    seq(floor((n^2+4*n+2)/2), n=0..20) # Gary Detlefs, Feb 10 2010
  • Mathematica
    Table[Floor[n^2/2] - 1, {n, 2, 60}] (* Robert G. Wilson v, Aug 31 2006 *)
    LinearRecurrence[{2, 0, -2, 1}, {1, 3, 7, 11}, 60] (* Harvey P. Dale, Jan 16 2015 *)
    Floor[Range[2, 20]^2/2] - 1 (* Eric W. Weisstein, Mar 27 2018 *)
    Table[((-1)^n + 2 n^2 - 5)/4, {n, 2, 20}] (* Eric W. Weisstein, Mar 27 2018 *)
    CoefficientList[Series[(-1 - x - x^2 + x^3)/((-1 + x)^3 (1 + x)), {x, 0, 20}], x] (* Eric W. Weisstein, Mar 27 2018 *)
  • PARI
    a(n) = n^2\2 - 1
    

Formula

a(2)=1; for n > 2, a(n) = a(n-1) + n - 1 + (n-1 mod 2). - Benoit Cloitre, Jan 12 2003
a(n) = T(n-1) + floor(n/2) - 1 = T(n) - floor((n+3)/2), where T(n) is the n-th triangular number (A000217). - Robert G. Wilson v, Aug 31 2006
Equals (n-1)-th row sums of triangles A134151 and A135152. Also, = binomial transform of [1, 2, 2, -2, 4, -8, 16, -32, ...]. - Gary W. Adamson, Nov 21 2007
G.f.: x^2*(1+x+x^2-x^3)/((1-x)^3*(1+x)). - R. J. Mathar, Sep 09 2008
a(n) = floor((n^2 + 4*n + 2)/2). - Gary Detlefs, Feb 10 2010
a(n) = abs(A188653(n)). - Reinhard Zumkeller, Apr 13 2011
a(n) = (2*n^2 + (-1)^n - 5)/4. - Bruno Berselli, Sep 14 2011
a(n) = a(-n) = A007590(n) - 1.
a(n) = A080827(n) - 2. - Kevin Ryde, Aug 24 2013
a(n) = 2*a(n-1) - 2*a(n-3) + a(n-4), n > 4. - Wesley Ivan Hurt, Aug 06 2015
a(n) = A000217(n-1) + A004526(n-2), for n > 1. - J. Stauduhar, Oct 20 2017
From Guenther Schrack, May 12 2018: (Start)
Set a(0) = a(1) = -1, a(n) = a(n-2) + 2*n - 2 for n > 1.
a(n) = A000982(n-1) + n - 2 for n > 1.
a(n) = 2*A033683(n) - 3 for n > 1.
a(n) = A061925(n-1) + n - 3 for n > 1.
a(n) = A074148(n) - n - 1 for n > 1.
a(n) = A105343(n-1) + n - 4 for n > 1.
a(n) = A116940(n-1) - n for n > 1.
a(n) = A179207(n) - n + 1 for n > 1.
a(n) = A183575(n-2) + 1 for n > 2.
a(n) = A265284(n-1) - 2*n + 1 for n > 1.
a(n) = 2*A290743(n) - 5 for n > 1. (End)
E.g.f.: 1 + x + ((x^2 + x - 2)*cosh(x) + (x^2 + x - 3)*sinh(x))/2. - Stefano Spezia, May 06 2021
Sum_{n>=2} 1/a(n) = 3/2 + tan(sqrt(3)*Pi/2)*Pi/(2*sqrt(3)) - cot(Pi/sqrt(2))*Pi/(2*sqrt(2)). - Amiram Eldar, Sep 15 2022

Extensions

Edited by Charles R Greathouse IV, Apr 23 2010

A064455 a(2n) = 3n, a(2n-1) = n.

Original entry on oeis.org

1, 3, 2, 6, 3, 9, 4, 12, 5, 15, 6, 18, 7, 21, 8, 24, 9, 27, 10, 30, 11, 33, 12, 36, 13, 39, 14, 42, 15, 45, 16, 48, 17, 51, 18, 54, 19, 57, 20, 60, 21, 63, 22, 66, 23, 69, 24, 72, 25, 75, 26, 78, 27, 81, 28, 84, 29, 87, 30, 90, 31, 93, 32, 96, 33, 99, 34, 102, 35, 105, 36, 108
Offset: 1

Views

Author

N. J. A. Sloane, Oct 02 2001

Keywords

Comments

Also number of 1's in n-th row of triangle in A071030. - Hans Havermann, May 26 2002
Number of ON cells at generation n of 1-D CA defined by Rule 54. - N. J. A. Sloane, Aug 09 2014
a(n)*A098557(n) equals the second right hand column of A167556. - Johannes W. Meijer, Nov 12 2009
Given a(1) = 1, for all n > 1, a(n) is the least positive integer not equal to a(n-1) such that the arithmetic mean of the first n terms is an integer. The sequence of arithmetic means of the first 1, 2, 3, ..., terms is 1, 2, 2, 3, 3, 4, 4, ... (A004526 disregarding its first three terms). - Rick L. Shepherd, Aug 20 2013

Examples

			a(13) = a(2*7 - 1) = 7, a(14) = a(2*7) = 21.
a(8) = 8-9+10-11+12-13+14-15+16 = 12. - _Bruno Berselli_, Jun 05 2013
		

Crossrefs

Interleaving of A000027 and A008585 (without first term).

Programs

  • ARIBAS
    maxarg := 75; for n := 1 to maxarg do if n mod 2 = 1 then write((n+1) div 2, " ") else write((n div 2)*3," "); end; end;
    
  • GAP
    a:=[];;  for n in [1..75] do if n mod 2 = 0 then Add(a,3*n/2); else Add(a,(n+1)/2); fi; od; a; # Muniru A Asiru, Oct 28 2018
    
  • Haskell
    import Data.List (transpose)
    a064455 n = n + if m == 0 then n' else - n'  where (n',m) = divMod n 2
    a064455_list = concat $ transpose [[1 ..], [3, 6 ..]]
    -- Reinhard Zumkeller, Oct 12 2013
    
  • Magma
    [(1/2)*n*(-1)^n+n+(1/4)*(1-(-1)^n): n in [1..80]]; // Vincenzo Librandi, Aug 10 2014
    
  • Maple
    A064455 := proc(n)
        if type(n,'even') then
            3*n/2 ;
        else
            (n+1)/2 ;
        end if;
    end proc: # R. J. Mathar, Aug 03 2015
  • Mathematica
    Table[ If[ EvenQ[n], 3n/2, (n + 1)/2], {n, 1, 70} ]
  • PARI
    a(n) = { if (n%2, (n + 1)/2, 3*n/2) } \\ Harry J. Smith, Sep 14 2009
    
  • PARI
    a(n)=if(n<3,2*n-1,((n-1)*(n-2))%(2*n-1)) \\ Jim Singh, Oct 14 2018
    
  • Python
    def A064455(n): return (3*n - (2*n-1)*(n%2))//2
    print([A064455(n) for n in range(1,81)]) # G. C. Greubel, Jan 30 2025

Formula

a(n) = (1/2)*n*(-1)^n + n + (1/4)*(1 - (-1)^n). - Stephen Crowley, Aug 10 2009
G.f.: x*(1+3*x) / ( (1-x)^2*(1+x)^2 ). - R. J. Mathar, Mar 30 2011
From Jaroslav Krizek, Mar 22 2011: (Start)
a(n) = n - A123684(n-1) for odd n.
a(n) = n + a(n-1) for even n.
a(n) = A123684(n) + A137501(n).
Abs( a(n) - A123684(n) ) = A052928(n). (End)
a(n) = Sum_{i=n..2*n} i*(-1)^i. - Bruno Berselli, Jun 05 2013
a(n) = n + floor(n/2)*(-1)^(n mod 2). - Bruno Berselli, Dec 14 2015
a(n) = (n^2-3*n+2) mod (2*n-1) for n>2. - Jim Singh, Oct 31 2018
E.g.f.: (1/2)*(x*cosh(x) + (1+3*x)*sinh(x)). - G. C. Greubel, Jan 30 2025

A072574 Triangle T(n,k) of number of compositions (ordered partitions) of n into exactly k distinct parts, 1<=k<=n.

Original entry on oeis.org

1, 1, 0, 1, 2, 0, 1, 2, 0, 0, 1, 4, 0, 0, 0, 1, 4, 6, 0, 0, 0, 1, 6, 6, 0, 0, 0, 0, 1, 6, 12, 0, 0, 0, 0, 0, 1, 8, 18, 0, 0, 0, 0, 0, 0, 1, 8, 24, 24, 0, 0, 0, 0, 0, 0, 1, 10, 30, 24, 0, 0, 0, 0, 0, 0, 0, 1, 10, 42, 48, 0, 0, 0, 0, 0, 0, 0, 0, 1, 12, 48, 72, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 12, 60, 120, 0
Offset: 1

Views

Author

Henry Bottomley, Jun 21 2002

Keywords

Comments

If terms in the compositions did not need to be distinct then the triangle would have values C(n-1,k-1), essentially A007318 offset.

Examples

			T(6,2)=4 since 6 can be written as 1+5=2+4=4+2=5+1.
Triangle starts (trailing zeros omitted for n>=10):
[ 1]  1;
[ 2]  1, 0;
[ 3]  1, 2, 0;
[ 4]  1, 2, 0, 0;
[ 5]  1, 4, 0, 0, 0;
[ 6]  1, 4, 6, 0, 0, 0;
[ 7]  1, 6, 6, 0, 0, 0, 0;
[ 8]  1, 6, 12, 0, 0, 0, 0, 0;
[ 9]  1, 8, 18, 0, 0, 0, 0, 0, 0;
[10]  1, 8, 24, 24, 0, 0, ...;
[11]  1, 10, 30, 24, 0, 0, ...;
[12]  1, 10, 42, 48, 0, 0, ...;
[13]  1, 12, 48, 72, 0, 0, ...;
[14]  1, 12, 60, 120, 0, 0, ...;
[15]  1, 14, 72, 144, 120, 0, 0, ...;
[16]  1, 14, 84, 216, 120, 0, 0, ...;
[17]  1, 16, 96, 264, 240, 0, 0, ...;
[18]  1, 16, 114, 360, 360, 0, 0, ...;
[19]  1, 18, 126, 432, 600, 0, 0, ...;
[20]  1, 18, 144, 552, 840, 0, 0, ...;
These rows (without the zeros) are shown in the Richmond/Knopfmacher reference.
From _Gus Wiseman_, Oct 17 2022: (Start)
Column n = 8 counts the following compositions.
  (8)  (1,7)  (1,2,5)
       (2,6)  (1,3,4)
       (3,5)  (1,4,3)
       (5,3)  (1,5,2)
       (6,2)  (2,1,5)
       (7,1)  (2,5,1)
              (3,1,4)
              (3,4,1)
              (4,1,3)
              (4,3,1)
              (5,1,2)
              (5,2,1)
(End)
		

Crossrefs

Columns (offset) include A057427 and A052928.
Row sums are A032020.
A008289 is the version for partitions (zeros removed).
A072575 counts strict compositions by maximum.
A097805 is the non-strict version, or A007318 (zeros removed).
A113704 is the constant instead of strict version.
A216652 is a condensed version (zeros removed).
A336131 counts splittings of partitions with distinct sums.
A336139 counts strict compositions of each part of a strict composition.

Programs

  • Mathematica
    Table[Length[Select[Join@@Permutations/@Select[IntegerPartitions[n],UnsameQ@@#&],Length[#]==k&]],{n,0,15},{k,1,n}] (* Gus Wiseman, Oct 17 2022 *)
  • PARI
    N=21;  q='q+O('q^N);
    gf=sum(n=0,N, n! * z^n * q^((n^2+n)/2) / prod(k=1,n, 1-q^k ) );
    /* print triangle: */
    gf -= 1; /* remove row zero */
    P=Pol(gf,'q);
    { for (n=1,N-1,
        p = Pol(polcoeff(P, n),'z);
        p += 'z^(n+1);  /* preserve trailing zeros */
        v = Vec(polrecip(p));
        v = vector(n,k,v[k]); /* trim to size n */
        print(v);
    ); }
    /* Joerg Arndt, Oct 20 2012 */

Formula

T(n, k) = T(n-k, k)+k*T(n-k, k-1) [with T(n, 0)=1 if n=0 and 0 otherwise] = A000142(k)*A060016(n, k).
G.f.: sum(n>=0, n! * z^n * q^((n^2+n)/2) / prod(k=1..n, 1-q^k ) ), rows by powers of q, columns by powers of z; includes row 0 (drop term for n=0 for this triangle, see PARI code); setting z=1 gives g.f. for A032020. [Joerg Arndt, Oct 20 2012]

A325677 Irregular triangle read by rows where T(n,k) is the number of Golomb rulers of length n with k + 1 marks, k > 0.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 4, 1, 4, 2, 1, 6, 6, 1, 6, 8, 1, 8, 18, 1, 8, 16, 1, 10, 30, 4, 1, 10, 34, 14, 1, 12, 48, 28, 1, 12, 48, 42, 1, 14, 72, 76, 1, 14, 72, 100, 1, 16, 96, 160, 8, 1, 16, 98, 190, 8, 1, 18, 126, 284, 40, 1, 18, 128, 316, 70
Offset: 1

Views

Author

Gus Wiseman, May 13 2019

Keywords

Comments

Also the number of length-k compositions of n such that every restriction to a subinterval has a different sum. A composition of n is a finite sequence of positive integers summing to n.

Examples

			Triangle begins:
   1
   1
   1   2
   1   2
   1   4
   1   4   2
   1   6   6
   1   6   8
   1   8  18
   1   8  16
   1  10  30   4
   1  10  34  14
   1  12  48  28
   1  12  48  42
   1  14  72  76
   1  14  72 100
   1  16  96 160   8
   1  16  98 190   8
   1  18 126 284  40
   1  18 128 316  70
Row n = 8 counts the following rulers:
  {0,8}  {0,1,8}  {0,1,3,8}
         {0,2,8}  {0,1,5,8}
         {0,3,8}  {0,1,6,8}
         {0,5,8}  {0,2,3,8}
         {0,6,8}  {0,2,7,8}
         {0,7,8}  {0,3,7,8}
                  {0,5,6,8}
                  {0,5,7,8}
and the following compositions:
  (8)  (17)  (125)
       (26)  (143)
       (35)  (152)
       (53)  (215)
       (62)  (251)
       (71)  (341)
             (512)
             (521)
		

Crossrefs

Row sums are A169942.
Row lengths are A325678(n) = A143824(n + 1) - 1.
Column k = 2 is A052928.
Column k = 3 is A325686.
Rightmost column is A325683.

Programs

  • Mathematica
    DeleteCases[Table[Length[Select[Join@@Permutations/@IntegerPartitions[n,{k}],UnsameQ@@ReplaceList[#,{_,s__,_}:>Plus[s]]&]],{n,15},{k,n}],0,{2}]

A055884 Euler transform of partition triangle A008284.

Original entry on oeis.org

1, 1, 2, 1, 2, 3, 1, 4, 4, 5, 1, 4, 8, 7, 7, 1, 6, 12, 16, 12, 11, 1, 6, 17, 25, 28, 19, 15, 1, 8, 22, 43, 49, 48, 30, 22, 1, 8, 30, 58, 87, 88, 77, 45, 30, 1, 10, 36, 87, 134, 167, 151, 122, 67, 42, 1, 10, 45, 113, 207, 270, 296, 247, 185, 97, 56, 1, 12, 54, 155, 295, 448, 510, 507, 394, 278, 139, 77
Offset: 1

Views

Author

Christian G. Bower, Jun 09 2000

Keywords

Comments

Number of multiset partitions of length-k integer partitions of n. - Gus Wiseman, Nov 09 2018

Examples

			From _Gus Wiseman_, Nov 09 2018: (Start)
Triangle begins:
   1
   1   2
   1   2   3
   1   4   4   5
   1   4   8   7   7
   1   6  12  16  12  11
   1   6  17  25  28  19  15
   1   8  22  43  49  48  30  22
   1   8  30  58  87  88  77  45  30
   ...
The fifth row {1, 4, 8, 7, 7} counts the following multiset partitions:
  {{5}}   {{1,4}}     {{1,1,3}}       {{1,1,1,2}}         {{1,1,1,1,1}}
          {{2,3}}     {{1,2,2}}      {{1},{1,1,2}}       {{1},{1,1,1,1}}
         {{1},{4}}   {{1},{1,3}}     {{1,1},{1,2}}       {{1,1},{1,1,1}}
         {{2},{3}}   {{1},{2,2}}     {{2},{1,1,1}}      {{1},{1},{1,1,1}}
                     {{2},{1,2}}    {{1},{1},{1,2}}     {{1},{1,1},{1,1}}
                     {{3},{1,1}}    {{1},{2},{1,1}}    {{1},{1},{1},{1,1}}
                    {{1},{1},{3}}  {{1},{1},{1},{2}}  {{1},{1},{1},{1},{1}}
                    {{1},{2},{2}}
(End)
		

Crossrefs

Row sums give A001970.
Main diagonal gives A000041.
Columns k=1-2 give: A057427, A052928.
T(n+2,n+1) gives A000070.
T(2n,n) gives A360468.

Programs

  • Maple
    h:= proc(n, i) option remember; expand(`if`(n=0, 1,
          `if`(i<1, 0, h(n, i-1)+x*h(n-i, min(n-i, i)))))
        end:
    g:= proc(n, i, j) option remember; expand(`if`(j=0, 1, `if`(i<0, 0, add(
          g(n, i-1, j-k)*x^(i*k)*binomial(coeff(h(n$2), x, i)+k-1, k), k=0..j))))
        end:
    b:= proc(n, i) option remember; expand(`if`(n=0, 1,
         `if`(i<1, 0, add(b(n-i*j, i-1)*g(i$2, j), j=0..n/i))))
        end:
    T:= (n, k)-> coeff(b(n$2), x, k):
    seq(seq(T(n,k), k=1..n), n=1..12);  # Alois P. Heinz, Feb 17 2023
  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    Table[Length[Join@@mps/@IntegerPartitions[n,{k}]],{n,5},{k,n}] (* Gus Wiseman, Nov 09 2018 *)

A099392 a(n) = floor((n^2 - 2*n + 3)/2).

Original entry on oeis.org

1, 1, 3, 5, 9, 13, 19, 25, 33, 41, 51, 61, 73, 85, 99, 113, 129, 145, 163, 181, 201, 221, 243, 265, 289, 313, 339, 365, 393, 421, 451, 481, 513, 545, 579, 613, 649, 685, 723, 761, 801, 841, 883, 925, 969, 1013, 1059, 1105, 1153, 1201, 1251, 1301, 1353, 1405
Offset: 1

Views

Author

Ralf Stephan following a suggestion from Luke Pebody, Oct 20 2004

Keywords

Crossrefs

Differs from A085913 at n = 61. Apart from leading term, identical to A080827.
Cf. A000217, A001844, A002522, A007494, A007590, A058331 (bisections).
From Guenther Schrack, Apr 17 2018: (Start)
First differences: A052928.
Partial sums: A212964(n) + n for n > 0.
Also A058331 and A001844 interleaved. (End)

Programs

  • Mathematica
    Array[Floor[(#^2 - 2 # + 3)/2] &, 54] (* or *)
    Rest@ CoefficientList[Series[x (-1 + x - x^2 - x^3)/((1 + x) (x - 1)^3), {x, 0, 54}], x] (* Michael De Vlieger, Apr 21 2018 *)
  • PARI
    a(n)=(n^2+3)\2-n \\ Charles R Greathouse IV, Aug 01 2013

Formula

a(n) = ceiling(n^2/2)-n+1. - Paul Barry, Jul 16 2006; index shifted by R. J. Mathar, Jul 29 2007
a(n) = ceiling(A002522(n-1)/2). - Branko Curgus, Sep 02 2007
From R. J. Mathar, Feb 20 2011: (Start)
G.f.: x *( -1+x-x^2-x^3 ) / ( (1+x)*(x-1)^3 ).
a(n) = 2*a(n-1) - 2*a(n-3) + a(n-4).
a(n+1) = (3 + 2*n^2 + (-1)^n)/4. (End)
a(n) = A007590(n-1) + 1 for n >= 2. - Richard R. Forberg, Aug 01 2013
a(n) = A000217(n) - A007494(n-1). - Bui Quang Tuan, Mar 27 2015
From Guenther Schrack, Apr 17 2018: (Start)
a(n) = (2*n^2 - 4*n + 5 -(-1)^n)/4.
a(n+2) = a(n) + 2*n for n > 0.
a(n) = 2*A033683(n-1) - 1 for n > 0.
a(n) = A047838(n-1) + 2 for n > 2.
a(n) = A074148(n-1) - n + 2 for n > 1.
a(n) = A183575(n-3) + 3 for n > 3.
a(n) = 2*A290743(n-1) - 3 for n > 0.
a(n) = 2*A290743(n-2) + A109613(n-5) for n > 4.
a(n) = A074148(n) - A014601(n-1) for n > 0. (End)
Sum_{n>=1} 1/a(n) = tanh(Pi/2)*Pi/2 + coth(Pi/sqrt(2))*Pi/(2*sqrt(2)) + 1/2. - Amiram Eldar, Sep 16 2022
E.g.f.: ((2 - x + x^2)*cosh(x) + (3 - x + x^2)*sinh(x) - 2)/2. - Stefano Spezia, Jan 28 2024

A091298 Triangle read by rows: T(n,k) is the number of plane partitions of n containing exactly k parts.

Original entry on oeis.org

1, 1, 2, 1, 2, 3, 1, 4, 3, 5, 1, 4, 7, 5, 7, 1, 6, 10, 13, 7, 11, 1, 6, 14, 20, 19, 11, 15, 1, 8, 18, 33, 32, 31, 15, 22, 1, 8, 25, 43, 56, 54, 43, 22, 30, 1, 10, 29, 66, 81, 99, 78, 64, 30, 42, 1, 10, 37, 83, 126, 150, 148, 118, 88, 42, 56, 1, 12, 44, 114, 174, 246, 235, 230, 166, 124, 56, 77
Offset: 1

Views

Author

Wouter Meeussen, Feb 24 2004

Keywords

Comments

First column is 1, representing the single-part {{n}}, last column is P(n), since the all-ones plane partitions form the Ferrers-Young plots of the (linear) partitions of n.
A plane partition of n is a two-dimensional table (or matrix) with nonnegative elements summing up to n, and nonincreasing rows and columns. (Zero rows and columns are ignored.) - M. F. Hasler, Sep 22 2018

Examples

			This plane partition of n=7: {{3,1,1},{2}} contains 4 parts: 3,1,1,2.
Triangle T(n,k) begins:
  1;
  1,  2;
  1,  2,  3;
  1,  4,  3,  5;
  1,  4,  7,  5,  7;
  1,  6, 10, 13,  7, 11;
  1,  6, 14, 20, 19, 11, 15;
  1,  8, 18, 33, 32, 31, 15, 22;
  1,  8, 25, 43, 56, 54, 43, 22, 30;
  1, 10, 29, 66, 81, 99, 78, 64, 30, 42;
  ...
		

Crossrefs

Row sums give A000219.
Column 1 is A000012. Column 2 is A052928. Diagonal and subdiagonal are A000041.

Programs

  • Mathematica
    (* see A089924 for "planepartition" *) Table[Length /@ Split[Sort[Length /@ Flatten /@ planepartitions[n]]], {n, 16}]
  • PARI
    A091298(n,k)=sum(i=1,#n=PlanePartitions(n),sum(j=1,#n[i],#n[i][j])==k)
    PlanePartitions(n,L=0,PP=List())={ n<2&&return([if(n,[[1]],[])]); for(N=1,n, my(P=apply(Vecrev, if(L, select(p->vecmin(L-Vecrev(p,#L))>=0, partitions(N,L[1],#L)), partitions(N)))); if(NM. F. Hasler, Sep 24 2018
Previous Showing 11-20 of 73 results. Next