login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

Sum of two powerful numbers greater than 1.
1

%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