%I #10 Oct 17 2012 08:55:11
%S 0,4,22,64,144,269,461,720,1072,1522,2092,2774,3626,4614,5776,7126,
%T 8694,10445,12461,14684,17204,19997,23077,26412,30156,34206,38600,
%U 43352,48532,54042,60072,66458,73338,80664,88450,96710,105638,114999
%N Sum {k + l*m <= n} (k + l*m), with 0 <= k,l,m <= n.
%H R. J. Mathar, <a href="/A106846/b106846.txt">Table of n, a(n) for n = 0..1000</a>
%p A106846 := proc(n)
%p local a,k,l,m ;
%p a := 0 ;
%p for k from 0 to n do
%p for l from 0 to n do
%p if l = 0 then
%p m := n ;
%p else
%p m := floor((n-k)/l) ;
%p end if;
%p if m >=0 then
%p m := min(m,n) ;
%p a := a+(m+1)*k+l*m*(m+1)/2 ;
%p end if;
%p end do:
%p end do:
%p a ;
%p end proc: # _R. J. Mathar_, Oct 17 2012
%Y Cf. A106633, A106634, A106847.
%K nonn
%O 0,2
%A _Ralf Stephan_, May 06 2005
|