|
|
|
|
2, 3, 5, 11, 17, 29, 37, 41, 61, 67, 73, 97, 101, 107, 109, 113, 131, 137, 149, 157, 181, 193, 197, 211, 227, 233, 241, 257, 269, 277, 281, 307, 317, 331, 337, 347, 353, 373, 389, 397, 401, 409, 421, 457, 461, 467, 491, 521, 541, 547, 557, 569, 577, 587, 601, 613, 617, 641
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
a(2) = 3 does not divide A000292(a(2)).
A000292(a(n)) = a(n) * k, n <> 2, k > 1.
|
|
LINKS
|
|
|
EXAMPLE
|
|
|
MATHEMATICA
|
t[n_] := n(n+1)(n+2)/6; seqQ[n_] := If[PrimeQ[n], Module[{tn = t[n], ans = True}, Do[If[Divisible[tn, t[i]], ans=False; Break[]], {i, 2, n-1}]; ans], False]; Select[Range[650], seqQ] (* Amiram Eldar, Nov 22 2018 *)
|
|
PROG
|
(PARI) t(n) = n*(n+1)*(n+2)/6;
isok(n) = if (isprime(n), my(tn = t(n)); for (i=2, n-1, if ((tn % t(i)) == 0, return (0))); return (1), return (0)); \\ Michel Marcus, Sep 23 2018
|
|
CROSSREFS
|
Cf. A000292 (tetrahedral numbers), A317986 (corresponding tetrahedral numbers).
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|