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.

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

This page as a plain text file.
%I A346808 #10 May 10 2024 02:17:58
%S A346808 61,64,66,67,69,70,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,
%T A346808 89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,
%U A346808 109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124
%N A346808 Numbers that are the sum of ten squares in ten or more ways.
%H A346808 Sean A. Irvine, <a href="/A346808/b346808.txt">Table of n, a(n) for n = 1..10000</a>
%e A346808 64 = 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 2^2 + 2^2 + 7^2
%e A346808    = 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 2^2 + 3^2 + 3^2 + 6^2
%e A346808    = 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 2^2 + 2^2 + 5^2 + 5^2
%e A346808    = 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 4^2 + 4^2 + 5^2
%e A346808    = 1^2 + 1^2 + 1^2 + 2^2 + 2^2 + 2^2 + 2^2 + 2^2 + 4^2 + 5^2
%e A346808    = 1^2 + 1^2 + 1^2 + 1^2 + 2^2 + 2^2 + 3^2 + 3^2 + 3^2 + 5^2
%e A346808    = 1^2 + 1^2 + 1^2 + 1^2 + 2^2 + 2^2 + 2^2 + 4^2 + 4^2 + 4^2
%e A346808    = 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 3^2 + 3^2 + 3^2 + 4^2 + 4^2
%e A346808    = 2^2 + 2^2 + 2^2 + 2^2 + 2^2 + 2^2 + 2^2 + 2^2 + 4^2 + 4^2
%e A346808    = 1^2 + 2^2 + 2^2 + 2^2 + 2^2 + 2^2 + 3^2 + 3^2 + 3^2 + 4^2
%e A346808    = 1^2 + 1^2 + 2^2 + 2^2 + 3^2 + 3^2 + 3^2 + 3^2 + 3^2 + 3^2
%e A346808 so 64 is a term.
%o A346808 (Python)
%o A346808 from itertools import combinations_with_replacement as cwr
%o A346808 from collections import defaultdict
%o A346808 keep = defaultdict(lambda: 0)
%o A346808 power_terms = [x**2 for x in range(1, 1000)]
%o A346808 for pos in cwr(power_terms, 10):
%o A346808     tot = sum(pos)
%o A346808     keep[tot] += 1
%o A346808     rets = sorted([k for k, v in keep.items() if v >= 10])
%o A346808     for x in range(len(rets)):
%o A346808         print(rets[x])
%Y A346808 Cf. A345558, A346803. Subsequence of A346807.
%K A346808 nonn
%O A346808 1,1
%A A346808 _David Consiglio, Jr._, Aug 04 2021