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-4 of 4 results.

A002377 Least number of 4th powers needed to represent n.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 5, 1, 2, 3
Offset: 1

Views

Author

Keywords

Comments

No terms are greater than 19, see A002804. - Charles R Greathouse IV, Aug 01 2013
Seven values of n need the maximum of 19 fourth powers. These form the arithmetic progression {79, 159, 239, 319, 399, 479, 559} each term being congruent to 79 mod 80. For n < 625 the available fourth powers are congruent to 1 or 16 mod 80, requiring 4*16 + 15*1 to sum to 79. However, 625 = 5^4 is congruent to 65 and 1*65 + 14*1 = 79. So for n > 625 and congruent to 79, only 15 fourth powers are needed to satisfy the mod 80 arithmetic. - Peter Munn, Apr 12 2017

References

  • D. H. Lehmer, Guide to Tables in the Theory of Numbers. Bulletin No. 105, National Research Council, Washington, DC, 1941, p. 82.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    Cnt4[n_] := Module[{k = 1}, While[Length[PowersRepresentations[n, k, 4]] == 0, k++]; k]; Array[Cnt4, 100] (* T. D. Noe, Apr 01 2011 *)
    seq[n_] := Module[{v = Table[0, {n}], s, p}, s = Sum[x^(k^4), {k, 1, n^(1/4)}] + O[x]^(n+1); p=1; For[k=1, k <= 19, k++, p *= s; For[i=1, i <= n, i++, If[v[[i]]==0 && Coefficient[p, x, i] != 0, v[[i]] = k]]]; v];
    seq[100] (* Jean-François Alcover, Sep 28 2019, after Andrew Howroyd *)
  • PARI
    seq(n)={my(v=vector(n), s=sum(k=1, sqrtint(sqrtint(n)), x^(k^4)) + O(x*x^n), p=1); for(k=1, 19, p*=s; for(i=1, n, if(!v[i] && polcoeff(p,i), v[i]=k))); v} \\ Andrew Howroyd, Jul 06 2018
    
  • Python
    from itertools import count
    from sympy.solvers.diophantine.diophantine import power_representation
    def A002377(n):
        if n == 1: return 1
        for k in count(1):
            try:
                next(power_representation(n,4,k))
            except:
                continue
            return k # Chai Wah Wu, Jun 25 2024

Extensions

More terms from Arlin Anderson (starship1(AT)gmail.com)

A186682 Total number of positive integers below 10^n requiring 18 positive biquadrates in their representation as sum of biquadrates.

Original entry on oeis.org

0, 2, 19, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24
Offset: 1

Views

Author

Martin Renner, Feb 25 2011

Keywords

Comments

A114322(n) + A186649(n) + A186651(n) + A186653(n) + A186655(n) + A186657(n) + A186659(n) + A186661(n) + A186663(n) + A186665(n) + A186667(n) + A186669(n) + A186671(n) + A186673(n) + A186675(n) + A186677(n) + A186680(n) + a(n) + A186684(n) = A002283(n)
a(n) = 24 for n >= 4. - Nathaniel Johnston, May 09 2011
Continued fraction expansion of (185-sqrt(145))/355. - Bruno Berselli, May 10 2011

Crossrefs

Cf. A046049.

Programs

  • Mathematica
    PadRight[{0, 2, 19}, 100, 24] (* Paolo Xausa, Jul 30 2024 *)

Formula

G.f.: x^2*(2+17*x+5*x^2)/(1-x). - Bruno Berselli, May 10 2011

Extensions

a(5)-a(6) from Lars Blomberg, May 08 2011
a(7) from Charles R Greathouse IV, May 08 2011
Terms after a(7) from Nathaniel Johnston, May 09 2011

A186683 Total number of n-digit numbers requiring 18 positive biquadrates in their representation as sum of biquadrates.

Original entry on oeis.org

0, 2, 17, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 1

Views

Author

Martin Renner, Feb 25 2011

Keywords

Comments

A161905(n) + A186650(n) + A186652(n) + A186654(n) + A186656(n) + A186658(n) + A186660(n) + A186662(n) + A186664(n) + A186666(n) + A186668(n) + A186670(n) + A186672(n) + A186674(n) + A186676(n) + A186678(n) + A186681(n) + a(n) + A186685(n) = A052268(n)
a(n) = 0 for n >= 5. - Nathaniel Johnston, May 09 2011

Crossrefs

Programs

  • Mathematica
    PadRight[{0, 2, 17, 5}, 100] (* Paolo Xausa, Jul 30 2024 *)

Formula

a(n) = A186682(n) - A186682(n-1).

A046050 Sum of 19 but no fewer nonzero fourth powers.

Original entry on oeis.org

79, 159, 239, 319, 399, 479, 559
Offset: 1

Views

Author

Keywords

Comments

Dickson noted that this sequence is complete to 4100. Deshouillers, Hennecart and Landreau showed that this sequence is complete up to 10^245, and Kawada, Wooley and Deshouillers showed that it is complete beyond 10^220.

References

  • J.-M. Deshouillers, K. Kawada and T. D. Wooley, On sums of sixteen biquadrates, Mem. Soc. Math. Fr. 100 (2005), pp. 120.

Crossrefs

Programs

  • Mathematica
    Select[Range[1000], (pr = PowersRepresentations[#, 19, 4]; test = pr != {} && FreeQ[pr, r_List /; (Times @@ r) == 0]; If[test, Print[#]]; test) &] (* Jean-François Alcover, Oct 30 2012 *)
  • PARI
    is(n)=n%80==79 && n<600 && n>0 \\ Charles R Greathouse IV, Jan 23 2014

Extensions

More terms from Arlin Anderson (starship1(AT)gmail.com). Jud McCranie remarks that probably all terms are shown.
Showing 1-4 of 4 results.