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.

Showing 1-10 of 58 results. Next

A001651 Numbers not divisible by 3.

Original entry on oeis.org

1, 2, 4, 5, 7, 8, 10, 11, 13, 14, 16, 17, 19, 20, 22, 23, 25, 26, 28, 29, 31, 32, 34, 35, 37, 38, 40, 41, 43, 44, 46, 47, 49, 50, 52, 53, 55, 56, 58, 59, 61, 62, 64, 65, 67, 68, 70, 71, 73, 74, 76, 77, 79, 80, 82, 83, 85, 86, 88, 89, 91, 92, 94, 95, 97, 98, 100, 101, 103, 104
Offset: 1

Views

Author

Keywords

Comments

Inverse binomial transform of A084858. - Benoit Cloitre, Jun 12 2003
Earliest monotonic sequence starting with (1,2) and satisfying the condition: "a(n)+a(n-1) is not in the sequence." - Benoit Cloitre, Mar 25 2004. [The numbers of the form a(n)+a(n-1) form precisely the complement with respect to the positive integers. - David W. Wilson, Feb 18 2012]
a(1) = 1; a(n) is least number which is relatively prime to the sum of all the previous terms. - Amarnath Murthy, Jun 18 2001
For n > 3, numbers having 3 as an anti-divisor. - Alexandre Wajnberg, Oct 02 2005
Also numbers n such that (n+1)*(n+2)/6 = A000292(n)/n is an integer. - Ctibor O. Zizka, Oct 15 2010
Notice the property described by Gary Detlefs in A113801: more generally, these numbers are of the form (2*h*n + (h-4)*(-1)^n-h)/4 (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 3). - Bruno Berselli, Nov 17 2010
A001651 mod 9 gives A141425. - Paul Curtz, Dec 31 2010. (Correct for the modified offset 1. - M. F. Hasler, Apr 07 2015)
The set of natural numbers (1, 2, 3, ...), sequence A000027; represents the numbers of ordered compositions of n using terms in the signed set: (1, 2, -4, -5, 7, 8, -10, -11, 13, 14, ...). This follows from (1, 2, 3, ...) being the INVERT transform of A011655, signed and beginning: (1, 1, 0, -1, -1, 0, 1, 1, 0, ...). - Gary W. Adamson, Apr 28 2013
Union of A047239 and A047257. - Wesley Ivan Hurt, Dec 19 2013
Numbers whose sum of digits (and digital root) is != 0 (mod 3). - Joerg Arndt, Aug 29 2014
The number of partitions of 3*(n-1) into at most 2 parts. - Colin Barker, Apr 22 2015
a(n) is the number of partitions of 3*n into two distinct parts. - L. Edson Jeffery, Jan 14 2017
Conjectured (and like even easily proved) to be the graph bandwidth of the complete bipartite graph K_{n,n}. - Eric W. Weisstein, Apr 24 2017
Numbers k such that Fibonacci(k) mod 4 = 1 or 3. Equivalently, sequence lists the indices of the odd Fibonacci numbers (see A014437). - Bruno Berselli, Oct 17 2017
Minimum value of n_3 such that the "rectangular spiral pattern" is the optimal solution for Ripà's n_1 X n_2 x n_3 Dots Problem, for any n_1 = n_2. For example, if n_1 = n_2 = 5, n_3 = floor((3/2)*(n_1 - 1)) + 1 = a(5). - Marco Ripà, Jul 23 2018
For n >= 54, a(n) = sat(n, P_n), the minimum number of edges in a P_n-saturated graph on n vertices, where P_n is the n-vertex path (see Dudek, Katona, and Wojda, 2003; Frick and Singleton, 2005). - Danny Rorabaugh, Nov 07 2017
From Roger Ford, May 09 2021: (Start)
a(n) is the smallest sum of arch lengths for the top arches of a semi-meander with n arches. An arch length is the number of arches covered + 1.
/\ The top arch has a length of 3. /\ The top arch has a length of 3.
/ \ Both bottom arches have a //\\ The middle arch has a length of 2.
//\/\\ length of 1. ///\\\ The bottom arch has a length of 1.
Example: a(6) = 8 /\ /\
//\\ /\ //\\ /\ 2 + 1 + 1 + 2 + 1 + 1 = 8. (End)
This is the lexicographically earliest increasing sequence of positive integers such that no polynomial of degree d can be fitted to d+2 consecutive terms (equivalently, such that no iterated difference is zero). - Pontus von Brömssen, Dec 26 2021

Examples

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

References

  • 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).

Crossrefs

