|
%I #39 Aug 18 2015 12:39:29
%S 4,8,9,25,27,32,36,49,64,100,121,125,128,144,169,196,216,225,243,256,
%T 289,324,343,361,400,441,484,529,576,676,729,784,841,900,961,1000,
%U 1024,1089,1156,1225,1331,1369,1444,1521,1600,1681,1728,1764
%N Complement of A259444.
%H Charles R Greathouse IV, <a href="/A259183/b259183.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) ~ n^2. - _Charles R Greathouse IV_, Aug 18 2015
%o (PARI) main(s)=my(v=vector(s), v2=[], x, r, n); v[1]=2; for(n=2, s, v[n]=v[n-1]+1; until(x==1, for(r=1, n-1, for(s=1, n-1, if((v[r]^v[s])===v[n], v2=concat(v2, v[n]); print1(v[n], ", ");v[n]++; x=0; break(2), x=1))))); v2;
%K nonn
%O 1,1
%A _Anders Hellström_, Jul 19 2015 following suggestion from _N. J. A. Sloane_
|
