|
%I #6 Oct 15 2013 03:36:55
%S 10,13,16,17,20,23,24,26,29,30,32,35,38,41,44,47,50,57,64,65,68,70,73,
%T 80,82,83,88,89,94,95,105,110,117,118,125,133,140,143,148,149,154,157,
%U 167,176,177,178,182,191,192,200,203,208,209,215,221,227,236,242
%N Numbers in A134651 which are the sum of two terms from A001043 in a unique way.
%C Numbers that are the sum of two terms of A001043, but not sum of a different pair of two terms.
%e Recall that A001043 begins with 5,8,12,18,24,30,36,42,52,60,68,78,84,90, ...
%e So 48 is not in the sequence because 48 = 18+30 = 24+24.
%e But 88 is in the sequence as it is only = 52+36.
%e And 89 is there too because = 84+5.
%o (PARI) issum(i, vss) = {nb = 0; for (j = 1, #vss, if (vss[j] > i, break); for (k = 1, j, sv = vss[j] + vss[k]; if (sv == i, nb++); if (sv > i, break););); return (nb == 1);}
%o lista(nn) = {vec = vector(nn, i, i); vss = select(i->((precprime((i-1)/2) + nextprime(i/2) == i) && (i>2)), vec); for (i = 1, nn, if (issum(i, vss), print1(i, ", ")););} \\ _Michel Marcus_, Oct 14 2013
%Y Cf. A001043, A134651, A135045.
%K nonn
%O 1,1
%A Zak Seidov, Feb 15 2008
%E Example, corrected and extended by _Michel Marcus_, Oct 14 2013
|
