|
%I #13 Jan 03 2016 17:35:03
%S 2,11,23,29,41,43,47,61,67,89,149,163,167,211,223,227,229,233,251,257,
%T 263,269,281,349,367,383,419,431,433,439,463,491,521,523,569,587,613,
%U 617,631,643,659,661,673,761,769,821,827,857,883,887,929,941,967,1123
%N Primes p the product of whose digits increased by one is prime.
%H Harry J. Smith, <a href="/A066725/b066725.txt">Table of n, a(n) for n = 1..1000</a>
%F A084354 MINUS A011540. - _R. J. Mathar_, Aug 26 2007
%e 2 * 2 * 3 + 1 = 13, a prime; since 223 is also prime, it is a term of the sequence.
%t Select[Prime[Range[200]],PrimeQ[Times@@IntegerDigits[#]+1]&] (* _Harvey P. Dale_, May 09 2012 *)
%o (PARI) ProdD(x)= { local(p=1); while (x>9 && p>0, p*=x%10; x\=10); return(p*x) } { n=0; for (m=1, 10^10, if (isprime(ProdD(p=prime(m)) + 1), write("b066725.txt", n++, " ", p); if (n==1000, return)) ) } \\ _Harry J. Smith_, Mar 19 2010
%K easy,nonn,base
%O 1,1
%A _Joseph L. Pe_, Jan 15 2002
%E Terms a(35) - a(54) from _Harry J. Smith_, Mar 19 2010
|
