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.

Showing 1-10 of 15 results. Next

A228552 Square root of the absolute value of A069191(n).

Original entry on oeis.org

1, 1, 1, 0, 1, 1, 2, 2, 0, 0, 0, 0, 0, 0, 0, 2, 2, 3, 5, 11, 8, 24, 48, 60, 56, 16, 12, 31, 155, 217, 588, 1148, 328, 164, 176, 132, 176, 395, 277, 697, 692, 191, 915, 76, 22742, 125664, 128079, 213885, 7371, 171654, 89678, 114902, 149465, 353497, 144573, 388325, 198676, 1738118, 1311164, 222898
Offset: 1

Views

Author

Zhi-Wei Sun, Aug 25 2013

Keywords

Comments

According to the comments of A069191, a(n) should be always integral. Note that a(2*n) is the absolute value of A228616(n) by the comments of A228591. We conjecture that a(n) > 0 for all n > 15.

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=Sqrt[Abs[Det[Table[If[PrimeQ[i+j]==True,1,0],{i,1,n},{j,1,n}]]]]
    Table[a[n],{n,1,20}]

A070223 Square roots of A069191, sorted.

Original entry on oeis.org

0, 1, 2, 3, 5, 8, 11, 12, 16, 24, 31, 48, 56, 60, 76, 132, 155, 164, 176, 191, 217, 277, 328, 395, 588, 692, 697, 915, 1148, 7371, 21100, 22742, 89678, 114902, 125664, 128079, 144573, 149465, 171654, 198676, 213885, 222898, 248633, 353497, 382075
Offset: 1

Views

Author

Robert G. Wilson v, May 25 2002

Keywords

Crossrefs

Cf. A069191.

A073364 Number of permutations p of (1,2,3,...,n) such that k+p(k) is prime for 1<=k<=n.

Original entry on oeis.org

1, 1, 1, 4, 1, 9, 4, 36, 36, 676, 400, 9216, 3600, 44100, 36100, 1223236, 583696, 14130081, 5461569, 158180929, 96275344, 5486661184, 2454013444, 179677645456, 108938283364, 5446753133584, 4551557699844, 280114147765321, 125264064932449, 9967796169000201
Offset: 1

Views

Author

Benoit Cloitre, Aug 23 2002

Keywords

Comments

a(n)=permanent(m), where the n X n matrix m is defined by m(i,j) = 1 or 0, depending on whether i+j is prime or composite respectively. - T. D. Noe, Oct 16 2007

Crossrefs

