%function [e,f]=ISPanneal(L,s,t,dt,lamda,w,N)
%
%獨(dú)立集問題(Independent Set Problem)的退火算法
%設(shè)有圖G=(V,E),要找V的最大獨(dú)立集，即找最大的V'屬于V,
%滿足其中的任意兩點(diǎn)都不關(guān)聯(lián)
%
%N為問題規(guī)模，即節(jié)點(diǎn)個(gè)數(shù)；w為權(quán)矩陣
%e為以0、1為元素的數(shù)組，來表示兩個(gè)子集
%lamda是一個(gè)大于1的罰函數(shù)因子
%
%L可取較大值，如500、1000；
%s取1、2等；t為初始溫度，參考范圍為0.5--2；
%dt為衰減因子，一般不小于0.9;
%L、s、t、dt應(yīng)通過多次試驗(yàn)來確定，以獲得優(yōu)化的結(jié)果
%參考《非數(shù)值并行算法--模擬退火算法》科學(xué)出版社

function [e,f]=ISPanneal(L,s,t,dt,lamda,w,N)
f=0;
e=zeros(1,N);

s0=0;
while 1
    a=0;
    for k=1:L
        [u,df]=ISPgen(e,N,w,lamda);
        if ISPacc(df,t)
              e(u)=1-e(u);
              f=f+df;
              a=1;
        end
    end
    disp(e);disp(f);
    t=t*dt
    if a==0
        s0=s0+1;
    else 
        s0=0;
    end
    if s0==s
        break;
    end
end