login
A249329
First row of spectral array W(log(Pi)).
0
1, 7, 8, 55, 62, 435, 497, 3440, 3937, 27208, 31145, 215199, 246344, 1702099, 1948443, 13462620
OFFSET
1,2
COMMENTS
log(Pi) = 1.144729885849400174143427351353058711647294812915311571513623...
The sequence is generated from the Beatty sequence (A059561) and from the complement of the Beatty sequence (A059562) for log(Pi).
LINKS
A. Fraenkel and C. Kimberling, Generalized Wythoff arrays, shuffles and interspersions, Discrete Mathematics 126 (1994) 137-149.
PROG
(PARI)
\\ Row i of the generalized Wythoff array W(h),
\\ where h is an irrational number between 1 and 2,
\\ and m is the number of terms in the vectors b and c.
row(h, i, m) = {
if(h<=1 || h>=2, print("Invalid value for h"); return);
my(
b=vector(m, n, floor(n*h)), \\ Beatty sequence for h
c=vector(m, n, floor(n*h/(h-1))), \\ Complement of b
w=[b[b[i]], c[b[i]]],
j=3
);
while(1,
if(j%2==1,
if(w[j-1]<=#b, w=concat(w, b[w[j-1]]), return(w))
,
if(w[j-2]<=#c, w=concat(w, c[w[j-2]]), return(w))
);
j++
)
}
allocatemem(10^9)
default(realprecision, 100)
row(log(Pi), 1, 10^7)
CROSSREFS
Sequence in context: A201919 A249310 A094556 * A041023 A041108 A248276
KEYWORD
nonn,more
AUTHOR
Colin Barker, Dec 03 2014
STATUS
approved