Partition function P: Difference between revisions

Content added Content deleted

Inline

@@ Line 188: / Line 188: @@
             if pd>nn then exit end if
             integer sgn = iff(remainder(i-1,4)<2 ? 1 : -1)
-            mpz pnmpd = partitionsP(nn-pd)
+            mpz pnmpd = pn[max(1,nn-pd)]
             if sgn=-1 then
                 mpz_sub(psum,psum,pnmpd)
@@ Line 205: / Line 205: @@
 {{out}}
 <pre>
-p(6666) = 193655306161707661080005073394486091998480950338405932486880600467114423441282418165863 (2.1s)
+p(6666) = 193655306161707661080005073394486091998480950338405932486880600467114423441282418165863 (0.8s)
 </pre>