|
%I #10 Aug 04 2018 14:34:52
%S 5,127,163,179,191,193,223,239,251,269,311,337,389,419,431,457,491,
%T 547,557,569,599,613,653,659,673,683,719,739,787,821,839,853,883,911,
%U 929,953,967,977,1117,1123,1201,1229,1249,1283,1289,1297,1303,1327,1381,1409,1423,1439,1451,1471,1481,1499
%N Primes p which have a unique representation as p = 2^i + q^j where i >= 0, j >= 1, q = odd prime.
%e 5 = 2+3 belongs to the sequence, but 23 = 2^2+19^1 = 2^4+7^1 does not.
%t maxp = 1500; Clear[cnt]; cnt[_] = 0;
%t pp = Prime[Range[PrimePi[maxp]]];
%t Do[p = 2^i + q^j; If[p <= maxp && PrimeQ[p], cnt[p] = cnt[p] + 1], {i, 0, Log[2, maxp] // Ceiling}, {j, 1, Log[3, maxp] // Ceiling}, {q, Rest[pp]}
%t ];
%t Select[pp, cnt[#] == 1&] (* _Jean-François Alcover_, Aug 04 2018 *)
%Y Subsequence of A115232. Cf. A115230, A115231.
%Y Cf. A000079, A061345.
%K nonn
%O 1,1
%A _Reinhard Zumkeller_, Jan 17 2006
%E Recomputed (based on recomputation of A115230) by _R. J. Mathar_ and _Reinhard Zumkeller_, Apr 29 2010.
%E Edited by _N. J. A. Sloane_, Apr 30 2010
%E Data corrected by _Jean-François Alcover_, Aug 04 2018
|
