%I #9 Sep 09 2019 10:58:13
%S 409,4409,4999,6869,44699,44909,46499,48409,64849,66889,68449,68909,
%T 86969,89689,98869,98899,480499,488689,490499,609809,806609,806999,
%U 866969,868669,868849,869489,869849,869899,869909,880949,6440809,6440999,6488969,6489449
%N Every digit of prime and its index contains a loop (only digits 0,4,6,8,9 in prime and its index).
%C Corresponding indices in A107624. Cf. A001744 Every digit contains a loop.
%H Harvey P. Dale, <a href="/A107625/b107625.txt">Table of n, a(n) for n = 1..376</a>
%t Do[id=Union[IntegerDigits[p=Prime[n]], IntegerDigits[n]];If[Count[id, 1]+Count[id, 2]+Count[id, 3]+Count[id, 5]+Count[id, 7]==0, Print[p]], {n, 10000}]
%t Module[{c={0,4,6,8,9}},Select[Prime/@(Rest[FromDigits/@Tuples[c,6]]), SubsetQ[ c,IntegerDigits[#]]&]] (* _Harvey P. Dale_, Sep 03 2015 *)
%o (PARI) is_a001744(n) = #setintersect(vecsort(digits(n), , 8), [1, 2, 3, 5, 7])==0
%o my(i=1); forprime(p=1, 65e5, if(is_a001744(p) && is_a001744(i), print1(p, ", ")); i++) \\ _Felix Fröhlich_, Sep 09 2019
%Y Cf. A001744, A107624.
%K nonn,base,dumb
%O 1,1
%A _Zak Seidov_, May 18 2005
%E More terms from _Harvey P. Dale_, Sep 03 2015
|