Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #10 Jun 08 2013 15:34:35
%S 1,49,169,36,1,1,2601,1089,1,8836,33489,44100,1,149769,128164,96721,1,
%T 156816,1225,40804,12321,831744,839056,1149184,1737124,3655744,407044,
%U 1890625,2208196,1089,1,1466521,6125625,2235025,2832489,1,3759721,6885376,8844676
%N Largest integer k in base n whose squared digits sum to sqrt(k).
%C Any d-digit number in base n meeting the criterion must also meet the condition d*(n-1)^2 < n^(d/2). Numerically, it can be shown this limits the candidate values to squares < 22*n^4. The larger values are statistically unlikely, and in fact the largest value of k in the first 1000 bases is ~9.96*n^4 in base 775.
%C a(n)=1 iff A226352(n)=1.
%H Christian N. K. Anderson, <a href="/A226353/b226353.txt">Table of n, a(n) for n = 2..1000</a>
%H Christian N. K. Anderson, <a href="/A226353/a226353.txt">Table of base, all solutions in base 10, and all solutions in base n</a> for bases 2 to 1000.
%e In base 8, the four solutions are the values {1,16,256,2601}, which are written as {1,20,400,5051} in base 8 and
%e sqrt(1) = 1 = 1^2
%e sqrt(16) = 4 = 2^2+0^2
%e sqrt(256) = 16 = 4^2+0^2+0^2
%e sqrt(2601)= 51 = 5^2+0^2+5^2+1^2
%o (R) inbase=function(n,b) { x=c(); while(n>=b) { x=c(n%%b,x); n=floor(n/b) }; c(n,x) }
%o for(n in 2:50) cat("Base",n,":",which(sapply((1:(4.7*n^2))^2,function(x) sum(inbase(x,n)^2)==sqrt(x)))^2,"\n")
%Y Cf. A226352, A226224.
%Y Cf. digital sums for digits at various powers: A007953, A003132, A055012, A055013, A055014, A055015.
%K nonn,base
%O 2,2
%A _Christian N. K. Anderson_ and _Kevin L. Schwartz_, Jun 04 2013