U L?h @s\ddlmZddlmZddlmZddlmZddlm Z ddl m Z ddl m Z dd lmZdd lmZmZmZmZmZdd lmZmZmZmZdd lmZdd lmZmZddl m!Z!m"Z"ddl#m$Z$ddl%m&Z&ddl'm(Z(ddl)m*Z*Gddde Z+e*dZ,e,-e.e.fe+ddl-m/Z/ddl0m1Z1m2Z2ddl3m4Z4m5Z5ddl6m7Z7m8Z8m9Z9dS)) annotations)Callable)product)_sympify)cacheit)S)Expr)PrecisionExhausted)expand_complexexpand_multinomial expand_mul_mexpand PoleError) fuzzy_bool fuzzy_not fuzzy_andfuzzy_or)global_parameters)is_gtis_lt) NumberKind UndefinedKind)sift)sympy_deprecation_warning)as_int) DispatchercseZdZUdZdZdZded<ded<edydd Zdzd d Z e d dddZ e d dddZ e ddZ eddZddZddZddZddZdd Zd!d"Zd#d$Zd%d&Zd'd(Zd)d*Zd+d,Zd-d.Zd/d0Zd1d2Zd3d4Zd5d6Zd7d8Z d9d:Z!d;d<Z"d=d>Z#d?d@Z$dAdBZ%dCdDZ&dEdFZ'dGdHZ(d{dIdJZ)dKdLZ*dMdNZ+dOdPZ,dQdRZ-dSdTZ.dUdVZ/dWdXZ0dYdZZ1d[d\Z2d]d^Z3d|d`daZ4d}dcddZ5d~dedfZ6edgdhZ7fdidjZ8dkdlZ9dmdnZ:dodpZ;dqdrZdwdxZ?Z@S)Powa% Defines the expression x**y as "x raised to a power y" .. deprecated:: 1.7 Using arguments that aren't subclasses of :class:`~.Expr` in core operators (:class:`~.Mul`, :class:`~.Add`, and :class:`~.Pow`) is deprecated. See :ref:`non-expr-args-deprecated` for details. Singleton definitions involving (0, 1, -1, oo, -oo, I, -I): +--------------+---------+-----------------------------------------------+ | expr | value | reason | +==============+=========+===============================================+ | z**0 | 1 | Although arguments over 0**0 exist, see [2]. | +--------------+---------+-----------------------------------------------+ | z**1 | z | | +--------------+---------+-----------------------------------------------+ | (-oo)**(-1) | 0 | | +--------------+---------+-----------------------------------------------+ | (-1)**-1 | -1 | | +--------------+---------+-----------------------------------------------+ | S.Zero**-1 | zoo | This is not strictly true, as 0**-1 may be | | | | undefined, but is convenient in some contexts | | | | where the base is assumed to be positive. | +--------------+---------+-----------------------------------------------+ | 1**-1 | 1 | | +--------------+---------+-----------------------------------------------+ | oo**-1 | 0 | | +--------------+---------+-----------------------------------------------+ | 0**oo | 0 | Because for all complex numbers z near | | | | 0, z**oo -> 0. | +--------------+---------+-----------------------------------------------+ | 0**-oo | zoo | This is not strictly true, as 0**oo may be | | | | oscillating between positive and negative | | | | values or rotating in the complex plane. | | | | It is convenient, however, when the base | | | | is positive. | +--------------+---------+-----------------------------------------------+ | 1**oo | nan | Because there are various cases where | | 1**-oo | | lim(x(t),t)=1, lim(y(t),t)=oo (or -oo), | | | | but lim( x(t)**y(t), t) != 1. See [3]. | +--------------+---------+-----------------------------------------------+ | b**zoo | nan | Because b**z has no limit as z -> zoo | +--------------+---------+-----------------------------------------------+ | (-1)**oo | nan | Because of oscillations in the limit. | | (-1)**(-oo) | | | +--------------+---------+-----------------------------------------------+ | oo**oo | oo | | +--------------+---------+-----------------------------------------------+ | oo**-oo | 0 | | +--------------+---------+-----------------------------------------------+ | (-oo)**oo | nan | | | (-oo)**-oo | | | +--------------+---------+-----------------------------------------------+ | oo**I | nan | oo**e could probably be best thought of as | | (-oo)**I | | the limit of x**e for real x as x tends to | | | | oo. If e is I, then the limit does not exist | | | | and nan is used to indicate that. | +--------------+---------+-----------------------------------------------+ | oo**(1+I) | zoo | If the real part of e is positive, then the | | (-oo)**(1+I) | | limit of abs(x**e) is oo. So the limit value | | | | is zoo. | +--------------+---------+-----------------------------------------------+ | oo**(-1+I) | 0 | If the real part of e is negative, then the | | -oo**(-1+I) | | limit is 0. | +--------------+---------+-----------------------------------------------+ Because symbolic computations are more flexible than floating point calculations and we prefer to never return an incorrect answer, we choose not to conform to all IEEE 754 conventions. This helps us avoid extra test-case code in the calculation of limits. See Also ======== sympy.core.numbers.Infinity sympy.core.numbers.NegativeInfinity sympy.core.numbers.NaN References ========== .. [1] https://en.wikipedia.org/wiki/Exponentiation .. [2] https://en.wikipedia.org/wiki/Zero_to_the_power_of_zero .. [3] https://en.wikipedia.org/wiki/Indeterminate_forms Tis_commutativeztuple[Expr, Expr]args_argsNcCs^|dkrtj}t|}t|}ddlm}t||s>t||rFtd||fD],}t|tsNtdt |j ddddd qN|r"|t j krt j S|t jkrt|t jrt jSt|t jrt|t jrt jSt|t jr|jrt j S|jd krt j S|t jkrt jS|t jkr|S|d kr,|s,t j S|jj d krn|t jkrd dlm}|t||jt||jSn`|jr~|js|jr|jr|j s|j!r|"r|j#r| }n|j$rt| | St j ||fkrt j S|t jkrt%|j&rt j St jSd dl'm(}|j)s |t jk r t||s ddl*m+}d dl'm,} d dl-m.} ||d d/\} } | | \} }t|| r|j0d |krt j| | S|j1r d dl2m3}m4}|||}|j!r |r || ||d d |t j5t j6kr t j| | S|7|}|dk r"|St8|||}|9|}t|tsJ|S|j:oV|j:|_:|S)Nr) Relationalz Relational cannot be used in Powzf Using non-Expr arguments in Pow is deprecated (in this case, one of the arguments is of type zf). If you really did intend to construct a power with this base, use the ** operator instead.z1.7znon-expr-args-deprecated)Zdeprecated_since_versionZactive_deprecations_target stacklevelFZAccumulationBoundsr AccumBounds) exp_polar) factor_termslog)fraction)sign)r-im);revaluater relationalr" isinstance TypeErrorr rtype__name__rComplexInfinityNaNInfinityrOne NegativeOnerZero is_finite __class__Exp1!sympy.calculus.accumulationboundsr'rminmax is_Symbol is_integer is_Integer is_numberis_Mul is_Numbercould_extract_minus_signis_evenis_oddabs is_infinite&sympy.functions.elementary.exponentialr(Zis_AtomZ exprtoolsr)r+Zsympy.simplify.radsimpr, as_coeff_Mulr is_Add$sympy.functions.elementary.complexesr-r. ImaginaryUnitPi _eval_power__new__Z _exec_constructor_postprocessorsr)clsber/r"argr'r(r)r+r,cexnumZdenr-r.sobjr]B/var/www/html/venv/lib/python3.8/site-packages/sympy/core/power.pyrSvs                           z Pow.__new__rcCs |jtjkrddlm}|SdSNrr*)baserr=rLr+)selfZargindexr+r]r]r^inverses  z Pow.inverser )returncCs |jdS)Nrr!rar]r]r^r`szPow.basecCs |jdSNrrdrer]r]r^expszPow.expcCs|jjtkr|jjStSdSN)rgkindrr`rrer]r]r^ris zPow.kindcCs dd|jfS)N)r4rTr]r]r^ class_keysz Pow.class_keycCsrddlm}m}|\}}||||rn|rn||||rPt| |S||||rnt| | SdS)Nr)askQ) Zsympy.assumptions.askrnro as_base_expintegerrGZevenrodd)raZ assumptionsrnrorUrVr]r]r^ _eval_refines  zPow._eval_refinecCs|\}}|tjkr"|||Sd}|jr4d}n.|jrBd}n |jdk rbddlm}m}m }m }ddl m } m } ddlm} dd} dd } |jr|d kr| |r|jd krtj|t| ||S|jd krt|| Sn0|jr|jrt|}|jrt||tj}t|dkd ks2|dkr8d}n|jrFd}n||jrjt|d kd krjd}nn| |rb| d tjtj|| tj|||d tj}|jr| |||dkr||}nd}nzl| d tjtj|| tj||| |d tj}|jr@| |||dkr@||}nd}Wntk r`d}YnX|dk r~|t|||SdS)Nrr)rWr.rer-rgr+)floorcSs6t|dddkrdS|\}}|jr2|dkr2dSdS)zZReturn True if the exponent has a literal 2 as the denominator, else None.qNrkT)getattras_numer_denomrB)rVndr]r]r^_half s  zPow._eval_power.._halfcSs8z|jddd}|jr|WSWntk r2YnXdS)zXReturn ``e`` evaluated to a Number with 2 significant digits, else None.rkTstrictN)evalfrFr )rVrvr]r]r^_n2s  zPow._eval_power.._n2r%TFrk)rprr6rBis_polaris_extended_realrOrWr.rtr-rLrgr+#sympy.functions.elementary.integersrv is_negativer9rrHrJ is_imaginaryrPis_extended_nonnegativerQHalfr )raotherrUrVr[rWr.rtr-rgr+rvr|rr]r]r^rRsd          "   zPow._eval_powerc Csr|j|j}}|jrn|jrn|jr6||dkr6tjSddlm}|jr|jr|jrt |t |t |}}}| }|dkr||kr| d|krt ||} t t || || |St t |||Sddl m} t|tr|jr|jr| ||}| t||dd|St|trn|jrn|jrnt | } | dkrn||} | | || }| t||dd|Sd S) aOA dispatched function to compute `b^e \bmod q`, dispatched by ``Mod``. Notes ===== Algorithms: 1. For unevaluated integer power, use built-in ``pow`` function with 3 arguments, if powers are not too large wrt base. 2. For very large powers, use totient reduction if $e \ge \log(m)$. Bound on m, is for safe factorization memory wise i.e. $m^{1/4}$. For pollard-rho to be faster than built-in pow $\log(e) > m^{1/4}$ check is added. 3. For any unevaluated power found in `b` or `e`, the step 2 will be recursed down to the base and the exponent such that the $b \bmod q$ becomes the new base and $\phi(q) + e \bmod \phi(q)$ becomes the new exponent, and then the computation for the reduced expression can be done. r)totientPr#r)ModFr/N)r`rgrB is_positiverr:Z%sympy.functions.combinatorial.numbersrrCint bit_lengthIntegerpowmodrr1rrD) rarwr`rgrrUrVmmbphirrr]r]r^ _eval_ModOs,       z Pow._eval_ModcCs|jjr|jjr|jjSdSrh)rgrBrr`rHrer]r]r^ _eval_is_evenszPow._eval_is_evencCst|}|dkr|jS|SNT)r_eval_is_extended_negativer;)raZext_negr]r]r^_eval_is_negatives zPow._eval_is_negativecCs|j|jkr|jjrdSn|jjr0|jjrdSn|jjrR|jjrDdS|jjrdSn|jjrl|jj r|jjSnn|jj r|jjrdSnX|jj r|jj r|jd}|jrdS|j r|jdkrdS|jj rddl m}||jj SdS)NTFr#rr*)r`rgrris_realis_extended_negativerHrIis_zeroris_extended_nonpositiverrBrLr+)rarr+r]r]r^_eval_is_extended_positives6    zPow._eval_is_extended_positivecCs|jtjkr |jjs|jjr dS|jjrJ|jjr<|jjrdS|jsJ|jrNdS|jr|jr|jsf|jrt |dj rt |dj rdS|j r|j r|j |j}|j S|jr|jr|djrdS|jr|jr|djrdSdS)NFTr)r is_rationalrBrrr9rrr;rrrFfuncrC)rarUrVcheckr]r]r^_eval_is_integers&      zPow._eval_is_integerc Cs|jtjkr8|jjrdS|jjr8dtj|jtjjSddl m }m}|jj}|dkr|jj |krv|jjjrv|jjS|jj t kr|jjtjkr|jjjr|jjSdS|jj}|dkrdS|r.|r.|jj rdS|jjr|jjrdS|jjr|jjrdS|jjr|jjrdS|jjr.|jjr.dS|r^|jjr^|jjdkr^t |j|j jS|jj}|jj}|r$|jjr|jjrdS|jjr$dSn|r||jjrdS|jjr|j\}}|r$|jr$t|j||j|ddjSn,|jtj tjfkr$|jdjdkr$dS|r|r|jtjkrBdS|jtj}|r|jjr|jr|jjr|jdjr|jrdS|||jtjj} | dk r| S|dkr|rt|jtr|jj dkrdSddl!m"} | |j|jtj} | j#r| jSdS) NTrkr)r+rgFrrrW)$r`rr=rgrrrPrQrHrLr+rrrrrBZis_extended_nonzerorr is_RationalrrIrN as_coeff_AddrCMulr9coeffr is_nonzeror1RationalprOrWr) rar+rgZreal_bZreal_eZim_bim_erXaokrWir]r]r^_eval_is_extended_reals $            zPow._eval_is_extended_realcCsD|jtjkr t|jj|jjgStdd|jDr@| r@dSdS)Ncss|] }|jVqdSrh)r).0rr]r]r^ ?sz'Pow._eval_is_complex..T) r`rr=rrgrrallr _eval_is_finiterer]r]r^_eval_is_complex:s zPow._eval_is_complexc Cs>|jjdkrdS|jjr8|jjr8|jj}|dk r4|SdS|jtjkrrd|jtjtj }|j rddS|jrndSdS|jjrddl m }||jj}|dk rdS|jj r|jj r|jjrdS|jj}|s|S|jjrdSd|jj}|r|jjS|S|jj dkr:ddlm}||j|jtj}d|j} | dk r:| SdS)NFrkTrr*r)r`rrrgrBrIrr=rQrPrHrLr+rrrrrOrW) rarrfr+ZimlogratZhalfrWrZisoddr]r]r^_eval_is_imaginaryBsL        zPow._eval_is_imaginarycCs@|jjr<|jjr|jjS|jjr,|jjr,dS|jtjkr=2)-th power cannot be a prime. rFN)r`rBrgrrer]r]r^_eval_is_primeszPow._eval_is_primecCsL|jjrH|jjrH|jdjr(|jdjsD|jdjrH|jjrH|jjrHdSdS)zS A power is composite if both base and exponent are greater than 1 rTN)r`rBrgrrrHrer]r]r^_eval_is_composites   zPow._eval_is_compositecCs|jjSrh)r`rrer]r]r^_eval_is_polarszPow._eval_is_polarcCsddlm}t|j|rT|j||}|j||}t||rH||S|||Sddlm}m }dd}||jks||kr|jt j kr|j rt|t r||j||S||j||St||jr|j|jkr||j|j} | jrt|| St||jrZ|j|jkrZ|jjdkr|jjtdd} |jjtdd} || | |\} } }| rZ||| }|dk r~t|t|j|}|Sn|j}g}g}|} |jjD]z}|||}|} || | |\} } }| r||| |dk r||qn|js|jsdS||q|rZt|}||dkrJt|j|dd n|jt|St||s||jr|jt j kr|jjr|jjr|jjtdd} |j||jjtdd} || | |\} } }| r||| }|dk rt|t|j|}|SdS) Nrr&ruc Ss$|\}}|\}}||kr |jr||}zt|ddd}Wn8tk rt|\} } | jrd| jpn| jon| j}YnX||dfSt|ts|f}t dd|DsdSzbt t|t|\}} |dkr| dkr|d 7}| t|8} | dkrd} nt | f|} d|| fWStk rYnXdS) a*Return (bool, pow, remainder_pow) where, if bool is True, then the exponent of Pow `old` will combine with `pow` so the substitution is valid, otherwise bool will be False. For noncommutative objects, `pow` will be an integer, and a factor `Pow(old.base, remainder_pow)` needs to be included. If there is no such factor, None is returned. For commutative objects, remainder_pow is always None. cti are the coefficient and terms of an exponent of self or old In this _eval_subs routine a change like (b**(2*x)).subs(b**x, y) will give y**2 since (b**x)**2 == b**(2*x); if that equality does not hold then the substitution should not occur so `bool` will be False. Fr}TNcss|] }|jVqdSrh)rB)rtermr]r]r^rsz1Pow._eval_subs.._check..)FNNrr) rr ValueErrorrprrrr1tuplerdivmodr) ct1ct2oldZcoeff1Zterms1Zcoeff2Zterms2rZcombinesrUrV remainder remainder_powr]r]r^_checks8       zPow._eval_subs.._checkF)Zas_Addrr)r>r'r1rgr`subs__rpow__rrLr+rr=Z is_Functionr_subsrFrrNZas_independentSymbolrZ as_coeff_mulr appendrrBAddis_Powrr)rarnewr'rUrVrgr+rlrrrrrresultZoargZnew_lZo_alrZnewaexpor]r]r^ _eval_subssv     :       &6  zPow._eval_subscCs<|j\}}|jr4|j|jkr4|jdkr4d|| fS||fS)aReturn base and exp of self. Explanation =========== If base a Rational less than 1, then return 1/Rational, -exp. If this extra processing is not needed, the base and exp properties will give the raw arguments. Examples ======== >>> from sympy import Pow, S >>> p = Pow(S.Half, 2, evaluate=False) >>> p.as_base_exp() (2, -2) >>> p.args (1/2, 2) >>> p.base, p.exp (1/2, 2) rr)r rrrw)rarUrVr]r]r^rp s zPow.as_base_expcCsrddlm}|jj|jj}}|r2||j|jS|rF|j||jS|dkrn|dkrnt|}||krn||SdS)Nr)adjointF)rOrrgrBr`rr )rarrrexpandedr]r]r^ _eval_adjoint=s zPow._eval_adjointcCs|ddlm}|jj|jj}}|r2||j|jS|rF|j||jS|dkrn|dkrnt|}||krn||S|jrx|SdS)Nr) conjugateF)rOrrgrBr`rr r)rarXrrrr]r]r^_eval_conjugateIs zPow._eval_conjugatecCsddlm}|jtjkr,|tj|jS|jj|jjp@|jj }}|rV|j|jS|rj||j|jS|dkr|dkrt |}||kr||SdS)Nr) transposeF) rOrr`rr=rrgrBrrKr )rarrrrr]r]r^_eval_transposeWs   zPow._eval_transposec sjj}tjkrXddlm}t||rX|jrXddlm }| |j f|j S|j r|dds|jdks||r|jrtfdd|jDSjrt|jdd d d \}}|rtfd d|Dt|SS) za**(n + m) -> a**n*a**mr)Sum)ProductforceFcsg|]}|qSr]rrxrUrar]r^ qsz.Pow._eval_expand_power_exp..cSs|jSrhrrr]r]r^sz,Pow._eval_expand_power_exp..Tbinarycsg|]}|qSr]rrrr]r^rus)r`rgrr=Zsympy.concrete.summationsrr1rZsympy.concrete.productsrrfunctionZlimitsrNgetr_all_nonneg_or_nonpposrr rr _from_args)rahintsrVrrrXncr]rr^_eval_expand_power_expes*    zPow._eval_expand_power_expc sTdd}j}j|js"S|jdd\}}|rfdd|D}jrjrbt|}n"tdd|dddD }|r|t|9}|S|sjt|dd St|g}t |d d d d \}}dd} t || } | d } || d7}| d} | t j } | rt j }t | d}|dkr0n|dkrF| |n~|dkr| rx|  }|t jk r| |n | t jn>| r|  }|t jk r| |n | t j| |~ |sԈjr| | |}|}njrtt | dkrjt j}|s&| djr&|| d9}t | dr:| }| D]}| | q>|t jk r| |n\| r|r| djr| dt jk r| t j| | d n || n || ~ | }||7}t j}|r2jrt |dd d d \}}tfdd|D}|tfdd|D9}|rP|jt|dd 9}|S)z(a*b)**n -> a**n * b**nrF)Zsplit_1cs&g|]}t|dr|jfn|qS)_eval_expand_power_base)hasattrrrr)rr]r^rsz/Pow._eval_expand_power_base..cSsg|] }|dqS)r%r]rr]r]r^rsNr%rcSs |jdkSNF)rrr]r]r^rrz-Pow._eval_expand_power_base..TrcSs4|tjkrtjS|j}|rdS|dkr0t|jSdSr)rrPrrr)rZpolarr]r]r^preds z)Pow._eval_expand_power_base..predr#rrrkcSs|jo|jjo|jjSrh)rrgrr`rDrr]r]r^rscs g|]}|j|jqSr])rr rrUrVrar]r^rscsg|]}j|ddqS)Frrrrr]r^rs)rr`rgrEZargs_cncrCrrrrrrPlenrpopr8r9rBAssertionErrorrFextendr)rarrrUZcargsrrrZ maybe_realrZsiftedZnonnegnegimagIrZnonnorzZnpowr])rVrrar^rys  "                             zPow._eval_expand_power_basec s2|j\}|}|jrX|jdkrXjrX|jst|j|j}|sH|S|||g}}||}|jrx| }t |D]}| ||qt |St |}jrgg}} jD] } | jr| | q| | q|rLt | } t |} |dkrt| |dd|| | St| |ddd} t| | dd|| | Sjr\}} |jr| jr|js| js||j| j|}|j| j|j| j}} n ||j|}|j|j| }} n.| js|| j|}|| j| j}} nd}t |t | ddf\}} }}|r|d@rV||| || |||}}|d8}||| | d|| }} |d}qtj}|dkr|||St|||||S| }ddlm}ddlm}|t||}||f|S|dkrt fdd jDS|d jr>t fd d jDSt fd d jDSn|jr|jdkrjrt|j|jkrd||  S|jr*jr*|d dsЈjdks|r*gg}}|jD],}|jr| ||n | |qt ||t !|gS|Sd S)zA(a + b + ..)**n -> a**n + n*a**(n-1)*b + .., n is nonzero integerrrkFdeepr)multinomial_coefficients)basic_from_dictcs g|]}jD] }||qqSr]r rrg)r`r]r^rSsz0Pow._eval_expand_multinomial..cs g|]}jD] }||qqSr]r r multir]r^rWscsg|] }|qSr]r])rrr r]r^r[srN)"r rrrNrCrrwrr_eval_expand_multinomialr make_argsrrris_Orderr r rD as_real_imagrrPZsympy.ntheory.multinomialrZsympy.polys.polyutilsr rrJrFrrrrr)rarrgrrzradicalZexpanded_base_nrZ order_termsZ other_termsrUrrr rkrXr{rrrr Zexpansion_dictrtailr])r`rr^rs         "        zPow._eval_expand_multinomialc sb|jjrddlm}|j}|jj|d\}}|s<|tjfStdt d\|dkr|j r~|j r~t |j|}||kr~|S||}nj|d|d}||| |}}|j r|j rt ||tj | }||kr|S|| }dd| D} tfd d| D} d d| D} tfd d| D} d d| D} tfd d| D} | |tj |i| ||i| || ifSddlm} m}m}|jjr|jj|d\}}|jr*|jtjkr*|jr|tjfS|jr*tj|j |jfS|||d||dtj} | ||}|| |j||j}}||||||fS|jtjkrddlm}|j\}}|r|j|f|}|j|f|}||||}}||||||fSddlm}m}|rNd|d<|j|f|}| d|kr.cs(g|] \\}}}|||qSr]r]rZaaZbbccrrUr]r^rs cSs$g|]}|ddddkr|qS)rrr#r]rr]r]r^rscs(g|] \\}}}|||qSr]r]rrr]r^rs cSs$g|]}|ddddkr|qS)rrr#rjr]rr]r]r^rscs(g|] \\}}}|||qSr]r]rrr]r^rs )atan2cossinrg)r.rtFcomplexignore)!rgrCZsympy.polys.polytoolsrr`rrr:symbolsDummyrFr rPtermsrr(sympy.functions.elementary.trigonometricrrrrrrrrrr=rLexpandrOr.rtr)rarrrrgZre_erexprmagrZre_partZim_part1Zim_part3rrrtrptprXr[r.rtrr]rr^rnsv      $  "  zPow.as_real_imagcCsFddlm}|j|}|j|}||||j||j|jSr_)rLr+r`diffrg)rar[r+ZdbaseZdexpr]r]r^_eval_derivatives   zPow._eval_derivativecCs|\}}|tjkr6ddlm}||jdd|S||}|jsP||}|jr|j r|j dkr| || |}| }| || S| ||S)NrrFr)rprr=rLrg _eval_evalfZ_evalfrCrrDrrrr$)raprecr`rgZ exp_functionr]r]r^r-s     zPow._eval_evalfcCsF|jj|rdS|jj|r>t|j|o:|jjo:|jdkSdSdS)NFrT)rghasr`bool_eval_is_polynomialrCraZsymsr]r]r^r1s  zPow._eval_is_polynomialcCs|jjr,|jjr,tt|jj|jjgr,dS|j| }|j sF|jS| \}}|j rb|j rbdS|jr|jrt|js~|j rdS||krdSn |j r|jS|tjkr|jr|jrdSdS)NTF)rgrBr`rrrrrrrprrr is_irrationalrr=r)rarrUrVr]r]r^_eval_is_rationals*    zPow._eval_is_rationalcCs6dd}|jjs||jrdS|jtjkr|j|j}|j|jkr|jjr|jjrVdS|jtj j rhdS|jtj tj j rdSn|jSn|jj r|jjdkr|jjS|jjdkr|jjr|jjS|jjrdS|jj r|jjSnR|jjr2|jjr2t |jjrt ||js*|jjdks*|jjr2|jj SdS)NcSs(z |djWStk r"YdSXdS)NrF)rr)r%r]r]r^_is_ones z'Pow._eval_is_algebraic.._is_oneTF)r`rrr=rr rgrZ is_algebraicrQrrPrrrBr3)rar5r[r]r]r^_eval_is_algebraicsB        zPow._eval_is_algebraiccCs8|jj|rdS|jj|r0|j|o.|jjSdSdSr)rgr/r`_eval_is_rational_functionrCr2r]r]r^r7"s   zPow._eval_is_rational_functionc Cs|j||}|jj}|r|S|j||}|dkr@|rs.    z Pow.matchesrc3 s$ddlm}m}ddlm}ddlm}ddlm} |j t j kr |jj |||d} | j rbd| S|| |d} | t jkr||||S| t jkr|S| | } || } }td|D](}|| |9}|j |||d}| |7} q| || ||7} ddlm}|| d d d Sdd lm}dd lm}||d d}|\}}|j|rRt||r|||||j||||dS|dk r||rddlm}td||gd\}}|||||||||}||}|}z>ddl m!}||t j"r|dk rt#|$|\}}Wn^t#t%tfk r|j|t&d|||d}|t j't j(rpt%|$|\}}YnX||rddl)m*}||+}|j,s|j-r|j.s||j/|||dkr||||j||||d}|||||kr|S|S|j0||d}||t j1j+dd}|j-sDt%|||t2rd| |j3r~|||||S|j,r||}||kr|||||7}|Sdd} fdd}!z|j$||d\}}"WnVt#t%fk r.||||ddkr$||||||YSt%YnX|j4rf|"t j5krfddl)m6}#|#|j$||d\}}"|"j7s|8}|j,r||S|j$||d\}}"|"j7s|||9}|j$||d\}}"|"j7st%ddl:m;}$|j||$||d}%i}&t<=|%D]*}| ||\}'}(|&>|(t j5|'|&|(<qt j1})t j5t j1i}*|&}+dd l?m@},mA}-|)|"j3r|-||)|,|)}.|+D]$}|*>|t j5|.|+||*|<q||!|+|&}+|)t j17})qVdd!lBmC}/|jDsb|j,rb|j3rb||E||}0|/|0j3r| ||d"d#||\}1}2nF|/|0j,rN| ||||j0|||d|\}1}2n| |||\}1}2n| |||\}1}2t j5}|*D]&}(|(|2}||*|(|1||7}q~|jDr|j7r||"|jFr|tG|ks z|||||7}Wn4t%k r||||j||||dYSX|S)$Nrru)limit)Ordersympify)rzlogxr)powsimpTrg)rcombine) powdenest)_illegal)r)rzrFcdir)Wildzc, ex)rTexclude) polygammark) logcombinerFrKrFF)r$c Sstjtj}}t|D]`}||rp|\}}||krxz||WStk rl|tjfYSXq||9}q||fSrh) rr8r:rrr/rpleadtermr)rrrrgfactorr`r]r]r^ coeff_exps   z$Pow._eval_nseries..coeff_expcsNi}t||D]:\}}||}|kr||tj||||||<q|Srh)rrrr:)Zd1Zd2rese1e2rYZmaxpowr]r^muls $zPow._eval_nseries..mul) nsimplify)ceiling) factorialffr.r%)HrLrgr+Zsympy.series.limitsrBZsympy.series.orderrCZsympy.core.sympifyrEr`rr=ZnseriesrZremoveOrr7rangeZsympy.simplify.powsimprGrInumbersrJZtrigsimprpr/r _eval_nseriessymbolrLr replaceZ'sympy.functions.special.gamma_functionsrNZ EulerGammarrRNotImplementedErrorr@r6r5Zsympy.simplify.simplifyrOcancelrrDr_eval_as_leading_termas_leading_termr8rrZis_Floatr:rZrsimplifyr$rr[rrr(sympy.functions.combinatorial.factorialsr\r]rOr.rBdirr=r)3rarrzrFrKrgr+rBrCrEZe_seriesZe0r(Z exp_seriesrrrGrIrJrUrVrLrXrYrN_rrOrUrr r'rTrYr{rZr[ZgpolyZgtermsZco1rVrr"Ztkr\r]rr.ndirZincoZinexr]rXr^rbs                $            "     (  $zPow._eval_nseriescCsxddlm}m}|j}|j}|jtjkr||j||d}||d} | tjkrX| |d} | j dkrltj| St d|n| |r||||} | j|||dSddl m} z|j|||d} Wnt k r|YSX|jsh| jrh| |sh|| ||} | | jr0|| |dd |S| | jrh||j|||d}|j dkrh|||S|| |SdS) NrrurQFzCannot expand %s around 0rPr^r%r_)rLrgr+r`rr=rhrr6rBrKrr/rOr.rBrrkrrrg)rarrFrKrgr+rVrUrWZarg0ltr.rrmZ log_leadtermr]r]r^rgos8             zPow._eval_as_leading_termcGs$ddlm}||j||||S)Nr)binomial)rjrorgr)rarzrprevious_termsror]r]r^ _taylor_terms zPow._taylor_termcs|jtjk r tj||f|S|dkr.tjS|dkr|\}}t|j||d}|j||d\}}|jr|\}}|jr|tjkr||} ||| } tj} | jst| j | j \} } ||| } | ||t||| || j fSt||}|jr,|j r,|j||d\}}||| \} } | \} }| tj ks||kr,| |t| ||fStj |||fS)aReturn the tuple (R, self/R) where R is the positive Rational extracted from self. Examples ======== >>> from sympy import sqrt >>> sqrt(4 + 4*sqrt(2)).as_content_primitive() (2, sqrt(1 + sqrt(2))) >>> sqrt(3 + 3*sqrt(2)).as_content_primitive() (1, sqrt(3)*sqrt(1 + sqrt(2))) >>> from sympy import expand_power_base, powsimp, Mul >>> from sympy.abc import x, y >>> ((2*x + 2)**2).as_content_primitive() (4, (x + 1)**2) >>> (4**((1 + y)/2)).as_content_primitive() (2, 4**(y/2)) >>> (3**((1 + y)/2)).as_content_primitive() (1, 3**((y + 1)/2)) >>> (3**((5 + y)/2)).as_content_primitive() (9, 3**((y + 1)/2)) >>> eq = 3**(2 + 2*x) >>> powsimp(eq) == eq True >>> eq.as_content_primitive() (9, 3**(2*x)) >>> powsimp(Mul(*_)) 3**(2*x + 2) >>> eq = (2 + 2*x)**y >>> s = expand_power_base(eq); s.is_Mul, s (False, (2*x + 2)**y) >>> eq.as_content_primitive() (1, (2*(x + 1))**y) >>> s = expand_power_base(_[1]); s.is_Mul, s (True, 2**y*(x + 1)**y) See docstring of Expr.as_content_primitive for more examples. )rclear)rp _keep_coeffas_content_primitiverrrr:rrrrwrErMr8)rarrzrUrVZceZpehr(ZcehrXr'Zicehrmer]r]r^r|s*+    $  zPow.as_content_primitivec Os|}|ddr|}|\}}|d}|rJ||}||krJ|S|j|}|j|} | r|rjdS|d}|dkrdSn | dkrdS|dS)NriTrF)rrirpequals is_constant) raZwrtflagsr%rUrVZbzrZeconZbconr]r]r^rs(      zPow.is_constantcCsF|j\}}||rB||sB||||}|||d|SdSrf)r r/r)rarzsteprUrVZnew_er]r]r^_eval_difference_delta(s zPow._eval_difference_delta)N)r)T)NF)r)Nr)FT)Ar4 __module__ __qualname____doc__r __slots____annotations__rrSrbpropertyr`rgri classmethodrmrsrRrrrrrrrrrrrrrrrrrprrrrrrrr,r-r1r4r6r7r8r9r<ryr>rbrgrqrsrurvrxryr|rr __classcell__r]r]rtr^rs X b    T8F1   {z S ' % # " < "   Qrpower)r)rr)rr{)rr!r N): __future__rtypingr itertoolsrrErcacherZ singletonrr%r rr rr r r rrZlogicrrrr parametersrr0rrrirrZsympy.utilities.iterablesrZsympy.utilities.exceptionsrZsympy.utilities.miscrZsympy.multipledispatchrrraddobjectrrarrrYrr{rcrr!r r]r]r]r^sJ             &