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

0, 1, 4, 9, 16, 25, 36, 49, 64, 81, 10, 11, 14, 19, 116, 125, 136, 149, 164, 181, 40, 41, 44, 49, 416, 425, 436, 449, 464, 481, 90, 91, 94, 99, 916, 925, 936, 949, 964, 981, 160, 161, 164, 169, 1616, 1625, 1636, 1649, 1664, 1681, 250, 251, 254, 259, 2516, 2525

Offset: 0

Author

Patrick De Geest, Mar 15 1999

Keywords

• nonn
• base

Links

• Alois P. Heinz, Table of n, a(n) for n = 0..10000

Crossrefs

See A068522 for another version.

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

Programs

Haskell

a048385 0 = 0
a048385 n = read (show (a048385 n') ++ show (m ^ 2)) :: Integer
            where (n', m) = divMod n 10
-- Reinhard Zumkeller, Jul 08 2014

MATLAB

m=1;
for u=0:200 digit=dec2base(u,10)-'0'; digitp=digit.^2;
    sol(m)=str2num(strrep(num2str(digitp), ' ', ''));m=m+1;
end
sol % Marius A. Burtea, May 17 2019

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

Maple

a:= n-> (s-> parse(cat(seq(parse(s[i])^2, i=1..length(s)))))(""||n):
seq(a(n), n=0..70);  # Alois P. Heinz, Jul 04 2014

Mathematica

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

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

Python

def digits(n):
    d=[]
    while n>0:
        d.append(n%10)
        n=n//10
    return d
def sqdig(n):
    new=0
    num=digits(n)
    spacing=0
    while num:
        k=num.pop(0)
        new+=(10**(spacing))*(k**2)
        if k>3:
            spacing+=1
        spacing+=1
    return new
# David Nacin, Aug 19 2012

Formula

a(n) >= n with equality iff n belongs to A007088. - Rémy Sigrist, May 17 2019