1

我是 Prolog 的新手，很难弄清楚我的代码有什么问题。我正在尝试查找循环素数 到目前为止，我有以下代码：

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 中的位数。
  • rotate 使 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.
4

1 回答 1

3

我认为问题出在你的numDigits/2谓词上。

它提供了各种各样的答案，所以程序一次又一次地回溯：

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

在程序中添加剪辑后numDigits/2似乎行为正确：

isCircularPrime(19).
false.
于 2013-04-26T07:13:09.977 回答