A061595 Product of digits + 1 is prime, sum of digits + 1 is prime and sum of digits - 1 is prime.

4, 6, 22, 66, 112, 114, 121, 123, 129, 132, 141, 147, 156, 165, 174, 189, 192, 198, 211, 213, 219, 231, 237, 273, 279, 291, 297, 312, 321, 327, 345, 354, 369, 372, 396, 411, 417, 435, 453, 459, 468, 471, 477, 486, 495, 516, 534, 543, 549, 561, 567, 576, 594

Offset

1,1

Links

Harry J. Smith, Table of n, a(n) for n = 1..1003

Example

For 147 we have (1*4*7) + 1 = 29, (1+4+7) + 1 = 13, (1+4+7) - 1 = 11.

Maple

a := proc (n) local nn: nn := convert(n, base, 10): if isprime(1+product(nn[j], j = 1 .. nops(nn))) = true and isprime(1+sum(nn[j], j = 1 .. nops(nn))) = true and isprime(-1+sum(nn[j], j = 1 .. nops(nn))) = true then n else end if end proc: seq(a(n), n = 1 .. 615); # Emeric Deutsch, Aug 02 2009

Prog

(PARI) SumD(x)= { s=0; while (x>9, s=s+x-10*(x\10); x=x\10); return(s + x) } ProdD(x)= { p=1; while (x>9, p=p*(x-10*(x\10)); x=x\10); return(p*x) } { x=9; for (n=0, 1000, until(isprime(s+1) && isprime(s-1) && isprime(ProdD(x)+1), x++; s=SumD(x)); write("b061595.txt", n, " ", x) ) } \\ Harry J. Smith, Jul 25 2009

Crossrefs

Subsequence of A167711.

Sequence in context: A192154 A289385 A151519 * A062940 A061596 A061597

Adjacent sequences: A061592 A061593 A061594 * A061596 A061597 A061598

Keyword

nonn,base

Author

Felice Russo, May 22 2001

Extensions

a(1)=4 and a(2)=6 added by Emeric Deutsch, Aug 02 2009

Status

approved