OFFSET
1,1
COMMENTS
Numbers k such that A004186(k) is prime. - Robert Israel, Apr 01 2024
If N is a term then all numbers with the same digits as N are terms too.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
142 is a term because its digits in decreasing order form 421 and this is prime.
MAPLE
dd:= proc(n) local L, i;
L:= sort(convert(n, base, 10));
add(L[i]*10^(i-1), i=1..nops(L))
end proc:
select(isprime @ dd, [$1..1000]); # Robert Israel, Apr 01 2024
MATHEMATICA
Select[Range[500], PrimeQ[FromDigits[ReverseSort[IntegerDigits[#]]]] &]
PROG
(Python)
from sympy import isprime
def ok(n): return isprime(int("".join(sorted(str(n), reverse=True))))
print([k for k in range(200) if ok(k)]) # Michael S. Branicky, Apr 01 2024
(Python)
from itertools import count, islice, combinations_with_replacement
from sympy import isprime
from sympy.utilities.iterables import multiset_permutations
def A371653_gen(): # generator of terms
for l in count(1):
xlist = []
for p in combinations_with_replacement('987654321', l):
if isprime(int(''.join(p))):
xlist.extend(int(''.join(d)) for d in multiset_permutations(p))
yield from sorted(xlist)
CROSSREFS
KEYWORD
nonn,base,easy
AUTHOR
César Eliud Lozada, Apr 01 2024
STATUS
approved