A198255 Number of ways to write n as the sum of two coprime squarefree semiprimes.
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 2, 0, 0, 0, 2, 0, 2, 0, 0, 0, 2, 0, 1, 0, 1, 0, 2, 0, 5, 0, 0, 0, 1, 0, 3, 1, 2, 0, 3, 0, 4, 1, 0, 1, 3, 0, 5, 0, 2, 0, 4, 0, 1, 2, 2, 0, 3, 0, 4, 2, 2, 1, 3, 0, 6, 1, 2, 0, 6
Offset: 1
Examples
Same as A197629 until a(97) since 97=(2)(3)(7)+(5)(11) and thus the value of the 97th term of A197629 is one greater. The 248th term of A197629 is the first even term which is one greater since 248=(3)(5)(7)+(11)(13). From _Antti Karttunen_, Sep 23 2018: (Start) For n = 97 there are following six solutions: 97 = 6+91 = 10+87 = 15+82 = 35+62 = 39+58 = 46+51, thus a(97) = 6. For n = 248 there are following seven solutions: 248 = 33+215 = 35+213 = 39+209 = 65+183 = 87+161 = 115+133 = 119+129, thus a(248) = 7. (End)
Links
- Antti Karttunen, Table of n, a(n) for n = 1..65537
Crossrefs
Cf. A197629.
Programs
-
MATLAB
function [asubn] = sps(n) % Returns the number of coprime, squarefree, semiprime, partitions a+b of n. r = 0; % r is the number of sps's of n k=6; while k < n/2, if gcd(k,n-k)==1 if length(factor(k)) == 2 if length(factor(n-k)) == 2 if prod(diff(factor(k)))*prod(diff(factor(n-k))) > 0 r = r + 1; end end end end k = k + 1; end asubn = r; end
-
PARI
A198255(n) = sum(k=4, (n-1)\2, gcd(k, n-k)==1&&(2==bigomega(k))&&(2==bigomega(n-k))&&issquarefree(k)&&issquarefree(n-k)); \\ Antti Karttunen, Sep 23 2018, after Charles R Greathouse IV's program for A197629
Comments