|
%I #48 Oct 24 2023 07:04:30
%S 1,16,17,81,82,97,98,256,257,272,273,337,338,353,354,625,626,641,642,
%T 706,707,722,723,881,882,897,898,962,963,978,979,1296,1297,1312,1313,
%U 1377,1378,1393,1394,1552,1553,1568,1569,1633,1634,1649,1650,1921,1922
%N Sums of distinct nonzero 4th powers.
%C 5134240 is the largest positive integer not in this sequence. - _Jud McCranie_
%C If we tightened the sequence requirement so that the sum was of more than one 4th power, we would remove exactly 32 4th powers from the terms: row 4 of A332065 indicates which 4th powers would remain. After a(1) = 1, the next number in this sequence that is in the analogous sequences for cubes and squares is a(24) = 881 = A364637(4). - _Peter Munn_, Aug 01 2023
%D The Penguin Dictionary of Curious and Interesting Numbers, David Wells, entry 5134240.
%H Donovan Johnson, <a href="/A003999/b003999.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (2,-1).
%F For n > 4244664, a(n) = n + 889576. - _Charles R Greathouse IV_, Sep 02 2011
%p (1+x)*(1+x^16)*(1+x^81)*(1+x^256)*(1+x^625)*(1+x^1296)*(1+x^2401)*(1+x^4096)*(1+x^6561)*(1+x^10000)
%t max = 2000; f[x_] := Product[1 + x^(k^4), {k, 1, 10}]; Exponent[#, x]& /@ List @@ Normal[Series[f[x], {x, 0, max}]] // Rest (* _Jean-François Alcover_, Nov 09 2012, after _Charles R Greathouse IV_ *)
%o (PARI) upto(lim)={
%o lim\=1;
%o my(v=List(),P=prod(n=1,lim^(1/4),1+x^(n^4),1+O(x^(lim+1))));
%o for(n=1,lim,if(polcoeff(P,n),listput(v,n)));
%o Vec(v)
%o }; \\ _Charles R Greathouse IV_, Sep 02 2011
%Y Cf. A046039 (complement).
%Y Cf. A003995, A003997, A194768, A194769 (analogs for squares, cubes, 5th and 6th powers).
%Y Cf. A001661, A332065, A364637.
%Y A217844 is a subsequence.
%K nonn,easy
%O 1,2
%A _N. J. A. Sloane_
|
