数独求解器无限递归

时间:2011-02-11 20:26:39

标签: recursion prolog backtracking

我正在写一个数独求解器。自从我触及prolog以来已经很长时间了,因此我不记得有关统一,回溯等的一切。我认为我会让系统永远回溯(但我没有得到任何堆栈异常 - 至少没有立即)。这是我到目前为止(这个谜题可以在http://en.wikipedia.org/wiki/File:Sudoku-by-L2G-20050714.svg找到):

% representation of the example puzzle
puzzle([5, 3, _, _, 7, _, _, _, _],
       [6, _, _, 1, 9, 5, _, _, _],
       [_, 9, 8, _, _, _, _, 6, _],
       [8, _, _, _, 6, _, _, _, 3],
       [4, _, _, 8, _, 3, _, _, 1],
       [7, _, _, _, 2, _, _, _, 6],
       [_, 6, _, _, _, _, 2, 8, _],
       [_, _, _, 4, 1, 9, _, _, 5],
       [_, _, _, _, 8, _, _, 7, 9]).

% solve(S)
% the starting point of the program
% saves the solution in the variable S
solve(R1, R2, C1) :-
    % save the rows into variables
    puzzle(R1, R2, R3, R4, R5, R6, R7, R8, R9),
    % solve for each row
    allunique(R1), allunique(R2), allunique(R3),
    allunique(R4), allunique(R5), allunique(R6),
    allunique(R7), allunique(R8), allunique(R9),
    % the columns must be created first
    nelement(R1, 1, C11), nelement(R2, 1, C21), nelement(R3, 1, C31),
    nelement(R4, 1, C41), nelement(R5, 1, C51), nelement(R6, 1, C61),
    nelement(R7, 1, C71), nelement(R8, 1, C81), nelement(R9, 1, C91),
    C1 = [C11, C21, C31, C41, C51, C61, C71, C81, C91],
    allunique(C1).

% allunique(List)
% Succeeds if all the numbers of List are between 1-9
% and each number exists only once
allunique([]). % Recursion stops when the list is empty

% A member should be between 1-9 and not a member of the tail
allunique([H|T]) :-
    allunique(T),
    member(H, [1, 2, 3, 4, 5, 6, 7, 8, 9]),
    not(member(H, T)).

% nelement(List, N-th, X)
% Saves the nth element of a list in variable X
nelement([H|_], 1, H). % The first element is the head

% All other elements will be found in the tail
nelement([_|T], N, X) :-
    N > 1,
    N1 is N-1,
    nelement(T, N1, X).

allunique(C1)会导致问题。似乎系统在第一列的第一个空盒子中放置了7并且从不改变它(或者至少在不久的将来不会改变它)。示例C1列表是[5, 6, 7, 8, 4, 7, 9, 8, 6]。我很想知道为什么会这样。

1 个答案:

答案 0 :(得分:4)

  1. 您的示例列表[5, 6, 7, 8, 4, 7, 9, 8, 6]不满足allunique,因为它包含8两次。
  2. solve/3不正确,因为它检查了所有行,但只检查了一列而没有“区域”(3x3正方形)。
  3. 评论中承诺的solve/1谓词没有出现,所以我无法测试你的代码; allunique/1nelement/3似乎很好。
  4. 即使你完成了这个程序,我也怀疑它是否会回复。我已经看到类似的Prolog程序运行了几个小时而没有找到解决方案。如果您想快速执行此操作,请查看CLP(fd)(链接适用于SWI,但SICStus,GNU和ECLiPSe具有类似的库。)