cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A162948 Numbers with a sum of digits equal to their smallest prime factor.

This page as a plain text file.
%I A162948 #8 Apr 17 2019 03:44:20
%S A162948 2,3,5,7,20,21,110,111,133,200,201,209,247,407,481,511,629,803,1010,
%T A162948 1011,1100,1101,1141,1387,1417,1651,1679,1853,2000,2001,2023,2119,
%U A162948 2159,2353,2401,2771,3031,3077,3097,3383,3439,3523,3749,3781,4577,4607,4913,5149,5161
%N A162948 Numbers with a sum of digits equal to their smallest prime factor.
%H A162948 Robert Israel, <a href="/A162948/b162948.txt">Table of n, a(n) for n = 1..10000</a>
%F A162948 {n: A007953(n) = A020639(n)}. - _R. J. Mathar_, Jul 19 2009
%e A162948 133 is in the sequence because 1 + 3 + 3 = 7 and 7 is the smallest prime factor of 133.
%p A162948 A007953 := proc(n) local d; add(d,d=convert(n,base,10)) ; end:
%p A162948 A020639 := proc(n) min(op(numtheory[factorset](n) )) ; end:
%p A162948 isA162948 := proc(n) RETURN( A007953(n) = A020639(n)) ; end:
%p A162948 for n from 1 to 6000 do if isA162948(n) then printf("%d,",n) ; fi; od: # _R. J. Mathar_, Jul 19 2009
%Y A162948 Cf: A007953, A020639.
%Y A162948 Even terms: A069537.
%K A162948 nonn,base
%O A162948 1,1
%A A162948 _Claudio Meller_, Jul 18 2009
%E A162948 Single-digit primes added by _R. J. Mathar_, Jul 19 2009