Programs

  • GAP
    Filtered([0..110],n->n mod 3<>0); # Muniru A Asiru, Jul 24 2018
    
  • Haskell
    a001651 = (`div` 2) . (subtract 1) . (* 3)
    a001651_list = filter ((/= 0) . (`mod` 3)) [1..]
    -- Reinhard Zumkeller, Jul 07 2012, Aug 23 2011
    
  • Magma
    [3*(2*n-1)/4-(-1)^n/4: n in [1..80]]; // Vincenzo Librandi, Jun 07 2011
    
  • Maple
    A001651 := n -> 3*floor(n/2) - (-1)^n; # Corrected by M. F. Hasler, Apr 07 2015
    A001651:=(1+z+z**2)/(z+1)/(z-1)**2; # Simon Plouffe in his 1992 dissertation
    a[1]:=1:a[2]:=2:for n from 3 to 100 do a[n]:=a[n-2]+3 od: seq(a[n], n=1..69); # Zerinvary Lajos, Mar 16 2008, offset corrected by M. F. Hasler, Apr 07 2015
  • Mathematica
    Select[Table[n,{n,200}],Mod[#,3]!=0&] (* Vladimir Joseph Stephan Orlovsky, Feb 18 2011 *)
    Drop[Range[200 + 1], {1, -1, 3}] - 1 (* József Konczer, May 24 2016 *)
    Floor[(3 Range[70] - 1)/2] (* Eric W. Weisstein, Apr 24 2017 *)
    CoefficientList[Series[(x^2 + x + 1)/((x - 1)^2 (x + 1)), {x, 0, 70}],
      x] (* or *)
    LinearRecurrence[{1, 1, -1}, {1, 2, 4}, 70] (* Robert G. Wilson v, Jul 25 2018 *)
  • PARI
    {a(n) = n + (n-1)\2}; /* Michael Somos, Jan 15 2011 */
    
  • PARI
    x='x+O('x^100); Vec(x*(1+x+x^2)/((1-x)*(1-x^2))) \\ Altug Alkan, Oct 22 2015
    
  • Python
    print([k for k in range(1, 105) if k%3]) # Michael S. Branicky, Sep 06 2021
    
  • Python
    def A001651(n): return (n<<1)-(n>>1)-1 # Chai Wah Wu, Mar 05 2024

Formula

a(n) = 3 + a(n-2) for n > 2.
a(n) = a(n-1) + a(n-2) - a(n-3) for n > 3.
a(2*n+1) = 3*n+1, a(2*n) = 3*n-1.
G.f.: x * (1 + x + x^2) / ((1 - x) * (1 - x^2)). - Michael Somos, Jun 08 2000
a(n) = (4-n)*a(n-1) + 2*a(n-2) + (n-3)*a(n-3) (from the Carlitz et al. article).
a(n) = floor((3*n-1)/2). [Corrected by Gary Detlefs]
a(1) = 1, a(n) = 2*a(n-1) - 3*floor(a(n-1)/3). - Benoit Cloitre, Aug 17 2002
a(n+1) = 1 + n - n mod 2 + (n + n mod 2)/2. - Reinhard Zumkeller, Dec 17 2002
a(1) = 1, a(n+1) = a(n) + (a(n) mod 3). - Reinhard Zumkeller, Mar 23 2003
a(1) = 1, a(n) = 3*(n-1) - a(n-1). - Benoit Cloitre, Apr 12 2003
a(n) = 3*(2*n-1)/4 - (-1)^n/4. - Benoit Cloitre, Jun 12 2003
Nearest integer to (Sum_{k>=n} 1/k^3)/(Sum_{k>=n} 1/k^4). - Benoit Cloitre, Jun 12 2003
Partial sums of A040001. a(n) = A032766(n-1)+1. - Paul Barry, Sep 02 2003
a(n) = T(n, 1) = T(n, n-1), where T is the array in A026386. - Emeric Deutsch, Feb 18 2004
a(n) = sqrt(3*A001082(n)+1). - Zak Seidov, Dec 12 2007
a(n) = A077043(n) - A077043(n-1). - Reinhard Zumkeller, Dec 28 2007
a(n) = A001477(n-1) + A008619(n-1). - Yosu Yurramendi, Aug 10 2008
Euler transform of length 3 sequence [2, 1, -1]. - Michael Somos, Sep 06 2008
A011655(a(n)) = 1. - Reinhard Zumkeller, Nov 30 2009
a(n) = n - 1 + ceiling(n/2). - Michael Somos, Jan 15 2011
a(n) = 3*A000217(n)+1 - 2*Sum_{i=1..n-1} a(i), for n>1. - Bruno Berselli, Nov 17 2010
a(n) = 3*floor(n/2) + (-1)^(n+1). - Gary Detlefs, Dec 29 2011
A215879(a(n)) > 0. - Reinhard Zumkeller, Dec 28 2012 [More precisely, A215879 is the characteristic function of A001651. - M. F. Hasler, Apr 07 2015]
a(n) = 2n - 1 - floor(n/2). - Wesley Ivan Hurt, Oct 25 2013
a(n) = (3n - 2 + (n mod 2)) / 2. - Wesley Ivan Hurt, Mar 31 2014
a(n) = A000217(n) - A000982(n-1). - Bui Quang Tuan, Mar 28 2015
1/1^3 - 1/2^3 + 1/4^3 - 1/5^3 + 1/7^3 - 1/8^3 + ... = 4 Pi^3/(3 sqrt(3)). - M. F. Hasler, Mar 29 2015
E.g.f.: (4 + sinh(x) - cosh(x) + 3*(2*x - 1)*exp(x))/4. - Ilya Gutkovskiy, May 24 2016
a(n) = a(n+k-1) + a(n-k) - a(n-1) for n > k >= 0. - Bob Selcoe, Feb 03 2017
a(n) = -a(1-n) for all n in Z. - Michael Somos, Jul 31 2018
a(n) = n + A004526(n-1). - David James Sycamore, Sep 06 2021
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(3*sqrt(3)) (A073010). - Amiram Eldar, Dec 04 2021
From Amiram Eldar, Nov 22 2024: (Start)
Product_{n>=1} (1 - (-1)^n/a(n)) = 1.
Product_{n>=2} (1 + (-1)^n/a(n)) = 2*Pi/(3*sqrt(3)) (A248897). (End)

Extensions

This is a list, so the offset should be 1. I corrected this and adjusted some of the comments and formulas. Other lines probably also need to be adjusted. - N. J. A. Sloane, Jan 01 2011
Offset of pre-2011 formulas verified or corrected by M. F. Hasler, Apr 07-18 2015 and by Danny Rorabaugh, Oct 23 2015

A072649 n occurs Fibonacci(n) times (cf. A000045).

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Jun 02 2002

Keywords

Comments

Number of digits in Zeckendorf-binary representation of n. E.g., the Zeckendorf representation of 12 is 8+3+1, which in binary notation is 10101, which consists of 5 digits. - Clark Kimberling, Jun 05 2004
First position where value n occurs is A000045(n+1), i.e., a(A000045(n)) = n-1, for n >= 2 and a(A000045(n)-1) = n-2, for n >= 3.
This is the number of distinct Fibonacci numbers greater than 0 which are less than or equal to n. - Robert G. Wilson v, Dec 10 2006
The smallest nondecreasing sequence a(n) such that a(Fibonacci(n-1)) = n. - Tanya Khovanova, Jun 20 2007

Examples

			1, 1, then F(2) 2's, then F(3) 3's, then F(4) 4's, ..., then F(k) k's, ...
		

Crossrefs

Cf. A001622 (golden ratio phi), A073010.
Used to construct A003714. Cf. also A002024, A072643, A072648, A072650.
Cf. A131234.
Partial sums: A256966, A256967.

Programs

  • Haskell
    a072649 n = a072649_list !! (n-1)
    a072649_list = f 1 where
       f n = (replicate (fromInteger $ a000045 n) n) ++ f (n+1)
    -- Reinhard Zumkeller, Jul 04 2011
    
  • Maple
    A072649 := proc(n)
        local j;
        for j from ilog[(1+sqrt(5))/2](n) while combinat[fibonacci](j+1)<=n do
        end do;
        j-1
    end proc:
    seq(A072649(n), n=1..120);  # Alois P. Heinz, Mar 18 2013
  • Mathematica
    Table[Table[n, {Fibonacci[n]}], {n, 10}] // Flatten (* Robert G. Wilson v, Jan 14 2007 *)
    a[n_] := Module[{j}, For[j = Floor@Log[GoldenRatio, n], Fibonacci[j+1] <= n, j++]; j-1];
    Table[a[n], {n, 1, 120}] (* Jean-François Alcover, Nov 17 2022, after Alois P. Heinz *)
  • PARI
    a(n) = -1+floor(log(((n+0.2)*sqrt(5)))/log((1+sqrt(5))/2))
    
  • PARI
    a(n)=local(m); if(n<1,0,m=0; until(fibonacci(m)>n,m++); m-2)
    
  • Python
    from sympy import fibonacci
    def a(n):
        if n<1: return 0
        m=0
        while fibonacci(m)<=n: m+=1
        return m-2
    print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 09 2017
    
  • Python
    def A072649(n):
        a, b, c = 0, 1, -2
        while a <= n:
            a, b = b, a+b
            c += 1
        return c # Chai Wah Wu, Nov 04 2024
    (MIT/GNU Scheme) (define (A072649 n) (let ((b (A072648 n))) (+ -1 b (floor->exact (/ n (A000045 (1+ b))))))) ;; (The implementation below is better)
    
  • Scheme
    (define (A072649 n) (if (<= n 3) n (let loop ((k 5)) (if (> (A000045 k) n) (- k 2) (loop (+ 1 k)))))) ;; (Use this with the memoized implementation of A000045 given under that entry. No floating point arithmetic is involved). - Antti Karttunen, Oct 06 2017

Formula

G.f.: (Sum_{n>1} x^Fibonacci(n))/(1-x). - Michael Somos, Apr 25 2003
From Hieronymus Fischer, May 02 2007: (Start)
a(n) = floor(log_phi((sqrt(5)*n + sqrt(5*n^2+4))/2)) - 1, where phi is A001622.
a(n) = floor(arcsinh(sqrt(5)*n/2)/log(phi)) - 1.
a(n) = A108852(n) - 2. (End)
a(n) = -1 + floor( log_phi( (n+0.2)*sqrt(5) ) ), where log_phi(x) is the logarithm to the base (1+sqrt(5))/2. - Ralf Stephan, May 14 2007
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(3*sqrt(3)) (A073010). - Amiram Eldar, Feb 18 2024

Extensions

Typo fixed by Charles R Greathouse IV, Oct 28 2009

A002324 Number of divisors of n == 1 (mod 3) minus number of divisors of n == 2 (mod 3).

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Coefficients of Dedekind zeta function for the quadratic number field of discriminant -3. See Formula section for the general expression. - N. J. A. Sloane, Mar 22 2022
Coefficients in expansion of Dirichlet series Product_p (1 - (Kronecker(m,p) + 1)*p^(-s) + Kronecker(m,p) * p^(-2s))^(-1) for m = -3.
(Number of points of norm n in hexagonal lattice) / 6, n>0.
The hexagonal lattice is the familiar 2-dimensional lattice (A_2) in which each point has 6 neighbors. This is sometimes called the triangular lattice.
The first occurrence of a(n) = 1, 2, 3, 4,... is at n= 1, 7, 49, 91, 2401, 637, ... as tabulated in A343771. - R. J. Mathar, Sep 21 2024

Examples

			G.f. = x + x^3 + x^4 + 2*x^7 + x^9 + x^12 + 2*x^13 + x^16 + 2*x^19 + 2*x^21 + ...
		

References

  • J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 112, first display.
  • J. W. L. Glaisher, Table of the excess of the number of (3k+1)-divisors of a number over the number of (3k+2)-divisors, Messenger Math., 31 (1901), 64-72.
  • D. H. Lehmer, Guide to Tables in the Theory of Numbers. Bulletin No. 105, National Research Council, Washington, DC, 1941, pp. 7-10.
  • 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).

