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 31-40 of 274 results. Next

A007310 Numbers congruent to 1 or 5 mod 6.

Original entry on oeis.org

1, 5, 7, 11, 13, 17, 19, 23, 25, 29, 31, 35, 37, 41, 43, 47, 49, 53, 55, 59, 61, 65, 67, 71, 73, 77, 79, 83, 85, 89, 91, 95, 97, 101, 103, 107, 109, 113, 115, 119, 121, 125, 127, 131, 133, 137, 139, 143, 145, 149, 151, 155, 157, 161, 163, 167, 169, 173, 175
Offset: 1

Views

Author

C. Christofferson (Magpie56(AT)aol.com)

Keywords

Comments

Numbers n such that phi(4n) = phi(3n). - Benoit Cloitre, Aug 06 2003
Or, numbers relatively prime to 2 and 3, or coprime to 6, or having only prime factors >= 5; also known as 5-rough numbers. (Edited by M. F. Hasler, Nov 01 2014: merged with comments from Zak Seidov, Apr 26 2007 and Michael B. Porter, Oct 09 2009)
Apart from initial term(s), dimension of the space of weight 2n cuspidal newforms for Gamma_0( 38 ).
Numbers k such that k mod 2 = 1 and (k+1) mod 3 <> 1. - Klaus Brockhaus, Jun 15 2004
Also numbers n such that the sum of the squares of the first n integers is divisible by n, or A000330(n) = n*(n+1)*(2*n+1)/6 is divisible by n. - Alexander Adamchuk, Jan 04 2007
Numbers n such that the sum of squares of n consecutive integers is divisible by n, because A000330(m+n) - A000330(m) = n*(n+1)*(2*n+1)/6 + n*(m^2+n*m+m) is divisible by n independent of m. - Kaupo Palo, Dec 10 2016
A126759(a(n)) = n + 1. - Reinhard Zumkeller, Jun 16 2008
Terms of this sequence (starting from the second term) are equal to the result of the expression sqrt(4!*(k+1) + 1) - but only when this expression yields integral values (that is when the parameter k takes values, which are terms of A144065). - Alexander R. Povolotsky, Sep 09 2008
For n > 1: a(n) is prime if and only if A075743(n-2) = 1; a(2*n-1) = A016969(n-1), a(2*n) = A016921(n-1). - Reinhard Zumkeller, Oct 02 2008
A156543 is a subsequence. - Reinhard Zumkeller, Feb 10 2009
Numbers n such that ChebyshevT(x, x/2) is not an integer (is integer/2). - Artur Jasinski, Feb 13 2010
If 12*k + 1 is a perfect square (k = 0, 2, 4, 10, 14, 24, 30, 44, ... = A152749) then the square root of 12*k + 1 = a(n). - Gary Detlefs, Feb 22 2010
A089128(a(n)) = 1. Complement of A047229(n+1) for n >= 1. See A164576 for corresponding values A175485(a(n)). - Jaroslav Krizek, May 28 2010
Cf. property described by Gary Detlefs in A113801 and in Comment: more generally, these numbers are of the form (2*h*n+(h-4)*(-1)^n-h)/4 (with h, n natural numbers), therefore ((2*h*n+(h-4)*(-1)^n-h)/4)^2-1 == 0 (mod h); in this case, a(n)^2 - 1 == 0 (mod 6). Also a(n)^2 - 1 == 0 (mod 12). - Bruno Berselli, Nov 05 2010 - Nov 17 2010
Numbers n such that ( Sum_{k = 1..n} k^14 ) mod n = 0. (Conjectured) - Gary Detlefs, Dec 27 2011
From Peter Bala, May 02 2018: (Start)
The above conjecture is true. Apply Ireland and Rosen, Proposition 15.2.2. with m = 14 to obtain the congruence 6*( Sum_{k = 1..n} k^14 )/n = 7 (mod n), true for all n >= 1. Suppose n is coprime to 6, then 6 is a unit in Z/nZ, and it follows from the congruence that ( Sum_{k = 1..n} k^14 )/n is an integer. On the other hand, if either 2 divides n or 3 divides n then the congruence shows that ( Sum_{k = 1..n} k^14 )/n cannot be integral. (End)
A126759(a(n)) = n and A126759(m) < n for m < a(n). - Reinhard Zumkeller, May 23 2013
(a(n-1)^2 - 1)/24 = A001318(n), the generalized pentagonal numbers. - Richard R. Forberg, May 30 2013
Numbers k for which A001580(k) is divisible by 3. - Bruno Berselli, Jun 18 2014
Numbers n such that sigma(n) + sigma(2n) = sigma(3n). - Jahangeer Kholdi and Farideh Firoozbakht, Aug 15 2014
a(n) are values of k such that Sum_{m = 1..k-1} m*(k-m)/k is an integer. Sums for those k are given by A062717. Also see Detlefs formula below based on A062717. - Richard R. Forberg, Feb 16 2015
a(n) are exactly those positive integers m such that the sequence b(n) = n*(n + m)*(n + 2*m)/6 is integral, and also such that the sequence c(n) = n*(n + m)*(n + 2*m)*(n + 3*m)/24 is integral. Cf. A007775. - Peter Bala, Nov 13 2015
Along with 2, these are the numbers k such that the k-th Fibonacci number is coprime to every Lucas number. - Clark Kimberling, Jun 21 2016
This sequence is the Engel expansion of 1F2(1; 5/6, 7/6; 1/36) + 1F2(1; 7/6, 11/6; 1/36)/5. - Benedict W. J. Irwin, Dec 16 2016
The sequence a(n), n >= 4 is generated by the successor of the pair of polygonal numbers {P_s(4) + 1, P_(2*s - 1)(3) + 1}, s >= 3. - Ralf Steiner, May 25 2018
The asymptotic density of this sequence is 1/3. - Amiram Eldar, Oct 18 2020
Also, the only vertices in the odd Collatz tree A088975 that are branch values to other odd nodes t == 1 (mod 2) of A005408. - Heinz Ebert, Apr 14 2021
From Flávio V. Fernandes, Aug 01 2021: (Start)
For any two terms j and k, the product j*k is also a term (the same property as p^n and smooth numbers).
From a(2) to a(phi(A033845(n))), or a((A033845(n))/3), the terms are the totatives of the A033845(n) itself. (End)
Also orders n for which cyclic and semicyclic diagonal Latin squares exist (see A123565 and A342990). - Eduard I. Vatutin, Jul 11 2023
If k is in the sequence, then k*2^m + 3 is also in the sequence, for all m > 0. - Jules Beauchamp, Aug 29 2024

Examples

			G.f. = x + 5*x^2 + 7*x^3 + 11*x^4 + 13*x^5 + 17*x^6 + 19*x^7 + 23*x^8 + ...
		

References

  • K. Ireland and M. Rosen, A Classical Introduction to Modern Number Theory, Springer-Verlag, 1980.

Crossrefs

A005408 \ A016945. Union of A016921 and A016969; union of A038509 and A140475. Essentially the same as A038179. Complement of A047229. Subsequence of A186422.
Cf. A000330, A001580, A002194, A019670, A032528 (partial sums), A038509 (subsequence of composites), A047209, A047336, A047522, A056020, A084967, A090771, A091998, A144065, A175885-A175887.
For k-rough numbers with other values of k, see A000027, A005408, A007775, A008364-A008366, A166061, A166063.
Cf. A126760 (a left inverse).
Row 3 of A260717 (without the initial 1).
Cf. A105397 (first differences).

Programs

Formula

a(n) = (6*n + (-1)^n - 3)/2. - Antonio Esposito, Jan 18 2002
a(n) = a(n-1) + a(n-2) - a(n-3), n >= 4. - Roger L. Bagula
a(n) = 3*n - 1 - (n mod 2). - Zak Seidov, Jan 18 2006
a(1) = 1 then alternatively add 4 and 2. a(1) = 1, a(n) = a(n-1) + 3 + (-1)^n. - Zak Seidov, Mar 25 2006
1 + 1/5^2 + 1/7^2 + 1/11^2 + ... = Pi^2/9 [Jolley]. - Gary W. Adamson, Dec 20 2006
For n >= 3 a(n) = a(n-2) + 6. - Zak Seidov, Apr 18 2007
From R. J. Mathar, May 23 2008: (Start)
Expand (x+x^5)/(1-x^6) = x + x^5 + x^7 + x^11 + x^13 + ...
O.g.f.: x*(1+4*x+x^2)/((1+x)*(1-x)^2). (End)
a(n) = 6*floor(n/2) - 1 + 2*(n mod 2). - Reinhard Zumkeller, Oct 02 2008
1 + 1/5 - 1/7 - 1/11 + + - - ... = Pi/3 = A019670 [Jolley eq (315)]. - Jaume Oliver Lafont, Oct 23 2009
a(n) = ( 6*A062717(n)+1 )^(1/2). - Gary Detlefs, Feb 22 2010
a(n) = 6*A000217(n-1) + 1 - 2*Sum_{i=1..n-1} a(i), with n > 1. - Bruno Berselli, Nov 05 2010
a(n) = 6*n - a(n-1) - 6 for n>1, a(1) = 1. - Vincenzo Librandi, Nov 18 2010
Sum_{n >= 1} (-1)^(n+1)/a(n) = A093766 [Jolley eq (84)]. - R. J. Mathar, Mar 24 2011
a(n) = 6*floor(n/2) + (-1)^(n+1). - Gary Detlefs, Dec 29 2011
a(n) = 3*n + ((n+1) mod 2) - 2. - Gary Detlefs, Jan 08 2012
a(n) = 2*n + 1 + 2*floor((n-2)/2) = 2*n - 1 + 2*floor(n/2), leading to the o.g.f. given by R. J. Mathar above. - Wolfdieter Lang, Jan 20 2012
1 - 1/5 + 1/7 - 1/11 + - ... = Pi*sqrt(3)/6 = A093766 (L. Euler). - Philippe Deléham, Mar 09 2013
1 - 1/5^3 + 1/7^3 - 1/11^3 + - ... = Pi^3*sqrt(3)/54 (L. Euler). - Philippe Deléham, Mar 09 2013
gcd(a(n), 6) = 1. - Reinhard Zumkeller, Nov 14 2013
a(n) = sqrt(6*n*(3*n + (-1)^n - 3)-3*(-1)^n + 5)/sqrt(2). - Alexander R. Povolotsky, May 16 2014
a(n) = 3*n + 6/(9*n mod 6 - 6). - Mikk Heidemaa, Feb 05 2016
From Mikk Heidemaa, Feb 11 2016: (Start)
a(n) = 2*floor(3*n/2) - 1.
a(n) = A047238(n+1) - 1. (suggested by Michel Marcus) (End)
E.g.f.: (2 + (6*x - 3)*exp(x) + exp(-x))/2. - Ilya Gutkovskiy, Jun 18 2016
From Bruno Berselli, Apr 27 2017: (Start)
a(k*n) = k*a(n) + (4*k + (-1)^k - 3)/2 for k>0 and odd n, a(k*n) = k*a(n) + k - 1 for even n. Some special cases:
k=2: a(2*n) = 2*a(n) + 3 for odd n, a(2*n) = 2*a(n) + 1 for even n;
k=3: a(3*n) = 3*a(n) + 4 for odd n, a(3*n) = 3*a(n) + 2 for even n;
k=4: a(4*n) = 4*a(n) + 7 for odd n, a(4*n) = 4*a(n) + 3 for even n;
k=5: a(5*n) = 5*a(n) + 8 for odd n, a(5*n) = 5*a(n) + 4 for even n, etc. (End)
From Antti Karttunen, May 20 2017: (Start)
a(A273669(n)) = 5*a(n) = A084967(n).
a((5*n)-3) = A255413(n).
A126760(a(n)) = n. (End)
a(2*m) = 6*m - 1, m >= 1; a(2*m + 1) = 6*m + 1, m >= 0. - Ralf Steiner, May 17 2018
From Amiram Eldar, Nov 22 2024: (Start)
Product_{n>=1} (1 - (-1)^n/a(n)) = sqrt(3) (A002194).
Product_{n>=2} (1 + (-1)^n/a(n)) = Pi/3 (A019670). (End)

