|
|
A100440
|
|
Number of distinct values of i*j + j*k + k*i with 1 <= i <= j <= k <= n.
|
|
7
|
|
|
1, 4, 10, 20, 33, 50, 68, 93, 123, 154, 193, 233, 276, 325, 377, 434, 500, 568, 643, 720, 804, 885, 979, 1068, 1168, 1274, 1381, 1495, 1615, 1746, 1876, 2005, 2148, 2285, 2437, 2596, 2748, 2908, 3077, 3241, 3425, 3608, 3796, 3979, 4181, 4388, 4585, 4804, 5015, 5237
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
|
|
LINKS
|
|
|
MAPLE
|
f:=proc(n) local i, j, k, t1; t1:={}; for i from 1 to n do for j from i to n do for k from j to n do t1:={op(t1), i*j+j*k+k*i}; od: od: od: t1:=convert(t1, list); nops(t1); end;
|
|
MATHEMATICA
|
f[n_] := Length[ Union[ Flatten[ Table[i*j + j*k + k*i, {i, n}, {j, i, n}, {k, j, n}] ]]]; Table[ f[n], {n, 48}] (* Robert G. Wilson v, Dec 14 2004 *)
|
|
PROG
|
(Haskell)
(PARI) first(n) = {my(v = vector(3*n^2, i, oo), res = vector(n)); forvec(x = vector(3, i, [1, n]), c = x[1]*x[2] + x[1]*x[3] + x[2]*x[3]; v[c] = min(x[3], v[c]); , 1); for(i = 1, #v, if(v[i] < oo, res[v[i]]++)); for(i = 2, #res, res[i] += res[i-1]); res } \\ David A. Corneth, Mar 23 2021
(Python)
from numba import njit
@njit()
def aupton(terms):
aset, alst = set(), []
for n in range(1, terms+1):
for i in range(1, n+1):
for j in range(i, n+1):
aset.add(i*j + j*n + n*i)
alst.append(len(aset))
return alst
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|