A116017 Numbers m such that m + sigma(m) is a repdigit.
1, 2, 3, 4, 5, 9, 34, 141, 198, 277, 297, 375, 499, 1420, 2651, 2777, 3554, 4999, 19050, 28660, 29128, 49999, 131061, 506311, 3844863, 3852517, 4761903, 4999999, 22222218, 37560831, 133878933, 506767303, 872011214, 1381799253, 1427435733, 2777777777, 3018915632, 3555555554
Offset: 1
Examples
22222218 + sigma(22222218) = 66666666.
Links
- Max Alekseyev, Table of n, a(n) for n = 1..61 (terms 1..45 from Hiroaki Yamanouchi, terms 46..51 from Giovanni Resta)
Programs
-
Mathematica
Do[If[Length[Union[IntegerDigits[n + DivisorSigma[1, n]]]]==1, Print[n]], {n, 60000000}] (* Farideh Firoozbakht, Aug 17 2006 *)
-
PARI
for(n=1, 10^7, d=digits(sigma(n)+n); c=0; for(i=1, #d-1, if(d[i]!=d[i+1], c++; break)); if(c==0, print1(n, ", "))) \\ Derek Orr, Aug 01 2014
-
Python
from sympy import divisors A116017 = [n for n in range(1,10**5) if len(set(str(n+sum(divisors(n))))) == 1] # Chai Wah Wu, Aug 11 2014
Extensions
More terms from Farideh Firoozbakht, Aug 17 2006, Dec 19 2007
a(36)-a(37) from Donovan Johnson, Feb 17 2013
a(38) from Farideh Firoozbakht, Aug 01 2014
Comments