|
|
A027637
|
|
a(n) = Product_{i=1..n} (4^i - 1).
|
|
20
|
|
|
1, 3, 45, 2835, 722925, 739552275, 3028466566125, 49615367752825875, 3251543125681443718125, 852369269595510700600441875, 893773106866112632882108339078125, 3748755223447856814435325652920396921875
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
Given probability p = 1/4^n that an outcome will occur at the n-th stage of an infinite process, then starting at n=1, 1-a(n)/A053763(n+1) is the probability that the outcome has occurred up to and including the n-th iteration. The limiting ratio is 1-A100221 ~ 0.3114625. - Bob Selcoe, Mar 01 2016
|
|
LINKS
|
|
|
FORMULA
|
a(n) ~ c * 2^(n*(n+1)), where c = Product_{k>=1} (1-1/4^k) = A100221 = 0.688537537120339715456514357293508184675549819378... . - Vaclav Kotesovec, Nov 21 2015
a(n) = 4^(binomial(n+1,2))*(1/4;1/4)_{n} = (4; 4)_{n}, where (a;q)_{n} is the q-Pochhammer symbol. - G. C. Greubel, Dec 24 2015
G.f.: Sum_{n>=0} 4^(n*(n+1)/2)*x^n / Product_{k=0..n} (1 + 4^k*x). - Ilya Gutkovskiy, May 22 2017
|
|
MAPLE
|
mul( 4^i-1, i=1..n) ;
end proc:
|
|
MATHEMATICA
|
a[ n_] := If[ n < 0, 0, Product[ (q^(2 k) - 1) / (q - 1), {k, n}] /. q -> 2]; (* Michael Somos, Sep 12 2014 *)
|
|
PROG
|
(Magma) [1] cat [&*[4^k-1: k in [1..n]]: n in [1..11]]; // Vincenzo Librandi, Dec 24 2015
(SageMath)
from sage.combinat.q_analogues import q_pochhammer
def A027637(n): return (-1)^n*q_pochhammer(n, 4, 4)
|
|
CROSSREFS
|
Sequences of the form q-Pochhammer(n, q, q): A005329 (q=2), A027871 (q=3), this sequence (q=4), A027872 (q=5), A027873 (q=6), A027875 (q=7), A027876 (q=8), A027877 (q=9), A027878 (q=10), A027879 (q=11), A027880 (q=12).
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|