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}]

A236375 Positive integers m with 2^(m-1)*phi(m) - 1 prime, where phi(.) is Euler's totient function.

Original entry on oeis.org

3, 7, 12, 15, 18, 31, 42, 108, 124, 140, 143, 155, 207, 327, 386, 463, 514, 823, 925, 1035, 1393, 1425, 2425, 3873, 5091, 5314, 5946, 12813, 14198, 15823, 19932, 22747, 37989, 38772
Offset: 1

Views

Author

Zhi-Wei Sun, Jan 24 2014

Keywords

Comments

According to the conjecture in A236374, this sequence should have infinitely many terms.
The prime 2^(a(34)-1)*phi(a(34)) - 1 = 2^(38771)*12888 - 1 has 11676 decimal digits.

Examples

			a(1) = 3 since neither 2^(1-1)*phi(1) - 1 = 0 nor 2^(2-1)*phi(2) - 1 = 1 is prime, but 2^(3-1)*phi(3) - 1 = 4*2 - 1 = 7 is prime.
		

Crossrefs

Programs

  • Mathematica
    q[m_]:=PrimeQ[2^(m-1)*EulerPhi[m]-1]
    n=0;Do[If[q[m],n=n+1;Print[n," ",m]],{m,1,10000}]
  • PARI
    s=[]; for(m=1, 1000, if(isprime(2^(m-1)*eulerphi(m)-1), s=concat(s, m))); s \\ Colin Barker, Jan 24 2014
Showing 1-2 of 2 results.