%function R=annealing(N,L,s,t,dt,C,R)
%N為問題規(guī)模，即節(jié)點個數(shù)；L可取較大值，如500、1000；
%s取1、2等；t為初始溫度，參考范圍為0.5--2；
%dt為衰減因子，一般不小于0.9;
%C為邊權(quán)矩陣，應(yīng)是一個強(qiáng)連通圖的邊權(quán)矩陣
%R為初始路徑，結(jié)果路徑也存放在R中
%L、s、t、dt應(yīng)通過多次試驗來確定，以獲得優(yōu)化的結(jié)果
%參考《非數(shù)值并行算法--模擬退火算法》科學(xué)出版社
function R=annealing(N,L,s,t,dt,C,R)
s0=0;
while 1
    a=0;
    for k=1:L
        [r,df]=calculate(R,C,N);
        if accept(t,df)
            R=r;a=1;
            disp(cost_sum(R,C,N));
        end
    end
    t=t*dt
    if a==0
        s0=s0+1;
    else 
        s0=0;
    end
    if s0==s
        break;
    end
end