生成约束执行和不执行的随机值

时间:2016-04-23 12:04:31

标签: prolog clpfd clpb

我有以下内容:

:-use_module(library(clpfd)).

list_index_value(List,Index,Value):-
  nth0(Index,List,Value).

length_conindexes_conrandomvector(Length,Conindexs,Randomvector):-
  length(Randomvector,Length),
  same_length(Conindexs,Ones),
  maplist(=(1),Ones),
  maplist(list_index_value(Randomvector),Conindexs,Ones),
  term_variables(Randomvector,Vars),
  maplist(random_between(0,1),Vars).

length_conindexes_notconrandomvector(Length,Conindexes,Randomvector):-
  length(Randomvector,Length),
  length(Conindexes,NumberOfCons),
  same_length(Conindexes,Values),
  sum(Values,#\=,NumberOfCons),
  maplist(list_index_value(Randomvector),Conindexes,Values),
  term_variables(Randomvector,Vars),
  repeat,
  maplist(random_between(0,1),Vars).

length_conindexes_conrandomvector/3用于生成1和0的随机向量,其中conindexes位置中的元素为1。

 ?-length_conindexes_conrandomvector(4,[0,1],R).
 R = [1, 1, 0, 1].

length_conindexes_notconrandomvector/3用于生成随机向量,其中并非所有的conindex都是1。

?- length_conindexes_notconrandomvector(3,[0,1,2],R).
R = [1, 0, 1] ;
R = [0, 1, 1] ;
R = [1, 1, 0] 

我觉得我已经用repeat命令“入侵”了。做这个的最好方式是什么?如果我使用标签,那么值不会是随机的?如果经常违反约束,则搜索效率非常低。这样做的最佳方式是什么?

1 个答案:

答案 0 :(得分:3)

在SWI-Prolog中,我会用 CLP(B)约束来完成所有这些。

例如 1

:- use_module(library(clpb)).

length_conindices_notconrandomvector(L, Cs, Rs):-
        L #> 0,
        LMax #= L - 1,
        numlist(0, LMax, Is),
        pairs_keys_values(Pairs, Is, _),
        list_to_assoc(Pairs, A),
        maplist(assoc_index_value(A), Cs, Vs),
        sat(~ *(Vs)),
        assoc_to_values(A, Rs).

assoc_index_value(A, I, V) :- get_assoc(I, A, V).

请注意,我还可以自由地使用 O (N 2 )方法将所需元素提取到 O ( N×log N)一。

示例查询:

?- length_conindices_notconrandomvector(4, [0,1], Rs).
Rs = [X1, X2, X3, X4],
sat(1#X1*X2).

始终建议将建模部分分离为自己的谓词,我们将其称为核心关系。要获得具体的解决方案,您可以使用random_labeling/2

?- length_conindices_notconrandomvector(4, [0,1], Rs),
   length(_, Seed),
   random_labeling(Seed, Rs).
Rs = [0, 1, 1, 1],
Seed = 0 ;
Rs = [1, 0, 0, 1],
Seed = 1 ;
Rs = [1, 0, 1, 1],
Seed = 2 ;
Rs = [1, 0, 0, 1],
Seed = 3 .

CLP(B)的random_labeling/2以这样的方式实施,即每个解决方案同样可能

<小时/> 1 我当然假设您:- use_module(library(clpfd)).已经~/.swiplrc