OFFSET
1,4
COMMENTS
Define a(1) = 0. To calculate a(n):
1. Expand (A + B)^n into 2^n words of length n consisting of letters A and B (i.e., use of the distributive and associative laws of multiplication but assume A and B do not commute).
2. To each of the 2^n words, associate a free binary necklace consisting of n "black and white pearls". Figuratively, all 2^n necklaces can be placed inside a treasure chest.
3. Remove all n-pearled necklaces which are found to have (at least) two adjacent white pearls from the chest.
4. If two necklaces are found to be equivalent, remove one of them from the chest. Continue until no two equivalent necklaces can be found in the chest.
5. Counting the total number of white pearls left in the chest gives a(n).
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..2000 (first 50 terms from Max Alekseyev)
Creighton Dement and Max Alekseyev, Notes on A113166
FORMULA
a(n) = Sum_{k=1..floor(n/2)} (k/(n-k))*Sum_{j=1..gcd(n,k)} binomial((n-k)*gcd(n,k,j)/gcd(n,k), k*gcd(n,k,j)/gcd(n,k)) (Alekseyev).
a(p) = Fibonacci(p-1) for all primes p. (Creighton Dement and Antti Karttunen, proved by Max Alekseyev).
a(n) = Sum_{d|n} phi(n/d)*Fibonacci(d-1), where phi=A000010. - Maxim Karimov and Vladislav Sulima, Aug 20 2021
MAPLE
with(numtheory): with(combinat):
a:= n-> add(phi(d)*fibonacci(n/d-1), d=divisors(n)):
seq(a(n), n=1..50); # Alois P. Heinz, Aug 21 2021
MATHEMATICA
a[n_] := Sum[EulerPhi[d]*Fibonacci[n/d - 1], {d, Divisors[n]}];
Array[a, 50] (* Jean-François Alcover, Jan 03 2022 *)
PROG
(PARI) A113166(n) = sum(k=1, n\2, k/(n-k) * sum(j=1, gcd(n, k), binomial((n-k)*gcd([n, k, j])/gcd(n, k), k*gcd([n, k, j])/gcd(n, k)) ))
(MATLAB)
function [res] = calcA113166(n)
d=divisors(n);
res=0;
for i=1:length(d)
res=res+eulerPhi(n/d(i))*fibonacci(d(i)-1);
end
end
% Maxim Karimov, Aug 21 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Creighton Dement, Jan 05 2006; Jan 08 2006; Jul 29 2006
EXTENSIONS
More terms from Max Alekseyev, Jun 20 2006
STATUS
approved