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.

A345487 Numbers that are the sum of seven squares in ten or more ways.

This page as a plain text file.
%I A345487 #14 May 10 2024 01:40:15
%S A345487 70,79,82,85,87,88,90,91,93,94,95,96,97,98,99,100,101,102,103,104,105,
%T A345487 106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,
%U A345487 123,124,125,126,127,128,129,130,131,132,133,134,135,136,137,138,139
%N A345487 Numbers that are the sum of seven squares in ten or more ways.
%H A345487 Sean A. Irvine, <a href="/A345487/b345487.txt">Table of n, a(n) for n = 1..1000</a>
%F A345487 Conjectures from _Chai Wah Wu_, Jan 05 2024: (Start)
%F A345487 a(n) = 2*a(n-1) - a(n-2) for n > 10.
%F A345487 G.f.: x*(-x^9 + x^8 - x^7 + x^6 - x^5 - x^4 - 6*x^2 - 61*x + 70)/(x - 1)^2. (End)
%e A345487 79 = 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 5^2 + 7^2
%e A345487    = 1^2 + 1^2 + 1^2 + 1^2 + 5^2 + 5^2 + 5^2
%e A345487    = 1^2 + 1^2 + 1^2 + 2^2 + 2^2 + 2^2 + 8^2
%e A345487    = 1^2 + 1^2 + 1^2 + 3^2 + 3^2 + 3^2 + 7^2
%e A345487    = 1^2 + 1^2 + 2^2 + 2^2 + 2^2 + 4^2 + 7^2
%e A345487    = 1^2 + 1^2 + 2^2 + 4^2 + 4^2 + 4^2 + 5^2
%e A345487    = 1^2 + 1^2 + 3^2 + 3^2 + 3^2 + 5^2 + 5^2
%e A345487    = 1^2 + 2^2 + 2^2 + 2^2 + 4^2 + 5^2 + 5^2
%e A345487    = 1^2 + 2^2 + 2^2 + 3^2 + 3^2 + 4^2 + 6^2
%e A345487    = 2^2 + 3^2 + 3^2 + 3^2 + 4^2 + 4^2 + 4^2
%e A345487    = 3^2 + 3^2 + 3^2 + 3^2 + 3^2 + 3^2 + 5^2
%e A345487 so 79 is a term.
%o A345487 (Python)
%o A345487 from itertools import combinations_with_replacement as cwr
%o A345487 from collections import defaultdict
%o A345487 keep = defaultdict(lambda: 0)
%o A345487 power_terms = [x**2 for x in range(1, 1000)]
%o A345487 for pos in cwr(power_terms, 7):
%o A345487     tot = sum(pos)
%o A345487     keep[tot] += 1
%o A345487     rets = sorted([k for k, v in keep.items() if v >= 10])
%o A345487     for x in range(len(rets)):
%o A345487         print(rets[x])
%Y A345487 Cf. A345477, A345486, A345497, A345506.
%K A345487 nonn
%O A345487 1,1
%A A345487 _David Consiglio, Jr._, Jun 20 2021