Question

鉴于参数makelist(A,N,K,M,L) :-应当扩展，并且如果A,N,K,M,L>0带有M>=2和L的整数编号应返回true：

[A^N mod M , A^(N+1) mod M , A^(N+2) mod M , ... , A*(N+K) mod M]

预期结果：

| ?- makelist(10,2,1,10000,L).
L = [100,1000]
| ?- makelist(2,6,10,100,L).
L = [64,28,56,12,24,48,96,92,84,68,36]
| ?- makelist(12345678,3,8,100,L).
L = [52,56,68,4,12,36,8,24,72]
| ?- makelist(2,3000,5,7,L).
L = [1,2,4,1,2,4]
| ?- makelist(2,555000,5,17,L).
L = [1,2,4,8,16,15]
| ?- makelist(2,3000000,5,21,L).
L = [1,2,4,8,16,11]
| ?- makelist(142857,98765432,9,100,L).
L = [1,57,49,93,1,57,49,93,1,57]

我认为从数字K到(N+K)的数字开始实现指数0的{{1}}元素

Answer 1

您可以从以下代码开始：

makelist(A,N,0,M,[X]) :- X is mod(A^N,M).
makelist(A,N,K,M,[X|L]) :-  
        Max is N + K,
        N <  Max,
        N1 is  N + 1,
        K1 is K-1,
        X is mod(A^N, M),
        makelist(A, N1, K1, M, L).

然后应用更多约束，例如M>= 2。

Answer 2

小指数的解决方案

众所周知的模幂的属性可以确保：

x ^ n mod m =（ x mod m ）^ n mod m

因此，对于指数较小的大碱基，可以使用以下代码：

makelist(A, N, K, M, L) :-
    A>0, N>0, K>0, M>0,
    K1 is N + K,
    findall(X, (between(N,K1,E), X is (A mod M)^I mod M), L).

以下是使用此代码找到的一些解决方案：

?- makelist(10,2,1,10000,L).
L = [100, 1000].

?- makelist(2,6,10,100,L).
L = [64, 28, 56, 12, 24, 48, 96, 92, 84|...].

?- makelist(12345678,3,8,100,L).
L = [52, 56, 68, 4, 12, 36, 8, 24, 72].

?- makelist(2,3000,5,7,L).
L = [1, 2, 4, 1, 2, 4].

?- makelist(2,555000,5,17,L).
L = [1, 2, 4, 8, 16, 15].

?- makelist(2,3000000,5,21,L).
L = [1, 2, 4, 8, 16, 11].

请注意，如果不使用此类属性，我们将无法计算出具有大指数的大型碱基的答案：

?- X is 142857^98765432 mod 100.
ERROR: Stack limit (0.5Gb) exceeded

?- X is (142857 mod 100)^98765432 mod 100.
X = 1.

大指数的解决方案

但是，对于大指数而言，此代码实际上效率很低。例如：

?- time(makelist(142857,98765432,9,100,L)).
% 47 inferences, 115.328 CPU in -544348879559065600.000 seconds (?% CPU, 0 Lips)
L = [1, 57, 49, 93, 1, 57, 49, 93, 1|...].

因此，对于大指数 n 的最佳方法是使用以下定义：

如果 n = 0，则 x ^ n = 1
如果 n = 1，则 x ^ n = x
如果 n 是偶数，则 x ^ n =（ x ^ 2）^（< em> n / 2）
如果 n 是奇数，则 x ^ n = x *（ x ^ 2）^ floor（ n / 2）

有了这个定义，我们可以仅使用O（lg n ）乘法来计算 x ^ n 。在此基础上，我们可以定义以下谓词：

fast_pow(X, N, M, P) :-
    (   N < 32
    ->  P is (X mod M)^N mod M
    ;   (   X1 is (X mod M)^2 mod M,
            N1 is N//2,
            fast_pow(X1, N1, M, P1),
            (   N mod 2 =:= 0
            ->  P is P1
            ;   P is (X mod M)*P1  mod M ) ) ).

fast_makelist(A, N, K, M, L) :-
    A>0, N>0, K>0, M>0,
    K1 is N + K,
    findall(P, (between(N,K1,E), fast_pow(A,E,M,P)), L).

现在，我们有：

?- time(fast_makelist(142857,98765432,9,100,L)).
% 1,704 inferences, 0.000 CPU in 0.000 seconds (?% CPU, Infinite Lips)
L = [1, 57, 49, 93, 1, 57, 49, 93, 1|...].

Prolog-制作清单

2 个答案:

小指数的解决方案

大指数的解决方案