OFFSET
0,3
COMMENTS
LINKS
Eric Weisstein's World of Mathematics, Erdős Squarefree Conjecture.
FORMULA
a(n) = binomial(2*n, n) / Product(p prime | p^2 > 2*n, floor(2*n/p) is odd).
EXAMPLE
binomial(10, 5) = 2^2*3^2*7. 2,3 <= sqrt(10), 7 > sqrt(10) so a(5) = 2^2*3^2 = 36.
PROG
(R)
library(primes)
n = 2*10^4
pp = generate_primes(max = sqrt(n))
npinb = list()
for(p in pp){
np = rep(0, n)
for(t in 1:(log(n)/log(p))) np[p*(1:(n/p))] = np[1:(n/p)] + 1 #multiplicities of prime factor p in the integers
npinb[[as.character(p)]] = cumsum(np)[2*(1:(n/2))] - 2*cumsum(np[1:(n/2)]) #multiplicities of prime factor p in the central binomial coefficients
}
res = rep(1, n/2)
for(p in pp){
select = p^2 <= 2*(1:(n/2))
res[select] = res[select]*p^npinb[[as.character(p)]][select]
}
#offset of res is 1
(PARI) a(n) = my(m=sqrtint(2*n), f=factor(binomial(2*n, n), m+1)); for (k=1, #f~, if (f[k, 1]>m, f[k, 1]=1)); factorback(f); \\ Michel Marcus, Mar 03 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Matthieu Pluntz, Mar 01 2023
STATUS
approved