OFFSET
1,1
COMMENTS
Similar to A387913.
When n+1 >= 5 is prime, a(n) seems to be divisible by n+1.
LINKS
Georg Fischer, Table of n, a(n) for n = 1..500
FORMULA
a(n) = 3*binomial(2*n+4,n+3)/(n+4)-(n+2)*(n+3)/2-1.
From Stefano Spezia, Sep 12 2025: (Start)
G.f.: 1 - (4 - 5*x + 2*x^2)/(1 - x)^3 + (1 - sqrt(1 - 4*x) + (sqrt(1 - 4*x) - 3)*x - 2*x^3)/(2*x^4).
E.g.f.: 1 - exp(x)*(8 + 6*x + x^2)/2 - 6*exp(2*x)*(1 - x)*BesselI(0, 2*x)/x^2 + 6*exp(2*x)*(1 - x + x^2)*BesselI(1, 2*x)/x^3. (End)
D-finite with recurrence (n+4)*a(n) +3*(-2*n-5)*a(n-1) +9*(n+1)*a(n-2) +2*(-2*n+1)*a(n-3) +3*(-4*n-1)=0. - R. J. Mathar, Sep 26 2025
EXAMPLE
The value a(3) = 74 is the Hodge number h^{1,2} of a Calabi-Yau hypersurface in a 4-dimensional toric variety.
PROG
(SageMath) [3*binomial(2*n+4, n+3)/(n+4)-(n+2)*(n+3)/2-1 for n in range(1, 33)]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
F. Chapoton, Sep 12 2025
EXTENSIONS
a(19) and following terms corrected by Georg Fischer, Sep 16 2025
STATUS
approved
