OFFSET
1,1
COMMENTS
Composite numbers n not in A176670 such that the sum of the 4th power of the digits of n equals the sum of the 4th power of the digits of the prime factors of n (with multiplicity). A176670 lists composite numbers having the same digits as their prime factors (with multiplicity), excluding zero digits.
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..1000
Patrick Costello, A new largest Smith number, Fibonacci Quarterly 40(4) (2002), 369-371.
Underwood Dudley, Smith numbers, Mathematics Magazine 67(1) (1994), 62-65.
S. S. Gupta, Smith Numbers, Mathematical Spectrum 37(1) (2004/5), 27-29.
S. S. Gupta, Smith Numbers.
Eric Weisstein's World of Mathematics, Smith number.
Wikipedia, Smith number.
A. Wilansky, Smith Numbers, Two-Year College Math. J. 13(1) (1982), p. 21.
Amin Witno, Another simple construction of Smith numbers, Missouri J. Math. Sci. 22(2) (2010), 97-101.
Amin Witno, Smith multiples of a class of primes with small digital sum, Thai Journal of Mathematics 14(2) (2016), 491-495.
EXAMPLE
3777 = 3*1259 is composite; sum of 4th power of the digits is 3^4 + 7^4 + 7^4 + 7^4 = 7284. Sum of 4th power of the digits of the prime factors 3, 1259 is 3^4 + 1^4 + 2^4 + 5^4 + 9^4 = 7284. The sums are equal, so 3777 is in the sequence.
17418 = 2*3*2903 is composite; sum of 4th power of the digits is 1^4 + 7^4 + 4^4 + 1^4 + 8^4 = 6755. Sum of 4th power of the digits of the prime factors 2, 3, 2903 is 2^4 + 3^4 + 2^4 + 9^4 + 0^4 + 3^4 = 6755. The sums are equal, so 17418 is in the sequence.
269636 = 2*2*67409 is composite; sum of 4th power of the digits is 2^4 + 6^4 + 9^4 + 6^4 + 3^4 + 6^4 = 10546. Sum of 4th power of the digits of the prime factors 2, 2, 67409 (with multiplicity) is 2^4 + 2^4 + 6^4 + 7^4 + 4^4 + 0^4 + 9^4 = 10546. The sums are equal, so 269636 is in the sequence.
MATHEMATICA
fQ[n_] := Block[{id = Sort@ IntegerDigits@ n, fid = Sort@ Flatten[ IntegerDigits@ Table[#[[1]], {#[[2]]}] & /@ FactorInteger@ n]}, While[ id[[1]] == 0, id = Drop[id, 1]]; While[ fid[[1]] == 0, fid = Drop[fid, 1]]; id != fid && Plus @@ (id^4) == Plus @@ (fid^4)]; k = 1; lst = {}; While[k < 10^6, If[f Q@ k, AppendTo[lst, k]; Print@ k]; k++]; lst
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Paul Weisenhorn, Dec 19 2010
STATUS
approved