OFFSET
9,3
COMMENTS
From Vladimir Shevelev, Apr 23 2011: (Start)
Also number of non-equivalent necklaces of 9 beads each of them painted by one of n colors.
The sequence solves the so-called Reis problem about convex k-gons in the case k=9 (see our comment to A032279). (End)
REFERENCES
N. Zagaglia Salvi, Ordered partitions and colourings of cycles and necklaces, Bull. Inst. Combin. Appl., 27 (1999), 37-40.
LINKS
Andrew Howroyd, Table of n, a(n) for n = 9..1000
Christian G. Bower, Transforms (2).
Hansraj Gupta, Enumeration of incongruent cyclic k-gons, Indian J. Pure and Appl. Math., 10 (1979), no. 8, 964-999.
F. Ruskey, Necklaces, Lyndon words, De Bruijn sequences, etc. [Cached copy, with permission, pdf format only]
Vladimir Shevelev, Necklaces and convex k-gons, Indian J. Pure and Appl. Math., 35 (2004), no. 5, 629-638.
Vladimir Shevelev, Necklaces and convex k-gons, Indian J. Pure and Appl. Math., 35 (2004), no. 5, 629-638.
Vladimir Shevelev, Spectrum of permanent's values and its extremal magnitudes in Lambda_n^3 and Lambda_n(alpha,beta,gamma), arXiv:1104.4051 [math.CO], 2011. (Cf. Section 5).
Index entries for linear recurrences with constant coefficients, signature (2,3,-6,-6,6,13,-2,-18,-1,11,3,0,0,-3,-11,1,18,2,-13,-6,6,6,-3,-2,1).
FORMULA
"DIK[ 9 ]" (necklace, indistinct, unlabeled, 9 parts) transform of 1, 1, 1, 1...
Put s(n,k,d) = 1, if n == k (mod d), and s(n,k,d) = 0, otherwise. Then a(n) =(1/3)*s(n,0,9) + (n-3)*(n-6)*s(n,0,3)/162 + (n-2)(n-4)*(n-6)*(n-8)*(945 + (n-1)*(n-3)*(n-5)*(n-7))/725760, if n is even; a(n) = (1/3)*s(n,0,9) + (n-3)*(n-6)*s(n,0,3)/162 +(n-1)*(n-3)*(n-5)*(n-7)*(945 + (n-2)*(n-4)*(n-6)*(n-8))/725760, if n is odd. - Vladimir Shevelev, Apr 23 2011
From Herbert Kociemba, Nov 05 2016: (Start)
G.f.: (1/2)*x^9*((1+x)/(1-x^2)^5 + 1/9*(1/(1-x)^9 - 2/(-1+x^3)^3 - 6/(-1+x^9))).
G.f.: k=9, x^k*((1/k)*Sum_{d|k} phi(d)*(1-x^d)^(-k/d) + (1+x)/(1-x^2)^floor((k+2)/2)/2. [edited by Petros Hadjicostas, Jul 18 2018] (End)
MATHEMATICA
k = 9; Table[(Apply[Plus, Map[EulerPhi[ # ]Binomial[n/#, k/# ] &, Divisors[GCD[n, k]]]]/n + Binomial[If[OddQ[n], n - 1, n - If[OddQ[k], 2, 0]]/2, If[OddQ[k], k - 1, k]/2])/2, {n, k, 50}] (* Robert A. Russell, Sep 27 2004 *)
k=9; CoefficientList[Series[x^k*(1/k Plus@@(EulerPhi[#] (1-x^#)^(-(k/#))&/@Divisors[k])+(1+x)/(1-x^2)^Floor[(k+2)/2])/2, {x, 0, 50}], x] (* Herbert Kociemba, Nov 04 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved