cp's OEIS Frontend

A358068 Numbers that share a (decimal) digit with the sum of their proper divisors.

%I A358068 #13 Oct 30 2022 23:05:40
%S A358068 6,11,12,13,14,16,17,18,19,20,21,26,28,31,32,35,40,41,42,44,46,51,56,
%T A358068 60,61,64,68,70,71,72,74,76,80,84,86,91,93,95,96,100,101,102,103,104,
%U A358068 106,107,108,109,110,111,112,113,114,120,121,124,125,126,127,128,130,131,132,135,136
%N A358068 Numbers that share a (decimal) digit with the sum of their proper divisors.
%H A358068 Wikipedia, <a href="https://en.wikipedia.org/wiki/Table_of_divisors">Table of divisors</a>
%e A358068 72 is in the sequence since 72 shares a digit with the sum of its proper divisors 123. Both numbers share the digit "2".
%o A358068 (Python)
%o A358068 from itertools import count, islice
%o A358068 from sympy import divisor_sigma
%o A358068 def A358068_gen(startvalue=2): # generator of terms >= startvalue
%o A358068     for n in count(max(startvalue,2)):
%o A358068         s = set(str(divisor_sigma(n)-n))
%o A358068         if any(d in s for d in set(str(n))):
%o A358068             yield n
%o A358068 A358068_list = list(islice(A358068_gen(),30)) # _Chai Wah Wu_, Oct 29 2022
%o A358068 (PARI) isok(n) = #setintersect(Set(digits(sigma(n)-n)), Set(digits(n))); \\ _Michel Marcus_, Oct 30 2022
%Y A358068 Cf. A001065, A357929.
%K A358068 nonn,base
%O A358068 1,1
%A A358068 _Wesley Ivan Hurt_, Oct 29 2022