解二元一次不定方程

lufeng369

使用了辗转相处法..

黑人类

用列主元高斯-约当消去法
可以解n元一次方程。

lufeng369

10 REM ax+by=c
20 DIM R(50),C(50)
21 PRINT "Ax+Bx=C"
30 INPUT "A=";A
31 INPUT "B=";B
32 INPUT "C=";C
40 R(0)=A:R(1)=B:C(0)=C:I=1
50 WHILE R(I):I=I+1
60 Q=INT(R(I-2)/R(I-1))
70 R(I)=R(I-2)-Q*R(I-1)
80 P=INT(C(I-2)/R(I-1))
90 C(I-1)=C(I-2)-P*R(I-1)
120 WEND
130 GCD=R(I-1)
131 REM PRINT GCD;",";R(I)","C(I-2):END
140 FOR J=I TO 2 STEP-1
150 M=C(J-2)-R(J-2)*Y
160 Y=M/R(J-1)
170 NEXT J
180 X=(C-B*Y)/A
190 TEXT:
200 PRINT A;"x+";B;"y=";C
210 PRINT "X=";X;"-";B/GCD;"t"
220 PRINT "Y=";Y;"+";A/GCD;"t"
230 PRINT "(t∈Ｚ)"