login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A048653 Numbers n such that the decimal digits of n^2 can be partitioned into two or more squares. 3

%I

%S 7,12,13,19,21,35,37,38,41,44,57,65,70,107,108,112,119,120,121,125,

%T 129,130,190,191,204,205,209,210,212,223,253,285,305,306,315,342,343,

%U 345,350,369,370,379,380,408,410,413,440,441,475,487,501,538,570,642,650

%N Numbers n such that the decimal digits of n^2 can be partitioned into two or more squares.

%H Reinhard Zumkeller, <a href="/A048653/b048653.txt">Table of n, a(n) for n = 1..10000</a>

%e 12 is present because 12^2=144 can be partitioned into three squares 1, 4 and 4.

%e 108^2 = 11664 = 1_16_64, 120^2 = 14400 = 1_4_400, so 108 and 120 are in the sequence.

%t (* This non-optimized program is not suitable to compute a large number of terms. *) split[digits_, pos_] := Module[{pos2}, pos2 = Transpose[{Join[ {1}, Most[pos+1]], pos}]; FromDigits[Take[digits, {#[[1]], #[[2]]}]]& /@ pos2]; sel[n_] := Module[{digits, ip, ip2, accu, nn}, digits = IntegerDigits[n^2]; ip = IntegerPartitions[Length[digits]]; ip2 = Flatten[ Permutations /@ ip, 1]; accu = Accumulate /@ ip2; nn = split[ digits, #]& /@ accu; SelectFirst[nn, Length[#]>1 && Flatten[ IntegerDigits[#] ] == digits && AllTrue[#, #>0 && IntegerQ[Sqrt[#]]&]&] ]; k = 1; Reap[Do[If[(s = sel[n]) != {}, Print["a(", k++, ") = ", n, " ", n^2, " ", s]; Sow[n]], {n, 1, 10^4}]][[2, 1]] (* _Jean-Fran├žois Alcover_, Sep 28 2016 *)

%o (Haskell)

%o a048653 n = a048653_list !! (n-1)

%o a048653_list = filter (f . show . (^ 2)) [1..] where

%o f zs = g (init $ tail $ inits zs) (tail $ init $ tails zs)

%o g (xs:xss) (ys:yss)

%o | h xs = h ys || f ys || g xss yss

%o | otherwise = g xss yss

%o where h ds = head ds /= '0' && a010052 (read ds) == 1

%o g _ _ = False

%o -- _Reinhard Zumkeller_, Oct 11 2011

%Y Cf. A048646, A048375, A010052, A000290; subsequence of A128783.

%K base,nice,nonn

%O 1,1

%A _Felice Russo_

%E Corrected and extended by _Naohiro Nomoto_, Sep 01 2001

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 19 21:28 EDT 2019. Contains 328244 sequences. (Running on oeis4.)