(PARI) P(n) = vecprod(primes(n));
a(n) = my(k=n+1, sp=Str(P(n))); while (#strsplit(Str(P(k)), sp) < 2, k++); k; \\
Michel Marcus, Dec 18 2025
(Python)
from itertools import count
from gmpy2 import mpz, digits, next_prime
from sympy import primorial, prime
if n == 0: return 4
m = mpz(primorial(n))
p = prime(n+1)
s = digits(m)
for k in count(n+1):
m *= p
if s in digits(m):
return k