OFFSET
0,2
COMMENTS
Number of ordered triples [k,l,m] with n = k+l*m and k, l, m all in the range [0..n].
From R. J. Mathar, Jun 30 2013: (Start)
A010766 is the following array A read by antidiagonals:
1, 1, 1, 1, 1, 1, ...
2, 1, 1, 1, 1, 1, ...
3, 2, 1, 1, 1, 1, ...
4, 2, 2, 1, 1, 1, ...
5, 3, 2, 2, 1, 1, ...
6, 3, 2, 2, 2, 1, ...
and apparently a(n) is the hook sum Sum_{k=0..n} A(n,k) + Sum_{r=0..n-1} A(r,n). (End)
LINKS
R. J. Mathar, Table of n, a(n) for n = 0..1000
FORMULA
From Ridouane Oudra, Apr 22 2024: (Start)
a(n) = 2*n + 1 + Sum_{k=1..n} floor(n/k);
a(n) = 2*n + 1 + Sum_{k=1..n} tau(k);
EXAMPLE
0+1*2 = 0+2*1 = 1+1*1 = 2+0*0 = 2+0*1 = 2+0*2 = 2+1*0 = 2+2*0 = 2, so a(2)=8.
a(3)=12: 3+0*0, 3+0*m (6), 2+1*1, 1+2*1 (2), 0+3*1 (2).
MAPLE
A106633 := proc(n)
local a, k, l, m ;
a := 0 ;
for k from 0 to n do
for l from 0 to n do
if l = 0 then
if k = n then
a := a+n+1 ;
end if;
else
m := (n-k)/l ;
if type(m, 'integer') then
a := a+1 ;
end if;
end if;
end do:
end do:
a ;
end proc: # R. J. Mathar, Oct 17 2012
MATHEMATICA
A106633[n_] := Module[{a, m}, a = 0; Do[If[l == 0, If[k == n, a = a + n + 1], m = (n - k)/l; If[IntegerQ[m], a = a + 1]], {k, 0, n}, {l, 0, n}]; a];
PROG
(PARI) list(n)={
my(v=vector(n), t);
for(i=2, n, for(j=1, min(n\i, i-1), v[i*j]+=2));
for(i=1, sqrtint(n), v[i^2]++);
concat(1, vector(n, k, 2*k+1+t+=v[k]))
}; \\ Charles R Greathouse IV, Oct 17 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Ralf Stephan, May 06 2005
EXTENSIONS
Definition clarified by N. J. A. Sloane, Jul 07 2012
STATUS
approved