login
Determinant of the n X n matrix with (i,j)-entry equal to 1 or 0 according as i + j and i + j + 2 are twin primes or not.
9

%I #28 Aug 28 2013 02:57:29

%S 0,-1,0,1,0,-1,0,1,0,0,0,0,0,-1,0,1,0,-1,0,1,0,-1,0,0,0,0,0,0,0,-1,0,

%T 1,0,-1,0,1,0,-1,0,1,0,-4,0,64,0,-324,0,81,0,-1,0,16,0,0,0,0,0,0,0,0

%N Determinant of the n X n matrix with (i,j)-entry equal to 1 or 0 according as i + j and i + j + 2 are twin primes or not.

%C Clearly p is odd if p and p + 2 are twin primes. If sigma is a permutation of {1,...,n}, and i + sigma (i) and i + sigma(i) + 2 are twin primes for all i = 1,...,n, then we must have sum_{i=1}^n (i + sigma(i)) == n (mod 2) and hence n is even. Therefore a(n) = 0 if n is odd.

%C By the general result mentioned in A228591, (-1)^n*a(2*n) equals the square of A228615(n).

%C Zhi-Wei Sun made the following general conjecture:

%C Let d be any positive even integer, and let D(d,n) be the n X n determinant with (i,j)-entry eual to 1 or 0 according as i + j and i + j + d are both prime or not. Then D(d,2*n) is nonzero for large n.

%C Note that when n is odd we have D(d,n) = 0 (just like a(n) = 0). Also, the conjecture implies de Polignac's conjecture that there are infinitely many primes p such that p and p + d are both prime.

%H Zhi-Wei Sun, <a href="/A228557/b228557.txt">Table of n, a(n) for n = 1..300</a>

%e a(1) = 0 since {2, 4} is not a twin prime pair.

%t a[n_]:=a[n]=Det[Table[If[PrimeQ[i+j]==True&&PrimeQ[i+j+2]==True,1,0],{i,1,n},{j,1,n}]]

%t Table[a[n],{n,1,100}]

%Y Cf. A001359, A006512, A228615, A069191, A071524, A228591, A228552, A228548, A228549, A228559, A228561, A228574, A228578.

%K sign

%O 1,42

%A _Zhi-Wei Sun_, Aug 25 2013