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-6 of 6 results.

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

A005321 Upper triangular n X n (0,1)-matrices with no zero rows or columns.

Original entry on oeis.org

1, 1, 2, 10, 122, 3346, 196082, 23869210, 5939193962, 2992674197026, 3037348468846562, 6189980791404487210, 25285903982959247885402, 206838285372171652078912306, 3386147595754801373061066905042, 110909859519858523995273393471390010
Offset: 0

Views

Author

Keywords

References

  • T. L. Greenough, Enumeration of interval orders without duplicated holdings, Preprint, circa 1976.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Column sums of A137252.

Programs

  • Mathematica
    max = 14; f[x_] := Sum[ x^n*Product[ (2^i-1) / (1+(2^i-1)*x), {i, 1, n}], {n, 0, max}]; CoefficientList[ Series[ f[x], {x, 0, max}], x] (* Jean-François Alcover, Nov 23 2011, after Vladeta Jovovic *)
  • PARI
    a(n) = 1 + sum(k=2, n, binomial(n,k)*sum(i=2, k, (-1)^i*prod(j=i+1, k, 2^j - 1))); \\ Michel Marcus, Oct 13 2019

Formula

a(n) = Sum_{k=0..n} binomial(n,k)*A005327(k+1).
G.f.: Sum_{n >= 0} x^n*Product_{i = 1..n} ((2^i-1)/(1 + (2^i-1)*x)). - Vladeta Jovovic, Mar 10 2008
From Peter Bala, Jul 06 2017: (Start)
Two conjectural continued fractions for the o.g.f.:
1/(1 - x/(1 - x/(1 - 6*x/(1 - 9*x/(1 - 28*x/(1 - 49*x/(1 - ... - 2^(n-1)*(2^n - 1)*x/(1 - (2^n - 1)^2*x/(1 - ...)))))))));
1 + x/(1 - 2*x/(1 - 3*x/(1 - 12*x/(1 - 21*x/(1 - ... - 2^n*(2^n - 1)*x/(1 - (2^(n+1) - 1)*(2^n - 1)*x/(1 - ...))))))). Cf. A289314 and A289315. (End)
a(n) = (-1)^n*Sum_{k=0..n} qS2(n,k)*[k]!*(-1)^k, where qS2(n,k) is the triangle A139382, and [k]! is q-factorial, q=2. - Vladimir Kruchinin, Oct 10 2019
a(n) = 1 + Sum_{k=2..n} binomial(n,k)*Sum{i=2..k} (-1)^i*Product_{j=i+1..k} (2^j - 1). See Greenough. - Michel Marcus, Oct 13 2019

Extensions

More terms from Max Alekseyev, Apr 27 2010

A002820 Number of n X n invertible binary matrices A such that A+I is invertible.

Original entry on oeis.org

1, 0, 2, 48, 5824, 2887680, 5821595648, 47317927329792, 1544457148312846336, 202039706313624586813440, 105823549214125066767168438272, 221819704567105547916502447159246848, 1860304261534304703934696550224148083769344, 62413833036707798343389591015829588620560344023040
Offset: 0

Views

Author

Keywords

Comments

Also number of linear orthomorphisms of GF(2)^n.

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

Cf. A002884.
Column k=0 of A346201 and of A346381.

Programs

  • Maple
    # (Maple program based on Dai et al. from N. J. A. Sloane, Aug 10 2011)
    N:=proc(n,i) option remember; if i = 1 then 1 else (2^n-2^(i-1))*N(n,i-1); fi; end;
    Oh:=proc(n) option remember; local r; global N;
    if n=0 then RETURN(1) elif n=1 then RETURN(0) else
    add( 2^(r-2)*N(n,r)*2^(r*(n-r))*Oh(n-r), r=2..n); fi; end;
    [seq(Oh(n),n=0..15)];
  • Mathematica
    ni[n_, i_] := ni[n, i] = If[i == 1, 1, (2^n - 2^(i-1))*ni[n, i-1]]; a[0] = 1; a[1] = 0; a[n_] := a[n] = Sum[ 2^(r-2)*ni[n, r]*2^(r*(n-r))*a[n-r], {r, 2, n}]; Table[a[n], {n, 1, 11}] (* Jean-François Alcover, Jan 19 2012, after Maple *)

Formula

