OFFSET
1,1
COMMENTS
Let M(n) denote the n X n matrix m(i,j)=0 if n divides ij-1, m(i,j) = 1 otherwise then det(M(n))=+1 if and only if n is prime ==3 (mod 4).
a(n) * A151763(n) = - a(n).
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
FORMULA
PROG
(PARI) { a(n)=isprime(n)*if(n%4-3, 0, 1) }; vector(100, n, a(n))
(Haskell)
a079261 n = fromEnum $ n `mod` 4 == 3 && a010051 n == 1
-- Reinhard Zumkeller, Oct 06 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Feb 04 2003
STATUS
approved