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 18 results. Next

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

A005329 a(n) = Product_{i=1..n} (2^i - 1). Also called 2-factorial numbers.

Original entry on oeis.org

1, 1, 3, 21, 315, 9765, 615195, 78129765, 19923090075, 10180699028325, 10414855105976475, 21319208401933844325, 87302158405919092510875, 715091979502883286756577125, 11715351900195736886933003038875, 383876935713713710574133710574817125
Offset: 0

Views

Author

Keywords

Comments

Conjecture: this sequence is the inverse binomial transform of A075272 or, equivalently, the inverse binomial transform of the BinomialMean transform of A075271. - John W. Layman, Sep 12 2002
To win a game, you must flip n+1 heads in a row, where n is the total number of tails flipped so far. Then the probability of winning for the first time after n tails is A005329 / A006125. The probability of having won before n+1 tails is A114604 / A006125. - Joshua Zucker, Dec 14 2005
Number of upper triangular n X n (0,1)-matrices with no zero rows. - Vladeta Jovovic, Mar 10 2008
Equals the q-Fibonacci series for q = (-2), and the series prefaced with a 1: (1, 1, 1, 3, 21, ...) dot (1, -2, 4, -8, ...) if n is even, and (-1, 2, -4, 8, ...) if n is odd. For example, a(3) = 21 = (1, 1, 1, 3) dot (-1, 2, -4, 8) = (-1, 2, -4, 24) and a(4) = 315 = (1, 1, 1, 3, 21) dot (1, -2, 4, -8 16) = (1, -2, 4, -24, 336). - Gary W. Adamson, Apr 17 2009
Number of chambers in an A_n(K) building where K=GF(2) is the field of two elements. This is also the number of maximal flags in an n-dimensional vector space over a field of two elements. - Marcos Spreafico, Mar 22 2012
Given probability p = 1/2^n that an outcome will occur at the n-th stage of an infinite process, then starting at n=1, A114604(n)/A006125(n+2) = 1-a(n)/A006125(n+1) is the probability that the outcome has occurred up to and including the n-th iteration. The limiting ratio is 1-A048651 ~ 0.7112119. These observations are a more formal and generalized statement of Joshua Zucker's Dec 14, 2005 comment. - Bob Selcoe, Mar 02 2016
Also the number of dominating sets in the n-triangular honeycomb rook graph. - Eric W. Weisstein, Jul 14 2017
Empirical: Letting Q denote the Hall-Littlewood Q basis of the symmetric functions over the field of fractions of the univariate polynomial ring in t over the field of rational numbers, and letting h denote the complete homogeneous basis, a(n) is equal to the absolute value of 2^A000292(n) times the coefficient of h_{1^(n*(n+1)/2)} in Q_{(n, n-1, ..., 1)} with t evaluated at 1/2. - John M. Campbell, Apr 30 2018
The series f(x) = Sum_{n>=0} x^(2^n-1)/a(n) satisfies f'(x) = f(x^2), f(0) = 1. - Lucas Larsen, Jan 05 2022

Examples

			G.f. = 1 + x + 3*x^2 + 21*x^3 + 315*x^4 + 9765*x^5 + 615195*x^6 + 78129765*x^7 + ...
		

References

  • Annie Cuyt, Vigdis Brevik Petersen, Brigitte Verdonk, Haakon Waadeland, and William B. Jones, Handbook of continued fractions for special functions, Springer, New York, 2008. (see 19.2.1)
  • Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, p. 358.
  • Mark Ronan, Lectures on Buildings (Perspectives in Mathematics; Vol. 7), Academic Press Inc., 1989.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A048651, A079555, A152476 (inverse binomial transform).
Column q=2 of A069777.

