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.
%I A334557 #22 Jul 17 2021 04:29:33 %S A334557 1,13,41,61,125,212,281,613,1156,1424,2225,3232,3316,4113,11125,11281, %T A334557 11525,12816,14913,16317,16441,19125,21284,21625,24128,25216,27521, %U A334557 31525,53125,56116,61321,65161,71325,82116,82217,83521,84313,111812,113125,113625,115336,115681,117125,118372 %N A334557 Numbers m such that m = p^2 + k^2, with p > 0, where p = A007954(m) = the product of digits of m. %H A334557 Giovanni Resta, <a href="/A334557/b334557.txt">Table of n, a(n) for n = 1..10000</a> %e A334557 13 is a term as p = 1*3 = 3 and 13 = 3^2 + 2^2. %e A334557 281 is a term as p = 2*8*1 = 16 and 281 = 16^2 + 5^2. %e A334557 118372 is a term as p = 1*1*8*3*7*2 = 336 and 118372 = 336^2 + 74^2. %o A334557 (PARI) isok(m) = my(p=vecprod(digits(m))); p && issquare(m - p^2); \\ _Michel Marcus_, May 06 2020 %Y A334557 Cf. A007954, A000404, A078134, A334542, A334558. %K A334557 nonn,base %O A334557 1,2 %A A334557 _Scott R. Shannon_, May 06 2020