cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 21-26 of 26 results.

A034323 a(n) = n-th quintic factorial number divided by 2.

Original entry on oeis.org

1, 7, 84, 1428, 31416, 848232, 27143424, 1004306688, 42180880896, 1982501402112, 103090072909824, 5876134155859968, 364320317663318016, 24409461283442307072, 1757481212407846109184, 135326053355404150407168, 11096736375143140333387776, 965416064637453209004736512
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • GAP
    a:=[1];; for n in [2..20] do a[n]:=(5*n-3)*a[n-1]; od; a; # G. C. Greubel, Feb 10 2019
  • Magma
    [(&*[5*j-3: j in [1..n]])/2: n in [1..20]]; // G. C. Greubel, Feb 10 2019
    
  • Maple
    f:= gfun:-rectoproc({a(n)=(5*n-3)*a(n-1),a(1)=1},a(n),remember):
    map(f, [$1..40]); # Robert Israel, Feb 10 2019
  • Mathematica
    Table[Product[5j-3,{j,n}]/2,{n,20}] (* Harvey P. Dale, Nov 25 2013 *)
  • PARI
    vector(20, n, prod(j=1, n, 5*j-3)/2) \\ G. C. Greubel, Feb 10 2019
    
  • Sage
    [product(5*j-3 for j in (1..n))/2 for n in (1..20)] # G. C. Greubel, Feb 10 2019
    

Formula

2*a(n) = (5*n-3)(!^5) = Product_{j=1..n} (5*j-3).
E.g.f.: (-1 + (1-5*x)^(-2/5))/2, with a(0) = 0.
a(n) ~ sqrt(2*Pi) * 5/(2*Gamma(2/5)) * n^(9/10) * (5*n/e)^n * (1 + (109/300)/n - ...). - Joe Keane (jgk(AT)jgk.org), Nov 24 2001
D-finite with recurrence a(n) = (5*n-3)*a(n-1). - Robert Israel, Feb 10 2019
From Amiram Eldar, Dec 19 2022: (Start)
a(n) = A047055(n)/2.
Sum_{n>=1} 1/a(n) = 2*(e/5^3)^(1/5)*(Gamma(2/5) - Gamma(2/5, 1/5)). (End)

A142460 Triangle read by rows: T(n,k) (1<=k<=n) given by T(n, 1) = T(n,n) = 1, otherwise T(n, k) = (m*n-m*k+1)*T(n-1,k-1) + (m*k-m+1)*T(n-1,k), where m = 5.

Original entry on oeis.org

1, 1, 1, 1, 12, 1, 1, 83, 83, 1, 1, 514, 1826, 514, 1, 1, 3105, 28310, 28310, 3105, 1, 1, 18656, 376615, 905920, 376615, 18656, 1, 1, 111967, 4627821, 22403635, 22403635, 4627821, 111967, 1, 1, 671838, 54377008, 478781506, 940952670, 478781506, 54377008, 671838, 1
Offset: 1

Views

Author

Roger L. Bagula, Sep 19 2008

Keywords

Comments

One of a family of triangles. For m = ...,-2,-1,0,1,2,3,4,5,... we get ..., A225372, A144431, A007318, A008292, A060187, A142458, A142459, A142560, ...

Examples

			Triangle begins as:
  1;
  1,      1;
  1,     12,        1;
  1,     83,       83,         1;
  1,    514,     1826,       514,         1;
  1,   3105,    28310,     28310,      3105,         1;
  1,  18656,   376615,    905920,    376615,     18656,        1;
  1, 111967,  4627821,  22403635,  22403635,   4627821,   111967,      1;
  1, 671838, 54377008, 478781506, 940952670, 478781506, 54377008, 671838, 1;
		

Crossrefs

Cf. A225372 (m=-2), A144431 (m=-1), A007318 (m=0), A008292 (m=1), A060187 (m=2), A142458 (m=3), A142459 (m=4), this sequence (m=5), A142561 (m=6), A142562 (m=7), A167884 (m=8), A257608 (m=9).
Cf. A047055 (row sums).

Programs

  • Maple
    A142460 := proc(n, k) if n = k then 1; elif k > n or k < 1 then 0 ; else (5*n-5*k+1)*procname(n-1, k-1)+(5*k-4)*procname(n-1, k) ; end if; end proc:
    seq(seq(A142459(n, k), k=1..n), n=1..10) ; # R. J. Mathar, May 11 2013
  • Mathematica
    T[n_, k_, m_]:= T[n, k, m]= If[k==1 || k==n, 1, (m*n-m*k+1)*T[n-1, k-1, m] + (m*k -m+1)*T[n-1, k, m] ];
    Table[T[n, k, 5], {n, 1, 10}, {k, 1, n}]//Flatten (* modified by G. C. Greubel, Mar 14 2022 *)
  • Sage
    def T(n,k,m): # A142460
        if (k==1 or k==n): return 1
        else: return (m*(n-k)+1)*T(n-1,k-1,m) + (m*k-m+1)*T(n-1,k,m)
    flatten([[T(n,k,5) for k in (1..n)] for n in (1..10)]) # G. C. Greubel, Mar 14 2022

