OFFSET
1,1
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..10000
PROG
(Python)
from sympy import prime, primerange
def auptopn(lim):
alst = []
for k, pk in enumerate(primerange(2, prime(lim)+1), start=1):
sdk = sum(map(int, str(k)))
if str(sdk) in str(pk): alst.append(sdk)
return alst
print(auptopn(1010)) # Michael S. Branicky, Jul 03 2021
CROSSREFS
KEYWORD
easy,nonn,base
AUTHOR
Zak Seidov, Sep 26 2002
EXTENSIONS
Corrected and extended by Sascha Kurz, Jan 30 2003
STATUS
approved