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 A224977 #19 May 10 2017 11:20:30 %S A224977 0,0,0,0,9,18,27,36,54,72,99,117,135,153,180,216,243,270,315,351,396, %T A224977 432,468,513,558,612,657,711,765,828,891,945,1017,1071,1143,1215,1278, %U A224977 1350,1431,1512,1593,1665,1746,1827,1917,2016,2106,2196,2295,2394,2493 %N A224977 n^2 minus sum of digits of n^2. %C A224977 a(n) mod 9 = 0 for all n. %H A224977 Christian N. K. Anderson, <a href="/A224977/b224977.txt">Table of n, a(n) for n = 0..9999</a> %H A224977 Christian N. K. Anderson, <a href="/A224977/a224977.jpg">Ulam spiral</a> of a(n) in red, with other multiples of 9 in gray. %F A224977 a(n) = A000290(n) - A004159(n). %e A224977 a(0) = 0^2 - 0 = 0. %e A224977 a(12) = 144 - (1+4+4) = 135. %t A224977 Table[n^2 - Sum[DigitCount[n^2][[i]]i, {i, 9}], {n, 50}] (* _Alonso del Arte_, Apr 21 2013 *) %t A224977 #^2-Total[IntegerDigits[#^2]]&/@Range[0,50] (* _Harvey P. Dale_, May 10 2017 *) %o A224977 (R) library(gmp); digsum<-function(x) sum(as.numeric(unlist(strsplit(as.character(x),split="")))) %o A224977 ans=rep(0,100); n=0 %o A224977 while(n<=100) ans[(n=n+1)]=n^2-digsum(n^2); ans %Y A224977 Cf. A067552. %Y A224977 Derived from A000290 and A004159. %K A224977 nonn,base,easy %O A224977 0,5 %A A224977 _Kevin L. Schwartz_ and _Christian N. K. Anderson_, Apr 21 2013