检查Prolog中的循环素数

时间:2013-04-26 06:44:46

标签: prolog primes

我是Prolog的新手,我很难搞清楚我的代码有什么问题。 我想找circular prime numbers 到目前为止,我有这段代码:

isCircularPrime(N):-
    prime(N),
    numDigits(N, Y),
    Y2 is Y-1,
    rotate(N,Y,N2),
    circularPrime(N2, Y2, Y2).

circularPrime(_, 0, _).     
circularPrime(N, 1, _):-prime(N).
circularPrime(N, I, J):-
    I > 1,
    prime(N),
    I2 is I-1,
    rotate(N,J,N2),
    circularPrime(N2,I2,J).

其中:

  • numDigits使Y为N中的位数。
  • 旋转使N2旋转N(例如12→21)。
  • 如果N是素数,则
  • 素数为真,否则为

目前,只要数字为素数但不是圆形素数(例如19,23等),它就会无限地递归。如果有人能解释我做错了什么,我真的很感激。

以下是其他代码,以防有人想要自己测试一下:

prime(2).
prime(3).
prime(N):-
    N>3,
    N mod 2 =\= 0,
    \+ divisible(N,3).

divisible(N,I):- N mod I =:= 0.
divisible(N,I):-
    I*I < N,
    I2 is I + 2,
    divisible(N, I2).

numDigits(N, 1):- 
    N<10, N>0.
numDigits(N, X):- 
    N2 is N/10, 
    numDigits(N2, Y), 
    X is (Y+1).

rotate(N, L, R):-
    LastDigit is N mod 10,
    Base is div(N, 10),
    Exponent is L - 1,
    Num is 10**Exponent,
    NewBase is LastDigit*Num,
    R is Base + NewBase.

1 个答案:

答案 0 :(得分:3)

我认为问题出在您的numDigits/2谓词中。

它提供了各种答案,因此程序会一次又一次地回溯:

numDigits(1, X).
X = 1 ;
X = 2 ;
X = 3 ;
X = 4 ;
X = 5 etc...

numDigits/2添加剪切后,程序似乎表现正常:

isCircularPrime(19).
false.