OFFSET
0,1
COMMENTS
a(0) = 9; binomial transform of sequence gives successive digits of sequence.
LINKS
N. J. A. Sloane, Transforms
MAPLE
revert := proc(n) local Linv, i, L ; L := convert(n, base, 10) ; Linv := [] ; for i from 1 to nops(L) do Linv := [op(Linv), op(-i, L)] ; od ; RETURN(Linv) ; end: A125004 := proc(nmax) local ldigs, T, diag, row ; T := array(1..nmax, 1..nmax) ; ldigs := [9, 1, 0] ; T[1, 1] := ldigs[1] ; for diag from 2 to nmax do T[diag, 1] := ldigs[diag] ; for row from diag-1 to 1 by -1 do T[row, diag-row+1] := T[row, diag-row]+T[row+1, diag-row] ; od ; if diag > 2 then ldigs := [op(ldigs), op(revert(T[1, diag])) ] ; fi ; od ; RETURN(T) ; end : nmax := 50 : T := A125004(nmax) : for i from 1 to nmax do printf("%d, ", T[1, i]) ; od ; # R. J. Mathar, Jan 10 2007
CROSSREFS
KEYWORD
nonn,easy,base
AUTHOR
N. J. A. Sloane, Jan 08 2007
EXTENSIONS
More terms from R. J. Mathar, Jan 10 2007
STATUS
approved