cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A003336 Numbers that are the sum of 2 positive 4th powers.

Original entry on oeis.org

2, 17, 32, 82, 97, 162, 257, 272, 337, 512, 626, 641, 706, 881, 1250, 1297, 1312, 1377, 1552, 1921, 2402, 2417, 2482, 2592, 2657, 3026, 3697, 4097, 4112, 4177, 4352, 4721, 4802, 5392, 6497, 6562, 6577, 6642, 6817, 7186, 7857, 8192, 8962, 10001, 10016, 10081, 10256, 10625
Offset: 1

Views

Author

Keywords

Comments

Numbers k such that k = x^4 + y^4 has a solution in positive integers x, y.
There are no squares in this sequence. - Altug Alkan, Apr 08 2016
As the order of addition doesn't matter we can assume terms are in nondecreasing order. - David A. Corneth, Aug 01 2020

Examples

			From _David A. Corneth_, Aug 01 2020: (Start)
16378801 is in the sequence as 16378801 = 43^4 + 60^4.
39126977 is in the sequence as 39126977 = 49^4 + 76^4.
71769617 is in the sequence as 71769617 = 19^4 + 92^4. (End)
		

Crossrefs

5906 is the first term in A060387 but not in this sequence. Cf. A020897.
Cf. A088687 (2 distinct 4th powers).
A###### (x, y): Numbers that are the form of x nonzero y-th powers.
Cf. A000404 (2, 2), A000408 (3, 2), A000414 (4, 2), A003072 (3, 3), A003325 (3, 2), A003327 (4, 3), A003328 (5, 3), A003329 (6, 3), A003330 (7, 3), A003331 (8, 3), A003332 (9, 3), A003333 (10, 3), A003334 (11, 3), A003335 (12, 3), A003336 (2, 4), A003337 (3, 4), A003338 (4, 4), A003339 (5, 4), A003340 (6, 4), A003341 (7, 4), A003342 (8, 4), A003343 (9, 4), A003344 (10, 4), A003345 (11, 4), A003346 (12, 4), A003347 (2, 5), A003348 (3, 5), A003349 (4, 5), A003350 (5, 5), A003351 (6, 5), A003352 (7, 5), A003353 (8, 5), A003354 (9, 5), A003355 (10, 5), A003356 (11, 5), A003357 (12, 5), A003358 (2, 6), A003359 (3, 6), A003360 (4, 6), A003361 (5, 6), A003362 (6, 6), A003363 (7, 6), A003364 (8, 6), A003365 (9, 6), A003366 (10, 6), A003367 (11, 6), A003368 (12, 6), A003369 (2, 7), A003370 (3, 7), A003371 (4, 7), A003372 (5, 7), A003373 (6, 7), A003374 (7, 7), A003375 (8, 7), A003376 (9, 7), A003377 (10, 7), A003378 (11, 7), A003379 (12, 7), A003380 (2, 8), A003381 (3, 8), A003382 (4, 8), A003383 (5, 8), A003384 (6, 8), A003385 (7, 8), A003387 (9, 8), A003388 (10, 8), A003389 (11, 8), A003390 (12, 8), A003391 (2, 9), A003392 (3, 9), A003393 (4, 9), A003394 (5, 9), A003395 (6, 9), A003396 (7, 9), A003397 (8, 9), A003398 (9, 9), A003399 (10, 9), A004800 (11, 9), A004801 (12, 9), A004802 (2, 10), A004803 (3, 10), A004804 (4, 10), A004805 (5, 10), A004806 (6, 10), A004807 (7, 10), A004808 (8, 10), A004809 (9, 10), A004810 (10, 10), A004811 (11, 10), A004812 (12, 10), A004813 (2, 11), A004814 (3, 11), A004815 (4, 11), A004816 (5, 11), A004817 (6, 11), A004818 (7, 11), A004819 (8, 11), A004820 (9, 11), A004821 (10, 11), A004822 (11, 11), A004823 (12, 11), A047700 (5, 2).
Cf. A000583 (4th powers).

