MATLAB-最小集合覆盖-贪心近似算法-以美国快递仓库选址覆盖为例

num=[];
z=zeros(353,546);
N=sum(sum(s{54}==1));
U=s{54};
Q=z;
s0=s;
t=0;
for j=1:40
max=0;
for i=1:53
if sum(sum(s0{i}&U==1))/N>max
t=i;
max=sum(sum(s0{i}&U==1))/N;
end
end
T=Q;
U=U.*(s0{t}==0);%将U在s0值为1的地方变为0 即U-s0;
Q=Q|s0{(t)};
if T==Q
break;
end
pause(0.01)%速度较快 需要pause可以看清过程
imh=image(cat(3,U,Q,s0{t}));

num(j)=t;%存储选择的集合编号
s0{t}=z;
if Q==s{54}%Q已经等于全集 故停止循环
break;
end
end