%I #30 Jun 01 2022 16:20:06
%S 1,1,1,1,19,17,1,409,604,82,20951,229931,411012,39080794,4382914408,
%T 6345486566,45119290746,581075656330,8672770990,869561574799171,
%U 71853663603175593,25509154378676494,24040267482771436703,102403319155457392955,11302410854347731819765
%N a(n) = A118478(n)*(A118478(n)+1) divided by the product of the first n primes.
%C When floor(A214089(n) / 2) = A118478(n), a(n) = A215085(n).
%H Chai Wah Wu, <a href="/A215021/b215021.txt">Table of n, a(n) for n = 1..31</a>
%o (Python 3.8+)
%o from itertools import combinations
%o from math import prod
%o from sympy import sieve, prime, primorial
%o from sympy.ntheory.modular import crt
%o def A215021(n): return 1 if n == 1 else (s := int(min(min(crt((m, (k:=primorial(n))//m), (0, -1))[0], crt((k//m, m), (0, -1))[0]) for m in (prod(d) for l in range(1, n//2+1) for d in combinations(sieve.primerange(prime(n)+1), l)))))*(s+1)//k # _Chai Wah Wu_, May 31 2022
%Y Cf. A118478, A002110.
%K nonn
%O 1,5
%A _N. J. A. Sloane_, Aug 02 2012, following a comment in A118478 made by _Robert G. Wilson v_ on May 13 2006
%E a(18)-a(25) from _J. Stauduhar_, Aug 17 2012