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.

A224977 n^2 minus sum of digits of n^2.

This page as a plain text file.
%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