sgæpddlmZddlmZddlmZmZmZddlm Z ddl m Z ddl m Z ddlmZddlmZdd lmZdd lmZdd lmZmZdd lmZdd lmZddlmZmZm Z m!Z!m"Z"m#Z#m$Z$m%Z%m&Z&m'Z'm(Z(m)Z)m*Z*m+Z+m,Z,m-Z-m.Z.m/Z/m0Z0m1Z1m2Z2m3Z3m4Z4m5Z5m6Z6m7Z7m8Z8m9Z9m:Z:m;Z;mZ>m?Z?m@Z@mAZAmBZBmCZCmDZDmEZEmFZFmGZGmHZHmIZIddlJmKZKmLZLmMZMddlNmOZOe dd\ZPZQZRZSZTZUZVZWZXZYZZZ[Z\Z]Z^Z_e dd\Z`ZaZbZcZdZeZfe dd\ZgZhZie ddZjdZkdZldZmdZndZodZpdZqdZrdZsd Ztd!Zud"Zvd#ZweKd$Zxd%Zyd&Zzd'Z{d(Z|d)Z}d*Z~d+Zd,Zd-Zd.Zd/Zd0Zd1Zd2Zd3Zd4Zd5Zd6Zd7Zd8Zd9Zd:Zd;Zd<Zd=Zd>Zd?Zd@ZeMdAZdBZdCZdDZdEZdFZdGZeMdHZdIZdJZdKZdLZyM)N)Add)Mul)Rationaloopi)Eq)S)symbols)Matrix) factorint)powsimp)_mexpand)default_sort_keyordered)sin) diophantine),diop_DN diop_solvediop_ternary_quadratic_normaldiop_general_pythagoreandiop_ternary_quadratic diop_lineardiop_quadraticdiop_general_sum_of_squaresdiop_general_sum_of_even_powersdescent diop_bf_DN divisible equivalentfind_DNldescentlength reconstruct partitionpower_representationprime_as_sum_of_two_squares square_factorsum_of_four_squaressum_of_three_squarestransformation_to_DNtransformation_to_normal classify_diopbase_solution_linear cornacchia sqf_normalgaussian_reduceholzer check_paramparametrize_ternary_quadratic sum_of_powerssum_of_squares_diop_ternary_quadratic_normal_nint_or_floor_odd_even _remove_gcd_can_do_sum_of_squaresDiophantineSolutionSetGeneralPythagoreanBinaryQuadratic)slowraisesXFAIL)signed_permutationsz.a, b, c, d, p, q, x, y, z, w, t, u, v, X, Y, ZT)integerzt_:7zm1:4n1c<ttt|SN)rr eqs c/var/www/html/venv/lib/python3.12/site-packages/sympy/solvers/diophantine/tests/test_diophantine.py diop_simplifyrJ's GHRL) **cHttdttdy)Nc2tttSrF)rrxrKrIz#test_input_format..,sk#a&1rKc4tttz dz SN)rrNrrOrKrIrPz#test_input_format..-sk!B$(3rK)r@ TypeErrorrOrKrItest_input_formatrU+s 912 934rKcvtdtk(sJttdtk(sJyrR)rsetr rOrKrI test_nosolsrX0s/ q>SU "" " qt  %% %rKcttdz tdz dzzddhk(sJttdz tdz zddhk(sJy)NrZr[rrNrOrKrItest_univariater_6sM q1uq1uqj( )dD\ 99 9 q1uq1uo &4, 66 6rKc ttdttdttdttdt dt dzzdt zzdz t gd d t dt dzdid fk(sJt t t ztzt t tgt t zd td id fk(sJt t t ztztzt dzztt t tgt t zd td t dzd td id fk(sJt t t zt tzzt dzzd zt t tgt t zd t tzd t dzd d d id fk(sJt t t ztztzdztt t tgt t zd td d dtd id fk(sJt t t zttzztt t tgt t zd ttzd idfk(sJt t t dzzd zt t gt t dzzd d d idfk(sJt t dzt dzztdzzdz t t tgd dt dzd tdzd t dzd idfk(sJt t dzt dzztdzzt t tgt dzd t dzd tdzd idfk(sJy)Nc2ttdzdz dz S)Nr[rSrZ)r,rNrOrKrIrPz$test_classify_diop..<smAqDFQJ7rKctdSNrZr,rOrKrIrPz$test_classify_diop..=s }Q/rKcPtttztztzdz Src)r,wrNyzrOrKrIrPz$test_classify_diop..>s ac!eAgk(BrKcTttdztdzztdzzdz S)NrSZ)r,rNrgrhrOrKrIrPz$test_classify_diop..?s% adQTkAqD6H26M(NrKr[*rZi univariateinhomogeneous_ternary_quadraticinhomogeneous_general_quadratichomogeneous_general_quadratic cubic_thuerjbigeneral_sum_of_even_powers$homogeneous_ternary_quadratic_normal) r@rT ValueErrorNotImplementedErrorr,rNrgrhrfrOrKrItest_classify_diopry;s) 978 :/0  BC  NO AqD2a4", - aaQT2 & 26 66 6 1q ! Aq AaCAq>#D&F FF F 1q1q!t+ , Aq! qsAq!QT1a35V1X XX X 1qsQT)A- . Aq AaCAaCAqD!Q24U3W WW W 1q1r) * Aq! qsAq!QAq13T/V VV V 1qs # Aq! qsAqsA&(G(I II I 1a4! $ A1a4Aq!<)1 11 1 A1q!t+{; < Aq AsAqD!QT1adA68TAV VV V A1q!t+ , Aq AqD!QT1adA.0V1X XX XrKc` ttdk(sJtdtzdk(sJtdtzdk(sJttdzdk(sJtdtzdzdk(sJtdtzdzdk(sJtttztt fk(sJtttzd ztt fk(sJtttzd z tt fk(sJtd tztz d z d k(sJtdtzdtzzd z dtzd z d tzd zfk(sJtdtzdtzz d z dtzd z dtzd z fk(sJtd tzdtzz d z dtzd zd tzd z fk(sJtd tzd tzzd tzdtzfk(sJtdtzdtzzd tztfk(sJtdtzdtzzdz dtzdz d tzdzfk(sJtdtzdtzzdz dk(sJtd tzdtzzdtzz d zdtzd zdtzd zfk(sJtdtzdtzzdtzz d ztdtzdt zzd zd tzdt zzd zfk(sJtdtzdtzzdtzz d zddk(sJtdtzdtzzdtzzd z dk(sJtd tzd tzzdtzz dz tdtzdt zzdzdtzd t zzdzfk(sJtdtzdtzz dtzzdz dtzdztdt zzt fk(sJtdt zdtzzdtzztz tt t tzdtzdt zzdtzzfk(sJttdy)N)rrZrS)r[rFrj)r)r|)NN)rr}r}NNN c&ttdz SNr[r^rOrKrIrPztest_linear..osj1orK) rrNrgt_0rht_1rft_2r@rTrOrKrI test_linearrTsP a=D  ac?d "" " ac?d "" " a!e  %% % acAg ' )) ) acAg % '' ' a!e sd ++ + a!eai S3$K // / a!eai S3$K // / acAgk "e ++ + acAaCi!m $3BsFQJ(? ?? ? acAaCi!m $3AcEAI(> >> > bdQqSj1n %!C%!)RVaZ)@ @@ @ acAaCi QsUBsFO 33 3 acAaCi RVSM 11 1 acAaCi!m $3BsFQJ(? ?? ? acAaCi!m $ 44 4 acAaCi!A#o) *quqy!C%!).D DD D acAaCi!A#o) *sAcEAcEMA4EquqQTu}WXGX.Y YY Y acAaCi!A#o)4 0I == = acAaCi!A#o) *.@ @@ @ acAaCi!A#o) *sBsFQsUNQ4F3QRSVQVY[H[.\ \\ \ acAaCi"Q$&* +# 33;/L LL L acAaCi"Q$&* +S#)QsURPSV^VXY\V\E\/] ]] ] 9-.rKcVtdtztzdtzzdtzz dzddhk(sJtdtztzdtzzdtzzdztk(sJtd tztzdtzzd tzz d z d hk(sJtd tztzdtzz dtzz d z dhk(sJtdtztzdtzzdtzzdzhdk(sJtdtztzdtzzdtzzdztk(sJtttztztzdzdtftdfhk(sJt dtztzsJy)NrS"rrZ)i{)r}irr[rj6)r)r{r}8r>i_riiiiirir|r@rZrr{0)rrNrgrWtrrOrKrI%test_quadratic_simple_hyperbolic_caserrs ac!ebdlRT)A- .  " ## # ac!eackAaC'!+ , 55 5 c!eAg!mac)B. /G9 << < c!eAg1nr!t+b0 1i[ @@ @ ac!eackBqD(1, -2F FF F ac!eackAaC'!+ , 55 5 acAgkAo &B7QG*< << < r!tAv  rKctdtdzzdtztzzdtdzzzdtzzdtzzdz dhk(sJtd tdzzd tdzzzd tzzd tzz d ztk(sJttdztdzzdtzzdtzzdzdhk(sJtdtdzzdtztzz dtdzzzdtzz dtzz dz dhk(sJtdtdzzd tztzzd tdzzzdz hdk(sJy)Nrnr[rrmi3)rr{rjrSr} rr{r{rrlV)r rr{r[rZrZrZr|r)rrNrgrWrOrKrItest_quadratic_elliptical_casersW bAg!Ao1a4/"Q$6A=D E) TT T a1fqAvo!+bd2R7 8CE AA A adQTkAaC'!A#-1 2xj @@ @ bAg!Ao1a4/"Q$6A=D E( SS S bAg1Q&AqD025 6 / 00 0rKctdtdzzdtztzz dtdzzzdtzzdtzzdzsJtdtdzzdtztzz dtdzzzdtzzd tzzdz sJtdtdzzdtztzzdtdzzzd tzzdtzzdz sJtd tdzzd tztzztdzz dtzzd z sJttdzdtztzztdzzdtzzdtzzd zsJttdzdtztzz tdzzdtzzdtzzd zsJttdzdtzz dzsJy)Nrr[rr}rrrrjrSrZ)(check_solutionsrNrgrOrKrItest_quadratic_parabolic_casers 1QT6BqDF?R1W4qs:QqS@2E FF F 1QT6BqDF?R1W4qs:RTAAE FF F 1QT6BqDF?R1W4qs:QqS@1D EE E 2ad7QqSU?QT1AaC7!; << < 1a4!A#a% ?? ? 1a4!A#a% ?? ? 1a4"Q$;+ ,, ,rKctdtztzsJtdtdzzdtztzz tdzzdzsJtdtdzzdtztzz dtdzzzdtzz dtzz dzsJtd tdzzd tztzzdtdzzzdtzz dtzzd z sJtd tdzzd tztzzdtdzzzd tzz d tzz dz sJtdtdzzdtztzz dtzz d tzz dz sJtdtztzdtdzzz dtzz dtzz d z sJttdztdzz dtzz dtzz sJttdzdtdzzz dtzz dtzz sJtdtdzzdtdzzz dtzz d tzz dz sJy)Nrrjr[r}r|rSrr rr rrrrrrOrKrItest_quadratic_perfect_squarers% 2a46 "" " 1QT6AaCE>AqD014 55 5 2ad7QqSU?Qq!tV3QqS82a4?"D EE E 2ad7RT!V+a1f4qs:QqS@2E FF F 1QT6BqDF?Qq!tV3bd:RTABF GG G 1QT6AaCE>AaC/14q8 99 9 3q57Qq!tV+ac11Q36; << < 1a4!Q$;1,qs2 33 3 1a4!AqD&=1Q3.14 55 5 1QT6AadF?QqS02a47!; << solverOrKrI!test_quadratic_non_perfect_squarers 1a4!A#:!Q$. // / 1QT6AadF?QqS01Q36 77 7 1a4!A#:1,qs2 33 3 1a4!AqD&=1Q3.14 55 5 1a4!Q$;1,qs2Q6 7 = = ?H: MM MrKc ddtzdtzdz zz tdtzdz zz}ttf}t|D]1}t|j t t ||s1Jy)Nir[rSrZrNrgrrJxreplacedictziprHvsols rItest_issue_9106rsm qsAaC!G} q!A#'{ *B AA2A T#a+->!?@@@ArKc tdztz tdzz }ttf}t|D]1}t|j t t ||s1Jyrrrs rItest_issue_18138rsW AAqDB AA2A T#a+->!?@@@ArKc>tdtdzzdtztzzdtdzzz dtzz dtzz dz sJtdtdzzdtztzz dtdzzzd tzz dz sJtd tz d tdzzztz d tztzz d tdzzz sJtd dtzzdtdzzzdtzzttzzdtdzzz sJy)Nrr[rrrrrrr}rrjrSrZrrOrKrItest_quadratic_non_perfect_slowrs 1QT6BqDF?Qq!tV3bd:RTABF GG G 2ad7QqSU?Qq!tV3ac9A= >> > 26AadF?Q.1Q61a4? @@ @ 1qs7Qq!tV+ac1AaC7!AqD&@ AA ArKc"tdddgk(sJtddgk(sJtdddgk(sJtd d dgk(sJtd d gk(sJtdd gk(sJtdddtfgk(sJtdddtztfgk(sJtddgk(sJtdddgk(sJtdddgk(sJtdddgk(sJtdd ddgk(sJtddgdk(sJtd ddgk(sJtdddgk(sJtdddgk(sJtd dd!gk(sJtd"dd#gk(sJtd dd$gk(sJtd%dd&gk(sJtd dd'gk(sJtd"dgk(sJtd(dd)gk(sJtd*dd+gk(sJtd,dd-gk(sJtd.dgk(sJtd d/gd0k(sJtd1dd2d3d d/d4k(sJtd d5gd6k(sJttd7d8hd9k(sJtd d:d;gk(sJtdd?gk(sJtd@dAdBgk(sJtd@dCgk(sJtdDdCdEgk(sJtdFddGddDdCsJtdDdHdIdJgk(sJtddtdtfgk(sJtddtgk(sJy)KNrSrrrr~irr[rZrrrrr}rrr)rrjiL)rrrr{rZr[))rr})r}r)rr)rrrZ)iri)iixi)lINd?lj\i)ii"i)lWRlQ]>I2lAAJtWjU^ )!r)')rr}r)rr}i")rrZ Sl!K2luGoXrr)rSrZ)m$rrrrTr))=rr)rrSr>iXiipi"'Mi)iSl)fMOl?_l;>x<i|rrZr)i ir)rrZ)rrZr[rr|{)rrZrri)rZ)rrZ)rrrrWrOrKrItest_DNrs 1a=VH $$ $ 3 r !! ! 3 x '' ' 3 x '' ' 3 r !! ! 1a=B   1a=aVH $$ $ 1a=ac1XJ && & 2r?b  1c?wi '' ' 1d y (( ( 1a=VH $$ $ 2q>ff- -- - 2s C CC C 2q>j\ )) ) 3?}o -- - 3??@ @@ @ 3?}o -- - 3?^_ __ _ 2q>gY && & 2q>i[ (( ( 2r?wi '' ' 3 r !! ! 2r?wi '' ' 3 y (( ( 5" "?!@ @@ @ 2r?b  2r?< << < b!S#r2 .$ 66 6 2r?H HH H wsB %F FF F 2r?{m ++ + 3 r !! ! 2r?w/ // / 3?wi '' ' 3 r !! ! 3 y (( ( b!S!S" -- - 3 7 3 33 3 1a Ax '' ' 1b!  "" "rKc*tddgdk(sJtddgdk(sJtddgk(sJtdd d gk(sJtd d d dgk(sJtdddgk(sJtdddgk(sJtddtdgk(sJtddtdtfgk(sJtddtdtztfdtztfgk(sJtddtdgk(sJtd dtgk(sJy)Nrr)r)rrZrr)r)rrSr)irrr|rZ)lu L"lEYrrrZrrr{riT)irrrrjr[rS)rrrOrKrI test_bf_pellrs8 b" !< << < b" !I II I c2 " $$ $ dA #B"C CC C b" &'!2 22 2 eR %B$C CC C c2 9+ -- - b!Q F8 ++ + aA Aq6( ** * aA AaC8bdAY"7 77 7 aA 6( ** * aQ 2 %% %rKctddddk(sJtddddk(sJtddd dk(sJtddd d k(sJtd d d dk(sJtdd ddk(sJy)Nr[rZrr|rjr}rSr~rrrrrr)r"rOrKrI test_lengthr#s !Q?a   "a q  "a  !! ! !Q q  !R  "" " !Q?a  rKc t|\}}|tttgz|zd}|tttgz|zd}t |j t ttf||f}t|jDcgc]$}t|jttg&c}}ttzttfD]}||jvsytdztdzdfD]}||jvsd||<|tdzdk7r;t|tdz|tdzxrt|d|tdzSycc}w)a9 Test whether X*Y, X, or Y terms are present in the equation after transforming the equation using the transformation returned by transformation_to_pell(). If they are not present we are good. Moreover, coefficient of X**2 should be a divisor of coefficient of Y**2 and the constant term. rrZFr[T)r*r XYrJsubsrrNrgrargsreversedas_independentkeysr) rHABur simplifiedrcoeffterms rIis_pell_transformation_okr,sS  #DAq 61a&> A q!A 61a&> A q!ArwwsAq6Aq6':;1%Dr!t%Ka%O PP P $SAX1Q%6A%=!%Cc!e%Kc%Q RR R $QTAaCE\3q!t8%;ac%ABqD%H2%M NN N $R1Wqs1u_r!Q$w%>A%E1%Lt%S TT TrKcttdzdtzz tdzz dk(sJttdzdtdzzz dz dk(sJttdzdtztzz dtdzzz dz dk(sJtdtdzzd tztzz tdzz d z d k(sJtdtdzzdtztzz tdzz d z d k(sJtdtdzzdtztzztdzz dk(sJtdtdzzdtztzz tdzzdtzzdtzz dz dk(sJy)Nr[rrSr})rSr}rjr)r}rrr)rrr)rTr)rZrrrl)ei)r rNrgrOrKrI test_find_DNrUsi 1a4!A#:1$ % // / 1a4!AqD&=1$ % // / 1a4!A#a%AqD(1, - 99 9 1QT6AaCE>AqD(2- .' 99 9 2ad7QqSU?AqD( )V 33 3 3q!t8ac!e#ad*QqS01Q36: ; NN NrKcgd}|D].\}}t||\}}}||dzz||dzzz|dzk(r.JtddJtddJy)N) rrrSrri)rjr)rrj[rrrZr{)rjrrri)rZrir[r{r)r!r abrfrNrgs rI test_ldescentr&_sz 4A'11a.1aAv!Q$!Q$&&&' B  ## # Aq> !! !rKcZtdtdzzdtdzzz tdzz sJtdtdzzdtdzzztdzz sJtdtdzzdtdzzztdzz sJtdtdzzd tdzzzdtdzzz sJttdzdtdzzztdzz sJtdtdzzdtdzzztdzz sJttdztdzztdzz sJtd tdzztdzzd tdzzz sJtd tdzztdzz d tdzzzsJtd tdzzdtdzzzdtdzzz sJtdtdzzdtdzzz dtdzzz sJtdtdzzdtdzzz dtdzzz sJy)Nr[Arihr}rjrSrrrrrrrri-|r1rrNrgrhrOrKrI"test_diop_ternary_quadratic_normalr-js 3q!t8eAqDj01a47 88 8 2ad7SAX-14 55 5 1QT6AadF?QT1 22 2 1QT6AadF?Qq!tV3 44 4 1a4!AqD&=1a4/ 00 0 1QT6AadF?QT1 22 2 1a4!Q$;A- .. . 2ad7QT>Bq!tG3 44 4 1QT6AqD=2ad72 33 3 3q!t8bAg-!Q$6 77 7 2ad7Qq!tV+c!Q$h6 77 7 3q!t8bAg-QT 9 :: :rKc t|}|ttttgz\}}}t |j ttttf|||f}t|jDcgc]}t|j|||gc}}||z||z||zfD]}||jvsyycc}w)NFT) r+r rNrgrhrJrrrrrrr) rHr rrZr rr rs rIis_normal_transformation_okr0ys $A1ay!!GAq!rwwsAq!9q!Qi'@ABJ *//RQ(+1++aAY78R SE1ac1Q3 5::<   Ss6"Ccttdzdtdzzztdzzdtztzz dtztzz dtztzzsJttdzdtdzzzdtdzzz sJttdzdtztzzsJtdtdzzdtdzzz dtztzz sJttdzdtztzzdtztzz dtztzzsJttdzdtztzzdtztzz ttzzsJttdztdzztdzzttzz ttzz ttzz sJttdzdtztzzdtdzzzsJtttzdtztzzdtztzzsJtdtztzdtztzzsJy) Nr[rSrrrdrr)r0rNrgrhrOrKrItest_transformation_to_normalr3s &q!ta1f}q!t';bd1f'Dr!tAv'MPRSTPTUVPV'V WW W &q!ta1f}s1a4x'? @@ @ &q!tbd1f} 55 5 &qAvAqD'82a46'A BB B &q!tbd1f}r!tAv'=1Q'F GG G &q!tbd1f}r!tAv'=!'C DD D &q!tad{QT'9AaC'?!A#'E!'K LL L &q!tac!e|a1f'< == = &qsQqSU{QqSU': ;; ; &qs1uqs1u} 55 5rKc* tdtdzztdzztdzzdtztzz sJttdztdzz tdzz ttzz ttzz sJtdtdzzttzz ttzz ttzz sJttdzttzz ttzz sJtdtdzzdtztzz ttzz sJtdtdzzdtdzzz ttzz sJtdtdzzdtdzzztdzz dtztzz dtztzzdtztzz sJtdtdzzdtztzz sJtdtdzzdtdzzz dtztzz tdzz sJttdzd tdzzz tdzz d tztzzdtztzz sJtd tdzzdtdzzzdtztzzdtztzzdtztzzsJttdzdtdzzztdzzttzz d tztzz sJttdzdtdzzztdzzttzz d tztzz dtztzzsJttdzdtdzzztdzzd tztzz d tztzz dtztzzsJtttzdtztzz d tztzzsJt tdztdzztdzzdk(sJt tdztdzzJt t ddtztzdtztzz dtdzzzdtztzz dtdzzz}t|dk(sJt dtdzzdtdzzztdzz dk(sJtttzdtztzzddtfk(sJdtztzdtztzz dtztzz dtdzzz}t|dtdzzdtztzz dtztzdtdzzzdtztzfk(sJtttzdtztzzdtztzt tdzztdzzttzfk(sJy)Nr[rjrSr}rrrr1rrkrrrc ttttfttzdtdzdtdzdtdzdiS)NrZr[rSr)r6rNrgrhrOrKrIrPz-test_diop_ternary_quadratic..s;&1ay 1aAq!Q$1a4+ -rKr|r)rr[r)rZrr[rr~r) rrNrhrgrr@rwrrDr3pqrrGs rItest_diop_ternary_quadraticr9s+ 1QT6AqD=1a4/!A#a%7 88 8 1a4!Q$;A-!3ac9 :: : 1QT6AaCAaC/ 00 0 1QT6AadF?QqS0 11 1 1QT6AadF?QT1AaCE9AaCEAAaCEI JJ J 1QT6BqDF? ++ + 2ad7Qq!tV+ac!e3ad: ;; ; 1a4"QT'>AqD02a4691Q3q5@ AA A 2ad7Qq!tV+ac!e3ac!e;ac!eC DD D 1a4!AqD&=1a4/!A#51Q> ?? ? 1a4!AqD&=1a4/!A#51Q>AaG HH H 1a4!AqD&=1a4/"Q$q&82a46ABqDFJ KK K 1Q31Q;Aa/ 00 0 (A1q!t); <@R RR R (A1 5 == = :-. Aa!A#a%!AqD& 1Q3q5 (1QT6 1B !" % 22 2 (1a4!AqD&1a4)? @   !!A#!A+ . Q    Aa!A#a%!A#a% !AqD& (B ( , 1a4!A#a%Aa!AqD&!A#a%0 11 1 acAaCEk "r!tAvs1a4x!Q$!&D DD DrKctdtdcxk(rdk(sJJtddk(sJtdtdcxk(rdk(sJJtdtdcxk(rdk(sJJtd td cxk(rdk(sJJtd dk(sJtd d k(sJtddk(sJtddk(sJtddk(sJttd dk(sJy)NrZr{rr}r~rjrr[rrrrrS4r5rirl)r'r rOrKrItest_square_factorr<s  }R0 5A 55 55 5  q   }R0 5A 55 55 5  }R0 5A 55 55 5   c 2 7a 77 77 7  q    !! !   !! !   !! !   ## # 3 (A -- -rKc$ttdztdzztdzz sJttdzdtztzztdzzsJtdtdzzdtdzzz tdzz sJtdtdzzdtdzzztdzz dtztzz dtztzzdtztzz sJttdztdzz tdzz sJttdzdtdzzz tdzz dtztzzd tztzz sJtd tztztdzzsJtd tdzzd tdzzz d tdzzz sJtd tdzzdtdzzz dtztzz dtztzz dtztzz sJtdtdzzdtdzzzdtztzzdtztzzdtztzzsJtd tdzzd tdzzz d tdzzz sJy)Nr[r(r)rSr}rr5rrr*rr+rrrkr,rOrKrI"test_parametrize_ternary_quadraticr@s@ 1a4!Q$;A- .. . 1a4!A#a%AqD02a469AaCEA BB B 1Q3q51a4< (( ( 3q!t8bAg-QT 9 :: : 3q!t8c!Q$h.Aa7"Q$q&@2a46I JJ J 2ad7Qq!tV+ac!e3ac!e;ac!eC DD D 3q!t8bAg-QT 9 :: :rKc,tdtztzttzzdtztzz sJtdtztzdtztzz dtztzz sJtdtztzdtztzzsJtttzttzzttzzsJtdtztzdtztzzdtztzzsJy)Nr[rSiYrrrr,rOrKrI test_no_square_ternary_quadraticrCs 1Q3q51Q3;1Q. // / 3q57SU1W,r!tAv5 66 6 2a46BqDF? ++ + 1Q319qs? ++ + 2a46BqDF?RT!V3 44 4rKcgd}|D].\}}t||\}}}||dzz||dzzz|dzk(r.Jttdttdttdy)N) rrrrrrr r!r"r[ctddS)Nr{rrrOrKrIrPztest_descent..s gb"orKctddS)NrrSrFrOrKrIrPztest_descent..s gamrKctddS)NrjrSrFrOrKrIrPztest_descent..s gamrK)rr@rTZeroDivisionErrorr#s rI test_descentrJsu gA'1!Q-1aAv!Q$!Q$&&&'  9-. 34 9+,rKctttz ttz zttz zsJtttz tdztdzztdzz zsJttdtzz dtzztdztdzztdzz zsJttdzdtdzzz dz sJttdzdtztzzsJttdzdtztzz tdzzsJtttdztdzz dz zsJttdtzdtzz dzzsJttdzdtdzzz dz tdztdzz dz zsJttdzdtdzzz dz tdtzz zsJttdztdzztdzz tdtzz dtzz dtzzzsJttdzdtztzz dtztzzsJttdzdtzz dzsJt ttz t t ttk(sJtdztdzzdz }t |d t | d k(sJt dtztzdtztzz dtzdtzfhk(sJtdztdzztdzzd z }d h}t ||k(sJtt|j}t |d |k(sJt tdzttdd zzdz tk(sJd tdzzdtdzzz tdzz }|j}ttttf|dhk(sJt |dtdzzdt dzzzdtdzzdtzt zz dt dzzzdtdzzdtzt zz dt dzzz fhk(sJdtdzzdtdzzztdzz }|j}ttttf|dhk(sJt |dtdzzt dzz dtdzzdtzt zzt dzz dtdzzdtzt zz dt dzzzfhk(sJdtdzzdtdzzzdtdzzz }|j}ttttf|dhk(sJt |dtdzzdt dzzz dtdzzd tzt zzd!t dzzz d"tdzzd#tzt zz d$t dzzzfhk(sJd%tdzzd&tdzzzd'tdzzz }|j}ttttf|d(hk(sJt |d)tdzzd*t dzzz d+tdzzd,tzt zzd-t dzzz d.tdzzd/tzt zz d0t dzzzfhk(sJtdzdtdzzzd1tdzzz }|j}ttttf|d2hk(sJt |d3tzt zdtdzzd3t dzzz tdzd1t dzzzfhk(sJt#t$d4t dtz tk(sJt dtz dtz zt&j(z hd5k(sJt tdztdzzdtzzdz d hd6k(sJt tdztdzzd7z d8z ttfd 9hd:k(sJt tdztdzzd7z d8z ttfd 9hd:k(sJttdztz sJttdztz sJt tdztz t*t*t*dzfhk(sJt tdztz t*t*dzt*fhk(sJy);Nr[rSrrZrmr}rjaTpermuterlrZr[rS\c)rr3i{riKfiijrR)rZrZr[r|i9)iiU i iiiiiaiiڙiV i=i i)r5i i>iipEi*i',iir)rr[rZrc:tttdzzdzS)Nr[rZrrNrgrOrKrIrPz"test_diophantine..-s AadFQJ(?rK>rSrrjrjrrSr|rZr>rjrZrjr{rrZrr{r{rZrrrrQ)symsrN>r[rrSr|rr[rr|r|rSr|rr[rSrSr[)rrNrgrhrfrrrrrWrBpopras_coefficients_dictr6r7r8r@rxr Halfr)rHbase_sol complete_solnr s rItest_diophantinerrs AEAE?AE2 33 3 AEAqD1a4K!Q$$67 88 8 A!GacMAqD1a4K!Q$,>? @@ @ 1a4!AqD&=1, -- - 1a4!A#a%< (( ( 1a4!A#a% ** * q1u R1X!6 66 6 A1r B r4 (KT,J JJ J qs2v!B 'QsUAcEN+; ;; ; A1q!t b B{H r?h && &+HLLN;> > AqD2ad7 QT !B  # # %E )1a)U ; <   r? AqD1QT6419s1uQw.1a47 QT DF1Hr!Q$w& () )) ) 1a4!AqD&1a4 B  # # %E )1a)U ; ;   r? !Q$A r!Q$w1QA- !Q$1Q1a4 !" "" " QT"QT' #ad( "B  # # %E )1a)U ;     r? AqD[6!Q$; &1 wqy{(B AqD[) A+q 2VAqD[@ B C CC C QT#ad( 3q!t8 #B  # # %E )1a)U ;     r? !Q$Yq!t #VAqD[5719%> > AqD51QT6>AadF2Qq!tV;bAgE FF F 2ad7QT>AadF2QT9Bq!tGC DD D adQTkAqD01a47 8 > >1aQR) > T TAqD[1a4 1Q!Aq!tad{QT/A B C DD DrKc tddD].}ttdtd|zD|z r.Jt t dzt dzzdz Jt t dzt dzztdzzdztk(sJt dzt dzztdzzdz t dhk(sJtdztdzzt dzzt dzztdzzdz tt ddk(sJtd td d fz}t|Dcgc]}|dz c}d z hd}t|k(sJttd dk(sJtdt dzz t dzz tdzz dhk(sJtdztdzzt dzzt"dzzdz t%t&fdycc}w)NrSrc3&K|] }|dz ywr[NrO.0is rI z9test_diop_general_sum_of_squares_quick..Ts"DA1a4"Dz:%ir[rlrOi!z:56Tnegativep>rrr[r[r[rrrr[rrrrrrjrjrjrrrZrZrZrSrrrZrZr[r}rrrZrZr}rrrrZr[rSrrrrZr}r}r}rrr[r[r[rrrr[rSrSrSrrr[rSr}r}rrrSrSrSrrrZrZrZrSrrrZrZr[rSrjrrZrZrSrjrrrZr[rSrSr}rrZrSrjr}r}rr[r[rjrjrrr[rSrSrjr}rrMir)r[r[r[rjct SrFrdrGsrIrPz8test_diop_general_sum_of_squares_quick..os rc(:rK)rangersumr rrNrgrhrWr rlenrrr$r%rurvr@rx)rvar base_solnrHs @rI&test_diop_general_sum_of_squares_quickrRs 1b\Js"D1C"DDqHIIIJ 'q!tad{Q 7 ?? ? &q!tad{QT'9A'= >#% GG G A1q!t y )B &r * ;   A1q!t ad "QT )D 0B *2q1 2a 77 7 $-7368 8C S!q!t! "S (BPI r?i '' ' {2t, - 77 7 rAqDy1a4'!Q$. /I; >> > A1q!t ad "Q &B  :;!"s8 G'cltdztdzztdzzdz }hd}t||k(sJy)Nr[i@B>irrrrrX r`irrr@r`hirrihrr`rrrrrrrrrrrirriXrrrrrr)rNrgrhr)rHrs rItest_issue_23807rrs= A1q!t g %BBI r?i '' 'rKcdD]4}tddD]#}t||D]}t||k(rJ%6ttddgk(sJtdddDcgc] }t|c}gdgdgdgk(sJttd d gk(sJttdd d gk(sJtdDcgc] }t|c}gd dd gdggk(sJycc}wcc}w) N)rrrZrrSr})rrrrrS)rrrrZr[)rrrZrZrZrrO)rZrZrZr[)rr$rlist)nkr7rs rItest_diop_partitionr}s #q! #Aq!_ #1v{"{ # ##  !Q B && &&q!Q/ 0DG 0/5; ;; ;  ! " %% %  !Q RD (( (&q\ *DG *y1a&1#.F FF F 1 +s #C+C0cdD] }t|\}}|dz|dzz|k(r JtdJtd}|dk(rt|dtusJy)N) r}rrr%ri% i iiYr[ri5 )iir)r&typeint)rr$r%anss rI test_prime_as_sum_of_two_squaresrsv > *1-1!tad{a  'q ) 11 1 %f -C * c!f!44 4!4rKcdD].}t|\}}}|dz|dzz|dzz|k(sJ|dk\r.JttdtdJtdJtddk(sJtd d k(sJy) N)rrZr[rrlmyl85Jdl8l9_ ri!i"i#i$i%i&r[rctdSNr{)r)rOrKrIrPz+test_sum_of_three_squares..s 3B7rKri<r)rrr}rjrrr[)r)r@rw)rr$r%rus rItest_sum_of_three_squaresrs1&q)1a!tad{QT!Q&&&Av v  :78  " ** *  * 22 2  #y 00 0  "i // /rKcbddlm}|dd}tdt|D|k(sJt t dt dD]a}t|}t|dk(sJtd |DsJtd |D|k(sJt|t|k(raJy) Nr)randintrZl@ kc3&K|] }|dz ywr~rOrs rIrz+test_sum_of_four_squares..s4q!t4rctdSr)r(rOrKrIrPz*test_sum_of_four_squares..s 226rKrrjc3&K|] }|dk\ ywrNrOrrs rIrz+test_sum_of_four_squares..s*a16*rc3&K|] }|dz ywr~rOrs rIrz+test_sum_of_four_squares..s(A1a4(r) sympy.core.randomrrr(r@rwrrallrsorted)rrresults rItest_sum_of_four_squaresrs) ?#A 4034 4 99 9 :67 4[.$Q'6{a*6****(((A---F|vf~--- .rKc 8gd}|D]K}|\}}}t|||} t|}t||k(sJd}|D] }|||zz} ||k(sJ 6t tddddddgk(sJt t d t t d t t d t td ddgk(sJt td d d dgk(sJt tddd gk(sJt tddd dgk(sJt tddddddgk(sJt tddddgk(sJt tddddgk(sJt tdddgk(sJt tdddddgk(sJt tddddgk(sJtt tddddk(sJd} dD] } t t| d| | z gk(r Jy#t$rYwxYw)N))rrSr[)r(r[rj)r[rZr[)rSrZrS)r}r[r[)i@0r[rj)ir[rSTrrr[rj)rZrZrSrS)rrr[rjc.ttdddS)Ng333333?r[rr%rOrKrIrPz+test_power_representation..st$8a$CDrKc.ttdddSNr[rrrOrKrIrPz+test_power_representation..t$8Aq$ABrKc.ttdddSrrrOrKrIrPz+test_power_representation..rrKr{rZr\rSr]rcr)zeros)rZr[r[r[r[r[)rrrrrrSr}Fr|rripi@)rrrr}rjr[rZ)r%nextr StopIterationrr@rwr4) teststestrr7rflchk_suml_ibigrs rItest_power_representationrs E1a Aq ) G1v{"{/C%Q.G/!|#| " $RAt4 5 |$ %% % :DE :BC :BC $RA. /2 55 5 $Q1- .4& 88 8 $Q1- ." 44 4 $Q1- .4& 88 8 $T1at< = /0 11 1 $T1au= >" DD D $RA. /H: == = $RA. /2 55 5 $Q1d3 4 @@ @ $Q1e4 5 ;; ; t()91a@A Ba GG G C $:M#q#'23r999:/!  s3H  HHctddd\}}t|dz||zzdz }|hdk(sJtddd\}}t||zd|zzd |zzd z }|hd k(sJy ) z< Test whether diophantine respects the assumptions. zm nT)rCpositiver[i>irZr}rrr_r}yrjr[za bFrSr>rrrrrr~rrr~rrrN)r r)mrdiofr$r%s rItest_assumptionsrs 5$ 6DAq q!taczC' (D M MM M 5$ 7DAq qsQqSy1Q3* +D Q QQ QrKc &t|}tj|}t|j}|j t |rD|j}|D]+}t|jt||dk(s+ny|rDy)z Determines whether solutions returned by diophantine() satisfy the original equation. Hope to generalize this so we can remove functions like check_ternay_quadratic, check_solutions_normal, check_solutions() )keyrFT) rr make_argsr free_symbolssortrrmrJrr)rHsfactorsrsolutionrs rIrrs BAmmBG r CHH!H" 557 AQVVCX$678A=   rKc Ddtztzdzdz}t|tdtzdz fhk(sJdtdzzdtztzzdtzzdtdzzzdtzzdz}t|tt dz fdtzdz tfhk(sJt ttdzzdz tdz dztfhk(sJt ttzdz tdtz fk(sJtd ddd d k(sJdtzdz dtzdzf}tdd dt|k(sJtdd dd td |Dk(sJtdddtk(sJtddddhk(sJtddddhk(sJttdtddddk(sJtdz tdzz tdzz tdzz}t!|t!| cxk(r\t"dzt$dzzt&dzz dt"zt&zdt$zt&zt"dzt$dzzt&dzzfk(sJJt)t+t-dtdz zt-dtdz zt-dtgd k(sJt)t+t/ddt-dtzt-dtgd k(sJt)t+t-dtzt/ddt-dtgd k(sJt1dddk(sJt3dt5d cxk(rdk(sJJt3d t5d  cxk(rdk(sJJt7ddddk(sJtt8dt;ddddk(sJtt<dtt<d t tdztdzzdz d!z hd"k(sJy)#Nr[rZr|rrrjrrSrrrrc3HK|]}|jtdywr)rr)r_s rIrz$test_diopcoverage..s:UA166!Q<:Us "rr}rrrlctdddS)NrjrrZ)r#rOrKrIrPz#test_diopcoverage..s{1b!4rKrrrTFrOctdS)N)r[rjr)r:rOrKrIrPz#test_diopcoverage../s k)4rKrkni6)rrZr}c~ttdztdzzttzzdtztzzdz S)Nr[rrrNrgrhrOrKrIrPz#test_diopcoverage..4s3 AqD1a4K!A#4E!A4MPR4R(SrKc:ttdztdzzS)NrSr[rYrOrKrIrPz#test_diopcoverage..5s AqD1a4K(@rKrc>rZrrZrrr{rr{r{rr{rrr)rNrgrrrrrr-tupler.rWr@rwr0rfrhrm1m2m3rr2r rr7r8r9r:rTr/rx)rHrs rItest_diopcoverager% s A#'A+ B b>sBsFQJ/0 00 0 1a4!A#a%"Q$ 1a4 '"Q$ . 3B b>q1"q&kBqD1Ha=9 99 9 !ad(Q, 'adUQYN+; ;; ; q1uqy !c1s7^ 33 3 1a4 0F :: : Q37BqD1H C 1b! , 33 3 1bD 1U:UQT:U5U UU U aB 35 (( ( aA 6( ** * aB F8 ++ + :45 1a #v -- - Q$A1 q!t #B #B ' "% , URU]RU "AbDG bDGRURU]RU* , ,, ,, , {1Q4!A#:qtacz1Q4!= >! CC C {8Aq>1Q4!8QqTA3? @A EE E {1Q4!8Xa^QqTA3? @A EE E "b !Q && & 758| , ,, ,, , 758| - -- -- - q!Q 9 ,, , 945 h 2    ST  @A !Q$A+,t3 4 S TT TrKcztdddddddk(sJtdddd d d d k(sJttd y) Nr[rrrjOrr[rrrrZr)r[rr[c"tddddddS)Nr[rrlrjr'r)r1rOrKrIrPztest_holzer..CsvaB2r:rK)r1r@rwrOrKrI test_holzerr*:sK !QAr2 &* 44 4 !Q1a $ 11 1 ::;rKcB fd}d\ dx\}}}}dx\}}}}||dk(sJ||dk(sJt |dzz |dzz |dzz z zksJt |dzz |dzz |dzz z zksJt||| } | |k(sJy)Nc8|dzz|dzzz|dzzz SrrO)rNrgrhr$r%rus rIrPz"test_fail_holzer..Hs(1a4!AqD&1QT61rK)rjr'r)rZrr(rr[)maxr1) rHrNrgrhxyzrrr/rhr$r%rus @@@rItest_fail_holzerr1Fs 1BGAq!GAq!cGAq!c s8q== s8q== qAvqAvqAv &!A#a% // / qAvqAvqAv &!A#a% // /q!Q1a A 8O8rKctdtzdtzzdtzztz t t t tzdt zdt zzdtzzfhk(sJy)Nrrrr)rrfrgrNrhrrrrOrKrItest_issue_9539r3Ts[ qsQqSy2a4'!+ , #sSy!C%"S&.1S5"8 9 : ;; ;rKctdtdztdzztdzzzdttzttzzttzzzz dhk(sJy)NrSr[rlrrrrrOrKrItest_issue_8943r6YsY  1a4!Q$;A QqS1Q3Y1_!55 7 ;  rKctdztdzztdzzdz tddhk(sJt dddhk(sJt t fdtdd }tdztdzz|dzzdz t d hk(sJttdztdzzdztk(sJt tdztdzzdz d tk(sJy) Nrjiq )rSrr)r[rjrr[ct dSr)rrGsrIrPz.test_diop_sum_of_even_powers..cs(GQ(OrKnegTr)rrrr)limit) rNrgrhrrr@rxr rrW)r9rHs @rItest_diop_sum_of_even_powersr;_s A1q!t d "B b>i3 33 3 *2q 1i5K KK K  OP %$ 'C A1sAv  $B *2 .:, >> > q!tad{Q '35 00 0 *1a4!Q$;?! D MM MrKc hd}tdztdzztdzztdzztdzzdz }t |t }t|dk(sJ||k(sJttdttdddgk(sJttddgk(sJttdddd gk(sJttd dgk(sJttdd gk(sJttd d dd gk(sJttd d gk(sJttdddgk(sJttddgk(sJttddddgk(sJttdddddgk(sJttdddgk(sJtdDcgc]!}ttt|dd#c}gdk(sJtdDcgc] }ttt|d"c}gdk(sJtdD][ttdd}|rtfd|DsJttd}tfd|Dr[Jttdttdttd d dd!gk(sJttd ddgk(sJttddddgk(sJttddd dd"d#gk(sJttddddgd$k(sJttd%ddgk(sJttd&d dgk(sJttd'd dd(gk(rd'd'k(sJttdd)zddgk(sJycc}wcc}w)*N>rrrZrZrrrr}rrrrZrSrrrrZrjr}rrrSrjrrrrSr}r}rrZrZr[rrrZrZrrrrZr[rSrSrrZrSrjrjrrZr}r}rrr[r[rSr}rr[rSr}rrrSrSrjr}rr[rrlc,ttddS)Nrr{)rr5rOrKrIrPz,test_sum_of_squares_powers..tstN2r$:;rKrZr\T)rrZrrSrr5rjr]r}2)r}r})rZrr)rZrZrZr[r[)rrrZrZrSr)rZrZrZrZrZrZrZrZr)rZrZrZrZr[r[rZrZr[r[r[r[r[rSr[rZrSrSrSrSrjrSrSr[r[rjrjrjrjr})rrrrrrZrrrZrrrZrrZrZrrZrZrrZr[rZrZrZrZrZrZrZrZrSc3HK|]}td|Dk(yw)c3&K|] }|dz ywr~rOrjs rIrz7test_sum_of_squares_powers...s 1!A 1rNrrrrs rIrz-test_sum_of_squares_powers..s!C1S 1q 11Q6C"c3HK|]}td|Dk(yw)c3&K|] }|dz ywr~rOrOs rIrz7test_sum_of_squares_powers...s'q!t'rNrQrRs rIrz-test_sum_of_squares_powers..s!93'Q''1,9rSc.ttdddS)Nr[r{rZrr4rOrKrIrPz,test_sum_of_squares_powers..stM!R$;<rKc.ttdddS)Nr[rZr{rWrOrKrIrPz,test_sum_of_squares_powers..stM!Q$;<rKr|rr)rrZrZ))rr})rZrjrkri)rr)r rrNrgrhrrrr@rwrr5rrWrr4)trurHrrs1s2s ` rItest_sum_of_squares_powersr]ks -C A1q!t ad "QT )C /B %b" -C s8r>> #:: :;< q!$ %$ // / q!$ % ++ + q!T* +x 77 7 sA& '2 -- - q!$ % ++ + q!T* + { :: : q!$ % ++ + q!$ %$ // / q!$ % ++ + r1% &66*: :: : r1d+ ,1* ** * q!$ %*B)C CC C;@9 EaC^Aq$/0 1 EJ  6;2Y ?C^Aq)* + ?D  2Y: 1d+ ,CCCCC 1% &9b9999 :  :<= :<=  b!Q' (XJ 66 6  b!Q' (B .. .  aA& 'D6 11 1  aAt, -)Y1G GG G  aAt, -1I II I  aA& '2 -- -  dAq) *b 00 0  dAq) *tf 4$$,G G  agq!, - 33 3= F @s &N?%OctdddusJtdddusJtddsJtddsJtddsJtddsJtdddusJy) NrSr{FrrZrrjr[)r;rOrKrItest__can_do_sum_of_squaresr_s !!R (E 11 1 !"a (E 11 1 !!Q '' ' !!Q '' ' !!Q '' ' !!Q '' ' !!Q '5 00 0rKcddlm}m}m}m}m}|dz|dzzdz }dh}t ||k(sJtt|j}t |d|k(sJ|dz|dzz|dzz|dzz|dzzd z }tt |d k(sJtt |dd k(sJhd }t d tdzzdtztzzdtdzzzdz d|k(sJy)NrrtrjrLrkTrMr[r(#i0rrrr) rzr$r%rurvrwrrWrBrmrrNrg) r$r%rurvrwrHrprqsolns rItest_diophantine_permute_signrcs'' A1 $BxH r?h && &+HLLN;.sk"Aq6:rK)rNrgrrr@rTrGs@rItest_issue_9538risH QqS1B r1 &C3+;*< << < 9:;rKctdtdzztdzzdtdzzz }t t |j \}}}||dzd|dzzz d|dzzd|z|zzd|z|zz d|dzzz |dzd|z|zz d|dzzzd|z|zz fhk(sJttdzdtdzzzdtdzzz }|d|z|z|dzd|dzzz |dzd|dzzzfhk(sJtdtdzzdtdzzztdzz }|d|dzz|dzz d|dzzd|z|zz|dzz d|dzzd|z|zz d|dzzzfhk(sJtdtdzzdtdzzzdtdzzz }|d|dzzd|dzzz d |z|zd |dzzd|dzzzfhk(sJtdtdzzdtdzzztdzz dtztzz d tztzzd tztzz d|dzzd|z|zz |dzz d|dzzd|z|zz|dzz d|dzzd|z|zz d|dzzzfk(sJtd tdzzdtdzzz dtdzzz d|dzzd|dzzz d|dzzd|z|zzd|dzzz d|dzzd|z|zzd|dzzzfk(sJy)Nr[r|rjrSHrrrrrr}rr*rr+i~ii ixi: iiib<)rrNrgrhrr r r3)rr7r8rs rItest_ternary_quadraticrlsAadFQTMAadF*+Aad''(GAq!  1qAv  1a4!A#a%!A#a%!AqD&( 1qs1u qAv!A%'( (( ( AqD1QT6MAadF*+A !A#a%A!Q$1qAv 67 77 7AadFQq!tVOad*+A !AqD&1a4-AqD1Q3q51a4!7 !Q$1Q1a4!" "" " AadFR1W$r!Q$w./A "QT'AadF"AaCE1QT6AadF?;< << < ( !Q$1a4!Q$1Q&1Q.1Q6 8 !Q$1QAqAv!A~14a1f !A7!Q$7<   ) AqD2ad7T!Q$Y& ( ad VAqD[ $q!t)eAgai"7q!t #AXQq(5A:5,7 77 7rKc& tggttk(sJjdk(sJjdk(sJt t fdt jgk(sJtttgttgjttfk(sJjttfk(sJt t fdjdtdhk(sJjddtftddtfhk(sJt t fdt jtdttitdtdigk(sJttttgttgtjd k(sJjtd ztzttz d fttd ztzttz d fhk(sJj!td tdzd tz d fhk(sJd tdzd tz d fhk(sJj!td td id hk(sJd d d hk(sJj!tditdzdtz d fhk(sJdtdzdtz d fhk(sJj!tdtd zdz tdzd fhk(sJddtd zdz tdzd fhk(sJj!td td idhk(sJd d dhk(sJj!tdtdidhk(sJdddhk(sJt t fdt t fdt t fdt t fdt t fdt t fdt t"fdtttgttg}|jtdtzf|jt dtzf|dddhk(sJy)NrOc0jtfSrF)addrN)r[srIrPz/test_diophantine_solution_set..srvvqd|rKc&jdS)Nr\ro)r\srIrPz/test_diophantine_solution_set..srvvd|rK)rSrjr{c&jSrF)update)r[r\srIrPz/test_diophantine_solution_set..sryy}rKrSrjr[rZrr)rTr{rZr}rr)rrrZ)rrZc(jdS)NrZ)rNrs3srIrPz/test_diophantine_solution_set..srwww|rKc*jdddSNrZr[rSrvrwsrIrPz/test_diophantine_solution_set.. srwwq!Q/rKc&jdS)NrOrqrwsrIrPz/test_diophantine_solution_set.. srvvbzrKc&jdS)N)rZr[rSrjrqrwsrIrPz/test_diophantine_solution_set.. srvvl3rKc&jdS)NrrqrwsrIrPz/test_diophantine_solution_set.. srvvf~rKcdddSrzrOrwsrIrPz/test_diophantine_solution_set.. sr!Q{rKcdS)NrZrrOrwsrIrPz/test_diophantine_solution_set..s b1grKrrtrr)r<rWr ryr@rwr dict_iteratorrNrgrr rorsrhrrrT)s4r[r\rxs @@@rItest_diophantine_solution_setrs B 'B r7ce   ::   ==B   :+,   " #r )) ) AA /B ::!Q   ==QF "" " :+,FF6N r7vh  IIfr1g r7vAw' '' ' :,-   " #B1 1a|'D DD D Aq Aq6 2B r}}  "" "FFAqD1Ha!eQ  r71q!a%+, ,, , 771a=a!eQUA./ // / a5a!eQUA&' '' ' 77Aq!Q< [M 11 1 a8 } $$ $ 77Aq6?BAq12 22 2 a5a"fa!eQ'( (( ( 771b>q!taxQ23 33 3 dBrrr )rHrs rI test_quadratic_parameter_passingrsq Qq1QT6 Br"((QF(;H BqD As1u:. .. . Aq>fX %% %rKN)sympy.core.addrsympy.core.mulrsympy.core.numbersrrrsympy.core.relationalrsympy.core.singletonr sympy.core.symbolr sympy.matrices.denser sympy.ntheory.factor_r sympy.simplify.powsimpr sympy.core.functionrsympy.core.sortingrr(sympy.functions.elementary.trigonometricrsympy.solvers.diophantiner%sympy.solvers.diophantine.diophantinerrrrrrrrrrrrrr r!r"r#r$r%r&r'r(r)r*r+r,r-r.r/r0r1r2r3r4r5r6r7r8r9r:r;r<r=r>sympy.testing.pytestr?r@rAsympy.utilities.iterablesrBr$r%rurvr7r8rNrgrhrfrr rrrr/rrrt_3t_4t_5t_6r"r#r$rDrJrUrXr_ryrrrrrrrrrrrrrrrr&r-r0r3r9r<r@rCrJrrr{rrrrrrrrrr%r*r1r3r6r;r]r_rcrfrirlrrrOrKrIrs811$"%'+*(881            54294d2D.1aAq!Q1aAq!Q$+FD$A!S#sCc VT * B T4 +5 & 7 X2/< 0- =NAABBH#V & <UO" ;  6!EH . ;5 -V5r D<@( G50&.$*:Z R.,T^ <  ;  N64r1 N<< < 78. b&rK