function [result,x,y]=TS_BS1(num_1,num_2,c_i,d_j,M_i,T)
%初始参数设置
runcount = 0;
tabulength = 3 + num_1;
candidate.x = round(rand(1,num_1));
y = zeros(num_1,num_2);
sum_d = 0;
for i = 1:num_2
sum_d = sum_d + d_j(i);
end
alfa = 0.9;
%初始化禁忌表
for i = 1:tabulength
tabulist(i).change = [0,0];
tabulist(i).x = 0;
tabulist(i).value = 0;
end
candidate.value = c_i*candidate.x';%目前的解
best_so_far.x = candidate.x;
best_so_far.value = candidate.value;
while runcount < T
runcount = runcount+1;
%产生领域解
kk = 1;
for i = 1:num_1-2
for j = i:num_1-1
flag = 1;
for n = 1:tabulength
tempflag = sum(abs(tabulist(n).change-[