Programs

  • Mathematica
    nn=12; Select[Union[Plus@@@(Tuples[Range[nn],{2}]^4)], # <= nn^4&] (* Harvey P. Dale, Dec 29 2010 *)
    Select[Range@ 11000, Length[PowersRepresentations[#, 2, 4] /. {0, } -> Nothing] > 0 &] (* _Michael De Vlieger, Apr 08 2016 *)
  • PARI
    list(lim)=my(v=List()); for(x=1, sqrtnint(lim\=1,4), for(y=1, min(sqrtnint(lim-x^4,4), x), listput(v, x^4+y^4))); Set(v) \\ Charles R Greathouse IV, Apr 24 2012; updated July 13 2024
    
  • PARI
    T=thueinit('x^4+1,1);
    is(n)=#thue(T,n)>0 && !issquare(n) \\ Charles R Greathouse IV, Feb 26 2017
    
  • Python
    def aupto(lim):
      p1 = set(i**4 for i in range(1, int(lim**.25)+2) if i**4 <= lim)
      p2 = set(a+b for a in p1 for b in p1 if a+b <= lim)
      return sorted(p2)
    print(aupto(10625)) # Michael S. Branicky, Mar 18 2021

Formula

{i: A216284(i) > 0}. - R. J. Mathar, Jun 04 2021

A209431 Numbers n such that x^4 + y^4 = n * z^4 is solvable in nonzero integers x,y,z with z > 1 and gcd(x,y,n) = 1.

Original entry on oeis.org

5906, 469297, 926977, 952577, 1127857, 1298257, 1347361, 1647377, 2455361, 3342817, 4928977, 5268706, 5519537, 8588161, 8879537, 9339361, 9391537, 9846017, 11414017, 14543026, 15547297, 16502722, 16657217, 16672322, 16830017, 19730162, 23672002, 25030097, 27681937, 27979762
Offset: 1

Views

Author

Jean-François Alcover, Mar 09 2012

Keywords

Comments

Values of z (1, 17, 41, 73, 89, ...) are elements of sequence A004625 (divisible only by primes congruent to 1 mod 8). The first composite z is 697 = 17*41: 41^4 + 822091^4 = 1935300738962*697^4.
Proof (after Ms. Adina Calvo) that values of z are divisible only by primes congruent to 1 mod 8: Let {x,y,z} be a nontrivial solution and p an odd prime divisor of z. Reducing the equation mod p, one gets in Z/pZ: x^4 + y^4 = 0 mod p. Hence (x*y^-1)^4 = -1, then x*y^-1 is an order-8 element of the multiplicative group (Z/pZ)*, which has p-1 elements. Therefore p is congruent to 1 mod 8.

Examples

			5906 is in the sequence because a^4 + b^4 = 5906*c^4 has the solution (a,b,c) = (25,149,17).
		

Crossrefs

Programs

  • Mathematica
    BiquadraticFreePart[n_] := Times @@ Power @@@ ({#[[1]], Mod[#[[2]], 4]} & /@ FactorInteger[n]); max = 10000; Sort[ Reap[Do[nz4 = x^4 + y^4; z4 = nz4/BiquadraticFreePart[nz4]; z = z4^(1/4); n = nz4/z4; If[z4 > 1 && IntegerQ[z] && GCD[x, y, z] == 1, Print[{n, x, y, z}]; Sow[n]], {x, 1, max}, {y, x, max}]][[2, 1]]]

Formula

Numbers in A060387 but not in A003336.

Extensions

Definition corrected by Hugo Pfoertner, Nov 08 2016

A085304 Least number of 4th powers required to represent n!.

Original entry on oeis.org

1, 1, 2, 6, 9, 10, 15, 15, 9, 10, 15, 6, 12, 12
Offset: 0

Views

Author

Labos Elemer, Jun 30 2003

Keywords

Examples

			n=6: 6!=720=625+81+14,length-of-solution=16>=a(6)
but 6!=720=2.256+13.16 seems shortest solution a(6)=15
after, see also A046046
n=7: 7!=5040=3.1296+4.256+8.16 so a(7)<=15 (uncertain);
n=8: a(8)<=9 because 8!=4.10000+1.256+4.16.
		

Crossrefs

Formula

"Shortest" solutions to n!=Sum[x(j)^4], j=1, .., m[n] with minimal value of m[n]: a(n)=Min{m[n]}. Per analogiam A084355.

Extensions

a(7)-a(11) from John W. Layman, Aug 13 2004
a(12) from Sean A. Irvine, Feb 11 2010
a(13) from Sean A. Irvine, Feb 15 2010
Showing 1-3 of 3 results.