A154313 Numbers n such that abs(A007605(n) - A007953(n)) < 2.

1, 2, 7, 13, 16, 29, 32, 38, 53, 54, 56, 63, 66, 68, 69, 76, 88, 94, 126, 156, 175, 176, 182, 183, 191, 192, 212, 213, 218, 227, 248, 252, 255, 258, 259, 280, 282, 286, 291, 293, 294, 295, 298, 306, 307, 321, 323, 324, 325, 326, 331, 334, 335, 338, 345, 348

Offset

1,2

Comments

Also, numbers n such that abs(A090431(n)) < 2. - Omar E. Pol, Jan 12 2009

Links

G. C. Greubel and Charles R Greathouse IV, Table of n, a(n) for n = 1..10000 (first 1000 terms from Greubel)

Example

If n = 16 then prime(n) = 53 and abs((5+3)-(1+6)) = 8-7 = 1, so 16 is in the sequence.

Mathematica

Select[Range[500], Abs[Apply[Plus, RealDigits[Prime[#]][[1]]] - Sum[DigitCount[#][[i]]*i, {i, 9}]] < 2 &] (* G. C. Greubel, Sep 10 2016 *)

Prog

(PARI) list(lim)=my(v=List(), n); forprime(p=2, , if(n++>lim, break); if(abs(sumdigits(n)-sumdigits(p))<2, listput(v, n))); Vec(v) \\ Charles R Greathouse IV, Sep 10 2016

Crossrefs

Cf. A000040, A141468, A007605, A007953, A090431.

Sequence in context: A018322 A063222 A030592 * A196012 A065104 A138645

Adjacent sequences: A154310 A154311 A154312 * A154314 A154315 A154316

Keyword

nonn,base,easy

Author

Juri-Stepan Gerasimov, Jan 07 2009

Extensions

Corrected and edited by Omar E. Pol, Jan 12 2009

Status

approved