sg5# 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%y))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]:}t|dk(st||z||D]}t|dk(s|t|<<td|D]d}t|dk(r|t|<|dz t|<%t|}||z}||zdk(rt||zt|<Ot||dz zt|<fy)N順rr?)_factor_totientintrange)maxnlimijxys L/var/www/html/venv/lib/python3.12/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-HQK"1+q1u-HQK .c |dk(rtS|}d}t|}||zdk(r||z}|dz }||zdk(r||z}dd|zz }t|}|dk(r|dk(rd|z} | || zz}|td|dz | z| z}t |d|dzd} t t td| z||t| |} t t td |zt| dz | ztt|||t| ||S|dk(rd|z} | || zz}|dkDr|td |dzdz | z| z}t |d|dzd} t t td| z||t| |} t t tdd |dzzzt| dztt|dz||t| ||S|||zz}||zdk(r>|dk(r3t ttd|tt|||StSt||stS||dz z|dz z} |td | dz |z}t |||d} t t td| z||t||} t t tdtd|ztt|||t| ||S|dk7s|dk\rt|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|dk(ryd|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 Av B A A q&A+ q Q q&A+ BB BqDA B Qw 6A#Ca#g A3q!a%%%,A%aAE215A'1 r4)*23-?C7"qAq!112!d+T3T"D* * 6A#Ca#g A1us2q!tax--4%aAE215A'(1Q3-T: t%C7B!a%=!Q78!Q$.6T"D* * AI q5A:AvVAq\*Xa[$/77Lq!$L1q5z1q5! s3q!$ $ !!Q *1 - HQqSM2t , QKw Qva|^ $ Xa[$ '/ C & &  AvaRc"bkB BK!tAvQR' "R HRL1,b 12578!tAvQR' "R HRL1,b 12578r"b$'BD)94@@AvsQwC"!1266" <AERUQYN*RU3q88A =w RL r2t d$ r2t   Q37CHRL1,b1 1R 7B qBEAI>!Q& &! +B 2b"d#RB%5t <ADA 71a=$ /FBDq!2& .A 2At 4d;A 1a=r.r;returncd}d}tD]A}|d|zdzz }||z }|dkr|S||}||dzz}d|kr|||z }||dzr| n|z }C|S)zt Calculate the partition function P(n) Parameters ========== n : int nonnegative integer rr3rr!r)r;prevrBpentar(npss r,_partition_recr]s A E W   1q Y 6  H H a!e  7 bMA 1q5aRa   Hr.c |dkry|dkr|tjz dkstjdk(r|dkr t|Sdtvr tt t j d|zdz dzzt jd |zz t jd z d zt jd z}t |d zd zx}}dt j dzzdt jdzz t j t jdzdz t j t jdzfd}tdt j|dz}|||dksJ|dkDr#|||dkr|dz}|dkDr |||dkr|}|dz}||z d kDr)||zdz}|||x}dkr|}n|dk\r|}||z d kDr)|} | dkDr tdt} ttt!dd||t#||} tt%d|} t'd | D]U} t)|| |}t+|| || | }t-| t|||} t/t1t3|dz}Wt 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.._MsJyy$q'zBtAq1uI.tyyDGA/FFFFr.r2(rz Input too bigr12)r] cache_lengthglobalsr-r$rfrCloglog2rgmaxceil ValueErrorrr rrrrr%r:rVr rabsr r)r;pbitsr=r@rmbigsmallrierMr\rMrNqrOdrjrkrls @@@r, _partitionrs 1u W ^88::R?  ' ' )Q .1v:a   !   1R#  1  xx|$&'( ) "  E59s?##D1 DGGQJDIIaL( )B 1 b B 3 BG a1c6" #C a:   (r!SzC'   (r!SzC' E 'C +/ 5[1 Q(NBc !C 3YE +/ A 5y)) A XmAq!4a8&)Q GF Xa[! $E 1a[* q!QK 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_targetcddlm}||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, npartitionsrsBR ! r.)F)& mpmath.libmprrrrrrr r r r r rrrrrsympy.external.gmpyrrrresidue_ntheoryrrsympy.utilities.decoratorrsympy.utilities.memoizationrrf itertoolsrr-r:rVr$r]rrr.r,rs@@@@@65C07 .0T=l$!Q c C  4J0#J0#J0Z  tBD Dr.