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-6 of 6 results.

A234310 Primes of the form 4^k + 4^m - 1, where k and m are positive integers.

Original entry on oeis.org

7, 19, 31, 67, 79, 127, 271, 1039, 1087, 1279, 4099, 4111, 4159, 5119, 8191, 16447, 20479, 65539, 65551, 65599, 81919, 131071, 262147, 262399, 263167, 266239, 524287, 1049599, 1114111, 1310719, 4194319, 4194559, 4195327, 16842751, 17825791, 67108879
Offset: 1

Views

Author

Zhi-Wei Sun, Dec 23 2013

Keywords

Comments

Clearly each term is congruent to 1 modulo 6.
By the conjecture in A234309, this sequence should have infinitely many terms.
Note that any Mersenne prime greater than 3 has the form 2^{2*k+1} - 1 = 4^k + 4^k - 1, where k is a positive integer.

Examples

			a(1) = 7 since 7 = 4^1 + 4^1 - 1 is prime.
a(2) = 19 since 19 = 4^1 + 4^2 - 1 is prime.
a(3) = 31 since 31 = 4^2 + 4^2 - 1 is prime.
		

Crossrefs

Programs

  • Mathematica
    n=0;Do[If[PrimeQ[4^k+4^m-1],n=n+1;Print[n," ",4^m+4^k-1]],{m,1,250},{k,1,m}]
  • PARI
    for(k=1,30,for(m=1,k,if(ispseudoprime(t=4^k+4^m-1),print1(t", ")))) \\ Charles R Greathouse IV, Dec 23 2013

A234475 Number of ways to write n = k + m with 2 < k <= m such that q(phi(k)*phi(m)/4) + 1 is prime, where phi(.) is Euler's totient function and q(.) is the strict partition function (A000009).

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 1, 2, 2, 3, 3, 4, 4, 3, 4, 5, 5, 4, 7, 7, 6, 5, 5, 7, 3, 6, 7, 7, 5, 7, 4, 8, 4, 7, 7, 8, 7, 4, 5, 5, 4, 4, 5, 5, 6, 5, 4, 5, 3, 5, 4, 6, 6, 4, 6, 5, 4, 3, 6, 4, 9, 4, 8, 6, 7, 6, 8, 4, 7, 4, 7, 8, 9, 2, 3, 1, 8, 6, 9, 6, 6, 6, 6, 4, 7, 5, 8, 8, 4, 5, 5, 9, 7, 10, 4, 10, 3, 7, 8, 6
Offset: 1

Views

Author

Zhi-Wei Sun, Dec 26 2013

Keywords

Comments

Conjecture: a(n) > 0 for all n > 5.
This implies that there are infinitely many primes p with p - 1 a term of A000009.

Examples

			a(6) = 1 since 6 = 3 + 3 with q(phi(3)*phi(3)/4) + 1 = q(1) + 1 = 2 prime.
a(76) = 1 since 76 = 18 + 58 with q(phi(18)*phi(58)/4) + 1 = q(42) + 1 = 1427 prime.
a(197) = 1 since 197 = 4 + 193 with q(phi(4)*phi(193)/4) + 1 = q(96) + 1 = 317789.
a(356) = 1 since 356 = 88 + 268 with q(phi(88)*phi(268)/4) + 1 = q(1320) + 1 = 35940172290335689735986241 prime.
		

Crossrefs

Programs

  • Mathematica
    f[n_,k_]:=PartitionsQ[EulerPhi[k]*EulerPhi[n-k]/4]+1
    a[n_]:=Sum[If[PrimeQ[f[n,k]],1,0],{k,3,n/2}]
    Table[a[n],{n,1,100}]

A234530 Primes p with q(p) + 1 also prime, where q(.) is the strict partition function (A000009).

Original entry on oeis.org

