|
%I #27 Feb 15 2021 22:42:01
%S 2,13,113,149,1013,1039,1427,1123,1439,1579,1237,10271,10453,10139,
%T 10253,10243,10457,11579,10789,10273,11239,12457,10729,13249,12347,
%U 13687,12539,14759,13799,10739,12637,12893,23957,13597,100493,12379,14593,101383,13789
%N Least prime with exactly n prime anagrams not equal to itself.
%C 13 has only one prime anagram (31), and no smaller prime has a prime anagram other than itself, so a(1) = 13.
%C 113 has 2 prime anagrams (131 and 311), and no smaller prime has two prime anagrams other than itself, so a(2) = 113.
%C 149 has 3 prime anagrams (419, 491, and 941), and no smaller prime has three prime anagrams other than itself, so a(3) = 149.
%H Michael S. Branicky, <a href="/A166921/b166921.txt">Table of n, a(n) for n = 0..2780</a> (terms 83..223 from P. CAMI and Chai Wah Wu, and terms 1..82 from P. CAMI)
%H Michael S. Branicky, <a href="/A166921/a166921.py.txt">Python program</a>
%e a(7) = prime 1123 with 7 prime anagrams 1213, 1231, 1321, 2113, 2131, 2311, 3121.
%o (Python) # see link for faster version
%o from sympy import isprime
%o from itertools import permutations
%o def anagrams(n):
%o s = str(n)
%o return set(int("".join(p)) for p in permutations(s) if p[0] != '0')
%o def num_prime_anagrams(n): return sum(isprime(i) for i in anagrams(n))
%o def a(n):
%o if n == 0: return 2
%o k = 3
%o while not isprime(k) or num_prime_anagrams(k) != n+1: k += 2
%o return k
%o print([a(n) for n in range(39)]) # _Michael S. Branicky_, Feb 13 2021
%Y Cf. A039986, A046810.
%K nonn,base
%O 0,1
%A _Pierre CAMI_, Oct 23 2009
%E Definition edited and a(0) added by _Chai Wah Wu_, Dec 26 2016
|
