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.

A223035 Prime numbers whose digits squared sum to a square.

This page as a plain text file.
%I A223035 #21 May 13 2013 19:02:09
%S A223035 2,3,5,7,43,263,269,1153,1531,1933,2063,2069,2287,2609,3319,3391,3511,
%T A223035 3931,4003,4441,4801,4889,5113,5399,5939,6029,6067,6203,6469,6607,
%U A223035 8849,9133,9539,10111,10177,10513,10531,10771,11149,11213,11273,11321,11491,11503
%N A223035 Prime numbers whose digits squared sum to a square.
%H A223035 Christian N. K. Anderson, <a href="/A223035/b223035.txt">Table of n, a(n) for n = 1..10000</a>
%e A223035 269 is a prime number, and 2^2+6^2+9^2 = 121 = 11^2.
%t A223035 Select[Prime[Range[2000]], IntegerQ[Sqrt[Total[IntegerDigits[#]^2]]] &] (* _T. D. Noe_, Apr 05 2013 *)
%o A223035 (R) ssod<-function(i) sum(as.numeric(strsplit(as.character(i),"")[[1]])^2)
%o A223035 issquare<-function(x) as.integer(sqrt(x))==sqrt(x)
%o A223035 x=as.bigz(c()); i=2
%o A223035 while(length(x)<10000) {if(issquare(ssod(i))) x=c(x,i); i=nextprime(i)}
%Y A223035 Prime numbers from the sequence A175396.
%Y A223035 Cf. A000040, A061910, A107288.
%K A223035 nonn,base
%O A223035 1,1
%A A223035 _Kevin L. Schwartz_ and _Christian N. K. Anderson_, Apr 04 2013