By standard s-m-n arguments, there exists some total function f such that phif(x)(y)=theta(x,y). But note that phif(x) is either the constant function epsilon (when x not in not K), in which case f(x) not in S(c), or the totally undefined function (when x in not K), in which case f(x) in S(c). Thus we have a reduction from not K to S(c) such that x in not K <=> f(x) in S(c). (Incidentally, this reduction shows that S(c) is non-r.e., not surprisingly, since every productive set is by definition non-r.e.)
Now note that, if we have dom phii subset S(c), then we have
(We use f-1(A) to denote the inverse image of A under f, i.e., the set of all x such that f(x) in A; we do not use f-1 to denote a function of any kind!) Moreover, f-1(dom phii) is r.e., so that there exists some phiz such that dom phiz = f-1(dom phii); to see that, just note that phiz(y) = phii(f(y)) has the desired properties. But then dom phiz subset not K and so, by our previous remarks, z in not K-dom phiz. It follows that f(z) in S(c)-dom phii. So all we need is to compute the z, i.e., to find a function g such that z=g(i); then the productive function for S(c) is just f*g. We saw that
has the desired properties. Thus, if f is the (total) function phij, then we have z=c(i,j) (using the composition function c) and the productive function for S(c) is just
where j is a known constant.