Crossrefs

Dedekind zeta functions for imaginary quadratic number fields of discriminants -3, -4, -7, -8, -11, -15, -19, -20 are A002324, A002654, A035182, A002325, A035179, A035175, A035171, A035170, respectively.
Dedekind zeta functions for real quadratic number fields of discriminants 5, 8, 12, 13, 17, 21, 24, 28, 29, 33, 37, 40 are A035187, A035185, A035194, A035195, A035199, A035203, A035188, A035210, A035211, A035215, A035219, A035192, respectively.

Programs

  • Haskell
    a002324 n = a001817 n - a001822 n  -- Reinhard Zumkeller, Nov 26 2011
    
  • Maple
    A002324 := proc(n)
        local a,pe,p,e;
        a :=1 ;
        for pe in ifactors(n)[2] do
            p := op(1,pe) ;
            e := op(2,pe) ;
            if p = 3 then
                ;
            elif modp(p,3) = 1 then
                a := a*(e+1) ;
            else
                a := a*(1+(-1)^e)/2 ;
            end if;
        end do:
        a ;
    end proc:
    seq(A002324(n),n=1..100) ; # R. J. Mathar, Sep 21 2024
  • Mathematica
    dn12[n_]:=Module[{dn=Divisors[n]},Count[dn,?(Mod[#,3]==1&)]-Count[ dn,?(Mod[#,3]==2&)]]; dn12/@Range[120]  (* Harvey P. Dale, Apr 26 2011 *)
    a[ n_] := If[ n < 1, 0, DivisorSum[ n, KroneckerSymbol[ -3, #] &]]; (* Michael Somos, Aug 24 2014 *)
    Table[DirichletConvolve[DirichletCharacter[3,2,m],1,m,n],{n,1,30}] (* Steven Foster Clark, May 29 2019 *)
    f[3, p_] := 1; f[p_, e_] := If[Mod[p, 3] == 1, e+1, (1+(-1)^e)/2]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 17 2020 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( sum(k=1, n, x^k / (1 + x^k + x^(2*k)), x * O(x^n)), n))}; \\ Michael Somos
    
  • PARI
    {a(n) = if( n<1, 0, sumdiv(n, d, (d%3==1) - (d%3==2)))};
    
  • PARI
    {a(n) = local(A, p, e); if( n<1, 0, A = factor(n); prod(k=1, matsize(A)[1], if( p=A[k,1], e=A[k,2]; if( p==3, 1, if( p%3==1, e+1, !(e%2))))))}; \\ Michael Somos, May 20 2005
    
  • PARI
    {a(n) = if( n<1, 0, qfrep([2,1; 1,2], n, 1)[n] / 3)}; \\ Michael Somos, Jun 05 2005
    
  • PARI
    {a(n) = if( n<1, 0, direuler(p=2, n, 1 / (1 - X) / (1 - kronecker(-3, p)*X))[n])}; \\ Michael Somos, Jun 05 2005
    
  • PARI
    my(B=bnfinit(x^2+x+1)); vector(100,n,#bnfisintnorm(B,n)) \\ Joerg Arndt, Jun 01 2024
    
  • Python
    from math import prod
    from sympy import factorint
    def A002324(n): return prod(e+1 if p%3==1 else int(not e&1) for p, e in factorint(n).items() if p != 3) # Chai Wah Wu, Nov 17 2022

Formula

From N. J. A. Sloane, Mar 22 2022 (Start):
The Dedekind zeta function DZ_K(s) for a quadratic field K of discriminant D is as follows.
Here m is defined by K = Q(sqrt(m)) (so m=D/4 if D is a multiple of 4, otherwise m=D).
DZ_K(s) is the product of three terms:
(a) Product_{odd primes p | D} 1/(1-1/p^s)
(b) Product_{odd primes p such that (D|p) = -1} 1/(1-1/p^(2s))
(c) Product_{odd primes p such that (D|p) = 1} 1/(1-1/p^s)^2
and if m is
0,1,2,3,4,5,6,7 mod 8, the prime 2 is to be included in term
-,c,a,a,-,b,a,a, respectively.
For Maple (and PARI) implementations, see link. (End)
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w) = u^2 - 3*v^2 + 4*w^2 - 2*u*w + w - v. - Michael Somos, Jul 20 2004
Has a nice Dirichlet series expansion, see PARI line.
G.f.: Sum_{k>0} x^k/(1+x^k+x^(2*k)). - Vladeta Jovovic, Dec 16 2002
a(3*n + 2) = 0, a(3*n) = a(n), a(3*n + 1) = A033687(n). - Michael Somos, Apr 04 2003
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^3), A(x^6)) where f(u1, u2, u3, u6) = (u1 - u3)*(u3 - u6) - (u2 - u6)^2. - Michael Somos, May 20 2005
Multiplicative with a(3^e) = 1, a(p^e) = e+1 if p == 1 (mod 3), a(p^e) = (1+(-1)^e)/2 if p == 2 (mod 3). - Michael Somos, May 20 2005
G.f.: Sum_{k>0} x^(3*k - 2) / (1 - x^(3*k - 2)) - x^(3*k - 1) / (1 - x^(3*k - 1)). - Michael Somos, Nov 02 2005
G.f.: Sum_{n >= 1} q^(n^2)(1-q)(1-q^2)...(1-q^(n-1))/((1-q^(n+1))(1-q^(n+2))...(1-q^(2n))). - Jeremy Lovejoy, Jun 12 2009
a(n) = A001817(n) - A001822(n). - R. J. Mathar, Mar 31 2011
A004016(n) = 6*a(n) unless n=0.
Dirichlet g.f.: zeta(s)*L(chi_2(3),s), with chi_2(3) the nontrivial Dirichlet character modulo 3 (A102283). - Ralf Stephan, Mar 27 2015
From Andrey Zabolotskiy, May 07 2018: (Start)
a(n) = Sum_{ m: m^2|n } A000086(n/m^2).
a(A003136(m)) > 0, a(A034020(m)) = 0 for all m. (End)
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/(3*sqrt(3)) = 0.604599... (A073010). - Amiram Eldar, Oct 11 2022

Extensions

More terms from David Radcliffe
Somos D.g.f. replaced with correct version by Ralf Stephan, Mar 27 2015

A000407 a(n) = (2*n+1)! / n!.

Original entry on oeis.org

1, 6, 60, 840, 15120, 332640, 8648640, 259459200, 8821612800, 335221286400, 14079294028800, 647647525324800, 32382376266240000, 1748648318376960000, 101421602465863680000, 6288139352883548160000, 415017197290314178560000
Offset: 0

Views

Author

Keywords

Comments

The e.g.f. of 1/a(n) = n!/(2*n+1)! is (exp(sqrt(x)) - exp(-sqrt(x)))/(2*sqrt(x)). - Wolfdieter Lang, Jan 09 2012
Product of the larger parts of the partitions of 2n+2 into exactly two parts. - Wesley Ivan Hurt, Jun 15 2013
For n > 0, a(n-1) = (2n-1)!/(n-1)!, the number of ways n people can line up in n labeled queues. The derivation is straightforward. Person 1 has (2n-1) choices - be first in line in one of the queues or get behind one of the other people. Person 2 has (2n-2) choices - choose one of the n queues or get behind one of the remaining n-2 people. Continuing in this fashion, we finally find that person n has to choose one of the n queues. - Dennis P. Walsh, Mar 24 2016
For n > 0, a(n-1) is the number of functions f:[n]->[2n] that are acyclic and injective. Note that f is acyclic if, for all x in [n], x is not a member of the set {f(x),f(f(x)), f(f(f(x))), ...}. - Dennis P. Walsh, Mar 25 2016
a(n) is the number of labeled maximal outerplanar graphs with n-3 vertices. - Allan Bickle, Feb 19 2024

Examples

			G.f. = 1 + 6*x + 60*x^2 + 840*x^3 + 15120*x^4 + 332640*x^5 + 8648640*x^6 + ...
For n=1 the a(1)=6 ways for 2 people to line up in 2 queues are as follows: Q1<P1,P2> Q2<>, Q1<P2,P1> Q2<>, Q1<P1> Q2<P2>, Q1<P2> Q2<P1>, Q1<> Q2<P1,P2>, Q1<> Q2<P2,P1>. - _Dennis P. Walsh_, Mar 24 2016
For the unique maximal outerplanar graph with 4 vertices, there are C(4,2)=6 ways to label the two degree 3 vertices, and the other two labels are forced.  Thus a(1) = 6.
		

References

  • L. W. Beineke and R. E. Pippert, Enumerating labeled k-dimensional trees and ball dissections, pp. 12-26 of Proceedings of Second Chapel Hill Conference on Combinatorial Mathematics and Its Applications, University of North Carolina, Chapel Hill, 1970. Reprinted with a slightly different title in Math. Annalen, 191 (1971), 87-98.
  • L. B. W. Jolley, Summation of Series, Dover, 1961.
  • Loren C. Larson, The number of essentially different nonattacking rook arrangements, J. Recreat. Math., 7 (No. 3, 1974), circa pages 180-181.
  • 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).

Crossrefs

A100622 is the "Number of topologically distinct solutions to the clone ordering problem for n clones" without the restriction that they be in a single contig (see [Newberg] for definition of contig).
Column m=0 of A292219.

Programs

  • Magma
    [Factorial(2*n+1) / Factorial(n): n in [0..20]]; // Vincenzo Librandi, Jun 16 2015
  • Maple
    For Maple program see A000903.
    a := n -> pochhammer(n+1,n+1); (for n>=0) # Peter Luschny, Feb 14 2009
  • Mathematica
    Table[(2n + 1)!/n!, {n, 0, 30}] (* Stefan Steinerberger, Apr 08 2006 *)
    a[ n_] := If[ n < 0, 1/2, 1] Pochhammer[ n + 1, n + 1]; (* Michael Somos, Jan 03 2015 *)
    a[ n_] := Which[ n < -1, -(-1)^n / (4 a[-n - 2]), n == -1, 1/2, True, (2 n + 1)! / n!]; (* Michael Somos, Jan 03 2015 *)
  • Maxima
    A000407(n):=(2*n+1)!/n!$
    makelist(A000407(n),n,0,30); /* Martin Ettl, Nov 05 2012 */
    
  • PARI
    a(n)=(2*n+1)!/n! \\ Charles R Greathouse IV, Jan 12 2012
    
  • PARI
    {a(n) = if( n<-1, -(-1)^n / (4 * a(-n-2)), n==-1, 1/2, (2*n + 1)! / n!)}; /* Michael Somos, Jan 03 2015 */
    

Formula

E.g.f.: (1 - 4*x)^(-3/2). - Michael Somos, Jan 03 2015
E.g.f.: Sum_{k>=0} a(k+2) * x^k / k! = (1 - 2*x - sqrt(1 - 4*x)) / 4.
E.g.f. for a(n-1), n >= 0, with a(-1) := 0 is (-1+1/(1-4*x)^(1/2))/2. 2*a(n) = (4*n+2)(!^4) := Product_{j=0..n} (4*j + 2), (one half of 4-factorial numbers). - Wolfdieter Lang
a(n) = C(n+1)*(n+2)!/2 for all n>=0. - Paul Barry, Feb 16 2005
For n>1, a(n) = (1/2)*A001813(n+1). - Zerinvary Lajos, Jun 06 2007
For asymptotics see the Robinson paper.
Sum_{n >=0} n!/a(n) = 2*Pi/3^(3/2) = 1.2091995761... = A248897 [Jolley eq 261]
G.f.: 1 / (1 - 6*x / (1 - 4*x / (1 - 10*x / (1 - 8*x / (1 - 14*x / ... ))))). - Michael Somos, May 12 2012
G.f.: 1/Q(0), where Q(k) = 1 + 2*(2*k-1)*x - 4*x*(k+1)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 03 2013
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - 2*x/(2*x + 1/(2*k+3)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 02 2013
a(n) = -(-1)^n / (4 * a(-2-n)) = a(n-1) * (4*n+2) for all n in Z. - Michael Somos, Jan 03 2015
a(n) = A087299(2*n + 1). - Michael Somos, Jan 03 2015
From Peter Bala, Feb 16 2015: (Start)
Recurrence equation: a(n) = 4*a(n-1) + 4*(2*n - 1)^2*a(n-2) with a(0) = 1 and a(1) = 6.
The integer sequence b(n) := a(n)*Sum_{k = 0..n} (-1)^k/(2*k + 1), beginning [1, 4, 52, 608, 12624, ...], satisfies the same second-order recurrence equation. This leads to Brouncker's generalized continued fraction expansion Sum_{k >= 0} (-1)^k/(2*k + 1) = Pi/4 = 1/(1 + 1^2/(2 + 3^2/(2 + 5^2/(2 + ... )))). Note b(n) = 2^n*A024199(n+1).
Recurrence equation: a(n) = (5*n + 2)*a(n-1) - 2*n*(2*n - 1)^2*a(n-2) with a(0) = 1 and a(1) = 6.
The integer sequence c(n) := a(n)*Sum_{k = 0..n} k!^2/(2*k + 1)!, beginning [1, 7, 72, 1014, 18276, ... ], satisfies the same second-order recurrence equation. This leads to the generalized continued fraction expansion Sum_{k >= 0} k!^2/(2*k + 1)! = 2*Pi/sqrt(27) = 2*A073010 = 1/(1 - 1/(7 - 12/(12 - 30/(17 - ... - 2*n*(2*n - 1)/((5*n + 2) - ... ))))). (End)
a(n) = Product_{k=n+1..(2*n+1)} k. - Carlos Eduardo Olivieri, Jun 03 2015
From Ilya Gutkovskiy, Jan 17 2017: (Start)
a(n) ~ 2^(2*n+3/2)*n^(n+1)/exp(n).
Sum_{n>=0} 1/a(n) = exp(1/4)*sqrt(Pi)*erf(1/2) = 1.184593072938653151..., where erf() is the error function. (End)
Sum_{n>=0} (-1)^n/a(n) = exp(-1/4)*sqrt(Pi)*erfi(1/2), where erfi() is the imaginary error function. - Amiram Eldar, Jan 18 2021
It follows from the comments above that we have a(n) = a(n-1)*(4*n+2), with a(1) = 6, a(0) = 1.
a(n) = A081125(2*n+1). - R. J. Mathar, Jun 07 2025

A102283 Period 3: repeat [0, 1, -1].

Original entry on oeis.org

0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1
Offset: 0

Views

Author

N. J. A. Sloane, Nov 02 2008

Keywords

Comments

The sequence is the non-principal Dirichlet character of the reduced residue system mod 3. (The other is A011655.) Associated Dirichlet L-functions are L(1, chi) = Sum_{n >= 1} a(n)/n = A073010, L(2, chi)= Sum_{n >= 1} a(n)/n^2 = A086724, or L(3, chi)= Sum_{n >= 1} a(n)/n^3 = A129404. [Jolley eq 310] - R. J. Mathar, Jul 15 2010
a(n) = 2*D(n) - L(n), where L(n) denotes the n-th Lucas number and D(n) denotes the so-called n-th quadrapell number -- defined and discussed by Dursun Tasci in his paper (see References below). We have D(n) = D(n-2) + 2*D(n-3) + D(n-4), D(0) = D(1) = D(2) = 1, D(3) = 2. G.f. D(x) = (1+x-x^3)/((1-x-x^2)(1+x+x^2)). - Roman Witula, Jul 31 2012
This is a strong elliptic divisibility sequence t_n as given in [Kimberling, p. 16] where x = -1, y = 0, z = -1. - Michael Somos, Nov 27 2019

Examples

			G.f. = x - x^2 + x^4 - x^5 + x^7 - x^8 + x^10 - x^11 + ... - _Michael Somos_, Nov 27 2019
		

References

  • M. N. Huxley, Area, Lattice Points and Exponential Sums, Oxford, 1996; p. 236.
  • L. B. W. Jolley, Summation of Series, Dover Publications (1961).
  • Paulo Ribenboim, My Numbers, My Friends: Popular Lectures on Number Theory, Springer-Verlag, NY, 2000, p. 6.

Crossrefs

Cf. A011655, A049347, A073010, A086724, A129404, A002324 (Mobius transform).

Programs

Formula

a(n) = A049347(n-1).
a(n) = -a(n-1) - a(n-2); a(0) = 0, a(1) = 1. G.f.: x/(1+x+x^2). - Philippe Deléham, Nov 03 2008
a(n) = -2*sin(4*Pi*n/3)/sqrt(3) = 2*sin(8*Pi*n/3)/sqrt(3). - Jaume Oliver Lafont, Dec 05 2008
a(n) = 2*sin(2*Pi*n/3)/sqrt(3). - Roman Witula, Jul 31 2012
a(n) = Legendre(n, 3), the Legendre symbol for p = 3. - Alonso del Arte, Feb 06 2013
a(n) = (-3/n), where (k/n) is the Kronecker symbol. See the Eric Weisstein and Wikipedia links. - Wolfdieter Lang, May 29 2013
Dirichlet g.f.: L(chi_2(3),s), with chi_2(3) the nontrivial Dirichlet character modulo 3. - Ralf Stephan, Mar 27 2015
a(n) = a(n-3) for n > 2. - Wesley Ivan Hurt, Jul 02 2016
E.g.f.: 2*sin(sqrt(3)*x/2)*exp(-x/2)/sqrt(3). - Ilya Gutkovskiy, Jul 02 2016
a(n) = H(2*n, 1, 1/2) for n > 0 where H(n, a, b) = hypergeom([a - n/2, b - n/2], [1 - n], 4). - Peter Luschny, Sep 03 2019
Euler transform of length 3 sequence [-1, 0, 1]. - Michael Somos, Nov 27 2019
a(n) = n - 3*floor((n+1)/3). - Wolfdieter Lang, Oct 07 2021

A086463 Decimal expansion of Pi^2/18.

Original entry on oeis.org

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

Views

Author

Eric W. Weisstein, Jul 21 2003

Keywords

Comments

The sequence of repeating coefficients [1,-1,-2,-1,1,2] in the sum in the formula section, is equal to the 6th column in A191898. - Mats Granvik, Mar 19 2012

Examples

			0.548311355616075478824138388882008396406316633735...
		

References

  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Section 1.4.1, p. 20.
  • A. Holroyd, Sharp Metastability Threshold for Two-Dimensional Bootstrap Percolation, Prob. Th. and Related Fields 125, 195-224, 2003.

Crossrefs

Programs

Formula

Sum[1/n^2/Binomial[2n,n], {n,Infinity}].
Pi^2/18 = A013661/3 = Sum[1/(i+0)^2 - 1/(i+1)^2 - 2/(i+2)^2 - 1/(i+3)^2 + 1/(i+4)^2 + 2/(i+5)^2, {i =1, 7, 13, 19, 25,.. infinity, stride of 6}]. - Mats Granvik, Mar 19 2012
Equals Sum_{k>=1} (H(k) - 2*H(2k))/((-3^k)*k). See Liu. - Michel Marcus, Feb 11 2020
Equals Sum_{k>=1} A007814(k)/k^2. - Amiram Eldar, Jul 13 2020
Equals (2/9) * Sum_{k>=0} (-1)^k*(7*k+5)*k!^3/((2*k+1)*(3*k+2)!) [Gosper 1974] - R. J. Mathar, Feb 07 2024
Continued fraction expansion: 1/(2 - 2/(13 - 48/(34 - 270/(65 - ... - 2*(2*n - 1)*n^3/(5*n^2 + 6*n + 2 - ... ))))). See A130549. - Peter Bala, Feb 16 2024

A057353 a(n) = floor(3n/4).

Original entry on oeis.org

0, 0, 1, 2, 3, 3, 4, 5, 6, 6, 7, 8, 9, 9, 10, 11, 12, 12, 13, 14, 15, 15, 16, 17, 18, 18, 19, 20, 21, 21, 22, 23, 24, 24, 25, 26, 27, 27, 28, 29, 30, 30, 31, 32, 33, 33, 34, 35, 36, 36, 37, 38, 39, 39, 40, 41, 42, 42, 43, 44, 45, 45, 46, 47, 48, 48, 49, 50, 51, 51, 52, 53, 54
Offset: 0

Views

Author

Keywords

Comments

The cyclic pattern (and numerator of the gf) is computed using Euclid's algorithm for GCD.
For n >= 2, a(n) is the number of different integers that can be written as floor(k^2/n) for k = 1, 2, 3, ..., n-1. Generalization of the 1st problem proposed during the 15th Balkan Mathematical Olympiad in 1998 where the question was asked for n = 1998 with a(1998) = 1498. - Bernard Schott, Apr 22 2022
For n > 1, a(n) is also the Hadwiger number of the (n+1)-cycle complement graph (up to at least n = 16). - Eric W. Weisstein, Mar 10 2025

References

  • N. Dershowitz and E. M. Reingold, Calendrical Calculations, Cambridge University Press, 1997.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics, Addison-Wesley, NY, 1994.

Crossrefs

Programs

Formula

G.f.: (1+x+x^2)*x^2/((1-x)*(1-x^4)). - Bruce Corrigan (scentman(AT)myfamily.com), Jul 03 2002
For all m>=0 a(4m)=0 mod 3; a(4m+1)=0 mod 3; a(4m+2)= 1 mod 3; a(4m+3) = 2 mod 3
a(n) = A002378(n) - A173562(n). - Reinhard Zumkeller, Feb 21 2010
a(n+1) = A140201(n) - A002265(n+1). - Reinhard Zumkeller, Jan 26 2011
a(n) = n-1 - A002265(n-1) = ( A007310(n) + A057077(n+1) )/4 for n>0. a(n) = a(n-1)+a(n-4)-a(n-5) for n>4. - Bruno Berselli, Jan 28 2011
a(n) = 1/8*(6*n + 2*cos((Pi*n)/2) + cos(Pi*n) - 2*sin((Pi*n)/2) - 3). - Ilya Gutkovskiy, Sep 18 2015
a(4n) = a(4n+1). - Altug Alkan, Sep 26 2015
Sum_{n>=2} (-1)^n/a(n) = Pi/(3*sqrt(3)) (A073010). - Amiram Eldar, Sep 29 2022

A248897 Decimal expansion of Sum_{i >= 0} (i!)^2/(2*i+1)!.

Original entry on oeis.org

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

Views

Author

Bruno Berselli, Mar 06 2015

Keywords

Comments

Value of the Borwein-Borwein function I_3(a,b) for a = b = 1. - Stanislav Sykora, Apr 16 2015
The area of a circle circumscribing a unit-area regular hexagon. - Amiram Eldar, Nov 05 2020

Examples

			1.2091995761561452337293855050947704881893774987284937170465899569254...
		

References

  • George Boros and Victor H. Moll, Irresistible integrals, Cambridge University Press (2006), pp. 120-121.
  • L. B. W. Jolley, Summation of Series, Dover (1961), No. 261, pp. 48, 49, (and No. 275).

Crossrefs

Cf. A091682 (Sum_{i >= 0} (i!)^2/(2*i)!).

Programs

  • Mathematica
    RealDigits[2 Sqrt[3] Pi/9, 10, 100][[1]]
  • PARI
    a = 2*Pi/(3*sqrt(3)) \\ Stanislav Sykora, Apr 16 2015

Formula

Equals 2*sqrt(3)*Pi/9 = 1 + 1/6 + 1/30 + 1/140 + 1/630 + 1/2772 + 1/12012 + ...
Equals m*I_3(m,m) = m*Integral_{x>=0} (x/(m^3+x^3)), for any m>0. - Stanislav Sykora, Apr 16 2015
Equals Integral_{x>=0} (1/(1+x^3)) dx. - Robert FERREOL, Dec 23 2016
From Peter Bala, Oct 27 2019: (Start)
Equals 3/4*Sum_{n >= 0} (n+1)!*(n+2)!/(2*n+3)!.
Equals Sum_{n >= 1} 3^(n-1)/(n*binomial(2*n,n)).
Equals 2*Sum_{n >= 1} 1/(n*binomial(2*n,n)). See Boros and Moll, pp. 120-121.
Equals Integral_{x = 0..1} 1/(1 - x^3)^(1/3) dx = Sum_{n >= 0} (-1)^n*binomial(-1/3,n) /(3*n + 1).
Equals 2*Sum_{n >= 1} 1/((3*n-1)*(3*n-2)) = 2*(1 - 1/2 + 1/4 - 1/5 + 1/7 - 1/8 + ...) (added Oct 30 2019). (End)
Equals Product_{k>=1} 9*k^2/(9*k^2 - 1). - Amiram Eldar, Aug 04 2020
From Peter Bala, Dec 13 2021: (Start)
Equals (2/3)*A093602.
Conjecture: for k >= 0, 2*sqrt(3)*Pi/9 = (3/2)^k * k!*Sum_{n = -oo..oo} (-1)^n/ Product_{j = 0..k} (3*n + 3*j + 1). (End)
Equals (3/4)*S - 1, where S = A248682. - Peter Luschny, Jul 22 2022
Equals Integral_{x=0..Pi/2} tan(x)^(1/3)/(sin(2*x) + 1) dx. See MIT Link. - Joost de Winter, Aug 26 2023
Continued fraction: 1/(1 - 1/(7 - 12/(12 - 30/(17 - ... - 2*n*(2*n - 1)/((5*n + 2) - ... ))))). See A000407. - Peter Bala, Feb 20 2024
Equals Sum_{n>=2} 1/binomial(n, floor(n/2)); and trivially if "floor" is replaced by "ceiling". - Richard R. Forberg, Aug 30 2024
Equals Product_{k>=2} (1 + (-1)^k/A001651(k)). - Amiram Eldar, Nov 22 2024
Equals 2*A073010 = 1/A086089 = sqrt(A214549) = exp(A256923) = A275486/2. - Hugo Pfoertner, Nov 22 2024
Equals 1 - (1/6) * Sum_{n>=1} A010815(n)/n. - Friedjof Tellkamp, Apr 05 2025
Equals A248181 - 2. - Pontus von Brömssen, Apr 05 2025

A020784 Decimal expansion of 1/sqrt(27).

Original entry on oeis.org

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

Views

Author

Keywords

Comments

This is the minimum ripple factor for a third-order Chebyshev filter for which the generalized reflectionless topology needs no negative elements. - Matthew A. Morgan, Oct 18 2017

Examples

			0.1924500897298752548363829268339858185492005837567089586728674....
		

References

  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Sections 8.4.3 and 8.16, pp. 495, 527.

Crossrefs

Programs

Formula

Equals Sum_{k>=0} binomial(2*k,k) * k/16^k. - Amiram Eldar, Aug 02 2020
Equals sqrt(3)/9. - Stefano Spezia, Dec 24 2024
Equals 1/A010482 = A020760/3 = sqrt(A021031) = A073010/Pi = A212886/2. - Hugo Pfoertner, Dec 24 2024

A000917 a(n) = (2n+3)!/(n!*(n+2)!).

Original entry on oeis.org

3, 20, 105, 504, 2310, 10296, 45045, 194480, 831402, 3527160, 14872858, 62403600, 260757900, 1085822640, 4508102925, 18668849760, 77138650050, 318107374200, 1309542023790, 5382578744400, 22093039119060, 90567738003600, 370847442355650, 1516927277253024
Offset: 0

Views

Author

Keywords

Comments

G.f.: c(x)*(4-c(x))/(1-4*x)^(3/2), c(x) = g.f. for Catalan numbers A000108 (agrees with Hansen, 1975, p. 99, (5.27.9)). Convolution of A038679 with A000984 (central binomial coefficients); also convolution of A038665 with A000302 (powers of 4). - Wolfdieter Lang, Dec 11 1999
Appears as diagonal in A003506. - Zerinvary Lajos, Apr 12 2006
a(n) is the number of double rises in all Grand Dyck paths of semilength n+2. Example: a(0)=3 because in the 6 (=A000984(2)) Grand Dyck paths of semilength 2, namely udud, (uu)dd, uddu, d(uu)d, dudu, dd(uu), we have a total of 3 uu's (shown between parentheses). - Emeric Deutsch, Nov 29 2008

References

  • Eldon R. Hansen, A Table of Series and Products, Prentice-Hall, Englewood Cliffs, NJ, 1975, p. 99, (5.27.9).

Crossrefs

Programs

  • Magma
    [(n+1)*Binomial(2*n+3, n+1): n in [0..25]]; // Vincenzo Librandi, Jun 01 2016
  • Maple
    a := proc(n) (n+1)*binomial(2*n+3, n+2) end: seq(a(n), n=0..23); # Zerinvary Lajos, Nov 26 2006
    seq((n+1)*binomial(2*n+4, n+2)/2, n=0..23); # Zerinvary Lajos, Feb 28 2007
  • Mathematica
    Table[(2*n + 3)!/(n!*(n + 2)!), {n, 0, 25}] (* T. D. Noe, Jun 20 2012 *)

Formula

a(n) = (n+1)*binomial(2*n+3, n+1) = (n+1)*A001700(n+1). - Vincenzo Librandi, Jun 01 2016
a(n) = (2*n+3)*A001791(n+1). - R. J. Mathar, Nov 09 2021
D-finite with recurrence +(n+2)*a(n) +10*(-n-1)*a(n-1) +12*(2*n+1)*a(n-2)=0. - R. J. Mathar, Nov 09 2021
D-finite with recurrence n*(n+2)*a(n) -2*(2*n+3)*(n+1)*a(n-1)=0. - R. J. Mathar, Nov 09 2021
From Amiram Eldar, Jan 24 2022: (Start)
Sum_{n>=0} 1/a(n) = 1 - Pi/(3*sqrt(3)) = 1 - A073010.
Sum_{n>=0} (-1)^n/a(n) = 6*log(phi)/sqrt(5) - 1, where phi is the golden ratio (A001622). (End)
Showing 1-10 of 58 results. Next