|
%I #13 Sep 25 2012 16:19:04
%S 8,12,13,16,17,18,20,24,25,29,31,32,33,34,35,36,40,41,43,44,45,48,50,
%T 52,53,54,57,58,59,61,63,64,65,68,72,73,74,76,80,81,85,88,89,90,91,96,
%U 97,98,99,100,104,106,108,109,112,113,116,117,121,124,125,127
%N Sum of two powerful numbers greater than 1.
%C A076871 is the primary sequence.
%e 8 = 4 + 4, 12 = 4 + 8, 13 = 4 + 9.
%o (PARI) isP(n)={
%o n>3 && vecmin(factor(n)[,2])>1
%o };
%o sumset(a,b)={
%o my(c=vector(#a*#b));
%o for(i=1,#a,
%o for(j=1,#b,
%o c[(i-1)*#b+j]=a[i]+b[j]
%o )
%o );
%o vecsort(c,,8)
%o };
%o upto(lim)={
%o my(v=select(isP, vector(floor(lim),i,i)));
%o select(n->n<=lim, sumset(v,v))
%o };
%K nonn
%O 1,1
%A _Charles R Greathouse IV_, Jul 06 2011
|
