cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A100440 Number of distinct values of i*j + j*k + k*i with 1 <= i <= j <= k <= n.

This page as a plain text file.
%I A100440 #25 Mar 23 2021 15:50:53
%S A100440 1,4,10,20,33,50,68,93,123,154,193,233,276,325,377,434,500,568,643,
%T A100440 720,804,885,979,1068,1168,1274,1381,1495,1615,1746,1876,2005,2148,
%U A100440 2285,2437,2596,2748,2908,3077,3241,3425,3608,3796,3979,4181,4388,4585,4804,5015,5237
%N A100440 Number of distinct values of i*j + j*k + k*i with 1 <= i <= j <= k <= n.
%C A100440 a(n) <= A000292(n); a(n) = number of terms in n-th row of the triangle in A200741. - _Reinhard Zumkeller_, Nov 21 2011
%H A100440 David A. Corneth, <a href="/A100440/b100440.txt">Table of n, a(n) for n = 1..3000</a>
%p A100440 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;
%t A100440 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 *)
%o A100440 (Haskell)
%o A100440 a100440 = length . a200741_row  -- _Reinhard Zumkeller_, Nov 21 2011
%o A100440 (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
%o A100440 (Python)
%o A100440 from numba import njit
%o A100440 @njit()
%o A100440 def aupton(terms):
%o A100440   aset, alst = set(), []
%o A100440   for n in range(1, terms+1):
%o A100440     for i in range(1, n+1):
%o A100440       for j in range(i, n+1):
%o A100440         aset.add(i*j + j*n + n*i)
%o A100440     alst.append(len(aset))
%o A100440   return alst
%o A100440 print(aupton(50)) # _Michael S. Branicky_, Mar 23 2021
%Y A100440 Cf. A027430, A100439, A102533, A102534, A200737.
%K A100440 nonn
%O A100440 1,2
%A A100440 _N. J. A. Sloane_, Nov 21 2004
%E A100440 More terms from _Robert G. Wilson v_, Dec 14 2004