%I #2 Mar 31 2012 13:20:27
%S 3,7,147,301,33411,1748509,36718689,4198170109,709490748421,
%T 82402282638039,1345903949754637,1564158644309443,855594778437265321,
%U 5136411178193150947,3703352459477261832787,261798531558431048025481
%N Numerator of determinant of n X n matrix with elements M[i,j] = (p^2 - p + 1)/(p*(p-1)) if i=j and 1 otherwise, where p=Prime[i].
%C All square prime divisors of a(n) {7,13,43,139,19,31,61,37,607,523,67,79,1201,241,1171,157,109,...} belong to A002476[n] Primes of form 6n + 1.
%t Numerator[ Table[ Det[ DiagonalMatrix[ Table[1/(Prime[i]*(Prime[i]-1)), {i, 1, n} ] + 1 ]], {n, 1, 150}]]
%Y Cf. A002476, A036689, A002061.
%K nonn
%O 1,1
%A _Alexander Adamchuk_, Aug 02 2006
|