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

Author

Juri-Stepan Gerasimov, Jan 07 2009

Keywords

• nonn
• base
• easy

Comments

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

Examples

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

Links

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

Crossrefs

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

Programs

Mathematica

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

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

Extensions

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