夫妻过河问题源码.docVIP

function y=boat() m=input(输入夫妻对数: ); k=input(输入船的最大容量: ); for i=1:m+1 A=[m,i-1]; B{i}=A; %B表示可取状态 end for i=m+2:2*m A=[i-(m+1),i-(m+1)]; B{i}=A; end for i=2*m+1:3*m+1 A=[0,i-(2*m+1)]; B{i}=A; end %下面的for语句是根据参数k算出所有决策d,用C表示. t=1; for u=0:1:k for v=0:1:k if u+v=1u+v=k C{t}=[u,v]; t=t+1; end end end D{1,1}=[m,m]; % D{1,1}表示初始状态. w=1; % w控制是否能够迭代成功. h=1; % h控制元素的个数. j=2; % j表示迭代的次数. x=0; % x是用来剔除掉不在允许状态集合中的元素. y=0; % y是防止重复剔除同一个元素. l=0; % l使用来标记元素是否被剔除. z(1)=0; z(2)=1; % z(j)表示上一步迭代后所剩下的所有元素. while w0 % p表示可以进行下一步迭代的所有元素. for p=z(j-1)+1:z(j) % q表示每一个元素要经过k*(k+3)/2种决策迭代. for q=1:k*(k+3)/2 D{h+1,j}=D{p,j-1}+(-1)^(j-1).*C{q}; h=h+1; l=1; %下面的for语句是核查迭代结果是否在允许状态集合B中. for e=1:3*m+1 if D{h,j}==B{e} x=1; break; end end if x==0 h=h-1; y=1; l=0; end %下面的if语句是剔除掉本次迭代过程中出现的重复元素. if j2y==0 for i=z(j)+1:h-1 if D{h,j}==D{i,j}; h=h-1; l=0; break; end end end %下面的if语句是剔除掉以前重复的元素. if j2y==0 if j==3 for i=z(1)+1:z(2) if D{h,j}==D{i,j-2} h=h-1; l=0; break; end end