OFFSET
1,1
COMMENTS
Each solution (a,b,c,d) corresponds to a Hadamard matrix of quaternion type H = [[A, B, C, D], [-B, A, -D, C], [-C, D, A, -B], [-D, -C, B, A]], where A and D are circulant matrices formed by a and d, respectively, and B=fliplr(circulant(b)) and C=fliplr(circulant(c)). The converse is not always true. To see this, set a=(-1, -1, -1, 1), b=(-1, -1, -1, 1), c=(-1, 1, 1, 1) and d=(1, -1, -1, -1). Then H is Hadamard but |dft(a)^2 + dft(d)^2| + |dft(b)^2 + dft(c)^2| = (16, 0, 16, 0).
16 is a divisor of a(n), for all n. If (a,b,c,d) is a solution, then each of the 16 tuples ((+-)a, (+-)b, (+-)c, (+-)d) is also a solution.
LINKS
L. D. Baumert and M. Hall, Hadamard matrices of the Williamson type, Math. Comp. 19:91 (1965) 442-447.
W. H. Holzmann, H. Kharaghani and B. Tayfeh-Rezaie, Williamson matrices up to order 59, Des. Codes Cryptogr. 46 (2008), 343-352.
Jeffery Kline, A complete list of solutions (a,b,c,d), for 1<=n<=5.
Jeffery Kline, Geometric Search for Hadamard Matrices, Theoret. Comput. Sci. 778 (2019), 33-46.
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Jeffery Kline, Dec 20 2018
STATUS
approved