Programs

  • GAP
    List([0..15],n->Product([1..n],i->2^i-1)); # Muniru A Asiru, May 18 2018
  • Magma
    [1] cat [&*[ 2^k-1: k in [1..n] ]: n in [1..16]]; // Vincenzo Librandi, Dec 24 2015
    
  • Maple
    A005329 := proc(n) option remember; if n<=1 then 1 else (2^n-1)*procname(n-1); end if; end proc: seq(A005329(n), n=0..15);
  • Mathematica
    a[0] = 1; a[n_] := a[n] = (2^n-1)*a[n-1]; a /@ Range[0,14] (* Jean-François Alcover, Apr 22 2011 *)
    FoldList[Times, 1, 2^Range[15] - 1] (* Harvey P. Dale, Dec 21 2011 *)
    Table[QFactorial[n, 2], {n, 0, 14}] (* Arkadiusz Wesolowski, Oct 30 2012 *)
    QFactorial[Range[0, 10], 2] (* Eric W. Weisstein, Jul 14 2017 *)
    a[ n_] := If[ n < 0, 0, (-1)^n QPochhammer[ 2, 2, n]]; (* Michael Somos, Jan 28 2018 *)
  • PARI
    a(n)=polcoeff(sum(m=0,n,2^(m*(m+1)/2)*x^m/prod(k=0,m,1+2^k*x+x*O(x^n))),n) \\ Paul D. Hanna, Sep 17 2009
    
  • PARI
    Dx(n,F)=local(D=F);for(i=1,n,D=deriv(D));D
    a(n)=local(A=1+x+x*O(x^n));for(i=1,n,A=1+sum(k=1,n,x^k/k!*Dx(k,x*A+x*O(x^n) ))); polcoeff(A,n) \\ Paul D. Hanna, Apr 21 2012
    
  • PARI
    {a(n) = if( n<0, 0, prod(k=1, n, 2^k - 1))}; /* Michael Somos, Jan 28 2018 */
    
  • PARI
    {a(n) = if( n<0, 0, (-1)^n * sum(k=0, n+1, (-1)^k * 2^(k*(k+1)/2) * prod(j=1, k, (2^(n+1-j) - 1) / (2^j - 1))))}; /* Michael Somos, Jan 28 2018 */
    

Formula

a(n)/2^(n*(n+1)/2) -> c = 0.2887880950866024212788997219294585937270... (see A048651, A048652).
From Paul D. Hanna, Sep 17 2009: (Start)
G.f.: Sum_{n>=0} 2^(n*(n+1)/2) * x^n / (Product_{k=0..n} (1+2^k*x)).
Compare to: 1 = Sum_{n>=0} 2^(n*(n+1)/2) * x^n/(Product_{k=1..n+1} (1+2^k*x)). (End)
G.f. satisfies: A(x) = 1 + Sum_{n>=1} x^n/n! * d^n/dx^n x*A(x). - Paul D. Hanna, Apr 21 2012
a(n) = 2^(binomial(n+1,2))*(1/2; 1/2){n}, where (a;q){n} is the q-Pochhammer symbol. - G. C. Greubel, Dec 23 2015
a(n) = Product_{i=1..n} A000225(i). - Michel Marcus, Dec 27 2015
From Peter Bala, Nov 10 2017: (Start)
O.g.f. as a continued fraction of Stieltjes' type: A(x) = 1/(1 - x/(1 - 2*x/(1 - 6*x/(1 - 12*x/(1 - 28*x/(1 - 56*x/(1 - ... -(2^n - 2^floor(n/2))*x/(1 - ... )))))))) (follows from Heine's continued fraction for the ratio of two q-hypergeometric series at q = 2. See Cuyt et al. 19.2.1).
A(x) = 1/(1 + x - 2*x/(1 - (2 - 1)^2*x/(1 + x - 2^3*x/(1 - (2^2 - 1)^2*x/(1 + x - 2^5*x/(1 - (2^3 - 1)^2*x/(1 + x - 2^7*x/(1 - (2^4 - 1)^2*x/(1 + x - ... ))))))))). (End)
0 = a(n)*(a(n+1) - a(n+2)) + 2*a(n+1)^2 for all n>=0. - Michael Somos, Feb 23 2019
From Amiram Eldar, Feb 19 2022: (Start)
Sum_{n>=0} 1/a(n) = A079555.
Sum_{n>=0} (-1)^n/a(n) = A048651. (End)

Extensions

Better definition from Leslie Ann Goldberg (leslie(AT)dcs.warwick.ac.uk), Dec 11 1999

A028362 Total number of self-dual binary codes of length 2n. Totally isotropic spaces of index n in symplectic geometry of dimension 2n.

Original entry on oeis.org

1, 3, 15, 135, 2295, 75735, 4922775, 635037975, 163204759575, 83724041661975, 85817142703524375, 175839325399521444375, 720413716161839357604375, 5902349576513949856852644375, 96709997811181068404530578084375
Offset: 1

Views

Author

Keywords

Comments

