%解搜索函數

%初始化禁忌表
for i=1:ant_n
    for j=1:city_n
        tobu(i,j)=0;
    end
end
%隨機放置螞蟻初始位置
for i=1:ant_n
    city=round((city_n-1)*rand)+1;  %隨機產生初始位置城市
    tobu(i,1)=city;
end

for i=1:ant_n  %對每一只螞蟻進行解構造
    n=1;
        
    for j=1:city_n-1  %每只螞蟻進行city_n-1次探索
        
        for k=1:city_n      %臨時池供選擇用
            temp_pool(k)=(hu_table(tobu(i,n),k)^b)*(ph_table(tobu(i,n),k)^a);
        end
        for k=1:n           %將已經訪問過的城市的選擇概率置0
            temp_pool(tobu(i,k))=0;
        end

        ada_sum=0;
        for k=1:city_n
            ada_sum=ada_sum+temp_pool(k);
        end
        r=rand*ada_sum;  %隨機產生一個數
        ada_temp=0;      %初始化累加值為0
        h=0;
        while(ada_temp<r)
            h=h+1;
            ada_temp=ada_temp+temp_pool(h);        
        end
        %退出循環時的h值即為被選擇的個體序號
        n=n+1;
        tobu(i,n)=h;    %當前訪問的城市加入禁忌表
    end
end