OFFSET
0,3
LINKS
L. Seidel, Über eine einfache Entstehungsweise der Bernoullischen Zahlen und einiger verwandten Reihen, Sitzungsberichte der mathematisch-physikalischen Classe der königlich bayerischen Akademie der Wissenschaften zu München, Vol. 7 (1877), 157-187.
MAPLE
PROG
(Sage)
def SIB(dim, m): # Algorithm of L. Seidel (1877).
E = matrix(ZZ, dim)
def plow(n, dir):
if dir : # right to left backward
E[n, 0] = int(n == 1)
for k in range(n-1, -1, -1) :
E[k, n-k] = E[k+1, n-k-1] - E[k, n-k-1]
else: # left to right forward
E[0, n] = 0
for k in range(1, n+1, 1) :
E[k, n-k] = E[k-1, n-k+1] + E[k-1, n-k]
dir = [mod(n, m) == 1 for n in (0..dim-1)]
for n in (0..dim-1): plow(n, dir[n])
return [E[0, k] if dir[k] else E[k, 0] for k in range(dim)]
A281588_list = lambda len: SIB(len, 4)
A281588_list(26)
CROSSREFS
KEYWORD
sign
AUTHOR
Peter Luschny, Feb 01 2017
STATUS
approved