login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A306595
Determinant of the circulant matrix whose first column corresponds to the binary digits of n.
2
0, 1, 1, 0, 1, 2, 2, 0, 1, 0, 0, 3, 0, -3, 3, 0, 1, 2, 2, 3, 2, 3, 3, 4, 2, 3, 3, 4, 3, 4, 4, 0, 1, 0, 4, 0, 0, -9, 9, 0, 4, 9, 0, 8, 9, 0, 8, 5, 0, 0, 9, 0, -9, -8, 0, -5, 0, 0, 8, 5, 0, -5, 5, 0, 1, 2, 2, 3, 2, 24, 24, 4, 2, 3, 3, 32, 3, 4, 32, 5, 2, 24, 3
OFFSET
0,6
COMMENTS
This sequence is the binary variant of A177894.
From Robert Israel, Mar 05 2019: (Start)
a(n) is divisible by A000120(n).
If A070939(n) is even then n is divisible by A000120(n)*A065359(n). (End)
FORMULA
a(A121016(n)) = 0 for any n > 0.
a(2^k) = 1 for any k >= 0.
a(A219325(n)) = A219325(n) for any n > 0.
EXAMPLE
For n = 13:
- the binary representation of 13 is "1101",
- the corresponding circulant matrix is:
[1 1 0 1]
[1 1 1 0]
[0 1 1 1]
[1 0 1 1]
- its determinant is -3,
- hence a(13) = -3.
MAPLE
a:= n-> `if`(n=1, 1, (l-> LinearAlgebra[Determinant](Matrix(nops(l),
shape=Circulant[l[-i]$i=1..nops(l)])))(convert(n, base, 2))):
seq(a(n), n=0..100); # Alois P. Heinz, Mar 05 2019
PROG
(PARI) a(n) = my (d=if (n, binary(n), [0])); my (m=matrix(#d, #d, i, j, d[1+(i-j)%#d])); return (matdet(m))
KEYWORD
sign,base,look
AUTHOR
Rémy Sigrist, Feb 27 2019
STATUS
approved