Formula

T(n, k, m) = (m*n - m*k + 1)*T(n-1, k-1, m) + (m*k - (m-1))*T(n-1, k, m), with T(t,1,m) = T(n,n,m) = 1, and m = 5.
Sum_{k=1..n} T(n, k, 5) = A047055(n-1).

Extensions

Edited by N. J. A. Sloane, May 08 2013, May 11 2013

A129116 Multifactorial array: A(k,n) = k-tuple factorial of n, for positive n, read by ascending antidiagonals.

Original entry on oeis.org

1, 1, 2, 1, 2, 6, 1, 2, 3, 24, 1, 2, 3, 8, 120, 1, 2, 3, 4, 15, 720, 1, 2, 3, 4, 10, 48, 5040, 1, 2, 3, 4, 5, 18, 105, 40320, 1, 2, 3, 4, 5, 12, 28, 384, 362880, 1, 2, 3, 4, 5, 6, 21, 80, 945, 3628800, 1, 2, 3, 4, 5, 6, 14, 32, 162, 3840, 39916800, 1, 2, 3, 4, 5, 6, 7, 24, 45, 280, 10395, 479001600
Offset: 1

Views

Author

Jonathan Vos Post, May 24 2007

Keywords

Comments

The term "Quintuple factorial numbers" is also used for the sequences A008546, A008548, A052562, A047055, A047056 which have a different definition. The definition given here is the one commonly used. This problem exists for the other rows as well. "n!2" = n!!, "n!3" = n!!!, "n!4" = n!!!!, etcetera. Main diagonal is A[n,n] = n!n = n.
Similar to A114423 (with rows and columns exchanged). - Georg Fischer, Nov 02 2021

Examples

			Table begins:
  k / A(k,n)
  1 | 1 2 6 24 120 720 5040 40320 362880 3628800 ... = A000142.
  2 | 1 2 3  8  15  48  105   384    945    3840 ... = A006882.
  3 | 1 2 3  4  10  18   28    80    162     280 ... = A007661.
  4 | 1 2 3  4   5  12   21    32     45     120 ... = A007662.
  5 | 1 2 3  4   5   6   14    24     36      50 ... = A085157.
  6 | 1 2 3  4   5   6    7    16     27      40 ... = A085158.
		

Crossrefs

Cf. A000142 (n!), A006882 (n!!), A007661 (n!!!), A007662(n!4), A085157 (n!5), A085158 (n!6), A114799 (n!7), A114800 (n!8), A114806 (n!9), A288327 (n!10).
Cf. A114423 (transposed).

Programs

  • Maple
    A:= proc(k,n) option remember; if n >= 1 then n* A(k, n-k) elif n >= 1-k then 1 else 0 fi end: seq(seq(A(1+d-n, n), n=1..d), d=1..16); # Alois P. Heinz, Feb 02 2009
  • Mathematica
    A[k_, n_] := A[k, n] = If[n >= 1, n*A[k, n-k], If[n >= 1-k, 1, 0]]; Table[ A[1+d-n, n], {d, 1, 16}, {n, 1, d}] // Flatten (* Jean-François Alcover, May 27 2016, after Alois P. Heinz *)

Formula

A(k,n) = n!k.
A(k,n) = M(n,k) in A114423. - Georg Fischer, Nov 02 2021

Extensions

Corrected and extended by Alois P. Heinz, Feb 02 2009

A303488 a(n) = n! * [x^n] 1/(1 - 5*x)^(n/5).

Original entry on oeis.org

1, 1, 14, 312, 9576, 375000, 17873856, 1004306688, 65006637696, 4763494479744, 389812500000000, 35237024762075136, 3487065897634615296, 374960171943074285568, 43532820293400237735936, 5427359437500000000000000, 723181462895975365595529216, 102563963819340862347122245632
Offset: 0

Views

Author

Ilya Gutkovskiy, Apr 24 2018

Keywords

Examples

			a(1) = 1;
a(2) = 2*7 = 14;
a(3) = 3*8*13 = 312;
a(4) = 4*9*14*19 = 9576;
a(5) = 5*10*15*20*25 = 375000, etc.
		

Crossrefs

