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.
%I A371372 #14 Jun 23 2024 16:11:28 %S A371372 0,1,5,40,405,4626,56360,716430,9392085,126044248,1723083930, %T A371372 23910223514,335912566824,4768447532200,68291880722182, %U A371372 985538181002940,14317376105810133,209213540276280758,3073003751985537656,45346188478477675122,671920054584212646330,9993514798883508502188 %N A371372 a(n) = Sum_{d|2*n} binomial(4*n/d-1, 2*n/d)*phi(d)/(4*n) for n>0 with a(0)=0. %C A371372 a(n) is the number of subsets of {1, 2, ..., 4*n-1} of size 2*n that sum to 3*n mod 4*n [Donderwinkel/Kolesnik]. %H A371372 Serte Donderwinkel and Brett Kolesnik, <a href="https://arxiv.org/abs/2403.12941">Asymptotics for Sinaĭ excursions</a>, arXiv:2403.12941 [math.PR], 2024. See Table 1 p. 4. %t A371372 Join[{0}, Table[Sum[Binomial[4*n/d - 1, 2*n/d] * EulerPhi[d] / (4*n), {d, Divisors[2*n]}], {n, 1, 20}]] (* _Vaclav Kotesovec_, Mar 20 2024 *) %o A371372 (PARI) a(n) = if (n==0, 0, sumdiv(2*n, d, binomial(4*n/d-1, 2*n/d)*eulerphi(d))/(4*n)); %o A371372 (Python) %o A371372 from math import comb %o A371372 from sympy import totient, divisors %o A371372 def A371372(n): return sum(comb((d<<1)-1,d)*totient((n<<1)//d) for d in divisors(n<<1,generator=True))//n>>2 if n else 0 # _Chai Wah Wu_, Mar 20 2024 %Y A371372 Cf. A333682. %K A371372 nonn %O A371372 0,3 %A A371372 _Michel Marcus_, Mar 20 2024