OFFSET
3,1
COMMENTS
These sequences arise from a bijection between bracelet colorings and unicyclic graphs, where each unlabeled rooted tree of size n represents a distinct color. Hence the total number of available colors is N = A000081(n).
Let a(n) be the number of bracelets of length 4 using exactly 2 colors, where the colors are chosen from a set of size N = A000081(n) (unlabeled rooted trees with n nodes). For a fixed set of 2 colors, the possible color distributions are (3,1) and (2,2). The (3,1) pattern gives 2 distinct bracelets (the single bead can be adjacent or opposite), and the (2,2) pattern gives 2 distinct bracelets (adjacent pairs or opposite pairs). Hence there are 2 + 2 = 4 distinct bracelets per color set. Therefore a(n) = 4 * binomial(N,2). Offset 3,1 because n >= 3.
EXAMPLE
For n=3, N = A000081(3) = 2 distinct colors, say X and Y. The only unordered pair of colors is {X, Y}. The 4 distinct bracelets of length 4 using exactly these two colors are:
X X X Y (single Y adjacent)
X X Y X (single Y opposite)
X X Y Y (adjacent pairs)
X Y X Y (opposite pairs)
All rotations and reflections are considered identical. Hence a(3) = 4.
For n=4, N = A000081(4) = 4 colors (A, B, C, D). Choosing any unordered pair of distinct colors, say A and B, there are again 4 bracelets as listed above (with X,Y replaced by A,B). Since there are C(4,2) = 6 such color pairs, the total number is 6 * 4 = 24, matching a(4) = 4 * C(4,2) = 24.
For n=5, N = A000081(5) = 9 colors. The number of unordered color pairs is C(9,2) = 36, and each gives 4 bracelets, so a(5) = 36 * 4 = 144.
PROG
(Python)
from math import comb
from functools import lru_cache
@lru_cache(maxsize=None)
def A000081(n):
if n <= 1:
return 1
total = 0
for j in range(1, n):
s = 0
for d in range(1, j + 1):
if j % d == 0:
s += d * A000081(d)
total += s * A000081(n - j)
return total // (n - 1)
def a(n):
return 4 * comb(A000081(n), 2)
print([a(n) for n in range(3, 27)])
CROSSREFS
KEYWORD
nonn
AUTHOR
Chen Wei-Shi, May 29 2026
STATUS
approved
