OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..700
EXAMPLE
Nondivisors of 22 are {3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21} and their sum of cubes 3^3 + 4^3 + 5^3 + 6^3 + 7^3 +8^3 + 9^3 + 10^3 + 12^3 + 13^3 + 14^3 + 15^3 + 16^3 + 17^3 + 18^3 + 19^3 + 20^3 + 21^3 = 52021 is prime.
MATHEMATICA
Q1[n_]:=(n > 0) && PrimeQ[n]; Select[Range[3502], Q1[(#(#+1)/2)^2 - DivisorSigma[3, #]] &]
Select[Range[4000], PrimeQ[Total[Complement[Range[#], Divisors[#]]^3]]&] (* Harvey P. Dale, May 02 2026 *)
PROG
(PARI) isok(m) = isprime(sum(k=1, m-1, if (m%k, k^3))); \\ Michel Marcus, Feb 09 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Jan 27 2025
STATUS
approved
