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

A236412 a(n) = |{0 < k < n: m = phi(k)/2 + phi(n-k)/8 is an integer with p(m)^2 + q(m)^2 prime}|, where phi(.) is Euler's totient, p(.) is the partition function (A000041) and q(.) is the strict partition function (A000009).

Original entry on oeis.org

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

Views

Author

Zhi-Wei Sun, Jan 24 2014

Keywords

Comments

Conjecture: a(n) > 0 for all n > 17.
We have verified this for n up to 65000.
The conjecture implies that there are infinitely positive integers m with p(m)^2 + q(m)^2 prime. See A236413 for a list of such numbers m. See also A236414 for primes of the form p(m)^2 + q(m)^2.

Examples

			a(15) = 1 since phi(2)/2 + phi(13)/8 = 1/2 + 12/8 = 2 with p(2)^2 + q(2)^2 = 2^2 + 1^2 = 5 prime.
a(69) = 1 since phi(5)/2 + phi(64)/8 = 2 + 4 = 6 with p(6)^2 + q(6)^2 = 11^2 + 4^2 = 137 prime.
a(89) = 1 since phi(73)/2 + phi(16)/8 = 36 + 1 = 37 with p(37)^2 + q(37)^2 = 21637^2 + 760^2 = 468737369 prime.
		

Crossrefs

Programs

  • Mathematica
    p[n_]:=IntegerQ[n]&&PrimeQ[PartitionsP[n]^2+PartitionsQ[n]^2]
    f[n_,k_]:=EulerPhi[k]/2+EulerPhi[n-k]/8
    a[n_]:=Sum[If[p[f[n,k]],1,0],{k,1,n-1}]
    Table[a[n],{n,1,100}]

A236413 Positive integers m with p(m)^2 + q(m)^2 prime, where p(.) is the partition function (A000041) and q(.) is the strict partition function (A000009).

Original entry on oeis.org

1, 2, 3, 4, 6, 17, 24, 37, 44, 95, 121, 162, 165, 247, 263, 601, 714, 742, 762, 804, 1062, 1144, 1149, 1323, 1508, 1755, 1833, 1877, 2330, 2380, 2599, 3313, 3334, 3368, 3376, 3395, 3504, 3688, 3881, 4294, 4598, 4611, 5604, 5696, 5764, 5988, 6552, 7206, 7540, 7689
Offset: 1

Views

Author

Zhi-Wei Sun, Jan 24 2014

Keywords

Comments

According to the conjecture in A236412, this sequence should have infinitely many terms.
See A236414 for primes of the form p(m)^2 + q(m)^2.
See also A236440 for a similar sequence.

Examples

			a(1) = 1 since p(1)^2 + q(1)^2 = 1^2 + 1^2 = 2 is prime.
a(2) = 2 since p(2)^2 + q(2)^2 = 2^2 + 1^2 = 5 is prime.
a(3) = 3 since p(3)^2 + q(3)^2 = 3^2 + 2^2 = 13 is prime.
		

Crossrefs

Programs

  • Mathematica
    pq[n_]:=PrimeQ[PartitionsP[n]^2+PartitionsQ[n]^2]
    n=0;Do[If[pq[m],n=n+1;Print[n," ",m]],{m,1,10000}]
Showing 1-2 of 2 results.