login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A173689 Numbers m such that the sum of square of factorial of decimal digits is square. 2

%I #19 Feb 23 2023 11:53:14

%S 0,1,2,3,4,5,6,7,8,9,122,202,212,220,221,244,424,442,1000,1001,1010,

%T 1011,1100,1101,1110,1111,2222,3333,3444,4344,4434,4443,4444,5555,

%U 6666,6677,6767,6776,6888,7667,7676,7766,7777,8688,8868,8886,8888,9999

%N Numbers m such that the sum of square of factorial of decimal digits is square.

%C Let the decimal expansion of m = d(0)d(1)...d(p). Numbers such that Sum_{k=0..p} (d(k)!)^2 is square.

%H Jinyuan Wang, <a href="/A173689/b173689.txt">Table of n, a(n) for n = 1..1487</a>

%e a(16) = 244 is in the sequence because (2!)^2 + (4!)^2 + (4!)^2 = 1156 = 34^2.

%p with(numtheory):for n from 0 to 10000 do:l:=length(n):n0:=n:s:=0:for m from

%p 1 to l do:q:=n0:u:=irem(q,10):v:=iquo(q,10):n0:=v :s:=s+(u!)^2:od: q:=sqrt(s):if

%p floor(q)= q then printf(`%d, `,n):else fi:od:

%t Select[Range[0,10000],IntegerQ[Sqrt[Total[(IntegerDigits[#]!)^2]]]&] (* _Harvey P. Dale_, Dec 19 2011 *)

%o (Python)

%o from itertools import count, islice, combinations_with_replacement

%o from math import factorial

%o from sympy.ntheory.primetest import is_square

%o from sympy.utilities.iterables import multiset_permutations

%o def A173689_gen(): # generator of terms

%o yield 0

%o for l in count(0):

%o for i in range(1,10):

%o fi = factorial(i)**2

%o yield from sorted(int(str(i)+''.join(map(str,k))) for j in combinations_with_replacement(range(10), l) for k in multiset_permutations(j) if is_square(fi+sum(map(lambda n:factorial(n)**2,j))))

%o A173689_list = list(islice(A173689_gen(),50)) # _Chai Wah Wu_, Feb 23 2023

%Y Cf. A173687, A173688.

%K nonn,base

%O 1,3

%A _Michel Lagneau_, Nov 25 2010

%E Offset changed to 1 by _Jinyuan Wang_, Feb 26 2020

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 10:22 EDT 2024. Contains 371268 sequences. (Running on oeis4.)