login
In base-10 notation replace digits of n with their squared values (Version 1).
19

%I #32 Sep 08 2022 08:44:57

%S 0,1,4,9,16,25,36,49,64,81,10,11,14,19,116,125,136,149,164,181,40,41,

%T 44,49,416,425,436,449,464,481,90,91,94,99,916,925,936,949,964,981,

%U 160,161,164,169,1616,1625,1636,1649,1664,1681,250,251,254,259,2516,2525

%N In base-10 notation replace digits of n with their squared values (Version 1).

%H Alois P. Heinz, <a href="/A048385/b048385.txt">Table of n, a(n) for n = 0..10000</a>

%F a(n) >= n with equality iff n belongs to A007088. - _Rémy Sigrist_, May 17 2019

%p a:= n-> (s-> parse(cat(seq(parse(s[i])^2, i=1..length(s)))))(""||n):

%p seq(a(n), n=0..70); # _Alois P. Heinz_, Jul 04 2014

%t Table[FromDigits[Flatten[IntegerDigits/@(IntegerDigits[n]^2)]],{n,0,80}] (* _Harvey P. Dale_, May 06 2019 *)

%o (Magma) [0] cat [StringToInteger(&cat[IntegerToString(h): h in Reverse([i^2: i in Intseq(n)])]): n in [1..55]]; // _Bruno Berselli_, Jul 31 2012

%o (Python)

%o def digits(n):

%o d=[]

%o while n>0:

%o d.append(n%10)

%o n=n//10

%o return d

%o def sqdig(n):

%o new=0

%o num=digits(n)

%o spacing=0

%o while num:

%o k=num.pop(0)

%o new+=(10**(spacing))*(k**2)

%o if k>3:

%o spacing+=1

%o spacing+=1

%o return new

%o # _David Nacin_, Aug 19 2012

%o (Haskell)

%o a048385 0 = 0

%o a048385 n = read (show (a048385 n') ++ show (m ^ 2)) :: Integer

%o where (n', m) = divMod n 10

%o -- _Reinhard Zumkeller_, Jul 08 2014

%o (PARI) a(n) = if (n, fromdigits(concat(apply(d -> my (d2=d^2); if (d2, digits(d2), [0]), digits(n)))), 0) \\ _Rémy Sigrist_, May 17 2019

%o (MATLAB) m=1;

%o for u=0:200 digit=dec2base(u,10)-'0'; digitp=digit.^2;

%o sol(m)=str2num(strrep(num2str(digitp), ' ', ''));m=m+1;

%o end

%o sol % _Marius A. Burtea_, May 17 2019

%Y See A068522 for another version.

%Y Cf. A007088, A048386, A048387, A048388, A048389.

%K nonn,base

%O 0,3

%A _Patrick De Geest_, Mar 15 1999