login
Primes whose multiplicative digital root is 2.
2

%I #18 Nov 05 2020 19:22:11

%S 2,37,43,73,137,173,211,223,317,367,389,431,673,827,839,929,983,1223,

%T 1279,1297,1367,1447,1621,1637,1693,1999,2111,2161,2179,2213,2269,

%U 2339,2393,2663,2699,2719,2791,2917,2969,2971,3167,3169,3221,3329,3463,3499,3617

%N Primes whose multiplicative digital root is 2.

%C Complement of A199981 with respect to A034049, numbers whose multiplicative digital root is 2.

%H Robert Israel, <a href="/A199980/b199980.txt">Table of n, a(n) for n = 1..10000</a>

%e Prime 389 is in sequence because 3*8*9=216, 2*1*6 =12, 1*2=2.

%p mdr:= proc(n) option remember;

%p local t;

%p t:= convert(convert(n,base,10),`*`);

%p if t < 10 then t else procname(t) fi

%p end proc:

%p select(t -> mdr(t) = 2 and isprime(t), [2, seq(i,i=3..10000,2)]); # _Robert Israel_, Nov 05 2020

%t t = {}; n = 0; While[Length[t] < 100, n = NextPrime[n]; s = n; While[s >= 10, s = Times @@ IntegerDigits[s]]; If[s == 2, AppendTo[t, n]]]; t (* _T. D. Noe_, Nov 15 2011 *)

%t Select[Prime[Range[600]],FixedPoint[Times@@IntegerDigits[#]&,#]==2&] (* _Harvey P. Dale_, Mar 28 2012 *)

%Y Cf. A199981 (nonprime numbers whose multiplicative digital root is 2).

%Y Includes A107612.

%Y Cf. A031347, A034049.

%K nonn,base

%O 1,1

%A _Jaroslav Krizek_, Nov 13 2011