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.

A038475 Sums of 3 distinct powers of 5.

This page as a plain text file.
%I A038475 #13 Apr 04 2025 14:35:49
%S A038475 31,131,151,155,631,651,655,751,755,775,3131,3151,3155,3251,3255,3275,
%T A038475 3751,3755,3775,3875,15631,15651,15655,15751,15755,15775,16251,16255,
%U A038475 16275,16375,18751,18755,18775,18875,19375,78131,78151,78155,78251,78255,78275,78751
%N A038475 Sums of 3 distinct powers of 5.
%H A038475 Amiram Eldar, <a href="/A038475/b038475.txt">Table of n, a(n) for n = 1..10000</a>
%t A038475 Sort[Plus @@@ Subsets[5^Range[0, 7], {3}]] (* _Amiram Eldar_, Jul 13 2022 *)
%o A038475 (Python)
%o A038475 from math import isqrt, comb
%o A038475 from sympy import integer_nthroot
%o A038475 def A038475(n): return 5**((r:=n-1-comb((m:=integer_nthroot(6*n,3)[0])+(t:=(n>comb(m+2,3)))+1,3))-comb((k:=isqrt(b:=r+1<<1))+(b>k*(k+1)),2))+5**((a:=isqrt(s:=n-comb(m-(t^1)+2,3)<<1))+((s<<2)>(a<<2)*(a+1)+1))+5**(m+t+1) # _Chai Wah Wu_, Apr 04 2025
%Y A038475 Base 5 interpretation of A038445.
%Y A038475 Cf. A000351, A038474, A038476, A038477.
%K A038475 nonn,easy
%O A038475 1,1
%A A038475 _Olivier Gérard_
%E A038475 Offset corrected by _Amiram Eldar_, Jul 13 2022