Reference gives a recurrence.
a(n) = 2^(n(n-1)/2) * A005327(n+1).

Extensions

More terms from Vladeta Jovovic, Mar 17 2000
Entry revised by N. J. A. Sloane, Aug 10 2011
a(0)=1 prepended by Alois P. Heinz, Jan 10 2025

A005014 Certain subgraphs of a directed graph (inverse binomial transform of A005321).

Original entry on oeis.org

1, 1, 7, 97, 2911, 180481, 22740607, 5776114177, 2945818230271, 3010626231336961, 6159741269315422207, 25217980756577338515457, 206535262396368402441592831, 3383460668577307168798173757441
Offset: 1

Views

Author

Keywords

References

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

Crossrefs

Pairwise sums of A005327.

Programs

  • Maple
    p:=proc(n) if n=0 then 1 else product(2^i-1,i=1..n) fi end: a:=n->(-1)^n+(p(n)+p(n-1))*sum((-1)^j/p(j),j=0..n-1): seq(a(n),n=1..14); # Emeric Deutsch, Jan 23 2005
  • Mathematica
    a[1] = 1; a[n_] := a[n] = (2^n-2)*a[n-1]-(-1)^n; Table[a[n], {n, 1, 14}] (* Jean-François Alcover, Jan 17 2014, after Vladeta Jovovic *)

Formula

a(n) = (-1)^n + (p(n) + p(n-1))Sum_{j=0..n-1} (-1)^j/p(j), where p(0)=1, p(k) = Product_{i=1..k} (2^i - 1) for k > 0. - Emeric Deutsch, Jan 23 2005
a(n) = (2^n-2)*a(n-1) - (-1)^n. - Vladeta Jovovic, Aug 20 2006
G.f.: Sum_{n>=0} (x^n*Product_{i=1..n} (2^i - 1)/(1 + 2^i*x)). - Vladeta Jovovic, Mar 10 2008

Extensions

More terms from Vladeta Jovovic, Aug 20 2006

A259970 Triangle read by rows: coefficients eta(n,k) arising from the study of completely transitive graphs on n nodes.

Original entry on oeis.org

1, 0, 1, 1, 3, 3, 6, 28, 42, 21, 91, 510, 1050, 945, 315, 2820, 18631, 48360, 61845, 39060, 9765, 177661, 1351413, 4220433, 6942915, 6357015, 3075975, 615195
Offset: 1

Views

Author

N. J. A. Sloane, Jul 11 2015

Keywords

Examples

			Triangle begins:
1,
0,1,
1,3,3,
6,28,42,21,
91,510,1050,945,315,
2820,18631,48360,61845,39060,9765,
177661,1351413,4220433,6942915,6357015,3075975,615195,
...
		

Crossrefs

Diagonals include A005327, A005328, A005329.
Row sums are A005014.

A259876 Triangle of numbers S(n,k) (0 <= k <= n) arising in the enumeration of interval orders without duplicated holdings.

Original entry on oeis.org

1, 1, -1, 3, -3, 1, 21, -21, 7, -1, 315, -315, 105, -15, 1, 9765, -9765, 3255, -465, 31, -1, 615195, -615195, 205065, -29295, 1953, -63, 1, 78129765, -78129765, 26043255, -3720465, 248031, -8001, 127, -1, 19923090075, -19923090075, 6641030025, -948718575, 63247905, -2040255, 32385, -255, 1
Offset: 0

Views

Author

N. J. A. Sloane, Jul 09 2015

Keywords

Examples

			Triangle begins:
     1;
     1,    -1;
     3,    -3,    1;
    21,   -21,    7,   -1;
   315,  -315,  105,  -15,  1;
  9765, -9765, 3255, -465, 31, -1;
  ...
		

References

  • T. L. Greenough, Enumeration of interval orders without duplicated holdings, Preprint, circa 1976.

Crossrefs

Row sums give A005327.
Column k=0 gives A005329.
Main diagonal gives A033999.
T(n+1,n) gives A225883(n+1).

Formula

T(n,k) = qfactorial(n)/qfactorial(k)*(-1)^(k), n>=k, where qfactorial(n) is A005329. - Vladimir Kruchinin, Feb 17 2020

Extensions

More terms from Alois P. Heinz, Feb 17 2020
Showing 1-6 of 6 results.