OFFSET
1,5
COMMENTS
EXAMPLE
a(6) = 31, because 13 and 31 use the same set of digits, and no smaller prime distinct from 13 has this property.
PROG
(PARI) a(n) = if (n<=4, return(-1)); if (n==5, my(p=prime(n), s=Set(digits(p)), q=1); while (!((q!=p) && isprime(q)), q = eval(concat(q, "1"))); return(q)); my(p=prime(n), s=Set(digits(p)), q=2); while(!((Set(digits(q)) == s) && (p!=q)), q=nextprime(q+1)); q; \\ Michel Marcus, Feb 19 2026
(Python)
from sympy import isprime, prime
from itertools import count, product
def a(n):
if n < 5: return -1
pn = prime(n)
S = sorted(set(str(pn)))
for d in count(len(S)):
for p in product(S, repeat=d):
if p[0] == "0" or len(set(p)) != len(S): continue
t = int("".join(p))
if t != pn and isprime(t):
return t
print([a(n) for n in range(1, 56)]) # Michael S. Branicky, Feb 20 2026
CROSSREFS
KEYWORD
sign,base
AUTHOR
Jean-Marc Rebert, Feb 19 2026
STATUS
approved
