A233417 a(n) = |{0 < k <= n/2: q(k)*q(n-k) + 1 is prime}|, where q(.) is the strict partition function (A000009).
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
Keywords
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.
Links
- Zhi-Wei Sun, Table of n, a(n) for n = 1..10000
- Z.-W. Sun, On a^n+ bn modulo m, arXiv preprint arXiv:1312.1166 [math.NT], 2013-2014.
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}]
Comments