> restart: > read("BM_res.mpl");read("rev.mpl"): > shft:=331: > f:=2*(x-shft)^11+3*(x-shft)^5 + (x-shft)^7; 11 5 7 f := 2 (x - 331) + 3 (x - 331) + (x - 331) > expand(f); 2 -10450520637186142189617162426 - 5246197415416405798015211 x 3 4 + 47548617057818627599595 x - 287302822100353848410 x 5 6 + 1215178099515045228 x - 3671232929043841 x 7 8 9 10 + 7922384395861 x - 11967348030 x + 12051710 x - 7282 x 11 + 2 x + 347298268904248529738432804 x > for i from 0 to 20 do a[i]:=expand(subs(x=2^(i+1)-s,f)); od: > > for i from 0 to 20 do b[i]:=expand(subs(x=3^(i+1)-s,f)); od: > p1:=rev(bermass_sres(a,9),z); p1 := -163954452963290742261906657024118060460253282651827981845\ 4959136085167952640931008496735150750988885893643480935056\ 923171624719902310400 s - 13871599560109134082170388611776\ 2097993340783185526311588064907993904217204013903019783277\ 22 772800 s - 340886458222298078692870303889908447637476287\ 21 32370100835026419559382838484111614303676384935936000 s - 7678749592263420294841610336336052504732461495068975866656\ 20 616528714792082021746498379656172929024000 s - 158450817\ 2387036172329874920273143722969730923317095167028766833608\ 19 603921734737408816934603930468352000 s - 299220858852763\ 0936280377745556244138305916996654956042823617295282038100\ 18 13809932693390568516693996339200 s - 5164007630844330784\ 0119896509593324162752622809762216907758387679419527096165\ 17 937314876111869827180239257600 s - 813022805047768731229\ 7734433916035262176765330072682890017211383201298240494139\ 16 241121781400609477091550822400 s - 116509396688167870545\ 4483547125655179248200176348192736584331040641995059903932\ 15 263913559040831737723118996684800 s - 151549195267028447\ 8356377296353467103387718014220455127119857172095499740932\ 14 56774558657488987376038130167788339200 s - 1783226043762\ 3988747234374927372547980249707553084476376931000081252229\ 13 551153126108932096362025353862181238813491200 s - 189029\ 2631526262098726469491099492758841388841172896333081811322\ 12 433653720048675226401250507682189279084329253312921600 s - 1593992143738096521793338181756781856413857248214034539305\ 26 3108730894191430624031539200 s - 41873001268479346783302\ 7530848330899751909053435217477921455374536407711664308224\ 27 00 s - 1004829841744628486705195239515294017241496157027\ 28 35593038198579471333478472089600 s - 2198576061857872109\ 0104776519716622802931952288505965812160751751706496204800\ 29 0 s - 43758138366513359856933779926213728988082972478782\ 30 1003206988162909025075200 s - 78993744538254751193990644\ 31 7328661863360687658196707355812067959583539200 s - 12889\ 1684782525011818258642787737399595672588336190363183486913\ 32 7408000 s - 51756666294459506389390318553909654629842939\ 23 7701109441001275252430422528100252229028754227200 s - 17\ 7055165128077418519290955357228915821753352285088054069233\ 24 3911014612713882756592749772800 s - 55512030221993726159\ 1645159383657766641386972530734416820328212570209930855810\ 25 2508339200 s - 18928749622434291744027160932538243289858\ 33 50997297460080201957376000 s - 37 2723768469628819347433606949875836523310874196180992000 s - 3011392758488629631059891376197158049000130419535904768000 36 s - 2913258285266010408041689881518985661033856003680656949248000 35 s - 2489252046134895551267673187261167102922186573971068\ 34 597305344000 s 38 - 2130476318769723545528750932516739313285734596608000 s 39 - 1419440205325737406997615099672268320362463232000 s 40 - 789573748338525431305444194057001673687040000 s 41 - 356708241418297501782393699734724280320000 s 42 - 125713755569301531490503543383654400000 s 43 - 32418738871460239825432561582080000 s 45 - 445764330122896421683200000 s 44 - 5440197037685852509506109440000 s - 13338245843384678\ 2514587619894855301570392718965799571262898263879432594550\ 1227687493802639591730238872947774987735022659428574822400\ 2 00 s - 70684714083349652069760602696195325037041712494412\ 4318821770160283742978952450024735231783924856899795351798\ 3 7317415030511050214604800 s - 274480898968404761495381857\ 0545295301583976780760637903311609046301435983316400627057\ 4 91899983266388501321099721608308396325758566400 s - 83277\ 5448722857624903416577004515252097452044777917126203971255\ 3025485437619006552184843370227776991199905612552836606101\ 5 500723200 s - 2055309387655277277192756047290888883486289\ 6382516381295994450970308821086381985393344140583757225830\ 6 6548800596755164035232563200 s - 424162559666448566944747\ 2308630364063529137244435609220168019226100203170768542980\ 7 771991610879721903527747820433172272146022400 s - 7467269\ 9076629244702833716449895566488923788451210541887045068360\ 5724466830237884494805339161860009436224688669410138783744\ 8 00 s - 11384144700404423866357916129071032744221497693625\ 6051635957156740622032412413926712921742562564699858559699\ 9 1991064192614400 s - 152062155918607749183323217944688345\ 6568709540180568659327508665721747053763062050141394801418\ 10 6779418159053924652495667200 s - 17961768624014521515669\ 3768606845879769864360243814110555699152888002597720054057\ 11 295157940150111143341435824702174003200 s - 984238922509\ 3143306959609717371120766503238377415422741572737316909543\ 9528861569672596409346041240878782518468663177031043222327\ 00108800 > p2:=rev(bermass_sres(b,9),z); p2 := -291925559246763117071458504936953917391845889115272698614\ 1747647998691506771015996707717405404156357462675230766433\ 278334436067214255455513047994657271548352869171200 s - 10\ 1662931038391610262106168411111731404018428467780062899718\ 5520313739106669906532938836737283173403379089422134804480\ 22 0 s - 15717310460723931176649882019949211229875168031638\ 9301863993560667041642606870865199109226107807237092728921\ 21 46876573286400 s - 1486555525235481776177651160166815061\ 0926535876823642377431527298192877567315768491573286399565\ 20 795083189191271145321948774400 s - 105044345915268537759\ 1311407000973881962965126427035549139265757258457915953101\ 19 7531770574033902072856150299296889057514514022400 s - 59\ 4518947807700360244644792922601986897965872841911102440074\ 3546579077843256560494807085711924610977831824520569298845\ 18 745073356800 s - 278861346815678524687920800461189345242\ 1883845062316910361667682358186812249077851137320452591703\ 17 249166280614860393158046132469760000 s - 110538303476307\ 3816298389782710022607242895341220968473927284204515105810\ 9950694952021710138204203676451288657678668435143845652463\ 16 61600 s - 3747264698679180483904204125516568524867807302\ 6450924619199654163282318575595706472769747628998580593849\ 15 8571861955678130569544655595110400 s - 10943827047221257\ 0615893588523710409031511253641208946581433664875460905680\ 4322602620153604827853211771521516067868036085377878972992\ 14 32358400 s - 2764906043817759300895119249641041623330553\ 9213018137620486547679421230896355767830323527251704765594\ 13 749007965207259766090390844956047992422400 s - 605336013\ 2232639886396222904804743808253208018668213826335476874903\ 9868386025873819129532647095520359017450916116056615743522\ 12 04097511336443904000 s - 1140871702601500925049028091804\ 9482407378277141169381898558708384869278327721613453119640\ 26 55614085555814400 s - 5637239755599919402584246675307246\ 7346755010284222966757612485422334427961948994525888125326\ 27 64447795200 s - 1359993740833421507340575822581725423164\ 1189317988714293854521202408020185958292738554984835055616\ 28 00 s + 4014315760576980471824864520887003455364602477255\ 29 256148346955145100025652104055477050611335168000 s + 274\ 6636345438021451055971875532368634693319968645922284420487\ 30 230327308642240635888175625011200 s - 732901728187597754\ 7113212504192333506785315307693629998446580296067559117176\ 31 89953668300800 s + 4130872004760816834587883239446941502\ 32 94612246520041512474463675841453139456621019136000 s + 2\ 6532267108355911465773156636171033577323813448984123742213\ 23 6297796094554222219091921682325726696682161743645900800 s + 8317839722140215611942018921998912411383859804902780197\ 24 745906573109272922613806521049569650778020513324413747200 s + 7050436707931857194533777933364064291558252320106223442\ 25 873728589593785313113033859695664275861046655188992000 s + 3117916924019922742727678093813508240517135476111446465921\ 33 463044937861456017227776000 s + 490407803256968221646133\ 37 381799360407166692845982122563557516065636352000 s + 703\ 4852121313226814213505548054020235780133656722253176482281\ 36 88661219328000 s + 8504056234457388420656865121602092496\ 35 05262909455234830216966506026366926848000 s + 2484792379\ 4421508532840702245893930998576412942033846645278981466363\ 34 30386325504000 s - 7624685089132852951116690217541954876\ 38 0786267370645030492765683712000 s - 78991881228542512501\ 39 868170050595852809047254354038068575469568000 s + 323983\ 40 85548450335149930604558438809264156734083999895388160000 s - 5931281195038303854301442201580852835341997270735257600000 41 s 42 + 650393094524769015632683083642948308069971100958720000 s 43 - 45691200583206224392840723214870096940934103040000 s 45 - 42598482835860787430465892726394060800000 s 44 + 1984970024399172284054905296550329329909760000 s - 14\ 6224780471635281211546123867059894197218678818037516616674\ 1982695406760842746840775991043613856290107867702463318862\ 2 13747192072922164007316771813015834023323238400 s - 47059\ 9042030402143570538356122476505656015203588709984385149918\ 6159396960443314349012669721091998818565823392935468376661\ 3 956359211555255202436990163223105136230400 s - 1093362832\ 8138747736975827903483734591639273714397018231812556394740\ 2453085072085890333489939378533496681784649950962696905168\ 4 614672177682497009642877278696243200 s - 1953059747418546\ 8260392852883354273912118516616813725597114769666868032070\ 5171276308557873522493922474170134039937970747482937700537\ 5 1008667258466558619838054400 s - 278901312136186221832473\ 0219883083802922655733582681310414629201361997755823076577\ 2459549796643210655268914913227851698928519542952665708278\ 6 167738938661273600 s - 3268239724532323574666308008949193\ 8134129337271264465959047542381897556095497827415292477826\ 8944647032925077183636396111663096919354186478173985475395\ 7 584000 s - 3200778668134580426897214407873599720655780181\ 3270043786449730590497835201182049509479229052830913040267\ 8 06545842043222475454144873846958280404970871193600 s - 26\ 5447963287385492923441250933089753826582215404592345390678\ 8971613504070229365219977867903577772651613059964033411519\ 9 0861756774315038217240498405376000 s - 188186589028071356\ 2636299461949871776311570721503428801083566216039660137180\ 5668231365285170393524905229889986251650920172864016492393\ 10 0583951317401600 s - 11480647195156237107037209739794091\ 8916488993898619511853083684301518536964903276222593838918\ 11 2700476856185029167584353524797854569280543007126323200 s - 2810998641922004238295990680676134002061173265422976860\ 4628870622154060437321507262840249497562740421005908395090\ 472819177023603052592478615613764832141767748878336000 > gcd(p1,p2); 2 176120502681600 s + 116591772775219200 s + 19295938394298777600 > solve(%=0); -331, -331 > >