OFFSET
1,2
COMMENTS
LINKS
Rémy Sigrist, Table of n, a(n) for n = 1..1000
Rémy Sigrist, Illustration of the first terms
PROG
(PARI) a(n)=my(h=matrix(n, n)); fordiv(n, d, fordiv(n/d, dd, for(x=1, d, for(y=1, dd, h[x, y]=max(h[x, y], n/d/dd))))); return(sum(x=1, n, sum(y=1, n, h[x, y])))
CROSSREFS
KEYWORD
nonn
AUTHOR
Rémy Sigrist, Mar 12 2017
STATUS
approved