|
%I #28 Sep 08 2022 08:46:12
%S 0,5,7,8,12,15,16,18,19,20,21,22,26,29,31,33,35,36,39,40,43,44,46,47,
%T 50,51,53,54,56,57,59,60,61,64,65,68,69,70,71,72,73,75,76,77,78,79,82,
%U 83,84,85,87,89,92,94,96,97,98,99,100,101,102,103,104,106
%N Numbers n>=0 such that (n+1)^3 - n^3 = 3*n^2+3*n+1 is not prime.
%C Complement of A111251.
%C Includes all members of A047383 except 1. - _Robert Israel_, May 12 2015
%H Charles R Greathouse IV, <a href="/A257772/b257772.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) ~ n. - _Charles R Greathouse IV_, May 22 2015
%e 5 is a term since (5+1)^3 - 5^3 = 91 = 13*7 is not prime.
%p remove(t -> isprime((t+1)^3-t^3), [$0..300]); # _Robert Israel_, May 12 2015
%t Select[Range[0, 200], ! PrimeQ[(#+1)^3 - #^3] &] (* _Giovanni Resta_, May 08 2015 *)
%o (UBASIC)
%o 10 print 0
%o 20 for n=1 to 200
%o 30 s = (n+1)^3 - n^3
%o 40 if prmdiv(s)<>s then print n
%o 50 next n
%o (Magma) [n: n in [0..120] | not IsPrime(3*n^2+3*n+1)]; // _Vincenzo Librandi_, May 13 2015
%o (PARI) for(n=0,100,if(!isprime(3*n^2+3*n+1),print1(n,", "))) \\ _Derek Orr_, May 19 2015
%Y Cf. A003215, A047383, A111251.
%K nonn,easy
%O 1,2
%A _Miguel Angel Velilla Mula_, May 08 2015
|
