login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A229270
Numbers n for which n’-n is prime, n' being the arithmetic derivative of n.
4
18, 210, 315, 330, 390, 462, 510, 546, 690, 726, 798, 870, 930, 966, 1110, 1218, 1230, 1290, 1302, 1554, 1590, 1770, 2010, 2130, 2190, 2310, 2370, 2490, 2730, 2910, 3030, 3210, 3270, 3570, 3810, 4110, 4290, 4470, 4530, 4830, 4890, 5010, 5430, 5790, 5910, 5970
OFFSET
1,1
LINKS
EXAMPLE
315 is in the list because 315’ = 318 and 318 - 315 = 3 that is prime.
MAPLE
with(numtheory); P:=proc(q) local a, n, p; for n from 1 to q do
a:=n*add(op(2, p)/op(1, p), p=ifactors(n)[2]); if isprime(a-n) then print(n); fi; od; end: P(10^5);
PROG
(Python)
from sympy import isprime, factorint
A229270 = [n for n in range(1, 10**5) if isprime(sum([int(n*e/p) for p, e in factorint(n).items()])-n)] # Chai Wah Wu, Aug 21 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Sep 18 2013
STATUS
approved