Programs

  • Mathematica
    Table[n! SeriesCoefficient[1/(1 - 5 x)^(n/5), {x, 0, n}], {n, 0, 17}]
    Table[Product[5 k + n, {k, 0, n - 1}], {n, 0, 17}]
    Table[5^n Pochhammer[n/5, n], {n, 0, 17}]

Formula

a(n) = Product_{k=0..n-1} (5*k + n).
a(n) = 5^n*Gamma(6*n/5)/Gamma(n/5).
a(n) ~ 6^(6*n/5-1/2)*n^n/exp(n).

A262246 Decimal expansion of Sum_{k>=0} (-1)^k/(5k+2).

Original entry on oeis.org

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

Views

Author

Gheorghe Coserea, Oct 06 2015

Keywords

Examples

			0.4069016342...
		

Crossrefs

Programs

  • Mathematica
    N[(1/5)*((Sqrt[5]-1)*Log[2] + Sqrt[5]*Log[Sin[3*Pi/10]] + (Pi/2)*Sec[Pi/10]), 100] (* G. C. Greubel, Oct 07 2015 *) (* fixed by Vaclav Kotesovec, Dec 11 2017 *)
  • PARI
    default(realprecision, 87);
    eval(vecextract(Vec(Str(sumalt(n=0, (-1)^(n)/(5*n+2)))), "3..-2"))

Formula

Sum_{n>=0} (-1)^n/(5n+2) = Integral_{x=0..1} x/(1+x^5)dx.
From G. C. Greubel, Oct 07 2015: (Start)
Sum_{n>=0} (-1)^n/(5n+2) = (1/5)*(sqrt(5)*log(phi) - log(2) + Pi*(5*phi^2)^(-1/4)), where 2*phi=1+sqrt(5).
Sum_{n>=0} (-1)^n/(5n+2) = (1/5)*(sqrt(5)*log(2*sin(3*Pi/10)) - log(2) + (Pi/2)*sec(Pi/10)).
(End)
Sum_{n>=0} (-1)^n/(5n+2) = (Psi(1/5) - Psi(7/10))/10 , see A200135 and A354643. - Robert Israel, Oct 08 2015
From Peter Bala, Feb 19 2024: (Start)
Equals (1/2)*Sum_{n >= 0} n!*(5/2)^n/(Product_{k = 0..n} 5*k + 2) = (1/2)*Sum_{n >= 0} n!*(5/2)^n/A047055(n+1) (apply Euler's series transformation to Sum_{k >= 0} (-1)^k/(5*k + 2)).
Continued fraction: 1/(2 + 2^2/(5 + 7^2/(5 + 12^2/(5 + ... + (5*n + 2)^2/(5 + ... ))))).
The slowly converging series representation Sum_{n >= 0} (-1)^n/(5*n + 2) for the constant can be accelerated to give the following faster converging series
1/4 + (5/2)*Sum_{n >= 0} (-1)^n/((5*n + 2)*(5*n + 7)) and
19/56 + (5^2/2)*Sum_{n >= 0} (-1)^n/((5*n + 2)*(5*n + 7)*(5*n + 12)).
These two series are the cases r = 1 and r = 2 of the general result:
for r >= 0, the constant equals C(r) + ((5/2)^r)*r!*Sum_{n >= 0} (-1)^n/((5*n + 2)*(5*n + 7)*...*(5*n + 5*r + 2)), where C(r) is the rational number (1/2)*Sum_{k = 0..r-1} (5/2)^k*k!/(2*7*12*...*(5*k + 2)). The general result can be proved by the WZ method as described in Wilf. (End)
From Peter Bala, Mar 03 2024: (Start)
Equals (1/2)*hypergeom([2/5, 1], [7/5], -1).
Gauss's continued fraction: 1/(2 + 2^2/(7 + 5^2/(12 + 7^2/(17 + 10^2/(22 + 12^2/(27 + 15^2/(32 + 17^2/(37 + 20^2/(42 + 22^2/(47 + ... )))))))))). (End)

A020039 Nearest integer to Gamma(n + 2/5)/Gamma(2/5).

Original entry on oeis.org

1, 0, 1, 1, 5, 20, 109, 695, 5142, 43193, 406016, 4222569, 48137291, 596902408, 7998492273, 115178288736, 1773745646540, 29089428603255, 506156057696641, 9313271461618196, 180677466355392996, 3685820313650017111
Offset: 0

Views

Author

Keywords

Comments

Gamma(n + 2/5)/Gamma(2/5) = 1, 2/5, 14/25, 168/125, 2856/625, 62832/3125, 1696464/15625, 54286848/78125, ... - R. J. Mathar, Sep 04 2016

Crossrefs

Programs

  • Maple
    Digits := 64:f := proc(n,x) round(GAMMA(n+x)/GAMMA(x)); end;
Previous Showing 21-26 of 26 results.