These numbers appear in the second column of A155103. - Mats Granvik, Jan 20 2009
a(n) = n terms in the sequence (1, 2, 4, 8, 16, ...) dot n terms in the sequence (1, 1, 3, 15, 135). Example: a(5) = 2295 = (1, 2, 4, 8, 16) dot (1, 1, 3, 15, 135) = (1 + 2 + 12 + 120 + 2160). - Gary W. Adamson, Aug 02 2010

Examples

			G.f. = x + 3*x^2 + 15*x^3 + 135*x^4 + 2295*x^5 + 75735*x^6 + 4922775*x^7 + ...
		

References

  • F. J. MacWilliams and N. J. A. Sloane, The Theory of Error-Correcting Codes, Elsevier-North Holland, 1978, p. 630.

Crossrefs

Cf. A155103. - Mats Granvik, Jan 20 2009
Cf. A005329, A006088. - Paul D. Hanna, Sep 16 2009

Programs

  • Magma
    [1] cat [&*[ 2^k+1: k in [1..n] ]: n in [1..16]]; // Vincenzo Librandi, Dec 24 2015
    
  • Maple
    seq(mul(1 + 2^j, j = 1..n-1), n = 1..20); # G. C. Greubel, Jun 06 2020
  • Mathematica
    Table[Product[2^i+1,{i,n-1}],{n,15}] (* or *) FoldList[Times,1, 2^Range[15]+1] (* Harvey P. Dale, Nov 21 2011 *)
    Table[QPochhammer[-2, 2, n - 1], {n, 15}] (* Arkadiusz Wesolowski, Oct 29 2012 *)
  • PARI
    {a(n)=polcoeff(sum(m=0,n,2^(m*(m-1)/2)*x^m/prod(k=0,m-1,1-2^k*x+x*O(x^n))),n)} \\ Paul D. Hanna, Sep 16 2009
    
  • PARI
    {a(n) = if( n<1, 0 , prod(k=1, n-1, 2^k + 1))}; /* Michael Somos, Jan 28 2018 */
    
  • PARI
    {a(n) = sum(k=0, n-1, 2^(k*(k+1)/2) * prod(j=1, k, (2^(n-j) - 1) / (2^j - 1)))}; /* Michael Somos, Jan 28 2018 */
    
  • Python
    for n in range(2,40,2):
      product = 1
      for i in range(1,n//2-1 + 1):
        product *= (2**i+1)
      print(product)
    # Nathan J. Russell, Mar 01 2016
    
  • Python
    from math import prod
    def A028362(n): return prod((1<Chai Wah Wu, Jun 20 2022
    
  • Sage
    from ore_algebra import *
    R. = QQ['x']
    A. = OreAlgebra(R, 'Qx', q=2)
    print((Qx - x - 1).to_list([0,1], 10))  # Ralf Stephan, Apr 24 2014
    
  • Sage
    from sage.combinat.q_analogues import q_pochhammer
    [q_pochhammer(n-1,-2,2) for n in (1..20)] # G. C. Greubel, Jun 06 2020
    
  • Scheme
    ;; With memoization-macro definec.
    (define (A028362 n) (A028362off0 (- n 1)))
    (definec (A028362off0 n) (if (zero? n) 1 (+ (A028362off0 (- n 1)) (* (expt 2 n) (A028362off0 (- n 1))))))
    ;; Antti Karttunen, Apr 15 2017

Formula

a(n) = Product_{i=1..n-1} (2^i+1).
Letting a(0)=1, we have a(n) = Sum_{k=0..n-1} 2^k*a(k) for n>0. a(n) is asymptotic to c*sqrt(2)^(n^2-n) where c=2.384231029031371724149899288.... = A079555 = Product_{k>=1} (1 + 1/2^k). - Benoit Cloitre, Jan 25 2003
G.f.: Sum_{n>=1} 2^(n*(n-1)/2) * x^n/(Product_{k=0..n-1} (1-2^k*x)). - Paul D. Hanna, Sep 16 2009
a(n) = 2^(binomial(n,2) - 1)*(-1; 1/2){n}, where (a;q){n} is the q-Pochhammer symbol. - G. C. Greubel, Dec 23 2015
From Antti Karttunen, Apr 15 2017: (Start)
a(n) = A048675(A285101(n-1)).
a(n) = b(n-1), where b(0) = 1, and for n > 0, b(n) = b(n-1) + (2^n)*b(n-1).
a(n) = Sum_{i=1..A000124(n-1)} A053632(n-1,i-1)*(2^(i-1)) [where the indexing of both rows and columns of irregular table A053632(row,col) is considered to start from zero].
(End)
G.f. A(x) satisfies: A(x) = x * (1 + A(2*x)) / (1 - x). - Ilya Gutkovskiy, Jun 06 2020
Conjectural o.g.f. as a continued fraction of Stieltjes type (S-fraction):
1/(1 - 3*x/(1 - 2*x/(1 - 10*x/(1 - 12*x/(1 - 36*x/(1 - 56*x/(1 - 136*x/(1 - 240*x/(1 - ... - 2^(n-1)*(2^n + 1)*x/(1 - 2^n*(2^n - 1)*x/(1 - ... ))))))))))). - Peter Bala, Sep 27 2023

A081845 Decimal expansion of Product_{k>=0} (1 + 1/2^k).

Original entry on oeis.org

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

Views

Author

Benoit Cloitre, Apr 09 2003

Keywords

Comments

Twice the product in A079555.

Examples

			4.76846205806274344829979857....
		

Crossrefs

Programs

  • Mathematica
    digits = 105; NProduct[1 + 1/2^k, {k, 0, Infinity}, WorkingPrecision -> digits+5, NProductFactors -> digits] // RealDigits[#, 10, digits]& // First (* Jean-François Alcover, Mar 04 2013 *)
    N[QPochhammer[-1, 1/2], 100] (* Vaclav Kotesovec, Dec 13 2015 *)
    2*N[QPochhammer[-1/2, 1/2], 200] (* G. C. Greubel, Dec 20 2015 *)
  • PARI
    prodinf(k=0,1/2^k,1) \\ Hugo Pfoertner, Feb 21 2020

Formula

lim sup Product_{k=0..floor(log_2(n))} (1 + 1/floor(n/2^k)) for n-->oo. - Hieronymus Fischer, Aug 20 2007
lim sup A132369(n)/A098844(n) for n-->oo. - Hieronymus Fischer, Aug 20 2007
lim sup A132269(n)/n^((1+log_2(n))/2) for n-->oo. - Hieronymus Fischer, Aug 20 2007
lim sup A132270(n)/n^((log_2(n)-1)/2) for n-->oo. - Hieronymus Fischer, Aug 20 2007
2*exp(Sum_{n>0} 2^(-n)*Sum_{k|n} -(-1)^k/k) = 2*exp(Sum_{n>0} A000593(n)/(n*2^n)). - Hieronymus Fischer, Aug 20 2007
lim sup A132269(n+1)/A132269(n) = 4.76846205806274344... for n-->oo. - Hieronymus Fischer, Aug 20 2007
Sum_{k>=1} (-1)^(k+1) * 2^k / (k*(2^k-1)) = log(A081845) = 1.562023833218500307570359922772014353168080202860122... . - Vaclav Kotesovec, Dec 13 2015
Equals 2*(-1/2; 1/2){infinity}, where (a;q){infinity} is the q-Pochhammer symbol. - G. C. Greubel, Dec 20 2015
Equals 1 + Sum_{n>=1} 2^n/((2-1)*(2^2-1)*...*(2^n-1)). - Robert FERREOL, Feb 21 2020
From Peter Bala, Jan 18 2021: (Start)
Constant C = 3*Sum_{n >= 0} (1/2)^n/Product_{k = 1..n} (2^k - 1).
Faster converging series:
C = (2*3*5)/(2^3)*Sum_{n >= 0} (1/4)^n/Product_{k = 1..n} (2^k - 1),
C = (2*3*5*9)/(2^6)*Sum_{n >= 0} (1/8)^n/Product_{k = 1..n} (2^k - 1),
C = (2*3*5*9*17)/(2^10)*Sum_{n >= 0} (1/16)^n/Product_{k = 1..n} (2^k - 1), and so on. The sequence [2,3,5,9,17,...] is A000051. (End)
From Amiram Eldar, Mar 20 2022: (Start)
Equals sqrt(2) * exp(log(2)/24 + Pi^2/(12*log(2))) * Product_{k>=1} (1 - exp(-2*(2*k-1)*Pi^2/log(2))) (McIntosh, 1995).
Equals 1/A083864. (End)
Equals lim_{n->oo} A020696(2^n)/A006125(n+1) (Sándor, 2021). - Amiram Eldar, Jun 29 2022

A132326 Decimal expansion of Product_{k>=1} (1+1/10^k).

Original entry on oeis.org

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

Views

Author

Hieronymus Fischer, Aug 20 2007

Keywords

Comments

Half the constant A132325.

Examples

			1.1122345691370506321260780670944...
		

Crossrefs

Programs

  • Mathematica
    digits = 105; NProduct[1+1/3^k, {k, 0, Infinity}, NProductFactors -> 100, WorkingPrecision -> digits+3] // N[#, digits+3]& // RealDigits[#, 10, digits]& // First (* Jean-François Alcover, Feb 18 2014 *)
    N[QPochhammer[-1/10,1/10]] (* G. C. Greubel, Dec 01 2015 *)
  • PARI
    prodinf(k=1, 1+1/10^k) \\ Amiram Eldar, May 20 2023

Formula

Equals (1/2)*lim sup_{n->oo} Product_{0<=k<=floor(log_10(n))} (1+1/floor(n/10^k)).
Equals (1/2)*lim sup_{n->oo} A132271(n)/n^((1+log_10(n))/2).
Equals (1/2)*lim sup_{n->oo} A132272(n)/n^((log_10(n)-1)/2).
Equals exp(Sum_{n>0} 10^(-n)*Sum_{k|n} -(-1)^k/k) = exp(Sum_{n>0} A000593(n)/(n*10^n)).
Equals (1/2)*lim sup_{n->oo} A132271(n+1)/A132271(n).
Equals (-1/10; 1/10){infinity}, where (a;q){infinity} is the q-Pochhammer symbol. - G. C. Greubel, Dec 01 2015
Equals (sqrt(2)/2) * exp(log(10)/24 + Pi^2/(12*log(10))) * Product_{k>=1} (1 - exp(-2*(2*k-1)*Pi^2/log(10))) (McIntosh, 1995). - Amiram Eldar, May 20 2023

A132324 Decimal expansion of Product_{k>=1} (1+1/3^k).

Original entry on oeis.org

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

Views

Author

Hieronymus Fischer, Aug 20 2007

Keywords

Comments

Half the constant A132323.

Examples

			1.56493401856701153793884910...
		

Crossrefs

Programs

  • Mathematica
    digits = 105; NProduct[1+1/3^k, {k, 1, Infinity}, NProductFactors -> 100, WorkingPrecision -> digits+5] // N[#, digits+5]& // RealDigits[#, 10, digits]& // First (* Jean-François Alcover, Feb 18 2014 *)
    N[QPochhammer[-1/3,1/3]] (* G. C. Greubel, Dec 01 2015 *)

Formula

(1/2)*lim sup Product{k=0..floor(log_3(n))} (1+1/floor(n/3^k)) for n-->oo.
(1/2)*lim sup A132327(n)/A132027(n) for n-->oo.
(1/2)*lim sup A132327(n)/n^((1+log_3(n))/2) for n-->oo.
(1/2)*lim sup A132328(n)/n^((log_3(n)-1)/2) for n-->oo.
exp(Sum_{n>0} 3^(-n)*Sum_{k|n} -(-1)^k/k) = exp(Sum_{n>0} A000593(n)/(n*3^n)).
(1/2)*lim sup A132327(n+1)/A132327(n) = 1.56493401856701153793884910... for n-->oo.
Equals (-1/3; 1/3){infinity}, where (a;q){infinity} is the q-Pochhammer symbol. - G. C. Greubel, Dec 01 2015
From Amiram Eldar, Feb 19 2022: (Start)
Equals (sqrt(2)/2) * exp(log(3)/24 + Pi^2/(12*log(3))) * Product_{k>=1} (1 - exp(-2*(2*k-1)*Pi^2/log(3))) (McIntosh, 1995).
Equals Sum_{n>=0} 1/A027871(n). (End)

A264685 Expansion of Product_{k>=1} (1 + x^k)/(1 - 2*x^k).

Original entry on oeis.org

1, 3, 9, 24, 60, 141, 324, 717, 1560, 3330, 7020, 14622, 30225, 61998, 126522, 257007, 520326, 1050396, 2116116, 4255584, 8547330, 17149350, 34382295, 68889840, 137969466, 276220962, 552865365, 1106356314, 2213644548, 4428657402, 8859340926, 17721640698
Offset: 0

Views

Author

Vaclav Kotesovec, Nov 21 2015

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 40; CoefficientList[Series[Product[(1 + x^k)/(1 - 2*x^k), {k, 1, nmax}], {x, 0, nmax}], x]

Formula

a(n) ~ c * 2^n, where c = A079555 / A048651 = Product_{k>=1} (2^k+1)/(2^k-1) = 8.25598793577825006554414084943227312652...

A203480 a(n) = v(n+1)/v(n), where v = A203479.

Original entry on oeis.org

4, 80, 6336, 1901824, 2167925760, 9505110118400, 162323441859870720, 10902076148767162433536, 2898720791385603198124032000, 3064112360434477703904869089280000, 12909951234577776926559241120412860416000
Offset: 1

Views

Author

Clark Kimberling, Jan 02 2012

Keywords

Crossrefs

Programs

  • Magma
    [(&*[2^j +2^(n+1) -2: j in [1..n]]): n in [1..20]]; // G. C. Greubel, Aug 28 2023
    
  • Mathematica
    (* First program *)
    f[j_]:= 2^j - 1; z = 15;
    v[n_]:= Product[Product[f[k] + f[j], {j,k-1}], {k,2,n}]
    Table[v[n], {n,z}]               (* A203479 *)
    Table[v[n+1]/v[n], {n,z-1}]      (* A203480 *)
    Table[v[n+1]/(4*v[n]), {n,z-1}]  (* A203481 *)
    (* Second program *)
    Table[Product[2^(n+1) +2^k -2, {k,n}], {n,20}] (* G. C. Greubel, Aug 28 2023 *)
  • SageMath
    [product(2^j+2^(n+1)-2 for j in range(1,n+1)) for n in range(1,21)] # G. C. Greubel, Aug 28 2023

Formula

a(n) = Product_{k=1..n} (2^k + 2^(n+1) - 2). - G. C. Greubel, Aug 28 2023
a(n) ~ c * 2^(n*(n+1)), where c = 1/QPochhammer(1/2, 1/4) = A079555 = 2.3842310290313717... - Vaclav Kotesovec, Aug 09 2025

A203481 a(n) = v(n+1)/(4*v(n)), where v = A203479.

Original entry on oeis.org

1, 20, 1584, 475456, 541981440, 2376277529600, 40580860464967680, 2725519037191790608384, 724680197846400799531008000, 766028090108619425976217272320000, 3227487808644444231639810280103215104000
Offset: 1

Views

Author

Clark Kimberling, Jan 02 2012

Keywords

Crossrefs

Programs

  • Magma
    [(&*[2^j + 2^(n+1) - 2: j in [1..n]])/4: n in [1..20]]; // G. C. Greubel, Aug 28 2023
    
  • Mathematica
    (* First program *)
    f[j_]:= 2^j - 1; z = 15;
    v[n_]:= Product[Product[f[k] + f[j], {j, k-1}], {k, 2, n}]
    Table[v[n], {n, z}]               (* A203479 *)
    Table[v[n+1]/v[n], {n, z-1}]      (* A203480 *)
    Table[v[n+1]/(4*v[n]), {n, z-1}]  (* A203481 *)
    (* Second program *)
    Table[Product[2^(n+1) +2^k -2, {k,n}]/4, {n,20}] (* G. C. Greubel, Aug 28 2023 *)
  • SageMath
    [product(2^j+2^(n+1)-2 for j in range(1,n+1))/4 for n in range(1,21)] # G. C. Greubel, Aug 28 2023

Formula

a(n) = (1/4)*Product_{k=1..n} (2^k + 2^(n+1) - 2). - G. C. Greubel, Aug 28 2023
a(n) ~ c * 2^(n*(n+1)-2), where c = 1/QPochhammer(1/2, 1/4) = A079555 = 2.3842310290313717... - Vaclav Kotesovec, Aug 09 2025

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

Original entry on oeis.org

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

Views

Author

Ilya Gutkovskiy, Apr 28 2020

Keywords

Examples

			(1 + 1/2) * (1 - 1/2^2) * (1 + 1/2^3) * (1 - 1/2^4) * (1 + 1/2^5) * ... = 1.2107241303010591801361728561...
		

Crossrefs

Programs

  • Mathematica
    RealDigits[QPochhammer[-1/2, -1/2], 10, 111] [[1]]
    N[QPochhammer[-2, 1/4]*QPochhammer[1/4]/3, 120] (* Vaclav Kotesovec, Apr 28 2020 *)
  • PARI
    prodinf(k=1, 1 - 1/(-2)^k) \\ Michel Marcus, Apr 28 2020

Formula

Equals Product_{k>=1} (4^k - 1)*(4^k + 2)/4^(2*k).
Equals exp(-Sum_{k>=1} A000203(k)/(k*(-2)^k)).
Showing 1-10 of 18 results. Next