|
%I #8 Aug 03 2014 14:26:15
%S 61,71,661,881,991,6661,99991,9999991,6666666661,7777777777771,
%T 666666666666666661,8888888888888888881,77777777777777777771,
%U 666666666666666666661,6666666666666666666661,77777777777777777777771
%N Primes of the form aaaa...aa1 where a is 6, 7, 8 or 9.
%C Easy-to-remember large primes can be formed in this manner.
%H Vincenzo Librandi, <a href="/A109549/b109549.txt">Table of n, a(n) for n = 1..43</a>
%F d=6, 7, 8, 9 a(n) = if prime then Sum[d*10^i, {i, 1, m}] + 1
%t d[n_] = If[5 + Mod[n, 6] > 0, 5 + Mod[n, 6], 1] a = Flatten[Table[Sum[d[k]*10^i, {i, 1, m}] + 1, {m, 1, 50}, {k, 1, 4}]] b = Flatten[Table[If[PrimeQ[a[[i]]] == True, a[[i]], {}], {i, 1, Length[a]}]]
%t Select[FromDigits/@Flatten[Table[PadLeft[{1},i,#]&/@{6,7,8,9},{i,2,100}],1],PrimeQ[#]&] (* _Vincenzo Librandi_, Dec 12 2011 *)
%Y Cf. A051200, A004022, A093176, A093177, A089346, A093174, A092571, A089345, A089347.
%K nonn
%O 1,1
%A _Roger L. Bagula_, Jun 26 2005
|
