list - 序言中最长的序列

我必须从Prolog中的列表中删除最长的质数序列。
我是Prolog的新手，我找不到找到最长序列的方法...
到目前为止，这是我所做的:

divisible(X,Y):-
  0 is X mod Y.
divisible(X,Y):-
  X > Y + 1,
  divisible(X,Y+1).

is_prime(2).
is_prime(3).
is_prime(P):-
  integer(P),
  P>3,
  P mod 2 =\= 0,
  not(divisible(P,3)).

这将从列表中删除素数。

removeP([],[]).
removeP([H],[H]):-
  not(is_prime(H)).
removeP([H|T],[H|L]):-
  not(is_prime(H)),
  removeP(T,L).
removeP([H|T],L):-
  is_prime(H),
  removeP(T,L).

在这里，我试图找到最长的序列，但是我不知道下一步该怎么做。

longest([],[]).
longest([H],[H]):-
  is_prime(H).
longest([H],[]):-
  not(in_prime(H)).
longest([H|T],L):-
  ....

最佳答案

库，特别是aggregate，有助于利用不确定性:

remove_longest(Pred, L, R) :-
    aggregate(max(C,Xc/Yc), P^(append([Xc,P,Yc],L), maplist(Pred,P), length(P,C)), max(C,X/Y)),
    append(X, Y, R).

谓词(对于您的情况为is_prime)保持通用。在此示例运行中，我仅使用原子“a”的标识:

?- remove_longest(=(a), [1,2,3,a,a,4,5,a], R).
R = [1, 2, 3, 4, 5, a].