A107625 Every digit of prime and its index contains a loop (only digits 0,4,6,8,9 in prime and its index).
409, 4409, 4999, 6869, 44699, 44909, 46499, 48409, 64849, 66889, 68449, 68909, 86969, 89689, 98869, 98899, 480499, 488689, 490499, 609809, 806609, 806999, 866969, 868669, 868849, 869489, 869849, 869899, 869909, 880949, 6440809, 6440999, 6488969, 6489449
Offset: 1
Links
- Harvey P. Dale, Table of n, a(n) for n = 1..376
Programs
-
Mathematica
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}] Module[{c={0,4,6,8,9}},Select[Prime/@(Rest[FromDigits/@Tuples[c,6]]), SubsetQ[ c,IntegerDigits[#]]&]] (* Harvey P. Dale, Sep 03 2015 *)
-
PARI
is_a001744(n) = #setintersect(vecsort(digits(n), , 8), [1, 2, 3, 5, 7])==0 my(i=1); forprime(p=1, 65e5, if(is_a001744(p) && is_a001744(i), print1(p, ", ")); i++) \\ Felix Fröhlich, Sep 09 2019
Extensions
More terms from Harvey P. Dale, Sep 03 2015
Comments