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

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