login
Sum of digits of n is the largest prime factor of n.
7

%I #20 Nov 30 2015 09:39:10

%S 2,3,5,7,12,50,70,308,320,364,476,500,605,700,704,715,832,935,1088,

%T 1183,1547,1729,2401,2584,2618,2704,2926,3080,3200,3536,3640,3952,

%U 4225,4760,4784,4913,5000,5491,5525,5819,5831,6050,6175,6517,6647,7000,7040,7150

%N Sum of digits of n is the largest prime factor of n.

%H Reinhard Zumkeller, <a href="/A052021/b052021.txt">Table of n, a(n) for n = 1..1000</a>

%F {n: A007953(n) = A006530(n)}. - _R. J. Mathar_, May 30 2010

%e 13685 has sum of digits '23' and 13685 = 5*7*17*'23'.

%p A007953 := proc(n) add(d,d=convert(n,base,10)) ; end proc:

%p A006530 := proc(n) numtheory[factorset](n) ; max(op(%)) ; end proc:

%p for n from 1 to 8000 do if A007953(n) = A006530(n) then printf("%d,",n) ; end if; end do: # _R. J. Mathar_, May 30 2010

%t Select[Range[2,8000],FactorInteger[#][[-1,1]]==Total[IntegerDigits[#]]&] (* _Harvey P. Dale_, Oct 17 2012 *)

%o (Haskell)

%o a052021 n = a052021_list !! (n-1)

%o a052021_list = tail $ filter (\x -> a007953 x == a006530 x) [1..]

%o -- _Reinhard Zumkeller_, Nov 06 2011

%Y Cf. A052018, A052019, A052020, A052022, A007953, A005349, A028834.

%K nonn,base,nice

%O 1,1

%A _Patrick De Geest_, Nov 15 1999

%E Single-digit primes added by _R. J. Mathar_, May 30 2010

%E Offset corrected by _Reinhard Zumkeller_, Nov 05 2011