OFFSET
1,4
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..200
EXAMPLE
Case n=4: there are the following 7 bracelets:
1212,
1213, 1232, 1323,
1234, 1243, 1324.
PROG
(PARI) \\ here U(n, k) is A208544(n, k) for n > 1.
U(n, k) = (sumdiv(n, d, eulerphi(n/d)*(k-1)^d)/n + if(n%2, 1-k, k*(k-1)^(n/2)/2))/2;
a(n)={if(n<1, n==0, sum(j=1, n, U(n, j)*sum(k=j, n, (-1)^(k-j)*binomial(k, j))))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrew Howroyd, Dec 20 2019
STATUS
approved