g5#&ddlmZmZmZmZmZmZmZmZm Z m Z m Z m Z m Z mZmZmZddlmZmZmZddlmZmZddlmZddlmZddlZddlmZd Zd Z d Z!eddgd e"d e"fdZ#d e"d e"fdZ$edddddZ%dS))fzerofrom_int from_rationalfonefhalfbitcountto_intmpf_mulmpf_divmpf_submpf_addmpf_sqrtmpf_pi mpf_cosh_sinhmpf_cosmpf_sin)gcdlegendrejacobi)_sqrt_mod_prime_poweris_quad_residue) deprecated)recurrence_memoNcountcd}dg|zadg|zat|dzdz}td|D]E}t|dkr2t||z||D]}t|dkr |t|<Ftd|D]z}t|dkr|t|<|dz t|<+t|}||z}||zdkrt||zt|<_t||dz zt|<{dS)N順rr?)_factor_totientintrange)maxnlimijxys U/var/www/html/ai-engine/env/lib/python3.11/site-packages/sympy/ntheory/partitions_.py_prer- s Dc$hGs4xH dCi..1 C 1c]]## 1:??1Q3a(( # #1:??!"GAJ 1d^^ . . 1:??GAJA#HQK  AJ qD q5A::"1+a-HQKK"1+q1u-HQKK . .c 4 |dkrtS|}d}t|}||zdkr||z}|dz }||zdk||z}dd|zz }t|}|dkrM|dkrd|z} | || zz}|td|dz | z| z}t |d|dzd} t t td| z||t| |} t t td |zt| dz | ztt|||t| ||S|dkrd|z} | || zz}|dkr|td |dzdz | z| z}t |d|dzd} t t td| z||t| |} t t tdd |dzzzt| dztt|dz||t| ||S|||zz}||zdkrT|dkrGt ttd|tt|||StSt||stS||dz z|dz z} |td | dz |z}t |||d} t t td| z||t||} t t tdtd|ztt|||t| ||S|dks|dkrt|dt|d}} d| |zz}||z|z|dzdz | zzt||z|z|zt |dz |z|z}| |z|z|dzdz |zzt||z|z| zt |dz |z|z}t t#|||t#||||S|dkrd|zd ztd t |dz |z|z}d|dz |dzdz dzz |dzzdzzdz}t t td t#||||t#|||Sd|zdztdt |dz |z|z}d||dzdz dzz dzzdz}t t#|||t#||||S)z Compute the inner sum in HRR formula [1]_ References ========== .. [1] https://msp.org/pjm/1956/6-1/pjm-v6-n1-p18-p.pdf rrr! @i@r )rr"rpowrr r rrrrrrrrrr#_a)nkpreck1epk2vpimodmarg_phid1d2n1n2s r,r:r:$s Avv B A A q&A++ q Q q&A++ BB BqDA B Qww 66A#Ca#g A3q!a%%%%,A%aAE2215A'1 r4))*23--??C7"qAq!1!1122!d++T33T""D** * 66A#Ca#g A1uus2q!tax---4%aAE2215A'(1Q3--T:: t%%C7B!a%=!Q788!Q$..66T""D** * AI q5A::AvvVAq\\**Xa[[$//777Lq!$$ L1q5z1q5! s3q!$$ $ !!Q * *1 - HQqSMM2t , , QKKw Qva||^ $ $ Xa[[$ ' '// C  && &  AvvaRc"bkkB BK!tAvQR' "R HRL1,b 1 12578!tAvQR' "R HRL1,b 1 12578r"b$''BD)9)94@@@AvvsQwC"!12666" <AERUQYN*RU3q88A =w RLL r2t  d$$ r2t    Q37CHRL1,b11 1R 7B qBEAI>!Q& &! +B 2b"d##RB%5%5t < <>ADA 71a==$ / /FBDq!!2&& . .A 2At 4 4d;;A 1a==r.r;returncd}d}tD]H}|d|zdzz }||z }|dkrn/||}||dzz}d|kr |||z }||dzr| n|z }I|S)zt Calculate the partition function P(n) Parameters ========== n : int nonnegative integer rr3rr!r)r;prevrBpentar(npss r,_partition_recr]s A E WW    1q Y 66 E H a!e  77 bMA 1q5 aRRa Hr.c |dkrdS|dkr |tz dks#tdkr|dkrt|Sdtvrtt t jd|zdz dzzt jd |zz t jd z d zt jd z}t |d zd zx}}dt jdzzdt j dzz t jt j dzdz t jt j dzfd}tdt j |dz}|||dksJ|dkr+|||dkr|dz}|dkr|||dk|}|dz}||z d kr.||zdz}|||x}dkr|}n|dkr|}||z d k.|} | dkrtdt} ttt!dd||t#||} tt%d|} t'd | D]q} t)|| |}t+|| || | }t-| t|||} t/t1t3|dz}rt t3t-| t4|S)zX Calculate the partition function P(n) Parameters ========== n : int ri@ Fr!i@8r"g@rr4 rg?d,r3KgUUUUUU?ctj}||z |||dz z ztj||z|z zzS)Nr)mathsqrtsinh)r;Nrgc1c2c3s r,_Mz_partition.._MsUy$$q''zBttAq1uI.tyDDGGA/F/FFFFr.r2(rz Input too bigr12)r] cache_lengthglobalsr-r$rfrCloglog2rgmaxceil ValueErrorrr rrrrr%r:rVr rabsr r)r;pbitsr=r@rmbigsmallrierMr\rMrNqrOdrjrkrls @@@r, _partitionrs_ 1uuq W ^88:::R??  ' ' ) )Q . .1v::a    !!   1R#  1  x||$&'( ) ) "  E59s?###D1 DGQJDIaLL( )B 1 b B 3 BGGGGGGG a1c6"" # #C 2a::     ((rr!SzzC''   ((rr!SzzC'' E 'C +// 5[1 "Q((NBc ! !CC 3YYE +// A 5yy))) A XmAq!44a88&))Q G GF Xa[[! $ $E 1a[[** q!QKK q!Q & & Awq!}}d + + S^^ $ $r ) vga--.. / //r.zpThe `sympy.ntheory.partitions_.npartitions` has been moved to `sympy.functions.combinatorial.numbers.partition`.z1.13z%deprecated-ntheory-symbolic-functions)deprecated_since_versionactive_deprecations_targetFc$ddlm}||S)a  Calculate the partition function P(n), i.e. the number of ways that n can be written as a sum of positive integers. .. deprecated:: 1.13 The ``npartitions`` function is deprecated. Use :class:`sympy.functions.combinatorial.numbers.partition` instead. See its documentation for more information. See :ref:`deprecated-ntheory-symbolic-functions` for details. P(n) is computed using the Hardy-Ramanujan-Rademacher formula [1]_. The correctness of this implementation has been tested through $10^{10}$. Examples ======== >>> from sympy.functions.combinatorial.numbers import partition >>> partition(25) 1958 References ========== .. [1] https://mathworld.wolfram.com/PartitionFunctionP.html r) partition)%sympy.functions.combinatorial.numbersr)r;verbosefunc_partitions r, npartitionsrs(BRQQQQQ >!  r.)F)& mpmath.libmprrrrrrr r r r r rrrrrsympy.external.gmpyrrrresidue_ntheoryrrsympy.utilities.decoratorrsympy.utilities.memoizationrrf itertoolsrr-r:rVr$r]rrr.r,rs:@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@6555555555CCCCCCCC000000777777 ...0T=T=T=l$!Q c C    4J0#J0#J0J0J0J0Z  tBDDD DDr.