A026741 a(n) = n if n odd, n/2 if n even.

Original entry on oeis.org

0, 1, 1, 3, 2, 5, 3, 7, 4, 9, 5, 11, 6, 13, 7, 15, 8, 17, 9, 19, 10, 21, 11, 23, 12, 25, 13, 27, 14, 29, 15, 31, 16, 33, 17, 35, 18, 37, 19, 39, 20, 41, 21, 43, 22, 45, 23, 47, 24, 49, 25, 51, 26, 53, 27, 55, 28, 57, 29, 59, 30, 61, 31, 63, 32, 65, 33, 67, 34, 69, 35, 71, 36, 73, 37, 75, 38
Offset: 0

Views

Author

J. Carl Bellinger (carlb(AT)ctron.com)

Keywords

Comments

a(n) is the size of largest conjugacy class in D_2n, the dihedral group with 2n elements. - Sharon Sela (sharonsela(AT)hotmail.com), May 14 2002
a(n+1) is the composition length of the n-th symmetric power of the natural representation of a finite subgroup of SL(2,C) of type D_4 (quaternion group). - Paul Boddington, Oct 23 2003
For n > 1, a(n) is the greatest common divisor of all permutations of {0, 1, ..., n} treated as base n + 1 integers. - David Scambler, Nov 08 2006 (see the Mathematics Stack Exchange link below).
From Dimitrios Choussos (choussos(AT)yahoo.de), May 11 2009: (Start)
Sequence A075888 and the above sequence are fitting together.
First 2 entries of this sequence have to be taken out.
In some cases two three or more sequenced entries of this sequence have to be added together to get the next entry of A075888.
Example: Sequences begin with 1, 3, 2, 5, 3, 7, 4, 9 (4 + 9 = 13, the next entry in A075888).
But it works out well up to primes around 50000 (haven't tested higher ones).
As A075888 gives a very regular graph. There seems to be a regularity in the primes. (End)
Starting with 1 = triangle A115359 * [1, 2, 3, ...]. - Gary W. Adamson, Nov 27 2009
From Gary W. Adamson, Dec 11 2009: (Start)
Let M be an infinite lower triangular matrix with (1, 1, 1, 0, 0, 0, ...) in every column, shifted down twice. This sequence starting with 1 = M * (1, 2, 3, ...)
M =
1;
1, 0;
1, 1, 0;
0, 1, 0, 0;
0, 1, 1, 0, 0;
0, 0, 1, 0, 0, 0;
0, 0, 1, 1, 0, 0, 0;
...
A026741 = M * (1, 2, 3, ...); but A002487 = lim_{n->infinity} M^n, a left-shifted vector considered as a sequence. (End)
A particular case of sequence for which a(n+3) = (a(n+2) * a(n+1)+q)/a(n) for every n > n0. Here n0 = 1 and q = -1. - Richard Choulet, Mar 01 2010
For n >= 2, a(n+1) is the smallest m such that s_n(2*m*(n-1))/(n-1) is even, where s_b(c) is the sum of digits of c in base b. - Vladimir Shevelev, May 02 2011
A001477 and A005408 interleaved. - Omar E. Pol, Aug 22 2011
Numerator of n/((n-1)*(n-2)). - Michael B. Porter, Mar 18 2012
Number of odd terms of n-th row in the triangles A162610 and A209297. - Reinhard Zumkeller, Jan 19 2013
For n >= 3, a(n) is the periodic of integer of spiral length ratio of spiral that have (n-1) circle centers. See illustration in links. - Kival Ngaokrajang, Dec 28 2013
This is the sequence of Lehmer numbers u_n(sqrt(R), Q) with the parameters R = 4 and Q = 1. It is a strong divisibility sequence, that is, gcd(a(n), a(m)) = a(gcd(n, m)) for all natural numbers n and m. Cf. A005013 and A108412. - Peter Bala, Apr 18 2014
The sequence of convergents of the 2-periodic continued fraction [0; 1, -4, 1, -4, ...] = 1/(1 - 1/(4 - 1/(1 - 1/(4 - ...)))) = 2 begins [0/1, 1/1, 4/3, 3/2, 8/5, 5/3, 12/7, ...]. The present sequence is the sequence of denominators; the sequence of numerators of the continued fraction convergents [0, 1, 4, 3, 8, 5, 12, ...] is A022998, also a strong divisibility sequence. - Peter Bala, May 19 2014
For n >= 3, (a(n-2)/a(n))*Pi = vertex angle of a regular n-gon. See illustration in links. - Kival Ngaokrajang, Jul 17 2014
For n > 1, the numerator of the harmonic mean of the first n triangular numbers. - Colin Barker, Nov 13 2014
The difference sequence is a permutation of the integers. - Clark Kimberling, Apr 19 2015
From Timothy Hopper, Feb 26 2017: (Start)
Given the function a(n, p) = n/p if n mod p = 0, else n, then a possible formula is: a(n, p) = n*(1 + (p-1)*((n^(p-1)) mod p))/p, p prime, (n^(p-1)) mod p = 1, n not divisible by p. (Fermat's Little Theorem). Examples: p = 2; a(n), p = 3; A051176(n), p = 5; A060791(n), p = 7; A106608(n).
Conjecture: lcm(n, p) = p*a(n, p), gcd(n, p) = n/a(n, p). (End)
Let r(n) = (a(n+1) + 1)/a(n+1) if n mod 2 = 1, a(n+1)/(a(n+1) + 2) otherwise; then lim_{k->oo} 2^(k+2) * Product_{n=0..k} r(n)^(k-n) = Pi. - Dimitris Valianatos, Mar 22 2021
Number of integers k from 1 to n such that gcd(n,k) is odd. - Amiram Eldar, May 18 2025

Examples

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

References

  • David Wells, Prime Numbers: The Most Mysterious Figures in Math. Hoboken, New Jersey: John Wiley & Sons (2005), p. 53.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers, 2nd Ed. Penguin (1997), p. 79.

Crossrefs

Signed version is in A030640. Partial sums give A001318.
Cf. A051176, A060819, A060791, A060789 for n / gcd(n, k) with k = 3..6. See also A106608 thru A106612 (k = 7 thru 11), A051724 (k = 12), A106614 thru A106621 (k = 13 thru 20).
Cf. A013942.
Cf. A227042 (first column). Cf. A005013 and A108412.

Programs

  • Haskell
    import Data.List (transpose)
    a026741 n = a026741_list !! n
    a026741_list = concat $ transpose [[0..], [1,3..]]
    -- Reinhard Zumkeller, Dec 12 2011
    
  • Magma
    [2*n/(3+(-1)^n): n in [0..70]]; // Vincenzo Librandi, Aug 14 2011
    
  • Maple
    A026741 := proc(n) if type(n,'odd') then n; else n/2; end if; end proc: seq(A026741(n), n=0..76); # R. J. Mathar, Jan 22 2011
  • Mathematica
    Numerator[Abs[Table[Det[DiagonalMatrix[Table[1/i^2 - 1, {i, 1, n - 1}]] + 1], {n, 20}]]] (* Alexander Adamchuk, Jun 02 2006 *)
    halfMax = 40; Riffle[Range[0, halfMax], Range[1, 2halfMax + 1, 2]] (* Harvey P. Dale, Mar 27 2011 *)
    a[ n_] := Numerator[n / 2]; (* Michael Somos, Jan 20 2017 *)
    Array[If[EvenQ[#],#/2,#]&,80,0] (* Harvey P. Dale, Jul 08 2023 *)
  • PARI
    a(n) = numerator(n/2) \\ Rick L. Shepherd, Sep 12 2007
    
  • Python
    def A026741(n): return n if n % 2 else n//2 # Chai Wah Wu, Apr 02 2021
  • Sage
    [lcm(n, 2) / 2 for n in range(77)] # Zerinvary Lajos, Jun 07 2009
    

Formula

G.f.: x*(1 + x + x^2)/(1-x^2)^2. - Len Smiley, Apr 30 2001
a(n) = 2*a(n-2) - a*(n-4) for n >= 4.
a(n) = n * 2^((n mod 2) - 1). - Reinhard Zumkeller, Oct 16 2001
a(n) = 2*n/(3 + (-1)^n). - Benoit Cloitre, Mar 24 2002
Multiplicative with a(2^e) = 2^(e-1) and a(p^e) = p^e, p > 2. - Vladeta Jovovic, Apr 05 2002
a(n) = n / gcd(n, 2). a(n)/A045896(n) = n/((n+1)*(n+2)).
For n > 0, a(n) = denominator of Sum_{i=1..n-1} 2/(i*(i+1)), numerator=A022998. - Reinhard Zumkeller, Apr 21 2012, Jul 25 2002 [thanks to Phil Carmody who noticed an error]
For n > 1, a(n) = GCD of the n-th and (n-1)-th triangular numbers (A000217). - Ross La Haye, Sep 13 2003
Euler transform of finite sequence [1, 2, -1]. - Michael Somos, Jun 15 2005
G.f.: x * (1 - x^3) / ((1 - x) * (1 - x^2)^2) = Sum_{k>0} k * (x^k - x^(2*k)). - Michael Somos, Jun 15 2005
a(n+3) + a(n+2) = 3 + a(n+1) + a(n). a(n+3) * a(n) = - 1 + a(n+2) * a(n+1). a(n) = -a(-n) for all n in Z. - Michael Somos, Jun 15 2005
For n > 1, a(n) is the numerator of the average of 1, 2, ..., n - 1; i.e., numerator of A000217(n-1)/(n-1), with corresponding denominators [1, 2, 1, 2, ...] (A000034). - Rick L. Shepherd, Jun 05 2006
Equals A126988 * (1, -1, 0, 0, 0, ...). - Gary W. Adamson, Apr 17 2007
For n >= 1, a(n) = gcd(n,A000217(n)). - Rick L. Shepherd, Sep 12 2007
a(n) = numerator(n/(2*n-2)) for n >= 2; A022998(n-1) = denominator(n/(2*n-2)) for n >= 2. - Johannes W. Meijer, Jun 18 2009
a(n) = A167192(n+2, 2). - Reinhard Zumkeller, Oct 30 2009
a(n) = A106619(n) * A109012(n). - Paul Curtz, Apr 04 2011
From R. J. Mathar, Apr 18 2011: (Start)
a(n) = A109043(n)/2.
Dirichlet g.f.: zeta(s-1)*(1 - 1/2^s). (End)
a(n) = A001318(n) - A001318(n-1) for n > 0. - Jonathan Sondow, Jan 28 2013
a((2*n+1)*2^p - 1) = 2^p - 1 + n*A151821(p+1), p >= 0 and n >= 0. - Johannes W. Meijer, Feb 03 2013
a(n+1) = denominator(H(n, 1)), n >= 0, with H(n, 1) = 2*n/(n+1) the harmonic mean of n and 1. a(n+1) = A227042(n, 1). See the formula a(n) = n/gcd(n, 2) given above. - Wolfdieter Lang, Jul 04 2013
a(n) = numerator(n/2). - Wesley Ivan Hurt, Oct 02 2013
a(n) = numerator(1 - 2/(n+2)), n >= 0; a(n) = denominator(1 - 2/n), n >= 1. - Kival Ngaokrajang, Jul 17 2014
a(n) = Sum_{i = floor(n/2)..floor((n+1)/2)} i. - Wesley Ivan Hurt, Apr 27 2016
Euler transform of length 3 sequence [1, 2, -1]. - Michael Somos, Jan 20 2017
G.f.: x / (1 - x / (1 - 2*x / (1 + 7*x / (2 - 9*x / (7 - 4*x / (3 - 7*x / (2 + 3*x))))))). - Michael Somos, Jan 20 2017
From Peter Bala, Mar 24 2019: (Start)
a(n) = Sum_{d|n, n/d odd} phi(d), where phi(n) is the Euler totient function A000010.
O.g.f.: Sum_{n >= 1} phi(n)*x^n/(1 - x^(2*n)). (End)
a(n) = A256095(2*n,n). - Alois P. Heinz, Jan 21 2020
E.g.f.: x*(2*cosh(x) + sinh(x))/2. - Stefano Spezia, Apr 28 2023
From Ctibor O. Zizka, Oct 05 2023: (Start)
For k >= 0, a(k) = gcd(k + 1, k*(k + 1)/2).
If (k mod 4) = 0 or 2 then a(k) = (k + 1).
If (k mod 4) = 1 or 3 then a(k) = (k + 1)/2. (End)
Sum_{n=1..oo} 1/a(n)^2 = 7*Pi^2/24. - Stefano Spezia, Dec 02 2023
a(n)*a(n+1) = A000217(n). - Rémy Sigrist, Mar 19 2025

Extensions

Better description from Jud McCranie
Edited by Ralf Stephan, Jun 04 2003

A000712 Generating function = Product_{m>=1} 1/(1 - x^m)^2; a(n) = number of partitions of n into parts of 2 kinds.

Original entry on oeis.org

1, 2, 5, 10, 20, 36, 65, 110, 185, 300, 481, 752, 1165, 1770, 2665, 3956, 5822, 8470, 12230, 17490, 24842, 35002, 49010, 68150, 94235, 129512, 177087, 240840, 326015, 439190, 589128, 786814, 1046705, 1386930, 1831065, 2408658, 3157789, 4126070, 5374390
Offset: 0

Views

Author

Keywords

Comments

For n >= 1, a(n) is also the number of conjugacy classes in the automorphism group of the n-dimensional hypercube. This automorphism group is the wreath product of the cyclic group C_2 and the symmetric group S_n, its order is in sequence A000165. - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Nov 04 2001
Also, number of noncongruent matrices in GL_n(Z): each Jordan block can only have +1 or -1 on the diagonal. - Michele Dondi (blazar(AT)lcm.mi.infn.it), Jun 15 2004
a(n) = Sum (k(1)+1)*(k(2)+1)*...*(k(n)+1), where the sum is taken over all (k(1),k(2),...,k(n)) such that k(1)+2*k(2)+...+n*k(n) = n, k(i)>=0, i=1..n, cf. A104510, A077285. - Vladeta Jovovic, Apr 21 2005
Convolution of partition numbers (A000041) with itself. - Graeme McRae, Jun 07 2006
Number of one-to-one partial endofunctions on n unlabeled points. Connected components are either cycles or "lines", hence two for each size. - Franklin T. Adams-Watters, Dec 28 2006
Equals A000716: (1, 3, 9, 22, 561, 108, ...) convolved with A010815. A000716 = the number of partitions of n into parts of 3 kinds = the Euler transform of [3,3,3,...]. - Gary W. Adamson, Oct 26 2008
Paraphrasing the g.f.: 1 + 2x + 5x^2 + ... = s(x) * s(x^2) * s(x^3) * s(x^4) * ...; where s(x) = 1 + 2x + 3x^2 + 4x^3 + ... is (up to a factor x) the g.f. of A000027. - Gary W. Adamson, Apr 01 2010
Also equals number of partitions of 2n in which the odd parts appear as many times in even as in odd positions. - Wouter Meeussen, Apr 17 2013
Also number of ordered pairs (R,S) with R a partition of r, S a partition of s, and r+s=n; see example. This corresponds to the formula a(n) = sum(r+s==n, p(r)*p(s) ) = Sum_{k=0..n} p(k)*p(n-k). - Joerg Arndt, Apr 29 2013
Also the number of all multi-graphs with exactly n-edges and with vertex degrees 1 or 2. - Ebrahim Ghorbani, Dec 02 2013
If one decomposes k-permutations into cycles and so-called paths, the number of different type of decompositions equals to a(k); see the paper by Chen, Ghorbani, and Wong. - Ebrahim Ghorbani, Dec 02 2013
Let T(n,k) be the number of partitions of n having parts 1 through k of two kinds, with T(n,0) = A000041(n), the number of partitions of n. Then a(n) = T(n,0) + T(n-1,1) + T(n-2,2) + T(n-3,3) + ... - Gregory L. Simay, May 18 2019
Also the number of orbits of projections in the partition monoid P_n under conjugation by permutations. - James East, Jul 21 2020

Examples

			Assume there are integers of two kinds: k and k'; then a(3) = 10 since 3 has the following partitions into parts of two kinds: 111, 111', 11'1', 1'1'1', 12, 1'2, 12', 1'2', 3, and 3'. - _W. Edwin Clark_, Jun 24 2011
There are a(4)=20 partitions of 4 into 2 sorts of parts. Here p:s stands for "part p of sort s":
01:  [ 1:0  1:0  1:0  1:0  ]
02:  [ 1:0  1:0  1:0  1:1  ]
03:  [ 1:0  1:0  1:1  1:1  ]
04:  [ 1:0  1:1  1:1  1:1  ]
05:  [ 1:1  1:1  1:1  1:1  ]
06:  [ 2:0  1:0  1:0  ]
07:  [ 2:0  1:0  1:1  ]
08:  [ 2:0  1:1  1:1  ]
09:  [ 2:0  2:0  ]
10:  [ 2:0  2:1  ]
11:  [ 2:1  1:0  1:0  ]
12:  [ 2:1  1:0  1:1  ]
13:  [ 2:1  1:1  1:1  ]
14:  [ 2:1  2:1  ]
15:  [ 3:0  1:0  ]
16:  [ 3:0  1:1  ]
17:  [ 3:1  1:0  ]
18:  [ 3:1  1:1  ]
19:  [ 4:0  ]
20:  [ 4:1  ]
- _Joerg Arndt_, Apr 28 2013
The a(4)=20 ordered pairs (R,S) of partitions for n=4 are
  ([4], [])
  ([3, 1], [])
  ([2, 2], [])
  ([2, 1, 1], [])
  ([1, 1, 1, 1], [])
  ([3], [1])
  ([2, 1], [1])
  ([1, 1, 1], [1])
  ([2], [2])
  ([2], [1, 1])
  ([1, 1], [2])
  ([1, 1], [1, 1])
  ([1], [3])
  ([1], [2, 1])
  ([1], [1, 1, 1])
  ([], [4])
  ([], [3, 1])
  ([], [2, 2])
  ([], [2, 1, 1])
  ([], [1, 1, 1, 1])
This list was created with the Sage command
   for P in PartitionTuples(2,4) : print P;
- _Joerg Arndt_, Apr 29 2013
G.f. = 1 + 2*x + 5*x^2 + 10*x^3 + 20*x^4 + 36*x^5 + 65*x^6 + 110*x^7 + 185*x^8 + ...
		

References

  • H. Gupta et al., Tables of Partitions. Royal Society Mathematical Tables, Vol. 4, Cambridge Univ. Press, 1958, p. 90.
  • J. Riordan, Combinatorial Identities, Wiley, 1968, p. 199.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 1, 1999; see Proposition 2.5.2 on page 78.

Crossrefs

Cf. A000165, A000041, A002107 (reciprocal of g.f.).
Cf. A002720.
Cf. A000716, A010815. - Gary W. Adamson, Oct 26 2008
Row sums of A175012. - Gary W. Adamson, Apr 03 2010
Column k=2 of A144064.

Programs

  • Haskell
    a000712 = p a008619_list where
       p _          0 = 1
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Nov 06 2012
    
  • Julia
    # DedekindEta is defined in A000594.
    A000712List(len) = DedekindEta(len, -2)
    A000712List(39) |> println # Peter Luschny, Mar 09 2018
    
  • Maple
    with(combinat): A000712:= n-> add(numbpart(k)*numbpart(n-k), k=0..n): seq(A000712(n), n=0..40); # Emeric Deutsch
  • Mathematica
    CoefficientList[ Series[ Product[1/(1 - x^n)^2, {n, 40}], {x, 0, 37}], x]; (* Robert G. Wilson v, Feb 03 2005 *)
    Table[Count[Partitions[2*n], q_ /; Tr[(-1)^Mod[Flatten[Position[q, ?OddQ]], 2]] === 0], {n, 12}] (* _Wouter Meeussen, Apr 17 2013 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ x]^-2, {x, 0, n}]; (* Michael Somos, Oct 12 2015 *)
    Table[Length@IntegerPartitions[n, All, Range@n~Join~Range@n], {n, 0, 15}] (* Robert Price, Jun 15 2020 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( 1 / eta(x + A)^2, n))}; /* Michael Somos, Nov 14 2002 */
    
  • PARI
    Vec(1/eta('x+O('x^66))^2) /* Joerg Arndt, Jun 25 2011 */
    
  • Python
    from sympy import npartitions
    def A000712(n): return (sum(npartitions(k)*npartitions(n-k) for k in range(n+1>>1))<<1) + (0 if n&1 else npartitions(n>>1)**2) # Chai Wah Wu, Sep 25 2023
  • SageMath
    # uses[EulerTransform from A166861]
    a = BinaryRecurrenceSequence(0, 1, 2, 2)
    b = EulerTransform(a)
    print([b(n) for n in range(40)]) # Peter Luschny, Nov 11 2020
    

Formula

a(n) = Sum_{k=0..n} p(k)*p(n-k), where p(n) = A000041(n).
Euler transform of period 1 sequence [ 2, 2, 2, ...]. - Michael Somos, Jul 22 2003
a(n) = A006330(n) + A001523(n). - Michael Somos, Jul 22 2003
a(0) = 1, a(n) = (1/n)*Sum_{k=0..n-1} 2*a(k)*sigma_1(n-k). - Joerg Arndt, Feb 05 2011
a(n) ~ (1/12)*3^(1/4)*n^(-5/4)*exp((2/3)*sqrt(3)*Pi*sqrt(n)). - Joe Keane (jgk(AT)jgk.org), Sep 13 2002
G.f.: Product_{i>=1} (1 + x^i)^(2*A001511(i)) (see A000041). - Jon Perry, Jun 06 2004
More precise asymptotics: a(n) ~ exp(2*Pi*sqrt(n/3)) / (4*3^(3/4)*n^(5/4)) * (1 - (Pi/(12*sqrt(3)) + 15*sqrt(3)/(16*Pi)) / sqrt(n) + (Pi^2/864 + 315/(512*Pi^2) + 35/192)/n). - Vaclav Kotesovec, Jan 22 2017
From Peter Bala, Jan 26 2016: (Start)
a(n) is odd iff n = 2*m and p(m) is odd.
a(n) = (2/n)*Sum_{k = 0..n} k*p(k)*p(n-k) for n >= 1.
Conjecture: : a(n) is divisible by 5 when n is congruent to 2, 3 or 4 modulo 5. (End)
Conjecture is proved in Hammond and Lewis. - Yen-chi R. Lin, Jun 24 2024
G.f.: exp(2*Sum_{k>=1} x^k/(k*(1 - x^k))). - Ilya Gutkovskiy, Feb 06 2018
With the convention that a(n) = 0 for n < 0 we have the recurrence a(n) = g(n) + Sum_{k >= 1} (-1)^(k+1)*(2*k + 1)*a(n - k*(k + 1)/2), where g(n) = (-1)^m if n = m*(3*m - 1)/2 is a generalized pentagonal number (A001318) else g(n) = 0. For example, n = 7 = -2*(3*(-2) - 1)/2 is a pentagonal number, g(7) = 1, and so a(7) = 1 + 3*a(6) - 5*a(4) + 7*a(1) = 1 + 195 - 100 + 14 = 110. - Peter Bala, Apr 06 2022
a(n) = p(n/2) + Sum_{k \in Z, k != 0} (-1)^{k-1} a(n-k^2), here p(n) = A000041(n) and p(x) = 0 when x is not an integer. - Yen-chi R. Lin, Jun 24 2024
Conjecture: a(25*n + 23) is divisible by 25 (checked for n < 400). - Peter Bala, Jan 13 2025

Extensions

More terms from Joe Keane (jgk(AT)jgk.org), Nov 17 2001
More terms from Michele Dondi (blazar(AT)lcm.mi.infn.it), Jun 15 2004
Definition rewritten by N. J. A. Sloane, Apr 02 2022

A005449 Second pentagonal numbers: a(n) = n*(3*n + 1)/2.

Original entry on oeis.org

0, 2, 7, 15, 26, 40, 57, 77, 100, 126, 155, 187, 222, 260, 301, 345, 392, 442, 495, 551, 610, 672, 737, 805, 876, 950, 1027, 1107, 1190, 1276, 1365, 1457, 1552, 1650, 1751, 1855, 1962, 2072, 2185, 2301, 2420, 2542, 2667, 2795, 2926, 3060, 3197, 3337, 3480
Offset: 0

Views

Author

Keywords

Comments

Number of edges in the join of the complete graph and the cycle graph, both of order n, K_n * C_n. - Roberto E. Martinez II, Jan 07 2002
Also number of cards to build an n-tier house of cards. - Martin Wohlgemuth, Aug 11 2002
The modular form Delta(q) = q*Product_{n>=1} (1-q^n)^24 = q*(1 + Sum_{n>=1} (-1)^n*(q^(n*(3*n-1)/2)+q^(n*(3*n+1)/2)))^24 = q*(1 + Sum_{n>=1} A033999(n)*(q^A000326(n)+q^a(n)))^24. - Jonathan Vos Post, Mar 15 2006
Row sums of triangle A134403.
Bisection of A001318. - Omar E. Pol, Aug 22 2011
Sequence found by reading the line from 0 in the direction 0, 7, ... and the line from 2 in the direction 2, 15, ... in the square spiral whose vertices are the generalized pentagonal numbers, A001318. - Omar E. Pol, Sep 08 2011
A general formula for the n-th second k-gonal number is given by T(n, k) = n*((k-2)*n+k-4)/2, n>=0, k>=5. - Omar E. Pol, Aug 04 2012
Partial sums give A006002. - Denis Borris, Jan 07 2013
A002260 is the following array A read by antidiagonals:
0, 1, 2, 3, 4, 5, ...
0, 1, 2, 3, 4, 5, ...
0, 1, 2, 3, 4, 5, ...
0, 1, 2, 3, 4, 5, ...
0, 1, 2, 3, 4, 5, ...
0, 1, 2, 3, 4, 5, ...
and a(n) is the hook sum: Sum_{k=0..n} A(n,k) + Sum_{r=0..n-1} A(r,n). - R. J. Mathar, Jun 30 2013
From Klaus Purath, May 13 2021: (Start)
This sequence and A000326 provide all integers m such that 24*m + 1 is a square. The union of the two sequences is A001318.
If A is a sequence satisfying the recurrence t(n) = 3*t(n-1) - 2*t(n-2) with the initial values either A(0) = 1, A(1) = n + 2 or A(0) = -1, A(1) = n - 1, then a(n) = (A(i)^2 - A(i-1)*A(i+1))/2^i + n^2 for i>0. (End)
a(n+1) is the number of Dyck paths of size (3,3n+2), i.e., the number of NE lattice paths from (0,0) to (3,3n+2) which stay above the line connecting these points. - Harry Richman, Jul 13 2021
Binomial transform of [0, 2, 3, 0, 0, 0, ...], being a(n) = 2*binomial(n,1) + 3*binomial(n,2). a(3) = 15 = [0, 2, 3, 0] dot [1, 3, 3, 1] = [0 + 6 + 9 + 0]. - Gary W. Adamson, Dec 17 2022
a(n) is the sum of longest side length of all nondegenerate integer-sided triangles with shortest side length n and middle side length (n + 1), n > 0. - Torlach Rush, Feb 04 2024

Examples

			From _Omar E. Pol_, Aug 22 2011: (Start)
Illustration of initial terms:
                                               O
                                             O O
                                 O         O O O
                               O O       O O O O
                     O       O O O     O O O O O
                   O O     O O O O     O O O O O
           O     O O O     O O O O     O O O O O
         O O     O O O     O O O O     O O O O O
    O    O O     O O O     O O O O     O O O O O
    O    O O     O O O     O O O O     O O O O O
    -    ---     -----     -------     ---------
    2     7        15         26           40
(End)
		

References

  • Henri Cohen, A Course in Computational Algebraic Number Theory, vol. 138 of Graduate Texts in Mathematics, Springer-Verlag, 2000.

Crossrefs

Cf. A016789 (first differences), A006002 (partial sums).
The generalized pentagonal numbers b*n+3*n*(n-1)/2, for b = 1 through 12, form sequences A000326, this sequence, A045943, A115067, A140090, A140091, A059845, A140672-A140675, A151542.
Cf. numbers of the form n*(n*k-k+4)/2 listed in A226488 (this sequence is the case k=3).
Cf. numbers of the form n*((2*k+1)*n+1)/2 listed in A022289 (this sequence is the case k=1).

Programs

Formula

a(n) = A110449(n, 1) for n>0.
G.f.: x*(2+x)/(1-x)^3. E.g.f.: exp(x)*(2*x + 3*x^2/2). a(n) = n*(3*n + 1)/2. a(-n) = A000326(n). - Michael Somos, Jul 18 2003
a(n) = A001844(n) - A000217(n+1) = A101164(n+2,2) for n>0. - Reinhard Zumkeller, Dec 03 2004
a(n) = Sum_{j=1..n} n+j. - Zerinvary Lajos, Sep 12 2006
a(n) = A126890(n,n). - Reinhard Zumkeller, Dec 30 2006
a(n) = 2*C(3*n,4)/C(3*n,2), n>=1. - Zerinvary Lajos, Jan 02 2007
a(n) = A000217(n) + A000290(n). - Zak Seidov, Apr 06 2008
a(n) = a(n-1) + 3*n - 1 for n>0, a(0)=0. - Vincenzo Librandi, Nov 18 2010
a(n) = A129267(n+5,n). - Philippe Deléham, Dec 21 2011
a(n) = 2*A000217(n) + A000217(n-1). - Philippe Deléham, Mar 25 2013
a(n) = A130518(3*n+1). - Philippe Deléham, Mar 26 2013
a(n) = (12/(n+2)!)*Sum_{j=0..n} (-1)^(n-j)*binomial(n,j)*j^(n+2). - Vladimir Kruchinin, Jun 04 2013
a(n) = floor(n/(1-exp(-2/(3*n)))) for n>0. - Richard R. Forberg, Jun 22 2013
a(n) = Sum_{i=1..n} (3*i - 1) for n >= 1. - Wesley Ivan Hurt, Oct 11 2013 [Corrected by Rémi Guillaume, Oct 24 2024]
a(n) = (A000292(6*n+k+1)-A000292(k))/(6*n+1) - A000217(3*n+k+1), for any k >= 0. - Manfred Arens, Apr 26 2015
Sum_{n>=1} 1/a(n) = 6 - Pi/sqrt(3) - 3*log(3) = 0.89036376976145307522... . - Vaclav Kotesovec, Apr 27 2016
a(n) = A000217(2*n) - A000217(n). - Bruno Berselli, Sep 21 2016
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*Pi/sqrt(3) + 4*log(2) - 6. - Amiram Eldar, Jan 18 2021
From Klaus Purath, May 13 2021: (Start)
Partial sums of A016789 for n > 0.
a(n) = 3*n^2 - A000326(n).
a(n) = A000326(n) + n.
a(n) = A002378(n) + A000217(n-1) for n >= 1. [Corrected by Rémi Guillaume, Aug 14 2024] (End)
From Klaus Purath, Jul 14 2021: (Start)
b^2 = 24*a(n) + 1 we get by b^2 = (a(n+1) - a(n-1))^2 = (a(2*n)/n)^2.
a(2*n) = n*(a(n+1) - a(n-1)), n > 0.
a(2*n+1) = n*(a(n+1) - a(n)). (End)
A generalization of Lajos' formula, dated Sep 12 2006, follows. Let SP(k,n) = the n-th second k-gonal number. Then SP(2k+1,n) = Sum_{j=1..n} (k-1)*n+j+k-2. - Charlie Marion, Jul 13 2024
a(n) = Sum_{k = 0..3*n} (-1)^(n+k+1) * binomial(k, 2) * binomial(3*n+k, 2*k). - Peter Bala, Nov 03 2024
For integer m, (6*m + 1)^2*a(n) + a(m) = a((6*m+1)*n + m). - Peter Bala, Jan 09 2025

A045943 Triangular matchstick numbers: a(n) = 3*n*(n+1)/2.

Original entry on oeis.org

0, 3, 9, 18, 30, 45, 63, 84, 108, 135, 165, 198, 234, 273, 315, 360, 408, 459, 513, 570, 630, 693, 759, 828, 900, 975, 1053, 1134, 1218, 1305, 1395, 1488, 1584, 1683, 1785, 1890, 1998, 2109, 2223, 2340, 2460, 2583, 2709, 2838, 2970, 3105, 3243, 3384, 3528
Offset: 0

Views

Author

Keywords

Comments

Also, 3 times triangular numbers, a(n) = 3*A000217(n).
In the 24-bit RGB color cube, the number of color-lattice-points in r+g+b = n planes at n < 256 equals the triangular numbers. For n = 256, ..., 765 the number of legitimate color partitions is less than A000217(n) because {r,g,b} components cannot exceed 255. For n = 256, ..., 511, the number of non-color partitions are computable with A045943(n-255), while for n = 512, ..., 765, the number of color points in r+g+b planes equals A000217(765-n). - Labos Elemer, Jun 20 2005
If a 3-set Y and an (n-3)-set Z are disjoint subsets of an n-set X then a(n-3) is the number of 3-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 19 2007
a(n) is also the smallest number that may be written both as the sum of n-1 consecutive positive integers and n consecutive positive integers. - Claudio Meller, Oct 08 2010
For n >= 3, a(n) equals 4^(2+n)*Pi^(1 - n) times the coefficient of zeta(3) in the following integral with upper bound Pi/4 and lower bound 0: int x^(n+1) tan x dx. - John M. Campbell, Jul 17 2011
The difference a(n)-a(n-1) = 3*n, for n >= 1. - Stephen Balaban, Jul 25 2011 [Comment clarified by N. J. A. Sloane, Aug 01 2024]
Sequence found by reading the line from 0, in the direction 0, 3, ..., and the same line from 0, in the direction 0, 9, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. This is one of the orthogonal axes of the spiral; the other is A032528. - Omar E. Pol, Sep 08 2011
A005449(a(n)) = A000332(3n + 3) = C(3n + 3, 4), a second pentagonal number of triangular matchstick number index number. Additionally, a(n) - 2n is a pentagonal number (A000326). - Raphie Frank, Dec 31 2012
Sum of the numbers from n to 2n. - Wesley Ivan Hurt, Nov 24 2015
Number of orbits of Aut(Z^7) as function of the infinity norm (n+1) of the representative integer lattice point of the orbit, when the cardinality of the orbit is equal to 5376 or 17920 or 20160. - Philippe A.J.G. Chevalier, Dec 28 2015
Also the number of 4-cycles in the (n+4)-triangular honeycomb acute knight graph. - Eric W. Weisstein, Jul 27 2017
Number of terms less than 10^k, k=0,1,2,3,...: 1, 3, 8, 26, 82, 258, 816, 2582, 8165, 25820, 81650, 258199, 816497, 2581989, 8164966, ... - Muniru A Asiru, Jan 24 2018
Numbers of the form 3*m*(2*m + 1) for m = 0, -1, 1, -2, 2, -3, 3, ... - Bruno Berselli, Feb 26 2018
Partial sums of A008585. - Omar E. Pol, Jun 20 2018
Column 1 of A273464. (Number of ways to select a unit lozenge inside an isosceles triangle of side length n; all vertices on a hexagonal lattice.) - R. J. Mathar, Jul 10 2019
Total number of pips in the n-th suit of a double-n domino set. - Ivan N. Ianakiev, Aug 23 2020

Examples

			From _Stephen Balaban_, Jul 25 2011: (Start)
T(n), the triangular numbers = number of nodes,
a(n-1) = number of edges in the T(n) graph:
       o    (T(1) = 1, a(0) = 0)
       o
      / \   (T(2) = 3, a(1) = 3)
     o - o
       o
      / \
     o - o  (T(3) = 6, a(2) = 9)
    / \ / \
   o - o - o
... [Corrected by _N. J. A. Sloane_, Aug 01 2024] (End)
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 543.

Crossrefs

The generalized pentagonal numbers b*n+3*n*(n-1)/2, for b = 1 through 12, form sequences A000326, A005449, A045943, A115067, A140090, A140091, A059845, A140672, A140673, A140674, A140675, A151542.
A diagonal of A010027.
Orbits of Aut(Z^7) as function of the infinity norm A000579, A154286, A102860, A002412, A115067, A008585, A005843, A001477, A000217.
Cf. A027480 (partial sums).
Cf. A002378 (3-cycles in triangular honeycomb acute knight graph), A028896 (5-cycles), A152773 (6-cycles).
This sequence: Sum_{k = n..2*n} k.
Cf. A304993: Sum_{k = n..2*n} k*(k+1)/2.
Cf. A050409: Sum_{k = n..2*n} k^2.
Similar sequences are listed in A316466.

Programs

Formula

a(n) is the sum of n+1 integers starting from n, i.e., 1+2, 2+3+4, 3+4+5+6, 4+5+6+7+8, etc. - Jon Perry, Jan 15 2004
a(n) = A126890(n+1,n-1) for n>1. - Reinhard Zumkeller, Dec 30 2006
a(n) + A145919(3*n+3) = 0. - Matthew Vandermast, Oct 28 2008
a(n) = A000217(2*n) - A000217(n-1); A179213(n) <= a(n). - Reinhard Zumkeller, Jul 05 2010
a(n) = a(n-1)+3*n, n>0. - Vincenzo Librandi, Nov 18 2010
G.f.: 3*x/(1-x)^3. - Bruno Berselli, Jan 21 2011
a(n) = A005448(n+1) - 1. - Omar E. Pol, Oct 03 2011
a(n) = A001477(n)+A000290(n)+A000217(n). - J. M. Bergot, Dec 08 2012
a(n) = 3*a(n-1)-3*a(n-2)+a(n-3) for n>2. - Wesley Ivan Hurt, Nov 24 2015
a(n) = A027480(n)-A027480(n-1). - Peter M. Chema, Jan 18 2017.
2*a(n)+1 = A003215(n). - Miquel Cerda, Jan 22 2018
a(n) = T(2*n) - T(n-1), where T(n) = A000217(n). In general, T(k)*T(n) = Sum_{i=0..k-1} (-1)^i*T((k-i)*(n-i)). - Charlie Marion, Dec 06 2020
E.g.f.: 3*exp(x)*x*(2 + x)/2. - Stefano Spezia, May 19 2021
From Amiram Eldar, Jan 10 2022: (Start)
Sum_{n>=1} 1/a(n) = 2/3.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*(2*log(2)-1)/3. (End)
Product_{n>=1} (1 - 1/a(n)) = -(3/(2*Pi))*cos(sqrt(11/3)*Pi/2). - Amiram Eldar, Feb 21 2023

A001082 Generalized octagonal numbers: k*(3*k-2), k=0, +- 1, +- 2, +-3, ...

Original entry on oeis.org

0, 1, 5, 8, 16, 21, 33, 40, 56, 65, 85, 96, 120, 133, 161, 176, 208, 225, 261, 280, 320, 341, 385, 408, 456, 481, 533, 560, 616, 645, 705, 736, 800, 833, 901, 936, 1008, 1045, 1121, 1160, 1240, 1281, 1365, 1408, 1496, 1541, 1633, 1680, 1776, 1825, 1925, 1976
Offset: 1

Views

Author

Keywords

Comments

Numbers of the form 3*m^2+2*m, m an integer.
3*a(n) + 1 is a perfect square.
a(n) mod 10 belongs to a periodic sequence: 0, 1, 5, 8, 6, 1, 3, 0, 6, 5, 5, 6, 0, 3, 1, 6, 8, 5, 1, 0. - Mohamed Bouhamida, Sep 04 2009
A089801 is the characteristic function. - R. J. Mathar, Oct 07 2011
Exponents of powers of q in one form of the quintuple product identity. (-x^-2 + 1) * q^0 + (x^-3 - x) * q^1 + (-x^-5 + x^3) * q^5 + (x^-6 - x^4) * q^8 + ... = Sum_{n>=0} q^(3*n^2 + 2*n) * (x^(3*n) - x^(-3*n - 2)) = Product_{k>0} (1 - x * q^(2*k - 1)) * (1 - x^-1 * q^(2*k - 1)) * (1 - q^(2*k)) * (1 - x^2 * q^(4*k)) * (1 - x^-2 * q^(4*k - 4)). - Michael Somos, Dec 21 2011
The offset 0 would also be valid here, all other entries of generalized k-gonal numbers have offset 0 (see cross references). - Omar E. Pol, Jan 12 2013
Also, x values of the Diophantine equation x(x+3)+(x+1)(x+2) = (x+y)^2+(x-y)^2. - Bruno Berselli, Mar 29 2013
Numbers n such that Sum_{i=1..n} 2*i*(n-i)/n is an integer (the addend is the harmonic mean of i and n-i). - Wesley Ivan Hurt, Sep 14 2014
Equivalently, integers of the form m*(m+2)/3 (nonnegative values of m are listed in A032766). - Bruno Berselli, Jul 18 2016
Exponents of q in the expansion of Sum_{n >= 0} ( q^n * Product_{k = 1..n} (1 - q^(2*k-1)) ) = 1 + q - q^5 - q^8 + q^16 + q^21 - - + + .... - Peter Bala, Dec 03 2020
Exponents of q in the expansion of Product_{n >= 1} (1 - q^(6*n))*(1 + q^(6*n-1))*(1 + q^(6*n-5)) = 1 + q + q^5 + q^8 + q^16 + q^21 + .... - Peter Bala, Dec 09 2020
Exponents of q in the expansion of Product_{n >= 1} (1 - q^n)^2*(1 - q^(4*n))^2 /(1 - q^(2*n)) = 1 - 2*q + 4*q^5 - 5*q^8 + 7*q^16 - + ... (a consequence of the quintuple product identity). The series coefficients are a signed version of A001651. - Peter Bala, Feb 16 2021
From Peter Bala, Nov 26 2024: (Start)
Apart from the first two terms, the exponents of q in the expansion of Sum_{n >= 1} q^(3*n+2) * (Product_{k = 2..n} 1 - q^(2*k-1)) = q^5 + q^8 - q^16 - q^21 + + - - ... (in Andrews, equation 8, replace q with q^2 and set x = q).
Exponents of q^2 in the expansion of Sum_{n >= 0} q^n / (Product_{k = 1..n+1 } 1 + q^(2*k-1)) = 1 + (q^2)^1 - (q^2)^5 - (q^2)^8 + (q^2)^16 + (q^2)^21 - - + + ... (Chen, equation 22). (End)

Examples

			For the ninth comment: 65 is in the sequence because 65 = 13*(13+2)/3 or also 65 = -15*(-15+2)/3. - _Bruno Berselli_, Jul 18 2016
		

Crossrefs

Partial sums of A022998.
Column 4 of A195152. A045944.
Sequences of generalized k-gonal numbers: A001318 (k=5), A000217 (k=6), A085787 (k=7), this sequence (k=8), A118277 (k=9), A074377 (k=10), A195160 (k=11), A195162 (k=12), A195313 (k=13), A195818 (k=14), A277082 (k=15), A274978 (k=16), A303305 (k=17), A274979 (k=18), A303813 (k=19), A218864 (k=20), A303298 (k=21), A303299 (k=22), A303303 (k=23), A303814 (k=24), A303304 (k=25), A316724 (k=26), A316725 (k=27), A303812 (k=28), A303815 (k=29), A316729 (k=30).
Cf. sequences of the form m*(m+k)/(k+1) listed in A274978. [Bruno Berselli, Jul 25 2016]

Programs

  • Haskell
    a001082 n = a001082_list !! n
    a001082_list = scanl (+) 0 $ tail a022998_list
    -- Reinhard Zumkeller, Mar 31 2012
    
  • Magma
    [n^2 - n - Floor(n/2)^2 : n in [1..50]]; // Wesley Ivan Hurt, Sep 14 2014
  • Maple
    seq(n*(n-1)-floor(n/2)^2, n=1..51); # Gary Detlefs, Feb 23 2010
  • Mathematica
    Table[If[EvenQ[n], n*(3*n-4)/4, (n-1) (3*n+1)/4], {n, 100}]
    LinearRecurrence[{1,2,-2,-1,1},{0,1,5,8,16},60] (* Harvey P. Dale, Feb 03 2024 *)
  • PARI
    {a(n) = if( n%2, (n-1) * (3*n + 1) / 4, n * (3*n - 4) / 4)};
    

Formula

a(n) = n*(3*n-4)/4 if n even, (n-1)*(3*n+1)/4 if n odd.
a(n) = n^2 - n - floor(n/2)^2.
G.f.: Sum_{n>=0} (-1)^n*[x^(a(2n+1)) + x^(a(2n+2))] = 1/1 - (x-x^2)/1 - (x^2-x^4)/1 - (x^3-x^6)/1 - ... - (x^k - x^(2k))/1 - ... (continued fraction where k=1..inf). - Paul D. Hanna, Aug 16 2002
a(n+1) = ceiling(n/2)^2 + A046092(floor(n/2)).
a(2n) = n(3n-2) = A000567(n), a(2n+1) = n(3n+2) = A045944(n). - Mohamed Bouhamida, Nov 06 2007
O.g.f.: -x^2*(x^2+4*x+1)/((x-1)^3*(1+x)^2). - R. J. Mathar, Apr 15 2008
a(n) = n^2+n-ceiling(n/2)^2 with offset 0 and a(0)=0. - Gary Detlefs, Feb 23 2010
a(n) = (6*n^2-6*n-1-(2*n-1)*(-1)^n)/8. - Luce ETIENNE, Dec 11 2014
E.g.f.: (3*x^2*exp(x) + x*exp(-x) - sinh(x))/4. - Ilya Gutkovskiy, Jul 15 2016
Sum_{n>=2} 1/a(n) = (9 + 2*sqrt(3)*Pi)/12. - Vaclav Kotesovec, Oct 05 2016
Sum_{n>=2} (-1)^n/a(n) = 3*log(3)/2 - 3/4. - Amiram Eldar, Feb 28 2022

Extensions

New sequence name from Matthew Vandermast, Apr 10 2003
Editorial changes by N. J. A. Sloane, Feb 03 2012
Edited by Omar E. Pol, Jun 09 2012

A048651 Decimal expansion of Product_{k >= 1} (1 - 1/2^k).

Original entry on oeis.org

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

Views

Author

Keywords

Comments

This is the limiting probability that a large random binary matrix is nonsingular (cf. A002884).
This constant is very close to 2^(13/24) * sqrt(Pi/log(2)) / exp(Pi^2/(6*log(2))) = 0.288788095086602421278899775042039398383022429351580356839... - Vaclav Kotesovec, Aug 21 2018
This constant is irrational (see Penn link). - Paolo Xausa, Dec 09 2024

Examples

			(1/2)*(3/4)*(7/8)*(15/16)*... = 0.288788095086602421278899721929230780088911904840685784114741...
		

References

  • Steven R. Finch, Mathematical Constants, Cambridge, 2003, pp. 318, 354-361.

Crossrefs

Programs

  • Mathematica
    RealDigits[ Product[1 - 1/2^i, {i, 100}], 10, 111][[1]] (* Robert G. Wilson v, May 25 2011 *)
    RealDigits[QPochhammer[1/2], 10, 100][[1]] (* Jean-François Alcover, Nov 18 2015 *)
  • PARI
    default(realprecision, 20080); x=prodinf(k=1, -1/2^k, 1); x*=10; for (n=0, 20000, d=floor(x); x=(x-d)*10; write("b048651.txt", n, " ", d)); \\ Harry J. Smith, May 07 2009

Formula

exp(-Sum_{k>0} sigma_1(k)/k*2^(-k)) = exp(-Sum_{k>0} A000203(k)/k*2^(-k)). - Hieronymus Fischer, Jul 28 2007
From Hieronymus Fischer, Aug 13 2007: (Start)
Equals lim inf Product_{k=0..floor(log_2(n))} floor(n/2^k)*2^k/n for n->oo.
Equals lim inf A098844(n)/n^(1+floor(log_2(n)))*2^(1/2*(1+floor(log_2(n)))*floor(log_2(n))) for n->oo.
Equals lim inf A098844(n)/n^(1+floor(log_2(n)))*2^A000217(floor(log_2(n))) for n->oo.
Equals lim inf A098844(n)/(n+1)^((1+log_2(n+1))/2) for n->oo.
Equals (1/2)*exp(-Sum_{n>0} 2^(-n)*Sum_{k|n} 1/(k*2^k)). (End)
Limit of A177510(n)/A000079(n-1) as n->infinity (conjecture). - Mats Granvik, Mar 27 2011
Product_{k >= 1} (1-1/2^k) = (1/2; 1/2){infinity}, where (a;q){infinity} is the q-Pochhammer symbol. - G. C. Greubel, Nov 27 2015
exp(Sum_{n>=1}(1/n/(1 - 2^n))) (according to Mathematica). - Mats Granvik, Sep 07 2016
(Sum_{k>0} (4^k-1)/(Product_{i=1..k} ((4^i-1)*(2*4^i-1))))*2 = 2/7 + 2/(3*7*31) + 2/(3*7*15*31*127)+2/(3*7*15*31*63*127*511) + ... (conjecture). - Werner Schulte, Dec 22 2016
Equals Sum_{k=-oo..oo} (-1)^k/2^((3*k+1)*k/2) (by Euler's pentagonal number theorem). - Amiram Eldar, Aug 13 2020
From Peter Bala, Dec 15 2020: (Start)
Constant C = Sum_{n >= 0} (-1)^n/( Product_{k = 1..n} (2^k - 1) ). The above conjectural result by Schulte follows by adding terms of this series in pairs.
C = (1/2)*Sum_{n >= 0} (-1/2)^n/( Product_{k = 1..n} (2^k - 1) ).
C = (3/8)*Sum_{n >= 0} (-1/4)^n/( Product_{k = 1..n} (2^k - 1) ).
1/C = Sum_{n >= 0} 2^(n*(n-1)/2)/( Product_{k = 1..n} (2^k - 1) ).
C = 1 - Sum_{n >= 0} (1/2)^(n+1)*Product_{k = 1..n} (1 - 1/2^k).
This latter identity generalizes as:
C = Sum_{n >= 0} (1/4)^(n+1)*Product_{k = 1..n} (1 - 1/2^k),
3*C = 1 - Sum_{n >= 0} (1/8)^(n+1)*Product_{k = 1..n} (1 - 1/2^k),
3*7*C = 6 + Sum_{n >= 0} (1/16)^(n+1)*Product_{k = 1..n} (1 - 1/2^k),
3*7*15*C = 91 - Sum_{n >= 0} (1/32)^(n+1)*Product_{k = 1..n} (1 - 1/2^k),
and so on, where the sequence [1, 0, 1, 6, 91, ...] is A005327.
(End)
From Amiram Eldar, Feb 19 2022: (Start)
Equals sqrt(2*Pi/log(2)) * exp(log(2)/24 - Pi^2/(6*log(2))) * Product_{k>=1} (1 - exp(-4*k*Pi^2/log(2))) (McIntosh, 1995).
Equals Sum_{n>=0} (-1)^n/A005329(n).
Equals exp(-A335764). (End)
Equals 1/A065446. - Hugo Pfoertner, Nov 23 2024

Extensions

Corrected by Hieronymus Fischer, Jul 28 2007

A211422 Number of ordered triples (w,x,y) with all terms in {-n,...,0,...,n} and w^2 + x*y = 0.

Original entry on oeis.org

1, 9, 17, 25, 41, 49, 57, 65, 81, 105, 113, 121, 137, 145, 153, 161, 193, 201, 225, 233, 249, 257, 265, 273, 289, 329, 337, 361, 377, 385, 393, 401, 433, 441, 449, 457, 505, 513, 521, 529, 545, 553, 561, 569, 585, 609, 617, 625, 657, 713, 753, 761
Offset: 0

Views

Author

Clark Kimberling, Apr 10 2012

Keywords

Comments

Suppose that S={-n,...,0,...,n} and that f(w,x,y,n) is a function, where w,x,y are in S. The number of ordered triples (w,x,y) satisfying f(w,x,y,n)=0, regarded as a function of n, is a sequence t of nonnegative integers. Sequences such as t/4 may also be integer sequences for all except certain initial values of n. In the following guide, such sequences are indicated in the related sequences column and may be included in the corresponding Mathematica programs.
...
sequence... f(w,x,y,n) ..... related sequences
A211415 ... w^2+x*y-1 ...... t+2, t/4, (t/4-1)/4
A211422 ... w^2+x*y ........ (t-1)/8, A120486
A211423 ... w^2+2x*y ....... (t-1)/4
A211424 ... w^2+3x*y ....... (t-1)/4
A211425 ... w^2+4x*y ....... (t-1)/4
A211426 ... 2w^2+x*y ....... (t-1)/4
A211427 ... 3w^2+x*y ....... (t-1)/4
A211428 ... 2w^2+3x*y ...... (t-1)/4
A211429 ... w^3+x*y ........ (t-1)/4
A211430 ... w^2+x+y ........ (t-1)/2
A211431 ... w^3+(x+y)^2 .... (t-1)/2
A211432 ... w^2-x^2-y^2 .... (t-1)/8
A003215 ... w+x+y .......... (t-1)/2, A045943
A202253 ... w+2x+3y ........ (t-1)/2, A143978
A211433 ... w+2x+4y ........ (t-1)/2
A211434 ... w+2x+5y ........ (t-1)/4
A211435 ... w+4x+5y ........ (t-1)/2
A211436 ... 2w+3x+4y ....... (t-1)/2
A211435 ... 2w+3x+5y ....... (t-1)/2
A211438 ... w+2x+2y ....... (t-1)/2, A118277
A001844 ... w+x+2y ......... (t-1)/4, A000217
A211439 ... w+3x+3y ........ (t-1)/2
A211440 ... 2w+3x+3y ....... (t-1)/2
A028896 ... w+x+y-1 ........ t/6, A000217
A211441 ... w+x+y-2 ........ t/3, A028387
A182074 ... w^2+x*y-n ...... t/4, A028387
A000384 ... w+x+y-n
A000217 ... w+x+y-2n
A211437 ... w*x*y-n ........ t/4, A007425
A211480 ... w+2x+3y-1
A211481 ... w+2x+3y-n
A211482 ... w*x+w*y+x*y-w*x*y
A211483 ... (n+w)^2-x-y
A182112 ... (n+w)^2-x-y-w
...
For the following sequences, S={1,...,n}, rather than
{-n,...,0,...n}. If f(w,x,y,n) is linear in w,x,y,n, then the sequence is a linear recurrence sequence.
A132188 ... w^2-x*y
A211506 ... w^2-x*y-n
A211507 ... w^2-x*y+n
A211508 ... w^2+x*y-n
A211509 ... w^2+x*y-2n
A211510 ... w^2-x*y+2n
A211511 ... w^2-2x*y ....... t/2
A211512 ... w^2-3x*y ....... t/2
A211513 ... 2w^2-x*y ....... t/2
A211514 ... 3w^2-x*y ....... t/2
A211515 ... w^3-x*y
A211516 ... w^2-x-y
A211517 ... w^3-(x+y)^2
A063468 ... w^2-x^2-y^2 .... t/2
A000217 ... w+x-y
A001399 ... w-2x-3y
A211519 ... w-2x+3y
A008810 ... w+2x-3y
A001399 ... w-2x-3y
A008642 ... w-2x-4y
A211520 ... w-2x+4y
A211521 ... w+2x-4y
A000115 ... w-2x-5y
A211522 ... w-2x+5y
A211523 ... w+2x-5y
A211524 ... w-3x-5y
A211533 ... w-3x+5y
A211523 ... w+3x-5y
A211535 ... w-4x-5y
A211536 ... w-4x+5y
A008812 ... w+4x-5y
A055998 ... w+x+y-2n
A074148 ... 2w+x+y-2n
A211538 ... 2w+2x+y-2n
A211539 ... 2w+2x-y-2n
A211540 ... 2w-3x-4y
A211541 ... 2w-3x+4y
A211542 ... 2w+3x-4y
A211543 ... 2w-3x-5y
A211544 ... 2w-3x+5y
A008812 ... 2w+3x-5y
A008805 ... w-2x-2y (repeated triangular numbers)
A001318 ... w-2x+2y
A000982 ... w+x-2y
A211534 ... w-3x-3y
A211546 ... w-3x+3y (triply repeated triangular numbers)
A211547 ... 2w-3x-3y (triply repeated squares)
A082667 ... 2w-3x+3y
A055998 ... w-x-y+2
A001399 ... w-2x-3y+1
A108579 ... w-2x-3y+n
...
Next, S={-n,...-1,1,...,n}, and the sequence counts the cases (w,x,y) satisfying the indicated inequality. If f(w,x,y,n) is linear in w,x,y,n, then the sequence is a linear recurrence sequence.
A211545 ... w+x+y>0; recurrence degree: 4
A211612 ... w+x+y>=0
A211613 ... w+x+y>1
A211614 ... w+x+y>2
A211615 ... |w+x+y|<=1
A211616 ... |w+x+y|<=2
A211617 ... 2w+x+y>0; recurrence degree: 5
A211618 ... 2w+x+y>1
A211619 ... 2w+x+y>2
A211620 ... |2w+x+y|<=1
A211621 ... w+2x+3y>0
A211622 ... w+2x+3y>1
A211623 ... |w+2x+3y|<=1
A211624 ... w+2x+2y>0; recurrence degree: 6
A211625 ... w+3x+3y>0; recurrence degree: 8
A211626 ... w+4x+4y>0; recurrence degree: 10
A211627 ... w+5x+5y>0; recurrence degree: 12
A211628 ... 3w+x+y>0; recurrence degree: 6
A211629 ... 4w+x+y>0; recurrence degree: 7
A211630 ... 5w+x+y>0; recurrence degree: 8
A211631 ... w^2>x^2+y^2; all terms divisible by 8
A211632 ... 2w^2>x^2+y^2; all terms divisible by 8
A211633 ... w^2>2x^2+2y^2; all terms divisible by 8
...
Next, S={1,...,n}, and the sequence counts the cases (w,x,y) satisfying the indicated relation.
A211634 ... w^2<=x^2+y^2
A211635 ... w^2A211790
A211636 ... w^2>=x^2+y^2
A211637 ... w^2>x^2+y^2
A211638 ... w^2+x^2+y^2
A211639 ... w^2+x^2+y^2<=n
A211640 ... w^2+x^2+y^2>n
A211641 ... w^2+x^2+y^2>=n
A211642 ... w^2+x^2+y^2<2n
A211643 ... w^2+x^2+y^2<=2n
A211644 ... w^2+x^2+y^2>2n
A211645 ... w^2+x^2+y^2>=2n
A211646 ... w^2+x^2+y^2<3n
A211647 ... w^2+x^2+y^2<=3n
A063691 ... w^2+x^2+y^2=n
A211649 ... w^2+x^2+y^2=2n
A211648 ... w^2+x^2+y^2=3n
A211650 ... w^3A211790
A211651 ... w^3>x^3+y^3; see Comments at A211790
A211652 ... w^4A211790
A211653 ... w^4>x^4+y^4; see Comments at A211790

Examples

			a(1) counts these 9 triples: (-1,-1,1), (-1, 1,-1), (0, -1, 0), (0, 0, -1), (0,0,0), (0,0,1), (0,1,0), (1,-1,1), (1,1,-1).
		

Crossrefs

Cf. A120486.

Programs

  • Mathematica
    t[n_] := t[n] = Flatten[Table[w^2 + x*y, {w, -n, n}, {x, -n, n}, {y, -n, n}]]
    c[n_] := Count[t[n], 0]
    t = Table[c[n], {n, 0, 70}] (* A211422 *)
    (t - 1)/8                   (* A120486 *)

A085787 Generalized heptagonal numbers: m*(5*m - 3)/2, m = 0, +-1, +-2 +-3, ...

Original entry on oeis.org

0, 1, 4, 7, 13, 18, 27, 34, 46, 55, 70, 81, 99, 112, 133, 148, 172, 189, 216, 235, 265, 286, 319, 342, 378, 403, 442, 469, 511, 540, 585, 616, 664, 697, 748, 783, 837, 874, 931, 970, 1030, 1071, 1134, 1177, 1243, 1288, 1357, 1404, 1476, 1525, 1600, 1651, 1729
Offset: 0

Author

Jon Perry, Jul 23 2003

Keywords

Comments

Zero together with the partial sums of A080512. - Omar E. Pol, Sep 10 2011
Second heptagonal numbers (A147875) and positive terms of A000566 interleaved. - Omar E. Pol, Aug 04 2012
These numbers appear in a theta function identity. See the Hardy-Wright reference, Theorem 355 on p. 284. See the g.f. of A113429. - Wolfdieter Lang, Oct 28 2016
Characteristic function is A133100. - Michael Somos, Jan 30 2017
40*a(n) + 9 is a square. - Bruno Berselli, Apr 18 2018
Numbers k such that the concatenation k225 is a square. - Bruno Berselli, Nov 07 2018
The sequence terms occur as exponents in the expansion of Sum_{n >= 0} q^(n*(n+1)) * Product_{k >= n+1} 1 - q^k = 1 - q - q^4 + q^7 + q^13 - q^18 - q^27 + + - - ... (see Hardy and Wright, Theorem 363, p. 290). - Peter Bala, Dec 15 2024

Examples

			From the first formula: a(5) = A000217(5) + A000217(2) = 15 + 3 = 18.
		

References

  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Fifth ed., Clarendon Press, Oxford, 2003, p. 284.

Crossrefs

Column 3 of A195152.
Sequences of generalized k-gonal numbers: A001318 (k=5), A000217 (k=6), this sequence (k=7), A001082 (k=8), A118277 (k=9), A074377 (k=10), A195160 (k=11), A195162 (k=12), A195313 (k=13), A195818 (k=14), A277082 (k=15), A274978 (k=16), A303305 (k=17), A274979 (k=18), A303813 (k=19), A218864 (k=20), A303298 (k=21), A303299 (k=22), A303303 (k=23), A303814 (k=24), A303304 (k=25), A316724 (k=26), A316725 (k=27), A303812 (k=28), A303815 (k=29), A316729 (k=30).

Programs

  • Haskell
    a085787 n = a085787_list !! n
    a085787_list = scanl (+) 0 a080512_list
    -- Reinhard Zumkeller, Apr 06 2015
  • Magma
    [5*n*(n+1)/8-1/16+(-1)^n*(2*n+1)/16: n in [0..60]]; // Vincenzo Librandi, Sep 11 2011
    
  • Mathematica
    Select[Table[(n*(n+1)/2-1)/5,{n,500}],IntegerQ] (* Vladimir Joseph Stephan Orlovsky, Feb 06 2012 *)
  • PARI
    t(n)=n*(n+1)/2
    for(i=0,40,print1(t(i)+t(floor(i/2)), ", "))
    
  • PARI
    {a(n) = (5*(-n\2)^2 - (-n\2)*3*(-1)^n) / 2}; /* Michael Somos, Oct 17 2006 */
    

Formula

a(n) = A000217(n) + A000217(floor(n/2)).
a(2*n-1) = A000566(n).
a(2*n) = A147875(n). - Bruno Berselli, Apr 18 2018
G.f.: x * (1 + 3*x + x^2) / ((1 - x) * (1 - x^2)^2). a(n) = a(-1-n) for all n in Z. - Michael Somos, Oct 17 2006
a(n) = 5*n*(n + 1)/8 - 1/16 + (-1)^n*(2*n + 1)/16. - R. J. Mathar, Jun 29 2009
a(n) = (A000217(n) + A001082(n))/2 = (A001318(n) + A118277(n))/2. - Omar E. Pol, Jan 11 2013
a(n) = A002378(n) - A001318(n). - Omar E. Pol, Oct 23 2013
Sum_{n>=1} 1/a(n) = 10/9 + (2*sqrt(1 - 2/sqrt(5))*Pi)/3. - Vaclav Kotesovec, Oct 05 2016
E.g.f.: (x*(9 + 5*x)*exp(x) - (1 - 2*x)*sinh(x))/8. - Franck Maminirina Ramaharo, Nov 07 2018
Sum_{n>=1} (-1)^(n+1)/a(n) = 5*log(5)/3 - 10/9 - 2*sqrt(5)*log(phi)/3, where phi is the golden ratio (A001622). - Amiram Eldar, Feb 28 2022

Extensions

New name from T. D. Noe, Apr 21 2006
Formula in sequence name added by Omar E. Pol, May 28 2012

A074377 Generalized 10-gonal numbers: m*(4*m - 3) for m = 0, +- 1, +- 2, +- 3, ...

Original entry on oeis.org

0, 1, 7, 10, 22, 27, 45, 52, 76, 85, 115, 126, 162, 175, 217, 232, 280, 297, 351, 370, 430, 451, 517, 540, 612, 637, 715, 742, 826, 855, 945, 976, 1072, 1105, 1207, 1242, 1350, 1387, 1501, 1540, 1660, 1701, 1827, 1870, 2002, 2047, 2185, 2232, 2376, 2425
Offset: 0

Author

W. Neville Holmes, Sep 04 2002

Keywords

Comments

Also called generalized decagonal numbers.
Odd triangular numbers decremented and halved.
It appears that this is zero together with the partial sums of A165998. - Omar E. Pol, Sep 10 2011 [this is correct, see the g.f., Joerg Arndt, Sep 29 2013]
Also, A033954 and positive members of A001107 interleaved. - Omar E. Pol, Aug 04 2012
Also, numbers m such that 16*m+9 is a square. After 1, therefore, there are no squares in this sequence. - Bruno Berselli, Jan 07 2016
Convolution of the sequences A047522 and A059841. - Ilya Gutkovskiy, Mar 16 2017
Numbers k such that the concatenation k5625 is a square. - Bruno Berselli, Nov 07 2018
Exponents in expansion of Product_{n >= 1} (1 + x^(8*n-7))*(1 + x^(8*n-1))*(1 - x^(8*n)) = 1 + x + x^7 + x^10 + x^22 + .... - Peter Bala, Dec 10 2020

Crossrefs

Cf. A001107 (10-gonal numbers).
Column 6 of A195152.
Sequences of generalized k-gonal numbers: A001318 (k=5), A000217 (k=6), A085787 (k=7), A001082 (k=8), A118277 (k=9), this sequence (k=10), A195160 (k=11), A195162 (k=12), A195313 (k=13), A195818 (k=14), A277082 (k=15), A274978 (k=16), A303305 (k=17), A274979 (k=18), A303813 (k=19), A218864 (k=20), A303298 (k=21), A303299 (k=22), A303303 (k=23), A303814 (k=24), A303304 (k=25), A316724 (k=26), A316725 (k=27), A303812 (k=28), A303815 (k=29), A316729 (k=30).
Cf. sequences of the form m*(m+k)/(k+1) listed in A274978. [Bruno Berselli, Jul 25 2016]

Programs

  • Magma
    [n^2+n-1/4+(-1)^n/4+n*(-1)^n/2: n in [0..50]]; // Vincenzo Librandi, Sep 29 2013
    
  • Mathematica
    CoefficientList[Series[x(1 +6x +x^2)/((1-x)(1-x^2)^2), {x, 0, 50}], x] (* Vincenzo Librandi, Sep 29 2013 *)
    LinearRecurrence[{1,2,-2,-1,1}, {0,1,7,10,22}, 50] (* G. C. Greubel, Nov 07 2018 *)
  • PARI
    a(n)=(2*n+3-4*(n%2))*(n-n\2)
    
  • PARI
    concat([0],Vec(x*(1 + 6*x + x^2)/((1 - x)*(1 - x^2)^2) +O(x^50))) \\ Indranil Ghosh, Mar 16 2017
    
  • Python
    def A074377(n): return (n+1>>1)*((n<<1)+(-1 if n&1 else 3)) # Chai Wah Wu, Mar 11 2025

Formula

(n(n+1)-2)/4 where n(n+1)/2 is odd.
G.f.: x*(1+6*x+x^2)/((1-x)*(1-x^2)^2). - Michael Somos, Mar 04 2003
a(2*k) = k*(4*k+3); a(2*k+1) = (2*k+1)^2+k. - Benoit Jubin, Feb 05 2009
a(n) = n^2+n-1/4+(-1)^n/4+n*(-1)^n/2. - R. J. Mathar, Oct 08 2011
Sum_{n>=1} 1/a(n) = (4 + 3*Pi)/9. - Vaclav Kotesovec, Oct 05 2016
E.g.f.: exp(x)*x^2 + (2*exp(x) - exp(-x)/2)*x - sinh(x)/2. - Ilya Gutkovskiy, Mar 16 2017
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2) - 4/9. - Amiram Eldar, Feb 28 2022
a(n) = (n+1)*(2*n-1)/2 if n is odd and a(n) = n*(2*n+3)/2 if n is even. - Chai Wah Wu, Mar 11 2025

Extensions

New name from T. D. Noe, Apr 21 2006
Formula in sequence name from Omar E. Pol, May 28 2012
Previous Showing 31-40 of 274 results. Next