十進BASIC 第2掲示板 過去ログ

LET n=1987829
LET p=3^(n-1)
LET t=p-INT(p/n)*n
PRINT -1 !trace
PRINT t !←←←←←　ここで実行時内部エラーが発生する
LET t=MOD(p,n)
PRINT -2 !trace
PRINT t
END

実行結果
-1

LET n=1987829
LET p=3^(n-1)
LET t=MOD(p,n)
PRINT -1 !trace
PRINT t !←←←←←　表示される数値が正しくない
LET t=p-INT(p/n)*n
PRINT -2 !trace
PRINT t !←←←←←　ここで実行時内部エラーが発生する
END

実行結果
-1
 836640
-2

LET t0=TIME
LET p=1987829
PRINT modpow(3,p-1,p)
PRINT "計算時間=";TIME-t0

LET t0=TIME
LET s=1
FOR k=1 TO p-1
   LET s=MOD(s*3,p)
NEXT k
PRINT s
PRINT "計算時間=";TIME-t0

END

EXTERNAL FUNCTION modpow(a,n,b) !a^n≡x mod b のxを返す　※nは非負整数
IF n<0 OR n<>INT(n) THEN !非負整数以外なら
   PRINT "modpow関数でパラメータが不適当です。"
   STOP
ELSE
   LET S=1
   DO WHILE n>0 !べき乗nを２進展開する
      IF MOD(n,2)=1 THEN LET S=MOD(S*a,b) !ビットが１なら計算する
      LET a=MOD(a*a,b)
      LET n=INT(n/2)
   LOOP
   LET modpow=S
END IF
END FUNCTION