restart:with(numtheory):Hit:=0:prepro:=0:cnt:=0:An48:=0:An26:=0:Digits:=700:K:=9079: ##################################################################################### A3L:=23 * 6959 * 69267032772761 * 24905421493145578416249151437563 * 135460119964929578352833498048615612058741648790897452004065796599887904767736939851478088233541311243872291330177006542445431480080950696661534013139274950734051561305811761089:A3V:=A3L*27000+K: #A3L:=180776592240767009075325425036351865190999552955783362327500337504233:A3V:=A3L*27000+K: appendto(`C:\\Users/usr/Desktop/yaf1/Nmr232_Allb.txt`): ################################## pp:=3000:start:=180776592240767009075325425036259428093500440758476179260165132103480:Mxx:=start*100^60: for k1k2 from start to Mxx by 30 do if frac(k1k2/30)=0 then K1K2:=15/2*k1k2:czc:=ifactor(k1k2,easy): ############################################################################## S_T_F1F2:=proc() local C1,radF,m,m1,m2,ABC,pv,C,FM1,FM2,Mm1,Mm2,val11,val12,val21,val22,G,FA,FB,F1v,F2v,V: global cnt,k1,k2,K1,K2,K,F1,F2,N1,M1,xb,Hit,Loc,S_loc,XB,A3L,A3V,_Z1,BC,Grp,BCFL,BCFV,BCFLxa,BCFVxa,BCxa,xa,M: Loc:=(BC-K)/27000: for XB from -pp to pp do S_loc:=solve(Loc-A3L): if whattype(S_loc[1])=integer or whattype(S_loc[1])=fraction or whattype(S_loc[2])=integer or whattype(S_loc[2])=fraction then Hit:=Hit+1:G:=(K*(K1[w]/15)+K2[w]/15)^2+4*K1[w]*K2[w]/225*(A3V-K): if nops(factor(BC))=2 or nops(factor(BC))=3 then FM1:=op(1,factor(BC)):FM2:=op(2,factor(BC)):Mm1:=S_loc[1]:Mm2:=S_loc[2]: val11:=coeff(FM1,M,1)*Mm1+coeff(FM1,M,0):val12:=coeff(FM2,M,1)*Mm1+coeff(FM2,M,0): val21:=coeff(FM1,M,1)*Mm2+coeff(FM1,M,0):val22:=coeff(FM2,M,1)*Mm2+coeff(FM2,M,0): fi: print():lprint(`F1_F2`,Grp,`xb=`,xb,BCFL,BCFV,`xb=`,XB,factor((BC-K)/27000),factor(BC),`m1=`,S_loc[1],`m2=`,S_loc[2]): print():lprint(`n1=`,K1[w]/15,`n2=`,K2[w]/15,`C=`,ifactor(G,easy),`N1=`,N1,`M1=`,M1,`F11=`,val11,`F12=`,val12,`F21=`,val21,`F22=`,val22):M:=S_loc[1]: print():lprint(`para1`,`N1=`,N1,`M1=`,M1):M:=S_loc[2]:print():lprint(`para2`,`N1=`,N1,`M1=`,M1):M:='M': print():lprint(`cnt=`,cnt,`k1=`,k1[v],`k2=`,k2[v],`K1=`,K1[w],`K2=`,K2[w],`x=`,x[w]): print(): ################################################################################################################################################ fi:od:XB:='XB':Loc:='Loc':_Z1:='_Z1': end: ######################################################################################################################################################################## S_T_FF1FF2:=proc() local C2,radFF,m,m1,m2,ABC,pv,C,FM1,FM2,Mm1,Mm2,val11,val12,val21,val22,G,FA,FB,F1v,F2v,V: global cnt,k1,k2,K1,K2,K,FF1,FF2,N1,M1,xb,Hit,Loc,S_loc,XB,A3L,A3V,_Z1,BC,Grp,BCFFL,BCFFV,BCFFLxa,BCFFVxa,BCxa,xa,M: Loc:=(BC-K)/27000: for XB from -pp to pp do S_loc:=solve(Loc-A3L): if whattype(S_loc[1])=integer or whattype(S_loc[1])=fraction or whattype(S_loc[2])=integer or whattype(S_loc[2])=fraction then Hit:=Hit+1:G:=(K*(K1[w]/15)+K2[w]/15)^2+4*K1[w]*K2[w]/225*(A3V-K): if nops(factor(BC))=2 or nops(factor(BC))=3 then FM1:=op(1,factor(BC)):FM2:=op(2,factor(BC)):Mm1:=S_loc[1]:Mm2:=S_loc[2]: val11:=coeff(FM1,M,1)*Mm1+coeff(FM1,M,0):val12:=coeff(FM2,M,1)*Mm1+coeff(FM2,M,0): val21:=coeff(FM1,M,1)*Mm2+coeff(FM1,M,0):val22:=coeff(FM2,M,1)*Mm2+coeff(FM2,M,0): fi: print():lprint(`FF1_FF2`,Grp,`xb=`,xb,BCFFL,BCFFV,`xb=`,XB,factor((BC-K)/27000),factor(BC),`m1=`,S_loc[1],`m2=`,S_loc[2]): print():lprint(`n1=`,K1[w]/15,`n2=`,K2[w]/15,`C=`,ifactor(G,easy),`N1=`,N1,`M1=`,M1,`FF11=`,val11,`FF12=`,val12,`FF21=`,val21,`FF22=`,val22):M:=S_loc[1]: print():lprint(`para1`,`N1=`,N1,`M1=`,M1):M:=S_loc[2]:print():lprint(`para2`,`N1=`,N1,`M1=`,M1):M:='M': print():lprint(`cnt=`,cnt,`k1=`,k1[v],`k2=`,k2[v],`K1=`,K1[w],`K2=`,K2[w],`x=`,x[w]): print(): ########################################################################################################### fi:od:XB:='XB':Loc:='Loc':_Z1:='_Z1': end: ##############################################################################lprint(C):lprint(r,K1[r],K2[r],x[r]):print():print():lprint(p,k1[p],k2[p]):print():lprint(B):print(): # pre processing K1K2_eq:=divisors(K1K2):C:=sort(convert(K1K2_eq,'list')):s:=nops(C)+1:if s<2300 then for r from 1 to nops(C) do s:=s-1:K1[r]:=C[r]:K2[r]:=C[s]:x[r]:=K1[r]*K+K2[r]:od: for w from 1 to nops(C) while prepro<>1 do if whattype((K1[w]*K+K2[w])/450)=integer then lft1:=abs(450*F1*k1k2-450*F2*1):lft2:=(450*F2*1-450*F1*k1k2):Pt1:=isolve(lft1-x[w]):Pt2:=isolve(lft2-x[w]): if whattype(Pt1)=set or whattype(Pt2)=set then prepro:=1:fi: fi:end do: if prepro=1 then k1k2_eq:=divisors(k1k2):B:=sort(convert(k1k2_eq,'list')):q:=nops(B)+1:if q<2300 then for p from 1 to nops(B) do q:=q-1:k1[p]:=B[p]:k2[p]:=B[q]:od: lprint(`k1k2=`,k1k2,czc,`K1K2=`,K1K2,`termsk1k2=`,nops(B),`termsK1K2=`,nops(C),`Total=`,nops(B)*nops(C)):print(): ###################################################################################################### for v from 1 to nops(B) do for w from 1 to nops(C) do if whattype((K1[w]*K+K2[w])/450)=integer then lft1:=(450*F1*k2[v]-450*F2*k1[v]):lft2:=(450*F2*k1[v]-450*F1*k2[v]): PQ1:=isolve(lft1-x[w]):PQ2:=isolve(lft2-x[w]):if whattype(PQ1)=set and whattype(PQ2)=set then print():lprint(`cnt=`,cnt,`k1=`,k1[v],`k2=`,k2[v],`K1=`,K1[w],`K2=`,K2[w],`x=`,x[w]):print(): lprint(`PQ1=`,PQ1,`PQ2=`,PQ2):print(): lprint(`1xx`,lhs(PQ1[1]),rhs(PQ1[1]),lhs(PQ1[2]),rhs(PQ1[2])):print():lprint(`2xx`,lhs(PQ2[1]),rhs(PQ2[1]),lhs(PQ2[2]),rhs(PQ2[2])):print(): if lhs(PQ1[1])=F1 then F1:=rhs(PQ1[1]):F2:=rhs(PQ1[2]):fi:if lhs(PQ1[1])=F2 then F2:=rhs(PQ1[1]):F1:=rhs(PQ1[2]):fi:F1F2:=factor(F1*F2*27000+K): if lhs(PQ2[1])=F1 then FF1:=rhs(PQ2[1]):FF2:=rhs(PQ2[2]):fi:if lhs(PQ2[1])=F2 then FF2:=rhs(PQ2[1]):FF1:=rhs(PQ2[2]):fi:FF1FF2:=factor(FF1*FF2*27000+K): if nops(F1F2)=2 or nops(FF1FF2)=2 then PP:=F1F2:QQ:=FF1FF2: ############################################################################## a1:=PP/op(2,PP):a2:=PP/op(1,PP):CHP1a:=op(2,a1):CHP2a:=op(2,a2): CHP1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHP2b:=coeff(op(1,a2),_Z1)+op(2,a2):lprint(`pp=`,PP,CHP1a,CHP2a,CHP1b,CHP2b):print(): a1:=QQ/op(2,QQ):a2:=QQ/op(1,QQ):CHQ1a:=op(2,a1):CHQ2a:=op(2,a2): CHQ1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHQ2b:=coeff(op(1,a2),_Z1)+op(2,a2):lprint(`qq=`,QQ,CHQ1a,CHQ2a,CHQ1b,CHQ2b):print(): ############################################################################## ######## 0119_499V258+10_260 FFU if CHP1a<0 and CHP2a<0 then a1:=-PP/op(2,PP):a2:=-PP/op(1,PP):CHP1a:=op(2,a1):CHP2a:=op(2,a2): CHP1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHP2b:=coeff(op(1,a2),_Z1)+op(2,a2):fi: if CHP1a>0 and CHP2a>0 then if frac((CHP1a-19)/30)=0 and frac((CHP2a-1)/30)=0 then tmp1a:=CHP1a:tmp2a:=CHP2a:CHP1a:=tmp2a:CHP2a:=tmp1a: tmp1b:=CHP1b:tmp2b:=CHP2b:CHP1b:=tmp2b:CHP2b:=tmp1b:fi:fi: if frac((CHP1a-1)/30)=0 then Grp:=`0119`:xb:=frac((CHP1a-1)/30/30)*30:n1a:=trunc((CHP1a-1)/30/30): if frac((xb-1)/3)=0 then m1a:=solve((32479+27000*L1-900*(4-xb)^2-570*xb+9900*n1a)-CHP2a) else m1a:=solve((23479+27000*L1-900*(4-xb)^2-570*xb+9900*n1a)-CHP2a):fi: n1b:=trunc((CHP1b-1)/30/30): if frac((xb-1)/3)=0 then m1b:=solve((32479+27000*L1-900*(4-xb)^2-570*xb+9900*n1b)-CHP2b) else m1b:=solve((23479+27000*L1-900*(4-xb)^2-570*xb+9900*n1b)-CHP2b):fi: N1:=n1a+(n1b-n1a)*M:M1:=m1a+(m1b-m1a)*M:lprint(`0119_F1F2`,`xb=`,xb,`n1a=`,n1a,`n1b=`,n1b,`m1a=`,m1a,`m1b=`,m1b,`N1=`,N1,`M1=`,M1):print(): if frac((xb-1)/3)=0 then BC:=(900*N1+30*XB+1)*(32479+27000*M1-900*(4-XB)^2-570*XB+9900*N1):BCxa:=(900*N1+30*xa+1)*(32479+27000*M1-900*(4-xa)^2-570*xa+9900*N1) else BC:=(900*N1+30*XB+1)*(23479+27000*M1-900*(4-XB)^2-570*XB+9900*N1):BCxa:=(900*N1+30*xa+1)*(23479+27000*M1-900*(4-xa)^2-570*xa+9900*N1):fi: XB:=xb:BCFL:=factor((BC-K)/27000):BCFV:=factor(BC):BCFLxa:=factor((BCxa-K)/27000):BCFVxa:=factor(BCxa):XB:='XB': S_T_F1F2();fi: ######## 0119_499V258+10_260 FFU if CHQ1a<0 and CHQ2a<0 then a1:=-QQ/op(2,QQ):a2:=-QQ/op(1,QQ):CHQ1a:=op(2,a1):CHQ2a:=op(2,a2): CHQ1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHQ2b:=coeff(op(1,a2),_Z1)+op(2,a2):fi: if CHQ1a>0 and CHQ2a>0 then if frac((CHQ1a-19)/30)=0 and frac((CHQ2a-1)/30)=0 then tmp1a:=CHQ1a:tmp2a:=CHQ2a:CHQ1a:=tmp2a:CHQ2a:=tmp1a: tmp1b:=CHQ1b:tmp2b:=CHQ2b:CHQ1b:=tmp2b:CHQ2b:=tmp1b:fi:fi: if frac((CHQ1a-1)/30)=0 then Grp:=`0119`:xb:=frac((CHQ1a-1)/30/30)*30:n1a:=trunc((CHQ1a-1)/30/30): if frac((xb-1)/3)=0 then m1a:=solve((32479+27000*L1-900*(4-xb)^2-570*xb+9900*n1a)-CHQ2a) else m1a:=solve((23479+27000*L1-900*(4-xb)^2-570*xb+9900*n1a)-CHQ2a):fi: n1b:=trunc((CHQ1b-1)/30/30): if frac((xb-1)/3)=0 then m1b:=solve((32479+27000*L1-900*(4-xb)^2-570*xb+9900*n1b)-CHQ2b) else m1b:=solve((23479+27000*L1-900*(4-xb)^2-570*xb+9900*n1b)-CHQ2b):fi: N1:=n1a+(n1b-n1a)*M:M1:=m1a+(m1b-m1a)*M:lprint(`0119_FF1FF2`,`xb=`,xb,`n1a=`,n1a,`n1b=`,n1b,`m1a=`,m1a,`m1b=`,m1b,`N1=`,N1,`M1=`,M1):print(): if frac((xb-1)/3)=0 then BC:=(900*N1+30*XB+1)*(32479+27000*M1-900*(4-XB)^2-570*XB+9900*N1):BCxa:=(900*N1+30*xa+1)*(32479+27000*M1-900*(4-xa)^2-570*xa+9900*N1) else BC:=(900*N1+30*XB+1)*(23479+27000*M1-900*(4-XB)^2-570*XB+9900*N1):BCxa:=(900*N1+30*xa+1)*(23479+27000*M1-900*(4-xa)^2-570*xa+9900*N1):fi: XB:=xb:BCFFL:=factor((BC-K)/27000):BCFFV:=factor(BC):BCFFLxa:=factor((BCxa-K)/27000):BCFFVxa:=factor(BCxa):XB:='XB': S_T_FF1FF2();fi: ######## 0119_499V258+10_260 FFU ######## 0707_457V_260 PFU if CHP1a<0 and CHP2a<0 then a1:=-PP/op(2,PP):a2:=-PP/op(1,PP):CHP1a:=op(2,a1):CHP2a:=op(2,a2): CHP1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHP2b:=coeff(op(1,a2),_Z1)+op(2,a2):fi: if CHP1a>0 and CHP2a>0 then if frac((CHP1a-7)/30)=0 and frac((CHP2a-7)/30)=0 then tmp1a:=CHP1a:tmp2a:=CHP2a:CHP1a:=tmp2a:CHP2a:=tmp1a: tmp1b:=CHP1b:tmp2b:=CHP2b:CHP1b:=tmp2b:CHP2b:=tmp1b:fi:fi: if frac((CHP1a-7)/30)=0 then Grp:=`0707`:xb:=frac((CHP1a-7)/30/30)*30:n1a:=trunc((CHP1a-7)/30/30): m1a:=solve((62497+27000*L1-1800*(22-xb)^2-1830*xb-900*n1a)-CHP2a):n1b:=trunc((CHP1b-7)/30/30): m1b:=solve((62497+27000*L1-1800*(22-xb)^2-1830*xb-900*n1b)-CHP2b):N1:=n1a+(n1b-n1a)*M:M1:=m1a+(m1b-m1a)*M:lprint(`0707_F1F2`,`xb=`,xb,`n1a=`,n1a,`n1b=`,n1b,`m1a=`,m1a,`m1b=`,m1b,`N1=`,N1,`M1=`,M1):print(): BC:=(900*N1+30*XB+7)*(62497+27000*M1-1800*(22-XB)^2-1830*XB-900*N1):BCxa:=(900*N1+30*xa+7)*(62497+27000*M1-1800*(22-xa)^2-1830*xa-900*N1): XB:=xb:BCFL:=factor((BC-K)/27000):BCFV:=factor(BC):BCFLxa:=factor((BCxa-K)/27000):BCFVxa:=factor(BCxa):XB:='XB': S_T_F1F2();fi: ######## 0707_457V_260 PFU if CHQ1a<0 and CHQ2a<0 then a1:=-QQ/op(2,QQ):a2:=-QQ/op(1,QQ):CHQ1a:=op(2,a1):CHQ2a:=op(2,a2): CHQ1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHQ2b:=coeff(op(1,a2),_Z1)+op(2,a2):fi: if CHQ1a>0 and CHQ2a>0 then if frac((CHQ1a-7)/30)=0 and frac((CHQ2a-7)/30)=0 then tmp1a:=CHQ1a:tmp2a:=CHQ2a:CHQ1a:=tmp2a:CHQ2a:=tmp1a: tmp1b:=CHQ1b:tmp2b:=CHQ2b:CHQ1b:=tmp2b:CHQ2b:=tmp1b:fi:fi: if frac((CHQ1a-7)/30)=0 then Grp:=`0707`:xb:=frac((CHQ1a-7)/30/30)*30:n1a:=trunc((CHQ1a-7)/30/30): m1a:=solve((62497+27000*L1-1800*(22-xb)^2-1830*xb-900*n1a)-CHQ2a):n1b:=trunc((CHQ1b-7)/30/30): m1b:=solve((62497+27000*L1-1800*(22-xb)^2-1830*xb-900*n1b)-CHQ2b):N1:=n1a+(n1b-n1a)*M:M1:=m1a+(m1b-m1a)*M:lprint(`0707_FF1FF2`,`xb=`,xb,`n1a=`,n1a,`n1b=`,n1b,`m1a=`,m1a,`m1b=`,m1b,`N1=`,N1,`M1=`,M1):print(): BC:=(900*N1+30*XB+7)*(62497+27000*M1-1800*(22-XB)^2-1830*XB-900*N1):BCxa:=(900*N1+30*xa+7)*(62497+27000*M1-1800*(22-xa)^2-1830*xa-900*N1): XB:=xb:BCFFL:=factor((BC-K)/27000):BCFFV:=factor(BC):BCFFLxa:=factor((BCxa-K)/27000):BCFFVxa:=factor(BCxa):XB:='XB': S_T_FF1FF2();fi: ######## 0707_457V_260 PFU ######## 1129_209V_260 FFU if CHP1a<0 and CHP2a<0 then a1:=-PP/op(2,PP):a2:=-PP/op(1,PP):CHP1a:=op(2,a1):CHP2a:=op(2,a2): CHP1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHP2b:=coeff(op(1,a2),_Z1)+op(2,a2):fi: if CHP1a>0 and CHP2a>0 then if frac((CHP1a-29)/30)=0 and frac((CHP2a-11)/30)=0 then tmp1a:=CHP1a:tmp2a:=CHP2a:CHP1a:=tmp2a:CHP2a:=tmp1a: tmp1b:=CHP1b:tmp2b:=CHP2b:CHP1b:=tmp2b:CHP2b:=tmp1b:fi:fi: if frac((CHP1a-11)/30)=0 then Grp:=`1129`:xb:=frac((CHP1a-11)/30/30)*30:n1a:=trunc((CHP1a-11)/30/30): m1a:=solve((25289+27000*L1-900*(7-xb)^2-570*xb+9900*n1a)-CHP2a):n1b:=trunc((CHP1b-11)/30/30): m1b:=solve((25289+27000*L1-900*(7-xb)^2-570*xb+9900*n1b)-CHP2b):N1:=n1a+(n1b-n1a)*M:M1:=m1a+(m1b-m1a)*M:lprint(`1129_F1F2`,`xb=`,xb,`n1a=`,n1a,`n1b=`,n1b,`m1a=`,m1a,`m1b=`,m1b,`N1=`,N1,`M1=`,M1):print(): BC:=(900*N1+30*XB+11)*(25289+27000*M1-900*(7-XB)^2-570*XB+9900*N1):BCxa:=(900*N1+30*xa+11)*(25289+27000*M1-900*(7-xa)^2-570*xa+9900*N1): XB:=xb:BCFL:=factor((BC-K)/27000):BCFV:=factor(BC):BCFLxa:=factor((BCxa-K)/27000):BCFVxa:=factor(BCxa):XB:='XB': S_T_F1F2();fi: ######## 1129_209V_260 FFU if CHQ1a<0 and CHQ2a<0 then a1:=-QQ/op(2,QQ):a2:=-QQ/op(1,QQ):CHQ1a:=op(2,a1):CHQ2a:=op(2,a2): CHQ1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHQ2b:=coeff(op(1,a2),_Z1)+op(2,a2):fi: if CHQ1a>0 and CHQ2a>0 then if frac((CHQ1a-29)/30)=0 and frac((CHQ2a-11)/30)=0 then tmp1a:=CHQ1a:tmp2a:=CHQ2a:CHQ1a:=tmp2a:CHQ2a:=tmp1a: tmp1b:=CHQ1b:tmp2b:=CHQ2b:CHQ1b:=tmp2b:CHQ2b:=tmp1b:fi:fi: if frac((CHQ1a-11)/30)=0 then Grp:=`1129`:xb:=frac((CHQ1a-11)/30/30)*30:n1a:=trunc((CHQ1a-11)/30/30): m1a:=solve((25289+27000*L1-900*(7-xb)^2-570*xb+9900*n1a)-CHQ2a):n1b:=trunc((CHQ1b-11)/30/30): m1b:=solve((25289+27000*L1-900*(7-xb)^2-570*xb+9900*n1b)-CHQ2b):N1:=n1a+(n1b-n1a)*M:M1:=m1a+(m1b-m1a)*M:lprint(`1129_FF1FF2`,`xb=`,xb,`n1a=`,n1a,`n1b=`,n1b,`m1a=`,m1a,`m1b=`,m1b,`N1=`,N1,`M1=`,M1):print(): BC:=(900*N1+30*XB+11)*(25289+27000*M1-900*(7-XB)^2-570*XB+9900*N1):BCxa:=(900*N1+30*xa+11)*(25289+27000*M1-900*(7-xa)^2-570*xa+9900*N1): XB:=xb:BCFFL:=factor((BC-K)/27000):BCFFV:=factor(BC):BCFFLxa:=factor((BCxa-K)/27000):BCFFVxa:=factor(BCxa):XB:='XB': S_T_FF1FF2();fi: ######## 1129_209V_260 FFU ######## 1313_523V246+7_260 PU if CHP1a<0 and CHP2a<0 then a1:=-PP/op(2,PP):a2:=-PP/op(1,PP):CHP1a:=op(2,a1):CHP2a:=op(2,a2): CHP1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHP2b:=coeff(op(1,a2),_Z1)+op(2,a2):fi: if CHP1a>0 and CHP2a>0 then if frac((CHP1a-13)/30)=0 and frac((CHP2a-13)/30)=0 then tmp1a:=CHP1a:tmp2a:=CHP2a:CHP1a:=tmp2a:CHP2a:=tmp1a: tmp1b:=CHP1b:tmp2b:=CHP2b:CHP1b:=tmp2b:CHP2b:=tmp1b:fi:fi: if frac((CHP1a-13)/30)=0 then Grp:=`1313`:xb:=frac((CHP1a-13)/30/30)*30:n1a:=trunc((CHP1a-13)/30/30): if frac(xb/2)=0 then m1a:=solve((21883-450*(xb-12)^2+27000*L1-900*n1a-30*xb)-CHP2a) else m1a:=solve((4783-450*(xb-11)^2+870*xb+27000*L1-900*n1a)-CHP2a):fi: n1b:=trunc((CHP1b-13)/30/30): if frac(xb/2)=0 then m1b:=solve((21883-450*(xb-12)^2+27000*L1-900*n1b-30*xb)-CHP2b) else m1b:=solve((4783-450*(xb-11)^2+870*xb+27000*L1-900*n1b)-CHP2b):fi: N1:=n1a+(n1b-n1a)*M:M1:=m1a+(m1b-m1a)*M:lprint(`1313_F1F2`,`xb=`,xb,`n1a=`,n1a,`n1b=`,n1b,`m1a=`,m1a,`m1b=`,m1b,`N1=`,N1,`M1=`,M1):print(): if frac(xb/2)=0 then BC:=(900*N1+30*XB+13)*(21883-450*(XB-12)^2+27000*M1-900*N1-30*XB):BCxa:=(900*N1+30*xa+13)*(21883-450*(xa-12)^2+27000*M1-900*N1-30*xa) else BC:=(900*N1+30*XB+13)*(4783-450*(XB-11)^2+870*XB+27000*M1-900*N1):BCxa:=(900*N1+30*xa+13)*(4783-450*(xa-11)^2+870*xa+27000*M1-900*N1):fi: XB:=xb:BCFL:=factor((BC-K)/27000):BCFV:=factor(BC):BCFLxa:=factor((BCxa-K)/27000):BCFVxa:=factor(BCxa):XB:='XB': S_T_F1F2();fi: ######## 1313_523V246+7_260 PU if CHQ1a<0 and CHQ2a<0 then a1:=-QQ/op(2,QQ):a2:=-QQ/op(1,QQ):CHQ1a:=op(2,a1):CHQ2a:=op(2,a2): CHQ1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHQ2b:=coeff(op(1,a2),_Z1)+op(2,a2):fi: if CHQ1a>0 and CHQ2a>0 then if frac((CHQ1a-13)/30)=0 and frac((CHQ2a-13)/30)=0 then tmp1a:=CHQ1a:tmp2a:=CHQ2a:CHQ1a:=tmp2a:CHQ2a:=tmp1a: tmp1b:=CHQ1b:tmp2b:=CHQ2b:CHQ1b:=tmp2b:CHQ2b:=tmp1b:fi:fi: if frac((CHQ1a-13)/30)=0 then Grp:=`1313`:xb:=frac((CHQ1a-13)/30/30)*30:n1a:=trunc((CHQ1a-13)/30/30): if frac(xb/2)=0 then m1a:=solve((21883-450*(xb-12)^2+27000*L1-900*n1a-30*xb)-CHQ2a) else m1a:=solve((4783-450*(xb-11)^2+870*xb+27000*L1-900*n1a)-CHQ2a):fi: n1b:=trunc((CHQ1b-13)/30/30): if frac(xb/2)=0 then m1b:=solve((21883-450*(xb-12)^2+27000*L1-900*n1b-30*xb)-CHQ2b) else m1b:=solve((4783-450*(xb-11)^2+870*xb+27000*L1-900*n1b)-CHQ2b):fi: N1:=n1a+(n1b-n1a)*M:M1:=m1a+(m1b-m1a)*M:lprint(`1313_FF1FF2`,`xb=`,xb,`n1a=`,n1a,`n1b=`,n1b,`m1a=`,m1a,`m1b=`,m1b,`N1=`,N1,`M1=`,M1):print(): if frac(xb/2)=0 then BC:=(900*N1+30*XB+13)*(21883-450*(XB-12)^2+27000*M1-900*N1-30*XB):BCxa:=(900*N1+30*xa+13)*(21883-450*(xa-12)^2+27000*M1-900*N1-30*xa) else BC:=(900*N1+30*XB+13)*(4783-450*(XB-11)^2+870*XB+27000*M1-900*N1):BCxa:=(900*N1+30*xa+13)*(4783-450*(xa-11)^2+870*xa+27000*M1-900*N1):fi: XB:=xb:BCFFL:=factor((BC-K)/27000):BCFFV:=factor(BC):BCFFLxa:=factor((BCxa-K)/27000):BCFFVxa:=factor(BCxa):XB:='XB': S_T_FF1FF2();fi: ######## 1313_523V246+7_260 PU ######## 1717_347V135+7_260 PD if CHP1a<0 and CHP2a<0 then a1:=-PP/op(2,PP):a2:=-PP/op(1,PP):CHP1a:=op(2,a1):CHP2a:=op(2,a2): CHP1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHP2b:=coeff(op(1,a2),_Z1)+op(2,a2):fi: if CHP1a>0 and CHP2a>0 then if frac((CHP1a-17)/30)=0 and frac((CHP2a-17)/30)=0 then tmp1a:=CHP1a:tmp2a:=CHP2a:CHP1a:=tmp2a:CHP2a:=tmp1a: tmp1b:=CHP1b:tmp2b:=CHP2b:CHP1b:=tmp2b:CHP2b:=tmp1b:fi:fi: if frac((CHP1a-17)/30)=0 then Grp:=`1717`:xb:=frac((CHP1a-17)/30/30)*30:n1a:=trunc((CHP1a-17)/30/30): if frac(xb/2)=0 then m1a:=solve((6437+450*(xb-1)^2-930*xb+27000*L1-900*n1a)-CHP2a) else m1a:=solve((13187+450*(xb-1)^2-930*xb+27000*L1-900*n1a)-CHP2a):fi: n1b:=trunc((CHP1b-17)/30/30): if frac(xb/2)=0 then m1b:=solve((6437+450*(xb-1)^2-930*xb+27000*L1-900*n1b)-CHP2b) else m1b:=solve((13187+450*(xb-1)^2-930*xb+27000*L1-900*n1b)-CHP2b):fi: N1:=n1a+(n1b-n1a)*M:M1:=m1a+(m1b-m1a)*M:lprint(`1717_F1F2`,`xb=`,xb,`n1a=`,n1a,`n1b=`,n1b,`m1a=`,m1a,`m1b=`,m1b,`N1=`,N1,`M1=`,M1):print(): if frac(xb/2)=0 then BC:=(900*N1+30*XB+17)*(6437+450*(XB-1)^2-930*XB+27000*M1-900*N1):BCxa:=(900*N1+30*xa+17)*(6437+450*(xa-1)^2-930*xa+27000*M1-900*N1) else BC:=(900*N1+30*XB+17)*(13187+450*(XB-1)^2-930*XB+27000*M1-900*N1):BCxa:=(900*N1+30*xa+17)*(13187+450*(xa-1)^2-930*xa+27000*M1-900*N1):fi: XB:=xb:BCFL:=factor((BC-K)/27000):BCFV:=factor(BC):BCFLxa:=factor((BCxa-K)/27000):BCFVxa:=factor(BCxa):XB:='XB': S_T_F1F2();fi: ######## 1717_347V135+7_260 PD if CHQ1a<0 and CHQ2a<0 then a1:=-QQ/op(2,QQ):a2:=-QQ/op(1,QQ):CHQ1a:=op(2,a1):CHQ2a:=op(2,a2): CHQ1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHQ2b:=coeff(op(1,a2),_Z1)+op(2,a2):fi: if CHQ1a>0 and CHQ2a>0 then if frac((CHQ1a-17)/30)=0 and frac((CHQ2a-17)/30)=0 then tmp1a:=CHQ1a:tmp2a:=CHQ2a:CHQ1a:=tmp2a:CHQ2a:=tmp1a: tmp1b:=CHQ1b:tmp2b:=CHQ2b:CHQ1b:=tmp2b:CHQ2b:=tmp1b:fi:fi: if frac((CHQ1a-17)/30)=0 then Grp:=`1717`:xb:=frac((CHQ1a-17)/30/30)*30:n1a:=trunc((CHQ1a-17)/30/30): if frac(xb/2)=0 then m1a:=solve((6437+450*(xb-1)^2-930*xb+27000*L1-900*n1a)-CHQ2a) else m1a:=solve((13187+450*(xb-1)^2-930*xb+27000*L1-900*n1a)-CHQ2a):fi: n1b:=trunc((CHQ1b-17)/30/30): if frac(xb/2)=0 then m1b:=solve((6437+450*(xb-1)^2-930*xb+27000*L1-900*n1b)-CHQ2b) else m1b:=solve((13187+450*(xb-1)^2-930*xb+27000*L1-900*n1b)-CHQ2b):fi: N1:=n1a+(n1b-n1a)*M:M1:=m1a+(m1b-m1a)*M:lprint(`1717_FF1FF2`,`xb=`,xb,`n1a=`,n1a,`n1b=`,n1b,`m1a=`,m1a,`m1b=`,m1b,`N1=`,N1,`M1=`,M1):print(): if frac(xb/2)=0 then BC:=(900*N1+30*XB+17)*(6437+450*(XB-1)^2-930*XB+27000*M1-900*N1):BCxa:=(900*N1+30*xa+17)*(6437+450*(xa-1)^2-930*xa+27000*M1-900*N1) else BC:=(900*N1+30*XB+17)*(13187+450*(XB-1)^2-930*XB+27000*M1-900*N1):BCxa:=(900*N1+30*xa+17)*(13187+450*(xa-1)^2-930*xa+27000*M1-900*N1):fi: XB:=xb:BCFFL:=factor((BC-K)/27000):BCFFV:=factor(BC):BCFFLxa:=factor((BCxa-K)/27000):BCFFVxa:=factor(BCxa):XB:='XB': S_T_FF1FF2();fi: ######## 1717_347V135+7_260 PD ######## 2323_413V_260 PFD if CHP1a<0 and CHP2a<0 then a1:=-PP/op(2,PP):a2:=-PP/op(1,PP):CHP1a:=op(2,a1):CHP2a:=op(2,a2): CHP1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHP2b:=coeff(op(1,a2),_Z1)+op(2,a2):fi: if CHP1a>0 and CHP2a>0 then if frac((CHP1a-23)/30)=0 and frac((CHP2a-23)/30)=0 then tmp1a:=CHP1a:tmp2a:=CHP2a:CHP1a:=tmp2a:CHP2a:=tmp1a: tmp1b:=CHP1b:tmp2b:=CHP2b:CHP1b:=tmp2b:CHP2b:=tmp1b:fi:fi: if frac((CHP1a-23)/30)=0 then Grp:=`2323`:xb:=frac((CHP1a-23)/30/30)*30:n1a:=trunc((CHP1a-23)/30/30): m1a:=solve((-10327+27000*L1+1800*(8-xb)^2+1770*xb-900*n1a)-CHP2a):n1b:=trunc((CHP1b-23)/30/30): m1b:=solve((-10327+27000*L1+1800*(8-xb)^2+1770*xb-900*n1b)-CHP2b):N1:=n1a+(n1b-n1a)*M:M1:=m1a+(m1b-m1a)*M:lprint(`2323_F1F2`,`xb=`,xb,`n1a=`,n1a,`n1b=`,n1b,`m1a=`,m1a,`m1b=`,m1b,`N1=`,N1,`M1=`,M1):print(): BC:=(900*N1+30*XB+23)*(-10327+27000*M1+1800*(8-XB)^2+1770*XB-900*N1):BCxa:=(900*N1+30*xa+23)*(-10327+27000*M1+1800*(8-xa)^2+1770*xa-900*N1): XB:=xb:BCFL:=factor((BC-K)/27000):BCFV:=factor(BC):BCFLxa:=factor((BCxa-K)/27000):BCFVxa:=factor(BCxa):XB:='XB': S_T_F1F2();fi: ######## 2323_413V_260 PFD if CHQ1a<0 and CHQ2a<0 then a1:=-QQ/op(2,QQ):a2:=-QQ/op(1,QQ):CHQ1a:=op(2,a1):CHQ2a:=op(2,a2): CHQ1b:=coeff(op(1,a1),_Z1)+op(2,a1):CHQ2b:=coeff(op(1,a2),_Z1)+op(2,a2):fi: if CHQ1a>0 and CHQ2a>0 then if frac((CHQ1a-23)/30)=0 and frac((CHQ2a-23)/30)=0 then tmp1a:=CHQ1a:tmp2a:=CHQ2a:CHQ1a:=tmp2a:CHQ2a:=tmp1a: tmp1b:=CHQ1b:tmp2b:=CHQ2b:CHQ1b:=tmp2b:CHQ2b:=tmp1b:fi:fi: if frac((CHQ1a-23)/30)=0 then Grp:=`2323`:xb:=frac((CHQ1a-23)/30/30)*30:n1a:=trunc((CHQ1a-23)/30/30): m1a:=solve((-10327+27000*L1+1800*(8-xb)^2+1770*xb-900*n1a)-CHQ2a):n1b:=trunc((CHQ1b-23)/30/30): m1b:=solve((-10327+27000*L1+1800*(8-xb)^2+1770*xb-900*n1b)-CHQ2b):N1:=n1a+(n1b-n1a)*M:M1:=m1a+(m1b-m1a)*M:lprint(`2323_FF1FF2`,`xb=`,xb,`n1a=`,n1a,`n1b=`,n1b,`m1a=`,m1a,`m1b=`,m1b,`N1=`,N1,`M1=`,M1):print(): BC:=(900*N1+30*XB+23)*(-10327+27000*M1+1800*(8-XB)^2+1770*XB-900*N1):BCxa:=(900*N1+30*xa+23)*(-10327+27000*M1+1800*(8-xa)^2+1770*xa-900*N1): XB:=xb:BCFFL:=factor((BC-K)/27000):BCFFV:=factor(BC):BCFFLxa:=factor((BCxa-K)/27000):BCFFVxa:=factor(BCxa):XB:='XB': S_T_FF1FF2();fi: ######## 2323_413V_260 PFD ################################################################# #lprint(``): cnt:=cnt+1: fi:fi:fi: if k2[v]=k1[v-1] then v:=nops(B)+1:w:=nops(C)+1:cnt:=2*cnt:Hit:=2*Hit:fi: _Z1:='_Z1':F1:='F1':F2:='F2':FF1:='FF1':FF2:='FF2': od:od:if Hit<>0 then lprint(`Solutions Found`,`Found=`,cnt,`Hits=`,Hit):print():fi:print(): fi:fi:fi:fi:prepro:=0:cnt:=0:Hit:=0:od: ##################################################################;