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 A367554 #13 Nov 23 2023 15:15:48 %S A367554 1,1,13,100,578,2530,8866,25300,60115,120175,204347,297160,371516, %T A367554 400024,371516,297160,204347,120175,60115,25300,8866,2530,578,100,13, %U A367554 1,1 %N A367554 a(n) is the number of 2n-regular circulant graphs of order 53. %H A367554 Brian Alspach and Marni Mishna, <a href="https://doi.org/10.1016/S0012-365X(02)00319-9">Enumeration of Cayley graphs and digraphs</a>, Discr. Math., 256 (2002), 527-539. %H A367554 Marni Mishna, <a href="https://www.sfu.ca/~mmishna/">Home page</a>. %H A367554 Marni Mishna, <a href="https://www.sfu.ca/~mmishna/research.html">Publications</a>. %H A367554 Marni Mishna, <a href="https://www.sfu.ca/~mmishna/PUB/sfu-thesis.pdf">Cayley Graph Enumeration</a>, Master's Thesis, Simon Fraser University, 2000. See p. 16 (which is p. 24 in the pdf). %F A367554 Sum_n a(n) = 2581428 = A049287(53) = A285620(53) = A000031((53-1)/2). - _Andrey Zabolotskiy_, Nov 22 2023 %o A367554 (SageMath) %o A367554 def a(k, p): %o A367554 return (2/(p-1)) * sum(euler_phi(d) * binomial((p-1)/(2*d), k/(2*d)) for d in divisors(gcd(k, p-1)/2)) # see Mishna; beware the missing prefactor (2/(p-1)) %o A367554 print([a(2*n, 53) for n in range(27)]) # _Andrey Zabolotskiy_, Nov 22 2023 %o A367554 (Python) %o A367554 from math import gcd, comb %o A367554 from sympy import totient, divisors %o A367554 def A367554(n): return sum(totient(d)*comb(26//d,n//d) for d in divisors(gcd(n,26),generator=True))//26 # _Chai Wah Wu_, Nov 23 2023 %Y A367554 Cf. A000031, A049287, A285620. %K A367554 nonn,fini,full %O A367554 0,3 %A A367554 _Andrey Zabolotskiy_, Nov 22 2023