|
|
A260804
|
|
Number of ways to write n as n = x * y * z * t + x + y + z + t where 1 <= x <= y <= z <= t <= n.
|
|
6
|
|
|
0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 2, 0, 2, 1, 2, 0, 2, 1, 2, 1, 1, 1, 4, 1, 2, 1, 2, 1, 3, 0, 3, 2, 2, 1, 5, 0, 1, 2, 4, 1, 4, 0, 3, 3, 2, 1, 4, 2, 3, 2, 2, 0, 5, 1, 4, 1, 2, 3, 6, 1, 2, 2, 5, 1, 4, 0, 4, 3, 3, 1, 6, 1, 2, 4, 4, 2, 4, 1, 4, 2, 2, 1, 8, 2, 4, 2, 4, 2, 5, 1, 4, 2, 2, 3, 8, 1, 3, 4, 4, 0, 4, 1, 6, 4, 3, 0
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,14
|
|
COMMENTS
|
Is the set of n for which a(n)=0 finite?
Note that this set contains only numbers n of the form prime + 1. Indeed, if n-1>=4 is a composite number, then n = p*q + 1, p>=2, q>=2. If p <= q, then, for x=1, y=1, z = p-1, t = q-1, we have
x*y*z*t + x + y + z + t = 1*1*(p-1)*(q-1) + 1 + 1 + (p-1) + (q-1) = p*q + 1 = n; so a(n) >= 1. If p > q, then we set x=1, y=1, z = q-1, t = p-1, and again a(n) >= 1.
Note also that limsup_{n->infinity} (a(n)) = infinity. Indeed, this limit is realized, say, on n = primorials +1 (A002110), since, when m goes to infinity, the number of representations of n - 1 = A002110(m) of the form p*q tends to infinity. On primorials +1 > 2 we have a subsequence: 0,1,3,8,27,... .
A generalization. For k>=2, let b_k(n) be the number of ways to write n as n = x_1 * x_2 *...* x_k + x_1 + x_2 + ... + x_k, where 1 <= x_1 <= x_2 <= ... <= x_k <= n.
Then, for n >= k-1, b_k(n) = 0 yields that n - k + 3 is prime with similar other comments. In particular, only b_2(n) = 0 if and only if n+1 is 1 or prime (cf. A072670). (End)
|
|
LINKS
|
|
|
FORMULA
|
|
|
MATHEMATICA
|
xmax = 9; ymax = 21; zmax = 98; (* When extending data, terms where maxima for x, y or z are reached have to be checked one by one. *)
r[n_] := r[n] = Module[{r1, r2, r3, rn}, r1 = Reap[Do[rn = Reduce[n == x y z t + x + y + z + t && 1 <= x <= y <= z <= t <= n, t, Integers]; If[rn =!= False, Sow[{x, y, z, t} /. {ToRules[rn]}]], {x, 1, xmax}, {y, 1, ymax}, {z, 1, zmax}]]; If[r1 == {Null, {}} , {}, r2 = r1[[2, 1]]; r3 = Flatten[r2, 1]; If[Max[r3[[All, 1]]] == xmax, Print[ "xmax reached at n = ", n]]; If[Max[r3[[All, 2]]] == ymax, Print["ymax reached at n = ", n]]; If[Max[r3[[All, 3]]] == zmax, Print["zmax reached at n = ", n]]; r3]];
a[n_] := Length[r[n]];
Table[Print["a(", n, ") = ", a[n], " ", r[n]]; a[n], {n, 0, 109}] (* Jean-François Alcover, Nov 19 2018 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|