OFFSET
1,4
MAPLE
a[1] := 1:b[1] := 1:for n from 2 to 100 do a[n] := floor(evalf(-1/2+1/2*sqrt(1+8*(n-1)*a[n-1])+1/10^50, 300)):b[n] := a[n]:a[n] := a[n]*(a[n]+1)/2:od:seq(b[l], l=1..100);
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Oct 26 2002
EXTENSIONS
More terms from Sascha Kurz, Jan 26 2003
STATUS
approved