cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A118770 Determinant of n X n matrix containing the first n^2 semiprimes in increasing order.

This page as a plain text file.
%I A118770 #17 Nov 21 2018 13:45:58
%S A118770 4,-14,-196,480,696,-57901,-525364,-409579,18528507,-237549252,
%T A118770 -2119519900,6713972874,18262155072,-19072020914992,162234208372185,
%U A118770 1471912942112734,6828673030820538,-35126752028893500,729026655790306778,-15365360727898374618
%N A118770 Determinant of n X n matrix containing the first n^2 semiprimes in increasing order.
%C A118770 Semiprime analog of A067276 Determinant of n X n matrix containing the first n^2 primes in increasing order. The first column contains the first n semiprimes in increasing order, the second column contains the next n semiprimes in increasing order, etc. Equivalently, first row contains first n semiprimes in increasing order, second row contains next n semiprimes in increasing order, etc. See also: A118713 a(n) = determinant of n X n circulant matrix whose first row is A001358(1), A001358(2), ..., A001358(n) where A001358(n) = n-th semiprime.
%H A118770 Harvey P. Dale, <a href="/A118770/b118770.txt">Table of n, a(n) for n = 1..570</a>
%e A118770 a(2) = -14 because of the determinant -14 =
%e A118770 |4,6 |
%e A118770 |9,10|.
%e A118770 a(6) = -57901 = the determinant
%e A118770 |4, 6, 9, 10, 14, 15,|
%e A118770 |21, 22, 25, 26, 33, 34,|
%e A118770 |35, 38, 39, 46, 49, 51,|
%e A118770 |55, 57, 58, 62, 65, 69,|
%e A118770 |74, 77, 82, 85, 86, 87,|
%e A118770 |91, 93, 94, 95, 106, 111|.
%t A118770 SemiPrimePi[ n_ ] := Sum[ PrimePi[ n/Prime @ i ] - i + 1, {i, PrimePi @ Sqrt @ n} ]; SemiPrime[ n_ ] := Block[ {e = Floor[ Log[ 2, n ] + 1 ], a, b}, a = 2^e; Do[ b = 2^p; While[ SemiPrimePi[ a ] < n, a = a + b ]; a = a - b/2, {p, e, 0, -1} ]; a + b/2 ]; Table[ Det[ Partition[ Array[ SemiPrime, n^2 ], n ] ], {n, 20} ] (* _Robert G. Wilson v_, May 26 2006 *)
%t A118770 Module[{nn=5000,spr},spr=Select[Range[nn],PrimeOmega[#]==2&];Table[Det[ Partition[ Take[spr,n^2],n]],{n,Sqrt[Length[spr]]}]] (* _Harvey P. Dale_, Nov 21 2018 *)
%Y A118770 Cf. A001358, A067276, A118713, A118772, A118779.
%K A118770 easy,sign
%O A118770 1,1
%A A118770 _Jonathan Vos Post_, May 22 2006
%E A118770 More terms from _Robert G. Wilson v_, May 26 2006
%E A118770 Typos in Mma program corrected by _Giovanni Resta_, Jun 12 2016