PRINT
PRINT "素因数分解 から 最大公約数 と 最小公倍数 を求める"
PRINT
INPUT PROMPT "自然数その１？ ":n1
INPUT PROMPT "自然数その２？ ":n2
PRINT
LET NN1=n1
LET NN2=n2
DIM A1(20),A2(20),B1(20),B2(20),C1(20),C2(20),D1(20),D2(20)
CALL SOINS(n1,A1,A2,T1)
PRINT NN1;"＝ ";
FOR i=1 TO T1
   PRINT A1(i);
   IF A2(i)>1 THEN
      PRINT "^";A2(i);
   END IF
   IF i<T1 THEN
      PRINT " * ";
   END IF
NEXT i
PRINT
CALL SOINS(n2,B1,B2,T2)
PRINT NN2;"＝ ";
FOR i=1 TO T2
   PRINT B1(i);
   IF B2(i)>1 THEN
      PRINT "^";B2(i);
   END IF
   IF i<T2 THEN
      PRINT " * ";
   END IF
NEXT i
PRINT
PRINT

LET Hi=0
FOR i=1 TO T1
   FOR j=1 TO T2
      IF A1(i)=B1(j) THEN
         LET Hi=Hi+1
         LET C1(Hi)=A1(i)
         IF A2(i)<B2(j) THEN
            LET C2(Hi)=A2(i)
         ELSE
            LET C2(Hi)=B2(j)
         END IF
      END IF
   NEXT j
NEXT i

LET Li=1
PRINT "最大公約数 ＝ ";
FOR i=1 TO Hi
   LET Li=Li*C1(i)^C2(i)
   PRINT C1(i);
   IF C2(i)>1 THEN
      PRINT "^";C2(i);
   END IF
   IF i<Hi THEN
      PRINT " * ";
   END IF
NEXT i
PRINT " ＝ ";Li

CALL SOINS(NN1*NN2/Li,D1,D2,T3)  ! 最大公約数 × 最小公倍数 ＝ ２つの数の積 より
LET Ls=1
PRINT "最小公倍数 ＝ ";
FOR i=1 TO T3
   LET Ls=Ls*D1(i)^D2(i)
   PRINT D1(i);
   IF D2(i)>1 THEN
      PRINT "^";D2(i);
   END IF
   IF i<T3 THEN
      PRINT " * ";
   END IF
NEXT i
PRINT " ＝ ";Ls
PRINT
PRINT NN1;"÷";Li;"＝";NN1/Li
PRINT NN2;"÷";Li;"＝";NN2/Li
PRINT Ls;"÷";NN1;"＝";Ls/NN1
PRINT Ls;"÷";NN2;"＝";Ls/NN2 
PRINT
END

EXTERNAL SUB SOINS(x,P(),Q(),W)
LET f=2
LET W=0
DO UNTIL x=1
   LET C=0
   DO WHILE MOD(x,f)=0
      IF C=0 THEN
         LET W=W+1
         LET P(W)=f
      END IF
      LET x=x/f
      LET C=C+1
   LOOP
   IF C=0 THEN GOTO 10
   LET Q(W)=C
10    LET f=f+1
   LOOP
END SUB