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.
a(4) = 11: the 11 compositions of this type of 4 into 4 parts being (4,0,0,0); (3,1,0,0); (3,0,1,0); (3,0,0,1); (2,2,0,0); (2,0,2,0); (2,0,0,2); (2,1,1,0); (2,1,0,1); (2,0,1,1); (1,1,1,1)
b:= proc(n,i,m) option remember; if n<0 then 0 elif n=0 then 1 elif i=1 then `if`(n<=m, 1, 0) else add(b(n-k, i-1, m), k=0..m) fi end: A:= (n,k)-> add(b(n-m, k-1, m), m=ceil(n/k)..n): seq(A(n,n), n=1..30); # Alois P. Heinz, Jun 14 2009
b[n_, i_, m_] := b[n, i, m] = Which[n<0, 0, n==0, 1, i==1, If[n <= m, 1, 0], True, Sum[b[n-k, i-1, m], {k, 0, m}]]; A[n_, k_] := Sum[b[n-m, k-1, m], {m, Ceiling[n/k], n}]; Table[A[n, n], {n, 1, 30}] (* Jean-François Alcover, Jul 15 2015, after Alois P. Heinz *)
N=120;v=vector(N,i,0);for(d=1,N,A=matrix(N,N,i,j,0);A[1,1]=1; for(i=1,N-1,for(j=0,N-1,s=0;for(k=0,min(j,d), s+=A[i,j-k+1]);A[i+1,j+1]=s)); for(i=d,N,v[i]+=A[i,i-d+1]));for(i=1,N,print1(v[i]", ")) \\ Robert Gerbicz, Apr 06 2011
Comments