Programs

  • Haskell
    a073364 n = length $ filter (all isprime)
                         $ map (zipWith (+) [1..n]) (permutations [1..n])
       where isprime n = a010051 n == 1  -- cf. A010051
    -- Reinhard Zumkeller, Mar 19 2011
  • Mathematica
    am[n_] := Permanent[Array[Boole[PrimeQ[2 #1 + 2 #2 - 1]]&, {n, n}]];
    ap[n_] := Permanent[Array[Boole[PrimeQ[2 #1 + 2 #2 + 1]]&, {n, n}]];
    a[n_] := If[n == 1, 1, If[EvenQ[n], am[n/2]^2, ap[(n-1)/2]^2]];
    Array[a, 28] (* Jean-François Alcover, Nov 03 2018 *)
  • PARI
    a(n)=sum(k=1,n!,n==sum(i=1,n,isprime(i+component(numtoperm(n,k),i))))
    
  • PARI
    a(n)={matpermanent(matrix(n,n,i,j,isprime(i + j)))} \\ Andrew Howroyd, Nov 03 2018
    

Formula

a(2n) = A000341(n)^2 and a(2n+1) = A134293(n)^2. - T. D. Noe, Oct 16 2007

Extensions

a(10) from Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Mar 14 2004
a(11) from Rick L. Shepherd, Mar 17 2004
a(12)-a(17) from John W. Layman, Jul 21 2004
More terms from T. D. Noe, Oct 16 2007

A228591 Determinant of the n X n (0,1)-matrix with (i,j)-entry equal to 1 if and only if i + j is 2 or an odd composite number.

Original entry on oeis.org

1, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, -1, -9, 81, 9, -1225, -2500, 2500, 2500, -225, -121, 841, 19044, -29584, -355216, 1527696, 141376, -40000, -40000, 10000, 59536, -258064, -139876, 935089, 885481, -16384, -1876900, 1710864, 818875456, -22896531856, -23799232900, 66328911936, 158281561, -45320023225
Offset: 1

Views

Author

Zhi-Wei Sun, Aug 27 2013

Keywords

Comments

Conjecture: a(n) = 0 for no n > 15.
We observe that (-1)^{n*(n-1)/2}*a(n) is always a square. This is a special case of the following general result established by Zhi-Wei Sun.
Theorem: Let M = (m_{i,j}) be an n X n symmetric matrix over a commutative ring. Suppose that the (i,j)-entry m_{i,j} is zero whenever i + j is even and greater than 2. If n is even, then (-1)^{n/2}*det(M) = D(n)^2, where D(n) denotes the determinant |m_{2i,2j-1}|{i,j = 1,...,n/2}. If n is odd, then (-1)^{(n-1)/2}*det(M) = m{1,1}*D(n)^2, where D(n) is the determinant |m_{2i,2j+1}|_{i,j = 1,...,(n-1)/2}.
This theorem extends the result mentioned in A069191.

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=Det[Table[If[(i+j==2)||(Mod[i+j,2]==1&&PrimeQ[i+j]==False),1,0],{i,1,n},{j,1,n}]]
    Table[a[n],{n,1,50}]

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

Original entry on oeis.org

1, -1, -1, -1, 0, 0, -1, 1, 1, 1, -1, 0, 0, 0, -1, -1, 1, 1, 1, -1, 2, 8, -18, -9, -1, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, -4096, -4096, 64, -20, -125, 5, -6, -216, 24, 152, 54872, -106742, 14045, 125, -21125, -274625, -274625, 10985, -16731, -970299, 1804275, 1312200, 373248, -691488, -192080
Offset: 1

Views

Author

Zhi-Wei Sun, Aug 27 2013

Keywords

Comments

Conjecture: a(n) is nonzero for any n > 35.
Clearly this conjecture implies the twin prime conjecture.

Examples

			a(1) = 1 since 2*(1 + 1) - 1 = 3 and 2*(1 + 1) + 1 = 5 are twin primes.
		

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=Det[Table[If[PrimeQ[2(i+j)-1]&&PrimeQ[2(i+j)+1],1,0],{i,1,n},{j,1,n}]]
    Table[a[n],{n,1,20}]
    Table[Det[Table[If[AllTrue[2(i+j)+{1,-1},PrimeQ],1,0],{i,k},{j,k}]],{k,60}] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Sep 21 2019 *)

A228557 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.

Original entry on oeis.org

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, 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
Offset: 1

Views

Author

Zhi-Wei Sun, Aug 25 2013

Keywords

Comments

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.
By the general result mentioned in A228591, (-1)^n*a(2*n) equals the square of A228615(n).
Zhi-Wei Sun made the following general conjecture:
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.
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.

Examples

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

Crossrefs

Programs

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

A228616 Determinant of the n X n matrix with (i,j)-entry equal to 1 or 0 according as 2*(i + j) - 1 is prime or not.

Original entry on oeis.org

1, 0, -1, -2, 0, 0, 0, -2, 3, 11, 24, -60, -16, 31, -217, -1148, -164, 132, 395, 697, -191, -76, 125664, -213885, -171654, 114902, 353497, -388325, -1738118, -222898, 248633, 382075, -1637075, 21100, 4049068, -189147708, -279083383, 472023163, -19063401, -6718578823
Offset: 1

Views

Author

Zhi-Wei Sun, Aug 27 2013

Keywords

Comments

Conjecture: a(n) is nonzero for any n > 7.
Clearly this conjecture implies that there are infinitely many primes.

Examples

			 a(1) = 1 since 2*(1 + 1) - 1 = 3 is a prime.
		

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=Det[Table[If[PrimeQ[2(i+j)-1],1,0],{i,1,n},{j,1,n}]]
    Table[a[n],{n,1,40}]

A228559 Determinant of the n X n matrix with (i,j)-entry equal to 1 or 0 according as i + j is a Sophie Germain prime or not.

Original entry on oeis.org

1, -1, -1, 1, 0, -1, -1, 1, 1, -1, 0, 1, 1, -1, -4, 16, 0, -64, -64, 64, 0, 0, 0, 0, 0, -64, -64, 64, 0, -16, -4, 1, 1, -1, 0, 4, 16, -64, -144, 324, 0, -81, -9, 1, 4, -16, 0, 64, 64, -64, 0, 0, 0, 0, 0, 262144, 4194304, -67108864, 0, 1073741824
Offset: 1

Views

Author

Zhi-Wei Sun, Aug 25 2013

Keywords

Comments

If p > 3 and 2*p+1 are both prime, then p == -1 (mod 6). If tau is a permutation of {1,...,n}, and i + tau(i) is a Sophie Germain prime for each i = 1,...,n, then n*(n+1) = sum_{i=1}^n (i + tau(i)) is congruent to -n or 2 - (n - 1) or 3 - (n - 1) or 3 + 3 - (n - 2) modulo 6, which is impossible when n == -1 (mod 6). Therefore a(6*n-1) = 0 for all n > 0.
Note also that (-1)^{n*(n-1)/2}*a(n) is always a square in view of the comments in A228591.
Conjecture: a(n) is nonzero if n is greater than 55 and not congruent to 5 modulo 6.
Zhi-Wei Sun also had some similar conjectures. For example, if b(n) denotes the n X n determinant with (i,j)-entry equal to 1 or 0 according as i + j and 2*(i + j) - 1 are both prime or not, then b(n) is nonzero when n is greater than 125 and not congruent to 3 modulo 6. (Just like a(6*n-1) = 0, we can show that b(6*n-3) = 0.)

Examples

			a(1) = 1 since 1 + 1 is a Sophie Germain prime.
		

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=Det[Table[If[PrimeQ[i+j]==True&&PrimeQ[2(i+j)+1]==True,1,0],{i,1,n},{j,1,n}]]
    Table[a[n],{n,1,20}]
    Det/@Table[If[AllTrue[{i+j,2(i+j)+1},PrimeQ],1,0],{n,60},{i,n},{j,n}] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, May 26 2017 *)

A228561 Determinant of the n X n matrix with (i,j)-entry equal to 1 or 0 according as i + j and 4*(i + j)^2 + 1 are both prime or not.

Original entry on oeis.org

1, -1, -1, 0, 1, -1, -1, 0, 1, -1, -1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, -1, -1, 0, 4, -16, -9, 25, 4, -81, -81, 81, 841, -5929, -3969, 19600, 69169, -667489, -285156, 80656, 276676, -790321, -60025, 3136, 10816, -40000, -45369, 221841, 86436, -168100, -12100, 13225, 11881, -87616, -71289, 729
Offset: 1

Views

Author

Zhi-Wei Sun, Aug 25 2013

Keywords

Comments

Conjecture: a(n) is nonzero for each n > 28.
This implies that there are infinitely many primes p with 4*p^2 + 1 also prime. Note also that (-1)^{n*(n-1)/2}*a(n) is always a square in view of the comments in A228591.

Examples

			a(1) = 1 since  1 + 1 = 2 and 4*2^2 + 1 = 17 are both prime.
		

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=Det[Table[If[PrimeQ[i+j]==True&&PrimeQ[4(i+j)^2+1]==True,1,0],{i,1,n},{j,1,n}]]
    Table[a[n],{n,1,30}]

A228574 Determinant of the 2*n X 2*n matrix with (i,j)-entry equal to 1 or 0 according as i + j is a prime congruent to 1 mod 4 or not.

Original entry on oeis.org

0, 1, 0, 1, 0, 16, 0, 1, 0, 1, 0, 6561, 0, 0, 0, 0, 0, 0, 0, 6561, 0, 456976, 0, 65536, 0, 84934656, 0, 12745506816, 0, 335563778560000, 0, 1105346784523536, 0, 441194850625, 0, 986262467993856, 0, 80385880645971214336, 0, 6387622009837971841
Offset: 1

Views

Author

Zhi-Wei Sun, Aug 25 2013

Keywords

Comments

For the (2*n-1) X (2*n-1) determinant with (i,j)-entry equal to 1 or 0 according as i + j is a prime congruent to 1 mod 4 or not, it is easy to see that it vanishes since sum_{i=1}^{2*n-1} (i + tau(i) - 1) is not a multiple of 4 for any permutation tau of {1,...,2n-1}.
Conjecture: a(2*n-1) = 0 for all n > 0, and a(2*n) is nonzero when n > 9.
Zhi-Wei Sun could prove the following related result:
Let m be any positive even integer, and let D(m, n) denote the n X n determinant with (i,j)-entry equal to 1 or 0 according as i + j is a prime congruent to 1 mod m or not. Then (-1)^{n*(n-1)/2}*D(m,n) is always an m-th power. (It is easy to see that D(m,n) = 0 if m does not divide n^2.)

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=Det[Table[If[Mod[i+j,4]==1&&PrimeQ[i+j]==True,1,0],{i,1,2n},{j,1,2n}]]
    Table[a[n],{n,1,20}]
Showing 1-10 of 15 results. Next