|
|
A107625
|
|
Every digit of prime and its index contains a loop (only digits 0,4,6,8,9 in prime and its index).
|
|
4
|
|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Corresponding indices in A107624. Cf. A001744 Every digit contains a loop.
|
|
LINKS
|
|
|
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 *)
|
|
PROG
|
(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
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base,dumb
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|