%function [e,f]=MCPanneal(L,s,t,dt,w,N)
%最大截問題(Max Cut Problem)的退火算法
%由于算法的缺陷，當f值在較長時間內不再改變時，可用cl+c停止程序
%N為問題規模，即節點個數；w為權矩陣
%L可取較大值，如500、1000；
%s取1、2等；t為初始溫度，參考范圍為0.5--2；
%dt為衰減因子，一般不小于0.9;
%e為以0、1為元素的數組，來表示兩個子集
%f為最大截權重和
%L、s、t、dt應通過多次試驗來確定，以獲得優化的結果
%參考《非數值并行算法--模擬退火算法》科學出版社

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

s0=0;
while 1
    a=0;
    for k=1:L
        [u,df]=MCPgen(e,N,w);
        if MCPacc(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