login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A076161 Numbers n such that n + sum of squares of digits of n (A258881) is a prime. 4
1, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 31, 34, 57, 73, 74, 75, 78, 91, 94, 97, 100, 101, 102, 103, 105, 107, 108, 109, 121, 122, 123, 126, 127, 128, 140, 142, 146, 148, 160, 161, 165, 166, 168, 182, 183, 188, 213, 216, 217, 234, 251, 275, 277, 297, 301 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

LINKS

Robert Israel, Table of n, a(n) for n = 1..10000

EXAMPLE

12 is a term because 12+(1^2+2^2) = 17 is a prime.

MAPLE

filter:= proc(n) local t; isprime(n+add(t^2, t=convert(n, base, 10))) end proc:

select(filter, [$1..1000]); # Robert Israel, Jan 30 2021

PROG

(Python)

from sympy import isprime

def ssd(n): return sum(int(d)**2 for d in str(n))

def ok(n): return isprime(n + ssd(n))

def aupto(limit): return [m for m in range(1, limit+1) if ok(m)]

print(aupto(301)) # Michael S. Branicky, Jan 30 2021

(PARI) isok(n) = isprime(n + norml2(digits(n))); \\ Michel Marcus, Jan 31 2021

CROSSREFS

Cf. A258881, A259391, A259567.

Sequence in context: A007959 A102492 A169905 * A318699 A131835 A262390

Adjacent sequences:  A076158 A076159 A076160 * A076162 A076163 A076164

KEYWORD

nonn,base

AUTHOR

Zak Seidov, Nov 01 2002

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 7 11:39 EDT 2022. Contains 357270 sequences. (Running on oeis4.)