OFFSET
0,8
MAPLE
b:= proc(n) option remember; `if`(n=0, 1, add(b(n-j)*add(
`if`(d::odd, d, 0), d=numtheory[divisors](j)), j=1..n)/n)
end:
a:= n-> igcd(b(n), combinat[numbpart](n)):
seq(a(n), n=0..120); # Alois P. Heinz, May 23 2023
MATHEMATICA
Table[GCD[PartitionsP[n], PartitionsQ[n]], {n, 0, 100}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, May 23 2023, inspired by Zhi-Wei Sun
STATUS
approved