%I #7 Feb 06 2019 02:20:49
%S 2,5,7,17,23,29,71,101,113,127,179,191,197,223,239,251,257,281,283,
%T 311,331,337,359,401,449,463,491,521,541,569,587,617,647,659,673,701,
%U 733,739,743,757,761,827,853,877,881,883,907,911,929,953,997,1051,1063,1093
%N Primes not of the form prime + k^2*2^(k-1).
%p From _R. J. Mathar_, Apr 27 2010: (Start)
%p A014477 := proc(n) (n+1)^2*2^n ; end proc:
%p isA176676 := proc(p) local n,a77 ; if not isprime(p) then return false; end if; for n from 1 do a77 := A014477(n-1) ; if isprime(p-a77) then return false; elif p-a77<2 then return true; end if; end do: end proc:
%p for n from 1 to 300 do p := ithprime(n) ; if isA176676(p) then printf("%d,",p) ; end if; end do: (End)
%Y Cf. A014477.
%K nonn
%O 1,1
%A _Juri-Stepan Gerasimov_, Apr 23 2010
%E Corrected (149 replaced by 179, 269 removed, 643 removed, 743 and 911 and 1053 inserted) by _R. J. Mathar_, Apr 27 2010
|