我编写了一个程序,它在Prolog列表中找到了最长的互质子序列(它不是还完美):
longest_lcs([A, B | Tail],X) :- gcd(A,B,1),lcs([B | Tail],X,A,1).
longest_lcs([A, B | Tail],X) :- lcs([B | Tail],X,A,0).
lcs([],G,_,_) :- rev(G,G1),write(G1).
lcs([A, B | Tail],G,Q,_) :- gcd(B,Q,1), gcd(A,B,1), lcs([B | Tail], [Q | G], A, 1),!.
lcs([A, B | Tail],G,Q,_) :- gcd(B,Q,1);gcd(A,B,1), lcs([B | Tail], [Q | G], A, 0).
lcs([A, B | Tail],G,_,0) :- gcd(A,B,1), lcs([B | Tail], G, A, 1).
lcs([A, B | Tail],G,_,1) :- lcs([B | Tail], G, A, 1).
lcs([A, B | Tail],G,_,0) :- lcs([B | Tail], G, A, 0).
lcs([A],G,Q,_) :- gcd(Q,A,1),lcs([], [A, Q | G], _, _).
目前我使用write
谓词输出子序列,但我需要它以下列方式运行:
?- longest_lcs([1,2,3,4],X).
X = [1,2,3,4]
?- longest_lcs([2,4,8,16],X).
X = []
我需要做出哪些修改,所以这有效吗?
答案 0 :(得分:4)
为什么要使用write / 1?专注于对你想要的明确的陈述性描述,并且顶层将为你写作。最长的互质子序列的可能公式是:它是一个互质子序列,没有其他互质子序列更长。代码看起来与此类似:
list_lcpsubseq(Ls, Subseq) :-
list_subseq(Ls, Subseq),
coprimes(Subseq),
length(Subseq, L),
\+ ( list_subseq(Ls, Others), coprimes(Others), length(Others, O), O > L ).