cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A107625 Every digit of prime and its index contains a loop (only digits 0,4,6,8,9 in prime and its index).

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