cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A366938 a(n) = Sum_{k=1..n} (-1)^(k-1) * binomial(k+2,3) * floor(n/k).

Original entry on oeis.org

1, -2, 9, -14, 22, -27, 58, -85, 91, -97, 190, -243, 213, -266, 460, -499, 471, -553, 778, -970, 896, -845, 1456, -1697, 1264, -1560, 2270, -2289, 2207, -2307, 3150, -3793, 3049, -3125, 4765, -5079, 4061, -4492, 6634, -6714, 5628, -6370, 7821, -9120, 7986, -7013
Offset: 1

Views

Author

Seiichi Manyama, Oct 29 2023

Keywords

Crossrefs

Partial sums of A320901.

Programs

  • PARI
    a(n) = sum(k=1, n, (-1)^(k-1)*binomial(k+2, 3)*(n\k));
    
  • Python
    from math import isqrt
    def A366938(n): return (((s:=isqrt(m:=n>>1))*(s+1)**3*(s+2)<<4)-(t:=isqrt(n))*(t+1)**2*(t+2)*(t+3)-sum((((q:=m//w)+1)*(q*(q+1)*(q+2)+(w*(w+1)*((w<<1)+1)<<1))<<4) for w in range(1,s+1))+sum(((q:=n//w)+1)*(q*(q+2)*(q+3)+(w*(w+1)*(w+2)<<2)) for w in range(1,t+1)))//24 # Chai Wah Wu, Oct 29 2023

Formula

G.f.: 1/(1-x) * Sum_{k>=1} x^k/(1+x^k)^4 = -1/(1-x) * Sum_{k>=1} binomial(k+2,3) * (-x)^k/(1-x^k).

A366939 a(n) = Sum_{k=1..n} (-1)^(k-1) * binomial(k+3,4) * floor(n/k).

Original entry on oeis.org

1, -3, 13, -26, 45, -70, 141, -228, 283, -366, 636, -879, 942, -1232, 1914, -2331, 2515, -3090, 4226, -5313, 5539, -6114, 8837, -10558, 9988, -11947, 15969, -17705, 18256, -20364, 26013, -30592, 29330, -31874, 42222, -47034, 44357, -49602, 64164, -69115, 66637, -74017
Offset: 1

Views

Author

Seiichi Manyama, Oct 29 2023

Keywords

Crossrefs

Programs

  • PARI
    a(n) = sum(k=1, n, (-1)^(k-1)*binomial(k+3, 4)*(n\k));
    
  • Python
    from math import isqrt
    from sympy import rf
    def A366939(n): return ((rf(s:=isqrt(m:=n>>1),3)*(s+1)*((s**2<<2)+13*s+8)<<3)-rf(t:=isqrt(n),5)*(t+1)+sum((((q:=m//w)+1)*(-q*(q+2)*((q**2<<2)+13*q+8)-5*w*(w+1)*((r:=w<<1)+1)*(r+3))<<3) for w in range(1,s+1))+sum(rf(q:=n//w,5)+5*(q+1)*rf(w,4) for w in range(1,t+1)))//120 # Chai Wah Wu, Oct 29 2023

Formula

G.f.: 1/(1-x) * Sum_{k>=1} x^k/(1+x^k)^5 = -1/(1-x) * Sum_{k>=1} binomial(k+3,4) * (-x)^k/(1-x^k).
Showing 1-2 of 2 results.