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 21-30 of 263 results. Next

A034856 a(n) = binomial(n+1, 2) + n - 1 = n*(n+3)/2 - 1.

Original entry on oeis.org

1, 4, 8, 13, 19, 26, 34, 43, 53, 64, 76, 89, 103, 118, 134, 151, 169, 188, 208, 229, 251, 274, 298, 323, 349, 376, 404, 433, 463, 494, 526, 559, 593, 628, 664, 701, 739, 778, 818, 859, 901, 944, 988, 1033, 1079, 1126, 1174, 1223, 1273, 1324, 1376, 1429, 1483
Offset: 1

Views

Author

Keywords

Comments

Number of 1's in the n X n lower Hessenberg (0,1)-matrix (i.e., the matrix having 1's on or below the superdiagonal and 0's above the superdiagonal).
If a 2-set Y and 2-set Z, having one element in common, are subsets of an n-set X then a(n-2) is the number of 3-subsets of X intersecting both Y and Z. - Milan Janjic, Oct 03 2007
Number of binary operations which have to be added to Moisil's algebras to obtain algebraic counterparts of n-valued Łukasiewicz logics. See the Wójcicki and Malinowski book, page 31. - Artur Jasinski, Feb 25 2010
Also (n + 1)!(-1)^(n + 1) times the determinant of the n X n matrix given by m(i,j) = i/(i+1) if i=j and otherwise 1. For example, (5+1)! * ((-1)^(5+1)) * Det[{{1/2, 1, 1, 1, 1}, {1, 2/3, 1, 1, 1}, {1, 1, 3/4, 1, 1}, {1, 1, 1, 4/5, 1}, {1, 1, 1, 1, 5/6}}] = 19 = a(5), and (6+1)! * ((-1)^(6+1)) * Det[{{1/2, 1, 1, 1, 1, 1}, {1, 2/3, 1, 1, 1, 1}, {1, 1, 3/4, 1, 1, 1}, {1, 1, 1, 4/5, 1, 1}, {1, 1, 1, 1, 5/6, 1}, {1, 1, 1, 1, 1, 6/7}}] = 26 = a(6). - John M. Campbell, May 20 2011
2*a(n-1) = n*(n+1) - 4, n>=0, with a(-1) = -2 and a(0) = -1, gives the values for a*c of indefinite binary quadratic forms [a, b, c] of discriminant D = 17 for b = 2*n + 1. In general D = b^2 - 4*a*c > 0 and the form [a, b, c] is a*x^2 + b*x*y + c*y^2. - Wolfdieter Lang, Aug 15 2013
a(n) is not divisible by 3, 5, 7, or 11. - Vladimir Shevelev, Feb 03 2014
With a(0) = 1 and a(1) = 2, a(n-1) is the number of distinct values of 1 +- 2 +- 3 +- ... +- n, for n > 0. - Derek Orr, Mar 11 2015
Also, numbers m such that 8*m+17 is a square. - Bruno Berselli, Sep 16 2015
Omar E. Pol's formula from Apr 23 2008 can be interpreted as the position of an element located on the third diagonal of an triangular array (read by rows) provided n > 1. - Enrique Pérez Herrero, Aug 29 2016
a(n) is the sum of the numerator and denominator of the fraction that is the sum of 2/(n-1) + 2/n; all fractions are reduced and n > 2. - J. M. Bergot, Jun 14 2017
a(n) is also the number of maximal irredundant sets in the (n+2)-path complement graph for n > 1. - Eric W. Weisstein, Apr 12 2018
From Klaus Purath, Dec 07 2020: (Start)
a(n) is not divisible by primes listed in A038890. The prime factors are given in A038889 and the prime terms of the sequence are listed in A124199.
Each odd prime factor p divides exactly 2 out of any p consecutive terms with the exception of 17, which appears only once in such an interval of terms. If a(i) and a(k) form such a pair that are divisible by p, then i + k == -3 (mod p), see examples.
If A is a sequence satisfying the recurrence t(n) = 5*t(n-1) - 2*t(n-2) with the initial values either A(0) = 1, A(1) = n + 4 or A(0) = -1, A(1) = n-1, then a(n) = (A(i)^2 - A(i-1)*A(i+1))/2^i for i>0. (End)
Mark each point on a 4^n grid with the number of points that are visible from the point; for n > 1, a(n) is the number of distinct values in the grid. - Torlach Rush, Mar 23 2021
The sequence gives the number of "ON" cells in the cellular automaton on a quadrant of a square grid after the n-th stage, where the "ON" cells lie only on the external perimeter and the perimeter of inscribed squares having the cell (1,1) as a unique common vertex. See Spezia link. - Stefano Spezia, May 28 2025

Examples

			From _Bruno Berselli_, Mar 09 2015: (Start)
By the definition (first formula):
----------------------------------------------------------------------
  1       4         8           13            19              26
----------------------------------------------------------------------
                                                              X
                                              X              X X
                                X            X X            X X X
                    X          X X          X X X          X X X X
          X        X X        X X X        X X X X        X X X X X
  X      X X      X X X      X X X X      X X X X X      X X X X X X
          X        X X        X X X        X X X X        X X X X X
----------------------------------------------------------------------
(End)
From _Klaus Purath_, Dec 07 2020: (Start)
Assuming a(i) is divisible by p with 0 < i < p and a(k) is the next term divisible by p, then from i + k == -3 (mod p) follows that k = min(p*m - i - 3) != i for any integer m.
(1) 17|a(7) => k = min(17*m - 10) != 7 => m = 2, k = 24 == 7 (mod 17). Thus every a(17*m + 7) is divisible by 17.
(2) a(9) = 53 => k = min(53*m - 12) != 9 => m = 1, k = 41. Thus every a(53*m + 9) and a(53*m + 41) is divisible by 53.
(3) 101|a(273) => 229 == 71 (mod 101) => k = min(101*m - 74) != 71 => m = 1, k = 27. Thus every a(101*m + 27) and a(101*m + 71) is divisible by 101. (End)
From _Omar E. Pol_, Aug 08 2021: (Start)
Illustration of initial terms:                             _ _
.                                           _ _           |_|_|_
.                              _ _         |_|_|_         |_|_|_|_
.                   _ _       |_|_|_       |_|_|_|_       |_|_|_|_|_
.          _ _     |_|_|_     |_|_|_|_     |_|_|_|_|_     |_|_|_|_|_|_
.   _     |_|_|    |_|_|_|    |_|_|_|_|    |_|_|_|_|_|    |_|_|_|_|_|_|
.  |_|    |_|_|    |_|_|_|    |_|_|_|_|    |_|_|_|_|_|    |_|_|_|_|_|_|
.
.   1       4         8          13            19              26
------------------------------------------------------------------------ (End)
		