2, 3, 11, 13, 29, 37, 47, 71, 79, 89, 103, 127, 131, 179, 181, 197, 233, 271, 331, 379, 499, 677, 691, 757, 887, 911, 1019, 1063, 1123, 1279, 1429, 1531, 1559, 1637, 2251, 2719, 3571, 4007, 4201, 4211, 4297, 4447, 4651, 4967, 5953, 6131, 7937, 8233, 8599, 8819, 9013, 11003, 11093, 11813, 12251, 12889, 12953, 13487, 13687, 15259
Offset: 1

Views

Author

Zhi-Wei Sun, Dec 27 2013

Keywords

Comments

By the conjecture in A234514, this sequence should have infinitely many terms.
It seems that a(n+1) < a(n) + a(n-1) for all n > 4.
See A234366 for primes of the form q(p) + 1 with p prime.
See also A234644 for a similar sequence.

Examples

			a(1) = 2 since 2 and q(2) + 1 = 2 are both prime.
a(2) = 3 since 3 and q(3) + 1 = 3 are both prime.
a(3) = 11 since 11 and q(11) + 1 = 13 are both prime.
		

Crossrefs

Programs

  • Mathematica
    n=0;Do[If[PrimeQ[PartitionsQ[Prime[k]]+1],n=n+1;Print[n," ",Prime[k]]],{k,1,10^5}]
    Select[Prime[Range[2000]],PrimeQ[PartitionsQ[#]+1]&] (* Harvey P. Dale, Apr 23 2017 *)

A233390 a(n) = |{0 < k < n: 2^k - 1 + q(n-k) is prime}|, where q(.) is the strict partition function (A000009).

Original entry on oeis.org

0, 1, 1, 1, 2, 1, 2, 2, 2, 1, 4, 4, 4, 2, 3, 2, 3, 3, 7, 4, 4, 5, 3, 4, 5, 5, 5, 6, 7, 6, 5, 4, 4, 9, 3, 6, 6, 5, 4, 7, 1, 4, 5, 6, 9, 6, 8, 6, 8, 4, 5, 8, 7, 4, 3, 4, 7, 6, 6, 3, 6, 5, 6, 4, 6, 8, 7, 8, 4, 5, 3, 6, 7, 7, 3, 10, 7, 5, 6, 10, 4, 8, 4, 6, 7, 6, 8, 10, 4, 6, 8, 9, 5, 6, 5, 7, 13, 5, 5, 6
Offset: 1

Views

Author

Zhi-Wei Sun, Dec 08 2013

Keywords

Comments

Conjecture: a(n) > 0 for all n > 1.
We have verified this for n up to 150000. For n = 124669, the least positive integer k with 2^k - 1 + q(n-k) prime is 13413.

Examples

			a(6) = 1 since 2^2 - 1 + q(4) = 3 + 2 = 5 is prime.
a(10) = 1 since 2^4 - 1 + q(6) = 15 + 4 = 19 is prime.
a(41) = 1 since 2^{16} - 1 + q(25) = 65535 + 142 = 65677 is prime.
a(127) = 1 since 2^{21} - 1 + q(106) = 2097151 + 728260 = 2825411 is prime.
a(153) = 1 since 2^{70} - 1 + q(83) = 1180591620717411303423 + 101698 = 1180591620717411405121 is prime.
a(164) = 1 since 2^{26} - 1 + q(138) = 67108863 + 8334326 = 75443189 is prime.
		

Crossrefs

Programs

  • Mathematica
    a[n_]:=Sum[If[PrimeQ[2^k-1+PartitionsQ[n-k]],1,0],{k,1,n-1}]
    Table[a[n],{n,1,100}]

A233417 a(n) = |{0 < k <= n/2: q(k)*q(n-k) + 1 is prime}|, where q(.) is the strict partition function (A000009).

Original entry on oeis.org

0, 1, 1, 2, 2, 2, 2, 3, 2, 2, 4, 5, 3, 1, 5, 7, 1, 3, 4, 4, 3, 2, 5, 3, 6, 6, 1, 6, 8, 6, 6, 4, 7, 7, 3, 5, 5, 6, 6, 5, 5, 3, 7, 8, 7, 7, 8, 8, 6, 4, 8, 8, 5, 3, 8, 8, 5, 15, 6, 8, 3, 9, 5, 6, 7, 9, 4, 6, 8, 9, 5, 4, 7, 8, 7, 6, 10, 9, 9, 8, 6, 6, 9, 9, 7, 12, 5, 10, 7, 7, 5, 3, 8, 10, 7, 5, 9, 7, 4, 5
Offset: 1

Views

Author

Zhi-Wei Sun, Dec 09 2013

Keywords

Comments

Conjecture: (i) a(n) > 0 for all n > 1. Similarly, for any integer n > 5, there is a positive integer k < n with q(k)*q(n-k) - 1 prime.
(ii) Let n > 1 be an integer. Then p(k) + q(n-k)^2 is prime for some 0 < k < n, where p(.) is the partition function (A000041). If n is not equal to 8, then k^3 + q(n-k)^2 is prime for some 0 < k < n.

Examples

			a(14) = 1 since q(1)*q(13) + 1 = 1*18 + 1 = 19 is prime.
a(17) = 1 since q(4)*q(13) + 1 = 2*18 + 1 = 37 is prime.
a(27) = 1 since q(13)*q(14) + 1 = 18*22 + 1 = 397 is prime.
		

Crossrefs

Programs

  • Mathematica
    a[n_]:=Sum[If[PrimeQ[PartitionsQ[k]*PartitionsQ[n-k]+1],1,0],{k,1,n/2}]
    Table[a[n],{n,1,100}]

A229835 Number of ways to write n = (p - 1)/6 + q, where p is a prime, and q is a term of the sequence A000009.

Original entry on oeis.org

0, 1, 2, 3, 3, 4, 5, 5, 5, 4, 6, 5, 7, 6, 7, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 6, 7, 10, 9, 6, 8, 6, 10, 8, 9, 7, 7, 10, 10, 9, 8, 7, 10, 7, 10, 3, 7, 12, 8, 10, 6, 8, 9, 6, 10, 8, 11, 7, 11, 8, 7, 9, 8, 12, 10, 8, 12, 7, 9, 10, 10, 8, 11, 10, 7, 10, 9, 14, 9, 9, 9, 8, 10, 10, 9, 7, 8, 9, 9, 8, 10, 9, 10, 10, 9, 7, 8, 7, 12, 8
Offset: 1

Views

Author

Zhi-Wei Sun, Dec 19 2013

Keywords

Comments

Conjecture: a(n) > 0 for all n > 1. Also, any integer n > 1 can be written as (p + 1)/6 + q, where p is a prime and q is a term of A000009.
We have verified this for n up to 2*10^8. Note that 26128189 cannot be written as (p - 1)/4 + q with p a prime and q a term of A000009. Also, 65152682 cannot be written as (p + 1)/4 + q with p a prime and q a term of A000009.

Examples

			a(2) = 1 since 2 = (7 - 1)/ 6 + 1 with 7 prime, and 1 = A000009(i) for i = 0, 1, 2.
a(3) = 2 since 3 = (7 - 1 )/6 + 2 with 7 prime and 2 = A000009(3) = A000009(4), and 3 = (13 - 1 )/6 + 1 with 13 prime and 1 = A000009(i) for i = 0, 1, 2.
		

Crossrefs

Programs

  • Mathematica
    Do[m=0;Do[If[PartitionsQ[k]>=n,Goto[aa]];If[k>1&&PartitionsQ[k]==PartitionsQ[k-1],Goto[bb]];
    If[PrimeQ[6(n-PartitionsQ[k])+1],m=m+1];Label[bb];Continue,{k,1,2n}];
    Label[aa];Print[n," ",m];Continue,{n,1,100}]
Showing 1-6 of 6 results.