请教问题

摇摆兔

最近，我写了一个程序，不过这个程序到现在都没有完成。
那个程序是有关于使用扩展欧几里得算法求解多元线性丢番图方程，越写越伤脑筋，我想这里又没有范例可以参考。一下就是那个未完成的程序。

1. PROGRAM EX1;
   
2. USES SYSTEM;
   
3. TYPE
   
4. 
   
5.         RECORD1=RECORD
   
6.                 A:INTEGER;
   
7.                 NEXT:POINTER1;
   
8.         RECORD2=RECORD
   
9.                 CONSTANT:INTEGER;
   
10.                 MINIMUM_DIFF:INTEGER;
    
11.                 NEXT1:POINTER3;
    
12.                 NEXT2:POINTER2;
    
13.         RECORD3=RECORD
    
14.                 ITEM=SET OF INTEGER;
    
15.         POINTER1=^RECORD1;
    
16.         POINTER2=^RECORD2;
    
17.         POINTER3=^RECORD3;
    
18. END;
    
19. VAR CURRENT_VAR:POINTER1;{factor of lhs}
    
20. VAR EXISTENCE:BOOLEAN;
    
21. VAR RHS:INTEGER;
    
22. VAR SOLUTION:POINTER2;{head pointer of the linked list of the solution}
    
23. VAR P:POINTER2;
    
24. VAR MULT:INTEGER;
    
25. VAR STEP:INTEGER;
    
26. VAR TEMP:POINTER2;
    
27. PROCEDURE EUCLID(A,B: INTEGER;VAR X,Y,GCD:INTEGER);
    
28.         VAR R,Q:INTEGER;VAR C,D:INTEGER;
    
29.         VAR PTR:POINTER1;
    
30. BEGIN
    
31.         IF B=0 THEN BEGIN PTR:=CURRENT_VAR^.NEXT;CURRENT_VAR:=PTR; MVARIABLE(A) END
    
32.                 ELSE BEGIN R:=A MOD B; Q:=A DIV B;EUCLID(B,R,X,Y,GCD); END;
    
33.         IF EXISTENCE=TRUE THEN
    
34.         BEGIN
    
35.         C:=Y;D:=X-Q*Y;
    
36.         X:=C;Y:=D;
    
37.         END ELSE
    
38.         HALT();
    
39.         GCD:=A;
    
40. END;
    
41. 
    
42. PROCEDURE SUB1(A,B:INTEGER;VAR P1,P2:POINTER2);
    
43.         VAR X,Y:INTEGER; VAR D:INTEGER;
    
44. BEGIN
    
45.         X:=1;Y:=0;EUCLID(A,B,X,Y,D);
    
46.         MULT:=MULT/D;
    
47.         P2^.MINIMUM_DIFF:=A/D*MULT;P2^.CONSTANT:=Y*MULT;
    
48.         TEMP^.MINIMUM_DIFF:=-B/D*MULT;TEMP^.CONSTANT:=X*MULT;
    
49.         STEP:=STEP+2;
    
50. END;
    
51. 
    
52. PROCEDURE MVARIABLE(GCD:INTEGER);
    
53.         VAR P1,P2:POINTER2;
    
54. BEGIN
    
55.         IF CURRENT_VAR=NIL THEN BEGIN
    
56.                 EXISTENCE:=(RHS MOD GCD=0);
    
57.                 EXIT();
    
58.                 END ELSE BEGIN
    
59.                 P1:=P;TEMP:=P1;
    
60.                 P2:=GETMEM(SIZEOF (POINTER2));
    
61.                 P1^.NEXT:=P2;
    
62.                 P:=P2;
    
63.                 END;
    
64.         IF CURRENT_VAR^.A>GCD THEN SUB1(CURRENT_VAR^.A,GCD,P1,P2)
    
65.         ELSE SUB1(GCD,CURRENT_VAR^.A,P2,P1);
    
66. END;
    
67. 
    
68. BEGIN
    
69.         STEP:=1;
    
70.         P:=GETMEM(SIZEOF (POINTER2));SOLUTION:=P;
    
71.         MULT:=RHS;TEMP^.MINIMUM_DIFF:=0;TEMP^.CONSTANT:=1;
    
72. END.
    

复制代码

flyliyang

这个真心不会.看上去很厉害....希望有高手能帮帮这位兄弟!

摇摆兔

我发现我写程序越写越思维越乱，然后这个程序就没有完成。

clxer

不知道你啥问题

摇摆兔

帮我写程序

clxer

不懂解方程