References

  • A. S. Karpenko, Łukasiewicz's Logics and Prime Numbers, 2006 (English translation).
  • G. C. Moisil, Essais sur les logiques non-chrysippiennes, Ed. Academiei, Bucharest, 1972.
  • Wójcicki and Malinowski, eds., Łukasiewicz Sentential Calculi, Wrocław: Ossolineum, 1977.

Crossrefs

Subsequence of A165157.
Triangular numbers (A000217) minus two.
Third diagonal of triangle in A059317.

Programs

  • Haskell
    a034856 = subtract 1 . a000096 -- Reinhard Zumkeller, Feb 20 2015
    
  • Magma
    [Binomial(n + 1, 2) + n - 1: n in [1..60]]; // Vincenzo Librandi, May 21 2011
    
  • Maple
    a := n -> hypergeom([-2, n-1], [1], -1);
    seq(simplify(a(n)), n=1..53); # Peter Luschny, Aug 02 2014
  • Mathematica
    f[n_] := n (n + 3)/2 - 1; Array[f, 55] (* or *) k = 2; NestList[(k++; # + k) &, 1, 55] (* Robert G. Wilson v, Jun 11 2010 *)
    Table[Binomial[n + 1, 2] + n - 1, {n, 53}] (* or *)
    Rest@ CoefficientList[Series[x (1 + x - x^2)/(1 - x)^3, {x, 0, 53}], x] (* Michael De Vlieger, Aug 29 2016 *)
  • Maxima
    A034856(n) := block(
            n-1+(n+1)*n/2
    )$ /* R. J. Mathar, Mar 19 2012 */
    
  • PARI
    A034856(n)=(n+3)*n\2-1 \\ M. F. Hasler, Jan 21 2015
    
  • Python
    def A034856(n): return n*(n+3)//2 -1 # G. C. Greubel, Jun 15 2025

Formula

G.f.: A(x) = x*(1 + x - x^2)/(1 - x)^3.
a(n) = A049600(3, n-2).
a(n) = binomial(n+2, 2) - 2. - Paul Barry, Feb 27 2003
With offset 5, this is binomial(n, 0) - 2*binomial(n, 1) + binomial(n, 2), the binomial transform of (1, -2, 1, 0, 0, 0, ...). - Paul Barry, Jul 01 2003
Row sums of triangle A131818. - Gary W. Adamson, Jul 27 2007
Binomial transform of (1, 3, 1, 0, 0, 0, ...). Also equals A130296 * [1,2,3,...]. - Gary W. Adamson, Jul 27 2007
Row sums of triangle A134225. - Gary W. Adamson, Oct 14 2007
a(n) = A000217(n+1) - 2. - Omar E. Pol, Apr 23 2008
From Jaroslav Krizek, Sep 05 2009: (Start)
a(n) = a(n-1) + n + 1 for n >= 1.
a(n) = n*(n-1)/2 + 2*n - 1.
a(n) = A000217(n-1) + A005408(n-1) = A005843(n-1) + A000124(n-1). (End)
a(n) = Hyper2F1([-2, n-1], [1], -1). - Peter Luschny, Aug 02 2014
a(n) = floor[1/(-1 + Sum_{m >= n+1} 1/S2(m,n+1))], where S2 is A008277. - Richard R. Forberg, Jan 17 2015
a(n) = A101881(2*(n-1)). - Reinhard Zumkeller, Feb 20 2015
a(n) = A253909(n+3) - A000217(n+3). - David Neil McGrath, May 23 2015
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n>3. - David Neil McGrath, May 23 2015
For n > 1, a(n) = 4*binomial(n-1,1) + binomial(n-2,2), comprising the third column of A267633. - Tom Copeland, Jan 25 2016
From Klaus Purath, Dec 07 2020: (Start)
a(n) = A024206(n) + A024206(n+1).
a(2*n-1) = -A168244(n+1).
a(2*n) = A091823(n). (End)
Sum_{n>=1} 1/a(n) = 3/2 + 2*Pi*tan(sqrt(17)*Pi/2)/sqrt(17). - Amiram Eldar, Jan 06 2021
a(n) + a(n+1) = A028347(n+2). - R. J. Mathar, Mar 13 2021
a(n) = A000290(n) - A161680(n-1). - Omar E. Pol, Mar 26 2021
E.g.f.: 1 + exp(x)*(x^2 + 4*x - 2)/2. - Stefano Spezia, Jun 05 2021
a(n) = A024916(n) - A244049(n). - Omar E. Pol, Aug 01 2021
a(n) = A000290(n) - A000217(n-2). - Omar E. Pol, Aug 05 2021

Extensions

More terms from Zerinvary Lajos, May 12 2006

A023531 a(n) = 1 if n is of the form m(m+3)/2, otherwise 0.

Original entry on oeis.org

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, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 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, 1
Offset: 0

Views

Author

Clark Kimberling, Jun 14 1998

Keywords

Comments

Can be read as table: a(n,m) = 1 if n = m >= 0, else 0 (unit matrix).
a(n) = number of 1's between successive 0's (see also A005614, A003589 and A007538). - Eric Angelini, Jul 06 2005
Triangle T(n,k), 0 <= k <= n, read by rows, given by A000004 DELTA A000007 where DELTA is the operator defined in A084938. - Philippe Deléham, Jan 03 2009
Sequence B is called a reverse reluctant sequence of sequence A, if B is triangle array read by rows: row number k lists first k elements of the sequence A in reverse order.
A023531 is reverse reluctant sequence of sequence A000007. - Boris Putievskiy, Jan 11 2013
Also the Bell transform (and the inverse Bell transform) of 0^n (A000007). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 19 2016
This is the turn sequence of the triangle spiral. To form the spiral: go a unit step forward, turn left a(n)*120 degrees, and repeat. The triangle sides are the runs of a(n)=0 (no turn). The sequence can be generated by a morphism with a special symbol S for the start of the sequence: S -> S,1; 1 -> 0,1; 0->0. The expansion lengthens each existing side and inserts a new unit side at the start. See the Fractint L-system in the links to draw the spiral this way. - Kevin Ryde, Dec 06 2019

Examples

			As a triangle:
       1
      0 1
     0 0 1
    0 0 0 1
   0 0 0 0 1
  0 0 0 0 0 1
G.f. = 1 + x^2 + x^5 + x^9 + x^14 + x^20 + x^27 + x^35 + x^44 + x^54 + ...
From _Kevin Ryde_, Dec 06 2019: (Start)
.
              1            Triangular spiral: start at S;
             / \             go a unit step forward,
            0   0   .        turn left a(n)*120 degrees,
           /     \   .       repeat.
          0   1   0   .
         /   / \   \   \   Each side's length is 1 greater
        0   0   0   0   0    than that of the previous side.
       /   /     \   \   \
      0   0   S---1   0   0
     /   /             \   \
    0   1---0---0---0---1   0
   /                         \
  1---0---0---0---0---0---0---1
(End)
		

Crossrefs

Programs

  • Haskell
    a023531 n = a023531_list !! n
    a023531_list = concat $ iterate ([0,1] *) [1]
    instance Num a => Num [a] where
       fromInteger k = [fromInteger k]
       (p:ps) + (q:qs) = p + q : ps + qs
       ps + qs         = ps ++ qs
       (p:ps) * qs'@(q:qs) = p * q : ps * qs' + [p] * qs
        *                = []
    -- Reinhard Zumkeller, Apr 02 2011
    
  • Maple
    seq(op([0$m,1]),m=0..10); # Robert Israel, Jan 18 2015
    # alternative
    A023531 := proc(n)
        option remember ;
        local m,t ;
        for m from 0 do
            t := m*(m+3)/2 ;
            if t > n then
                return 0 ;
            elif t = n then
                return 1 ;
            end if;
        end do:
    end proc:
    seq(A023531(n),n=0..40) ; # R. J. Mathar, May 15 2025
  • Mathematica
    If[IntegerQ[(Sqrt[9+8#]-3)/2],1,0]&/@Range[0,100] (* Harvey P. Dale, Jul 27 2011 *)
    a[ n_] := If[ n < 0, 0, Boole @ IntegerQ @ Sqrt[ 8 n + 9]]; (* Michael Somos, May 17 2014 *)
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 2, 0, x^(1/2)] / (2 x^(1/8)) - 1) / x, {x, 0, n}]; (* Michael Somos, May 17 2014 *)
  • PARI
    {a(n) = if( n<0, 0, issquare(8*n + 9))}; /* Michael Somos, May 17 2014 */
    
  • PARI
    A023531(n)=issquare(8*n+9) \\ M. F. Hasler, Apr 12 2018
    
  • Python
    from math import isqrt
    def A023531(n): return int((k:=n+1<<1)==(m:=isqrt(k))*(m+1)) # Chai Wah Wu, Nov 09 2024
  • Sage
    def A023531_row(n) :
        if n == 0: return [1]
        return [0] + A023531_row(n-1)
    for n in (0..9): print(A023531_row(n))  # Peter Luschny, Jul 22 2012
    

Formula

If (floor(sqrt(2*n))-(2*n/(floor(sqrt(2*n)))) = -1, 1, 0). - Gerald Hillier, Sep 11 2005
a(n) = 1 - A023532(n); a(n) = 1 - mod(floor(((10^(n+2) - 10)/9)10^(n+1 - binomial(floor((1+sqrt(9+8n))/2), 2) - (1+floor(log((10^(n+2) - 10)/9, 10))))), 10). - Paul Barry, May 25 2004
a(n) = floor((sqrt(9+8n)-1)/2) - floor((sqrt(1+8n)-1)/2). - Paul Barry, May 25 2004
a(n) = round(sqrt(2n+3)) - round(sqrt(2n+2)). - Hieronymus Fischer, Aug 06 2007
a(n) = ceiling(2*sqrt(2n+3)) - floor(2*sqrt(2n+2)) - 1. - Hieronymus Fischer, Aug 06 2007
From Franklin T. Adams-Watters, Jun 29 2009: (Start)
G.f.: (1/2 x^{-1/8}theta_2(0,x^{1/2}) - 1)/x, where theta_2 is a Jacobi theta function.
G.f. for triangle: Sum T(n,k) x^n y^k = 1/(1-x*y). Sum T(n,k) x^n y^k / n! = Sum T(n,k) x^n y^k / k! = exp(x*y). Sum T(n,k) x^n y^k / (n! k!) = I_0(2*sqrt(x*y)), where I is the modified Bessel function of the first kind. (End)
a(n) = A000007(m), where m=(t*t+3*t+4)/2-n, t=floor((-1+sqrt(8*n-7))/2). - Boris Putievskiy, Jan 11 2013
The row polynomials are p(n,x) = x^n = (-1)^n n!Lag(n,-n,x), the normalized, associated Laguerre polynomials of order -n. As the prototypical Appell sequence with e.g.f. exp(x*y), its raising operator is R = x and lowering operator, L = d/dx, i.e., R p(n,x) = p(n+1,x), and L p(n,x) = n * p(n-1,x). - Tom Copeland, May 10 2014
a(n) = A010054(n+1) if n >= 0. - Michael Somos, May 17 2014
a(n) = floor(sqrt(2*(n+1)+1/2)-1/2) - floor(sqrt(2*n+1/2)-1/2). - Mikael Aaltonen, Jan 18 2015
a(n) = A003057(n+3) - A003057(n+2). - Robert Israel, Jan 18 2015
a(A000096(n)) = 1; a(A007701(n)) = 0. - Reinhard Zumkeller, Feb 14 2015
Characteristic function of A000096. - M. F. Hasler, Apr 12 2018
Sum_{k=1..n} a(k) ~ sqrt(2*n). - Amiram Eldar, Jan 13 2024

A055998 a(n) = n*(n+5)/2.

Original entry on oeis.org

0, 3, 7, 12, 18, 25, 33, 42, 52, 63, 75, 88, 102, 117, 133, 150, 168, 187, 207, 228, 250, 273, 297, 322, 348, 375, 403, 432, 462, 493, 525, 558, 592, 627, 663, 700, 738, 777, 817, 858, 900, 943, 987, 1032, 1078, 1125, 1173, 1222, 1272
Offset: 0

Views

Author

Barry E. Williams, Jun 14 2000

Keywords

Comments

If X is an n-set and Y a fixed (n-3)-subset of X then a(n-3) is equal to the number of 2-subsets of X intersecting Y. - Milan Janjic, Aug 15 2007
Bisection of A165157. - Jaroslav Krizek, Sep 05 2009
a(n) is the number of (w,x,y) having all terms in {0,...,n} and w=x+y-1. - Clark Kimberling, Jun 02 2012
Numbers m >= 0 such that 8m+25 is a square. - Bruce J. Nicholson, Jul 26 2017
a(n-1) = 3*(n-1) + (n-1)*(n-2)/2 is the number of connected, loopless, non-oriented, multi-edge vertex-labeled graphs with n edges and 3 vertices. Labeled multigraph analog of A253186. There are 3*(n-1) graphs with the 3 vertices on a chain (3 ways to label the middle graph, n-1 ways to pack edges on one of connections) and binomial(n-1,2) triangular graphs (one way to label the graphs, pack 1 or 2 or ...n-2 on the 1-2 edge, ...). - R. J. Mathar, Aug 10 2017
a(n) is also the number of vertices of the quiver for PGL_{n+1} (see Shen). - Stefano Spezia, Mar 24 2020
Starting from a(2) = 7, this is the 4th column of the array: natural numbers written by antidiagonals downwards. See the illustration by Kival Ngaokrajang and the cross-references. - Andrey Zabolotskiy, Dec 21 2021

References

  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, p. 193.

Crossrefs

a(n) = A095660(n+1, 2): third column of (1, 3)-Pascal triangle.
Row n=2 of A255961.

Programs

Formula

G.f.: x*(3-2*x)/(1-x)^3.
a(n) = A027379(n), n > 0.
a(n) = A126890(n,2) for n > 1. - Reinhard Zumkeller, Dec 30 2006
a(n) = A000217(n) + A005843(n). - Reinhard Zumkeller, Sep 24 2008
If we define f(n,i,m) = Sum_{k=0..n-i} binomial(n,k)*Stirling1(n-k,i)*Product_{j=0..k-1} (-m-j), then a(n) = -f(n,n-1,3), for n >= 1. - Milan Janjic, Dec 20 2008
a(n) = A167544(n+8). - Philippe Deléham, Nov 25 2009
a(n) = a(n-1) + n + 2 with a(0)=0. - Vincenzo Librandi, Aug 07 2010
a(n) = Sum_{k=1..n} (k+2). - Gary Detlefs, Aug 10 2010
a(n) = A034856(n+1) - 1 = A000217(n+2) - 3. - Jaroslav Krizek, Sep 05 2009
Sum_{n>=1} 1/a(n) = 137/150. - R. J. Mathar, Jul 14 2012
a(n) = 3*n + A000217(n-1) = 3*n - floor(n/2) + floor(n^2/2). - Wesley Ivan Hurt, Jun 15 2013
a(n) = Sum_{i=3..n+2} i. - Wesley Ivan Hurt, Jun 28 2013
a(n) = 3*A000217(n) - 2*A000217(n-1). - Bruno Berselli, Dec 17 2014
a(n) = A046691(n) + 1. Also, a(n) = A052905(n-1) + 2 = A055999(n-1) + 3 for n>0. - Andrey Zabolotskiy, May 18 2016
E.g.f.: x*(6+x)*exp(x)/2. - G. C. Greubel, Apr 05 2019
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*log(2)/5 - 47/150. - Amiram Eldar, Jan 10 2021
From Amiram Eldar, Feb 12 2024: (Start)
Product_{n>=1} (1 - 1/a(n)) = -5*cos(sqrt(33)*Pi/2)/(4*Pi).
Product_{n>=1} (1 + 1/a(n)) = 15*cos(sqrt(17)*Pi/2)/(2*Pi). (End)

A028552 a(n) = n*(n+3).

Original entry on oeis.org

0, 4, 10, 18, 28, 40, 54, 70, 88, 108, 130, 154, 180, 208, 238, 270, 304, 340, 378, 418, 460, 504, 550, 598, 648, 700, 754, 810, 868, 928, 990, 1054, 1120, 1188, 1258, 1330, 1404, 1480, 1558, 1638, 1720, 1804, 1890, 1978, 2068, 2160, 2254, 2350, 2448, 2548, 2650
Offset: 0

Views

Author

Keywords

Comments

n*(n-3), for n >= 3, is the number of [body] diagonals of an n-gonal prism. - Antreas P. Hatzipolakis (xpolakis(AT)otenet.gr)
a(n) = A028387(n)-1. Half of the difference between n(n+1)(n+2)(n+3) and the largest square less than it. Calling this difference "SquareMod": a(n) = (1/2)*SquareMod(n(n+1)(n+2)(n+3)). - Rainer Rosenthal, Sep 04 2004
n != -2 such that x^4 + x^3 - n*x^2 + x + 1 is reducible over the integers. Starting at 10: n such that x^4 + x^3 - n*x^2 + x + 1 is a product of irreducible quadratic factors over the integers. - James R. Buddenhagen, Apr 19 2005
If a 3-set Y and a 3-set Z, having two element in common, are subsets of an n-set X then a(n-4) is the number of 3-subsets of X intersecting both Y and Z. - Milan Janjic, Oct 03 2007
Starting with offset 1 = binomial transform of [4, 6, 2, 0, 0, 0, ...]. - Gary W. Adamson, Jan 09 2009
The sequence provides all nonnegative integers m such that 4*m + 9 is a square. - Vincenzo Librandi, Mar 03 2013
The second-order linear recurrence relations b(n)=3*b(n-1) + a(m-3)*b(n-2), n>=2, b(0)=0, b(1)=1, have closed form solutions involving only powers of m and 3-m where m>=4 is a positive integer; and lim_{n->infinity} b(n+1)/b(n) = 4. - Felix P. Muga II, Mar 18 2014
If a rook is placed at a corner of an n X n chessboard, the expected number of moves for it to reach the opposite corner is a(n-1). (See Mathematics Stack Exchange link.) - Eric M. Schmidt, Oct 29 2014
Partial sums of the even composites (which are A005843 without the 2). - R. J. Mathar, Sep 09 2015
a(n) is the number of segments necessary to represent n squares of area 1, 4, ..., n^2 having the upper and left sides overlapped:
__ __ __ __
|| || | || | | || | | |
|___| | __| | | __| | |
| __ | | __ | |
| __ __|
4 10 18 28 - Stefano Spezia, May 29 2023

Examples

			G.f. = 4*x + 10*x^2 + 18*x^3 + 28*x^4 + 40*x^5 + 54*x^6 + 70*x^7 + 88*x^8 + ...
		

Crossrefs

Programs

Formula

a(n) = 2*A000096(n).
a(A002522(n)) = A156798(n). - Reinhard Zumkeller, Feb 16 2009
a(n) = a(n-1) + 2*(n+1) for n>0, with a(0)=0. - Vincenzo Librandi, Aug 05 2010
Sum_{n>=1} 1/a(n) = 11/18 via Sum_{n>=0} 1/((n+x)*(n+y)) = (psi(x)-psi(y))/ (x-y). - R. J. Mathar, Mar 22 2011
G.f.: 2*x*(2 - x)/(1 - x)^3. - Arkadiusz Wesolowski, Dec 31 2011
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), with a(0)=0, a(1)=4, a(2)=10. - Harvey P. Dale, Feb 05 2012
a(n) = 4*C(n+1,2) - 2*C(n,2) for n>=0. - Felix P. Muga II, Mar 11 2014
a(-3 - n) = a(n) for all n in Z. - Michael Somos, Mar 18 2014
E.g.f.: (x^3 + 4*x)*exp(x). - G. C. Greubel, Jul 20 2017
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/3 - 5/18. - Amiram Eldar, Jan 15 2021
From Amiram Eldar, Feb 12 2024: (Start)
Product_{n>=1} (1 - 1/a(n)) = 2*cos(sqrt(13)*Pi/2)/Pi.
Product_{n>=1} (1 + 1/a(n)) = -6*cos(sqrt(5)*Pi/2)/Pi. (End)

A008865 a(n) = n^2 - 2.

Original entry on oeis.org

-1, 2, 7, 14, 23, 34, 47, 62, 79, 98, 119, 142, 167, 194, 223, 254, 287, 322, 359, 398, 439, 482, 527, 574, 623, 674, 727, 782, 839, 898, 959, 1022, 1087, 1154, 1223, 1294, 1367, 1442, 1519, 1598, 1679, 1762, 1847, 1934, 2023, 2114, 2207, 2302, 2399, 2498
Offset: 1

Views

Author

Keywords

Comments

For n >= 2, least m >= 1 such that f(m, n) = 0 where f(m,n) = Sum_{i=0..m} Sum_{k= 0..i} (-1)^k*(floor(i/n^k) - n*floor(i/n^(k+1))). - Benoit Cloitre, May 02 2004
For n >= 3, the a(n)-th row of Pascal's triangle always contains a triple forming an arithmetic progression. - Lekraj Beedassy, Jun 03 2004
Let C = 1 + sqrt(2) = 2.414213...; and 1/C = 0.414213... Then a(n) = (n + 1 + 1/C) * (n + 1 - C). Example: a(6) = 34 = (7 + 0.414...) * (7 - 2.414...). - Gary W. Adamson, Jul 29 2009
The sequence (n-4)^2-2, n = 7, 8, ... enumerates the number of non-isomorphic sequences of length n, with entries from {1, 2, 3} and no two adjacent entries the same, that minimally contain each of the thirteen rankings of three players (111, 121, 112, 211, 122, 212, 221, 123, 132, 213, 231, 312, 321) as embedded order isomorphic subsequences. By "minimally", we mean that the n-th symbol is necessary for complete inclusion of all thirteen words. See the arXiv paper below for proof. If n = 7, these sequences are 1213121, 1213212, 1231213, 1231231, 1231321, 1232123, and 1232132, and for each case, there are 3! = 6 isomorphs. - Anant Godbole, Feb 20 2013
a(n), n >= 0, with a(0) = -2, gives the values for a*c of indefinite binary quadratic forms [a, b, c] of discriminant D = 8 for b = 2*n. In general D = b^2 - 4*a*c > 0 and the form [a, b, c] is a*x^2 + b*x*y + c*y^2. - Wolfdieter Lang, Aug 15 2013
With a different offset, this is 2*n^2 - (n + 1)^2, which arises in one explanation of why Bertrand's postulate does not automatically prove Legendre's conjecture: as n gets larger, so does the range of numbers that can have primes that satisfy Bertrand's postulate yet do nothing for Legendre's conjecture. - Alonso del Arte, Nov 06 2013
x*(x + r*y)^2 + y*(y + r*x)^2 can be written as (x + y)*(x^2 + s*x*y + y^2). For r >= 0, the sequence gives the values of s: in fact, s = (r + 1)^2 - 2. - Bruno Berselli, Feb 20 2019
For n >= 2, the continued fraction expansion of sqrt(a(n)) is [n-1; {1, n-2, 1, 2n-2}]. For n=2, this collapses to [1; {2}]. - Magus K. Chu, Sep 06 2022

Examples

			G.f. = -x + 2*x^2 + 7*x^3 + 14*x^4 + 23*x^5 + 34*x^6 + 47*x^7 + 62*x^8 + 79*x^9 + ...
		

Crossrefs

Cf. A145067 (Zero followed by partial sums of A008865).
Cf. A028871 (primes).
Cf. A263766 (partial products).
Cf. A270109. [Bruno Berselli, Mar 17 2016]

Programs

  • Haskell
    a008865 = (subtract 2) . (^ 2) :: Integral t => t -> t
    a008865_list = scanl (+) (-1) [3, 5 ..]
    -- Reinhard Zumkeller, May 06 2013
    
  • Magma
    [n^2 - 2: n in [1..60]]; // Vincenzo Librandi, May 01 2014
  • Mathematica
    Range[50]^2 - 2 (* Harvey P. Dale, Mar 14 2011 *)
  • PARI
    {for(n=1, 47, print1(n^2-2, ","))} \\ Klaus Brockhaus, Oct 17 2008
    

Formula

For n > 1: a(n) = A143053(A000290(n)), A143054(a(n)) = A000290(n). - Reinhard Zumkeller, Jul 20 2008
G.f.: (x-5*x^2+2*x^3)/(-1+3*x-3*x^2+x^3). - Klaus Brockhaus, Oct 17 2008
E.g.f.: (x^2 + x -2)*exp(x) + 2. - G. C. Greubel, Aug 19 2017
a(n+1) = A101986(n) - A101986(n-1) = A160805(n) - A160805(n-1). - Reinhard Zumkeller, May 26 2009
For n > 1, a(n) = floor(n^5/(n^3 + n + 1)). - Gary Detlefs, Feb 10 2010
a(n) = a(n-1) + 2*n - 1 for n > 1, a(1) = -1. - Vincenzo Librandi, Nov 18 2010
Right edge of the triangle in A195437: a(n) = A195437(n-2, n-2). - Reinhard Zumkeller, Nov 23 2011
a(n)*a(n-1) + 2 = (a(n) - n)^2 = A028552(n-2)^2. - Bruno Berselli, Dec 07 2011
a(n+1) = A000096(n) + A000096(n-1) for all n in Z. - Michael Somos, Nov 11 2015
From Amiram Eldar, Jul 13 2020: (Start)
Sum_{n>=1} 1/a(n) = (1 - sqrt(2)*Pi*cot(sqrt(2)*Pi))/4.
Sum_{n>=1} (-1)^n/a(n) = (1 - sqrt(2)*Pi*cosec(sqrt(2)*Pi))/4. (End)
Assume offset 0. Then a(n) = 2*LaguerreL(2, 1 - n). - Peter Luschny, May 09 2021
From Amiram Eldar, Feb 05 2024: (Start)
Product_{n>=1} (1 - 1/a(n)) = sqrt(2/3)*sin(sqrt(3)*Pi)/sin(sqrt(2)*Pi).
Product_{n>=2} (1 + 1/a(n)) = -Pi/(sqrt(2)*sin(sqrt(2)*Pi)). (End)

A106566 Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, 1, 1, 1, 1, 1, 1, 1, ... ] DELTA [1, 0, 0, 0, 0, 0, 0, 0, ... ] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 2, 1, 0, 5, 5, 3, 1, 0, 14, 14, 9, 4, 1, 0, 42, 42, 28, 14, 5, 1, 0, 132, 132, 90, 48, 20, 6, 1, 0, 429, 429, 297, 165, 75, 27, 7, 1, 0, 1430, 1430, 1001, 572, 275, 110, 35, 8, 1, 0, 4862, 4862, 3432, 2002, 1001, 429, 154, 44, 9, 1
Offset: 0

Views

Author

Philippe Deléham, May 30 2005

Keywords

Comments

Catalan convolution triangle; g.f. for column k: (x*c(x))^k with c(x) g.f. for A000108 (Catalan numbers).
Riordan array (1, xc(x)), where c(x) the g.f. of A000108; inverse of Riordan array (1, x*(1-x)) (see A109466).
Diagonal sums give A132364. - Philippe Deléham, Nov 11 2007

Examples

			Triangle begins:
  1;
  0,   1;
  0,   1,   1;
  0,   2,   2,  1;
  0,   5,   5,  3,  1;
  0,  14,  14,  9,  4,  1;
  0,  42,  42, 28, 14,  5, 1;
  0, 132, 132, 90, 48, 20, 6, 1;
From _Paul Barry_, Sep 28 2009: (Start)
Production array is
  0, 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, 1, 1, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1, 1, 1, 1 (End)
		

Crossrefs

The three triangles A059365, A106566 and A099039 are the same except for signs and the leading term.
See also A009766, A033184, A059365 for other versions.
The following are all versions of (essentially) the same Catalan triangle: A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.

Programs

  • Magma
    A106566:= func< n,k | n eq 0 select 1 else (k/n)*Binomial(2*n-k-1, n-k) >;
    [A106566(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Sep 06 2021
    
  • Maple
    A106566 := proc(n,k)
        if n = 0 then
            1;
        elif k < 0 or k > n then
            0;
        else
            binomial(2*n-k-1,n-k)*k/n ;
        end if;
    end proc: # R. J. Mathar, Mar 01 2015
  • Mathematica
    T[n_, k_] := Binomial[2n-k-1, n-k]*k/n; T[0, 0] = 1; Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Feb 18 2017 *)
    (* The function RiordanArray is defined in A256893. *)
    RiordanArray[1&, #(1-Sqrt[1-4#])/(2#)&, 11] // Flatten (* Jean-François Alcover, Jul 16 2019 *)
  • PARI
    {T(n, k) = if( k<=0 || k>n, n==0 && k==0, binomial(2*n - k, n) * k/(2*n - k))}; /* Michael Somos, Oct 01 2022 */
  • Sage
    def A106566(n, k): return 1 if (n==0) else (k/n)*binomial(2*n-k-1, n-k)
    flatten([[A106566(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Sep 06 2021
    

Formula

T(n, k) = binomial(2n-k-1, n-k)*k/n for 0 <= k <= n with n > 0; T(0, 0) = 1; T(0, k) = 0 if k > 0.
T(0, 0) = 1; T(n, 0) = 0 if n > 0; T(0, k) = 0 if k > 0; for k > 0 and n > 0: T(n, k) = Sum_{j>=0} T(n-1, k-1+j).
Sum_{j>=0} T(n+j, 2j) = binomial(2n-1, n), n > 0.
Sum_{j>=0} T(n+j, 2j+1) = binomial(2n-2, n-1), n > 0.
Sum_{k>=0} (-1)^(n+k)*T(n, k) = A064310(n). T(n, k) = (-1)^(n+k)*A099039(n, k).
Sum_{k=0..n} T(n, k)*x^k = A000007(n), A000108(n), A000984(n), A007854(n), A076035(n), A076036(n), A127628(n), A126694(n), A115970(n) for x = 0,1,2,3,4,5,6,7,8 respectively.
Sum_{k>=0} T(n, k)*x^(n-k) = C(x, n); C(x, n) are the generalized Catalan numbers.
Sum_{j=0..n-k} T(n+k,2*k+j) = A039599(n,k).
Sum_{j>=0} T(n,j)*binomial(j,k) = A039599(n,k).
Sum_{k=0..n} T(n,k)*A000108(k) = A127632(n).
Sum_{k=0..n} T(n,k)*(x+1)^k*x^(n-k) = A000012(n), A000984(n), A089022(n), A035610(n), A130976(n), A130977(n), A130978(n), A130979(n), A130980(n), A131521(n) for x= 0,1,2,3,4,5,6,7,8,9 respectively. - Philippe Deléham, Aug 25 2007
Sum_{k=0..n} T(n,k)*A000108(k-1) = A121988(n), with A000108(-1)=0. - Philippe Deléham, Aug 27 2007
Sum_{k=0..n} T(n,k)*(-x)^k = A000007(n), A126983(n), A126984(n), A126982(n), A126986(n), A126987(n), A127017(n), A127016(n), A126985(n), A127053(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Oct 27 2007
T(n,k)*2^(n-k) = A110510(n,k); T(n,k)*3^(n-k) = A110518(n,k). - Philippe Deléham, Nov 11 2007
Sum_{k=0..n} T(n,k)*A000045(k) = A109262(n), A000045: Fibonacci numbers. - Philippe Deléham, Oct 28 2008
Sum_{k=0..n} T(n,k)*A000129(k) = A143464(n), A000129: Pell numbers. - Philippe Deléham, Oct 28 2008
Sum_{k=0..n} T(n,k)*A100335(k) = A002450(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A100334(k) = A001906(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A099322(k) = A015565(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A106233(k) = A003462(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A151821(k+1) = A100320(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A082505(k+1) = A144706(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A000045(2k+2) = A026671(n). - Philippe Deléham, Feb 11 2009
Sum_{k=0..n} T(n,k)*A122367(k) = A026726(n). - Philippe Deléham, Feb 11 2009
Sum_{k=0..n} T(n,k)*A008619(k) = A000958(n+1). - Philippe Deléham, Nov 15 2009
Sum_{k=0..n} T(n,k)*A027941(k+1) = A026674(n+1). - Philippe Deléham, Feb 01 2014
G.f.: Sum_{n>=0, k>=0} T(n, k)*x^k*z^n = 1/(1 - x*z*c(z)) where c(z) the g.f. of A000108. - Michael Somos, Oct 01 2022

Extensions

Formula corrected by Philippe Deléham, Oct 31 2008
Corrected by Philippe Deléham, Sep 17 2009
Corrected by Alois P. Heinz, Aug 02 2012

A014106 a(n) = n*(2*n + 3).

Original entry on oeis.org

0, 5, 14, 27, 44, 65, 90, 119, 152, 189, 230, 275, 324, 377, 434, 495, 560, 629, 702, 779, 860, 945, 1034, 1127, 1224, 1325, 1430, 1539, 1652, 1769, 1890, 2015, 2144, 2277, 2414, 2555, 2700, 2849, 3002, 3159, 3320, 3485, 3654, 3827, 4004, 4185, 4370
Offset: 0

Views

Author

Keywords

Comments

If Y is a 2-subset of a 2n-set X then, for n >= 1, a(n-1) is the number of (2n-2)-subsets of X intersecting Y. - Milan Janjic, Nov 18 2007
This sequence can also be derived from 1*(2+3)=5, 2*(3+4)=14, 3*(4+5)=27, and so forth. - J. M. Bergot, May 30 2011
Consider the partitions of 2n into exactly two parts. Then a(n) is the sum of all the parts in the partitions of 2n + the number of partitions of 2n + the total number of partition parts of 2n. - Wesley Ivan Hurt, Jul 02 2013
a(n) is the number of self-intersecting points of star polygon {(2*n+3)/(n+1)}. - Bui Quang Tuan, Mar 25 2015
Bisection of A000096. - Omar E. Pol, Dec 16 2016
a(n+1) is the number of function calls required to compute Ackermann's function ack(2,n). - Olivier Gérard, May 11 2018
a(n-1) is the least denominator d > n of the best rational approximation of sqrt(n^2-2) by x/d (see example and PARI code). - Hugo Pfoertner, Apr 30 2019
The number of cells in a loose n X n+1 rectangular spiral where n is even. See loose rectangular spiral image. - Jeff Bowermaster, Aug 05 2019
a(n-1) is the dimension of the second cohomology group of 2n+1-dimensional Heisenberg Lie algebra h_{2n+1}. - Rafik Khalfi, Jan 27 2025

Examples

			a(5-1) = 44: The best approximation of sqrt(5^2-2) = sqrt(23) by x/d with d <= k is 24/5 for all k < 44, but sqrt(23) ~= 211/44 is the first improvement. - _Hugo Pfoertner_, Apr 30 2019
		

References

  • Jolley, Summation of Series, Dover (1961).

Crossrefs

Cf. A091823. See A110325 for another version.

Programs

  • Magma
    [n*(2*n+3): n in [0..50]]; // Vincenzo Librandi, Apr 25 2011
  • Maple
    A014106 := proc(n) n*(2*n+3) ; end proc: # R. J. Mathar, Feb 13 2011
    seq(k*(2*k+3), k=1..100); # Wesley Ivan Hurt, Jul 02 2013
  • Mathematica
    Table[n (2 n + 3), {n, 0, 120}] (* Michael De Vlieger, Apr 02 2015 *)
    LinearRecurrence[{3,-3,1},{0,5,14},50] (* Harvey P. Dale, Jul 21 2023 *)
  • PARI
    a(n)=2*n^2+3*n
    
  • PARI
    \\ least denominator > n in best rational approximation of sqrt(n^2-2)
    for(n=2,47,for(k=n,oo,my(m=denominator(bestappr(sqrt(n^2-2),k)));if(m>n,print1(k,", ");break(1)))) \\ Hugo Pfoertner, Apr 30 2019
    

Formula

a(n) - 1 = A091823(n). - Howard A. Landman, Mar 28 2004
A014107(-n) = a(n), A000384(n+1) = a(n)+1. - Michael Somos, Nov 06 2005
G.f.: x*(5 - x)/(1 - x)^3. - Paul Barry, Feb 27 2003
E.g.f: x*(5 + 2*x)*exp(x). - Michael Somos, Nov 06 2005
a(n) = a(n-1) + 4*n + 1, n > 0. - Vincenzo Librandi, Nov 19 2010
a(n) = 4*A000217(n) + n. - Bruno Berselli, Feb 11 2011
Sum_{n>=1} 1/a(n) = 8/9 -2*log(2)/3 = 0.4267907685155920.. [Jolley eq. 265]
Sum_{n>=1} (-1)^(n+1)/a(n) = 4/9 + log(2)/3 - Pi/6. - Amiram Eldar, Jul 03 2020
From Leo Tavares, Jan 27 2022: (Start)
a(n) = A000384(n+1) - 1. See Hex-tangles illustration.
a(n) = A014105(n) + n*2. See Second Hex-tangles illustration.
a(n) = 2*A002378(n) + n. See Ob-tangles illustration.
a(n) = A005563(n) + 2*A000217(n). See Trap-tangles illustration. (End)

A023532 a(n) = 0 if n is of the form m*(m+3)/2, otherwise 1.

Original entry on oeis.org

0, 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, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 0

Views

Author

Keywords

Comments

From Stark: "alpha = 0.101101110111101111101111110 ... is irrational. For if alpha were rational, its decimal expansion would be periodic and have a period of length r starting with the k-th digit of the expansion.
"But by the very nature of alpha, there will be blocks of r digits, all 1, in this expansion after the k-th digit and the periodicity would then guarantee that everything after such a block of r digits would also be all ones.
"This contradicts the fact that there will always be zeros occurring after any given point in the expansion of alpha. Hence alpha is irrational."
a(A000096(n)) = 0; a(A007401(n)) = 1. - Reinhard Zumkeller, Dec 04 2012
Sequence B is called a reverse reluctant sequence of sequence A, if B is triangle array read by rows: row number k lists first k elements of the sequence A in reverse order. A023532 is reverse reluctant sequence of sequence A211666. - Boris Putievskiy, Jan 11 2013
An example of a sequence with infinite critical exponent [Vaslet]. - N. J. A. Sloane, May 05 2013

Examples

			From _Boris Putievskiy_, Jan 11 2013: (Start)
As a triangular array written by rows, the sequence begins:
  0;
  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;
  ...
(End)
		

References

  • Harold M. Stark, An Introduction to Number Theory, The MIT Press, Cambridge, Mass, eighth printing 1994, page 170.

Crossrefs

Essentially the same sequence as A114607 and A123110. - N. J. A. Sloane, Feb 07 2020

Programs

  • Haskell
    a023532 = (1 -) . a010052 . (+ 9) . (* 8)
    a023532_list = concat $ iterate (\rs -> 1 : rs) [0]
    -- Reinhard Zumkeller, Dec 04 2012
    
  • Maple
    A023532 := proc(n)
        option remember ;
        local m,t ;
        for m from 0 do
            t := m*(m+3)/2 ;
            if t > n then
                return 1 ;
            elif t = n then
                return 0 ;
            end if;
        end do:
    end proc:
    seq(A023532(n),n=0..40) ; # R. J. Mathar, May 15 2025
  • Mathematica
    a = {}; Do[a = Append[a, Join[ {0}, Table[1, {n} ] ] ], {n, 1, 13} ]; a = Flatten[a]
    Table[PadLeft[{0},n,1],{n,0,20}]//Flatten (* Harvey P. Dale, Jul 10 2019 *)
  • PARI
    for(n=1,9,print1("0, ");for(i=1,n,print1("1, "))) \\ Charles R Greathouse IV, Jun 16 2011
    
  • PARI
    a(n)=!issquare(8*n+9) \\ Charles R Greathouse IV, Jun 16 2011
    
  • Python
    from sympy.ntheory.primetest import is_square
    def A023532(n): return bool(is_square((n<<3)+9))^1 # Chai Wah Wu, Feb 10 2023

Formula

a(n) = 0 if and only if 8n+9 is a square. - Charles R Greathouse IV, Jun 16 2011
Blocks of lengths 1, 2, 3, 4, ... of ones separated by a single zero.
a(n) = 1 - floor((sqrt(9+8n)-1)/2) + floor((sqrt(1+8n)-1)/2). - Paul Barry, May 25 2004
a(n) = A211666(m), where m = (t^2 + 3*t + 4)/2n - n, t = floor((-1 + sqrt(8*n-7))/2). - Boris Putievskiy, Jan 11 2013
a(n) = [A002262(n) < A003056(n)]. - Yuchun Ji, May 18 2020
a(n) = 1-A023531(n). - R. J. Mathar, May 15 2025

Extensions

Additional comments from Robert G. Wilson v, Nov 06 2000

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

A007401 Add n-1 to n-th term of 'n appears n times' sequence (A002024).

Original entry on oeis.org

1, 3, 4, 6, 7, 8, 10, 11, 12, 13, 15, 16, 17, 18, 19, 21, 22, 23, 24, 25, 26, 28, 29, 30, 31, 32, 33, 34, 36, 37, 38, 39, 40, 41, 42, 43, 45, 46, 47, 48, 49, 50, 51, 52, 53, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76
Offset: 1

Views

Author

Keywords

Comments

Complement of A000096 = increasing sequence of positive integers excluding n*(n+3)/2. - Jonathan Vos Post, Aug 13 2005
As a triangle: (1; 3,4; 6,7,8; 10,11,12,13; ...), Row sums = A127736: (1, 7, 21, 46, 85, 141, 217, ...). - Gary W. Adamson, Oct 25 2007
Odd primes are a subsequence except 5, cf. A004139. - Reinhard Zumkeller, Jul 18 2011
A023532(a(n)) = 1. - Reinhard Zumkeller, Dec 04 2012
T(n,k) = ((n+k)^2+n-k)/2 - 1, n,k > 0, read by antidiagonals. - Boris Putievskiy, Jan 14 2013
A023531(a(n)) = 0. - Reinhard Zumkeller, Feb 14 2015

Examples

			From _Boris Putievskiy_, Jan 14 2013: (Start)
The start of the sequence as table:
   1,  3,  6, 10, 15, 21, 28, ...
   4,  7, 11, 16, 22, 29, 37, ...
   8, 12, 17, 23, 30, 38, 47, ...
  13, 18, 24, 31, 39, 48, 58, ...
  19, 25, 32, 40, 49, 59, 70, ...
  26, 33, 41, 50, 60, 71, 83, ...
  34, 42, 51, 61, 72, 84, 97, ...
  ...
The start of the sequence as triangle array read by rows:
   1;
   3,  4;
   6,  7,  8;
  10, 11, 12, 13;
  15, 16, 17, 18, 19;
  21, 22, 23, 24, 25, 26;
  28, 29, 30, 31, 32, 33, 34;
  ...
Row number r contains r numbers r*(r+1)/2, r*(r+1)/2+1, ..., r*(r+1)/2+r-1. (End)
		

References

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

Crossrefs

Programs

  • Haskell
    a007401 n = a007401_list !! n
    a007701_list = [x | x <- [0..], a023531 x == 0]
    -- Reinhard Zumkeller, Feb 14 2015, Dec 04 2012
    
  • Mathematica
    f[n_] := n + Floor[ Sqrt[2n] - 1/2]; Array[f, 66]; (* Robert G. Wilson v, Feb 13 2011 *)
  • PARI
    a(n)=n+floor(sqrt(n+n)-1/2) \\ Charles R Greathouse IV, Feb 13 2011
    
  • PARI
    for(m=1,9, for(n=m*(m+1)/2,(m^2+3*m-2)/2, print1(n", "))) \\ Charles R Greathouse IV, Feb 13 2011
    
  • Python
    from math import isqrt
    def A007401(n): return n-1+(isqrt(n<<3)+1>>1) # Chai Wah Wu, Oct 18 2022

Formula

From Boris Putievskiy, Jan 14 2013: (Start)
a(n) = A014132(n) - 1.
a(n) = A003057(n)^2 - A114327(n) - 1.
a(n) = ((t+2)^2 + i - j)/2-1, where
i = n-t*(t+1)/2,
j = (t*t+3*t+4)/2-n,
t = floor((-1+sqrt(8*n-7))/2). (End)
Previous Showing 21-30 of 263 results. Next