%function [e,c,f]=GCPanneal2(L,s,t,dt,lamda,d,b,n)
%
%圖著色問題(Graph Colouring Problem)的退火算法
%GCP問題可看為將頂點集劃分為最少個數獨立集的問題
%
%求解此問題有兩種算法，
%GCPanneal1適用于度數小于20的情形
%GCPanneal2適用于各種度數
%在GCPanneal1中，w(i)表示賦予顏色i的權值
%
%
%n為問題規模，即節點個數；b為關聯矩陣
%lamda是一個大于1的罰函數因子
%d為圖G的最大度數，最小著色上界為d+1
%e(u)表示u被著的顏色號%
%c(i)表示著以顏色i的頂點個數%
%
%L可取較大值，如500、1000；
%s取1、2等；t為初始溫度，參考范圍為0.5--2；
%dt為衰減因子，一般不小于0.9;
%L、s、t、dt應通過多次試驗來確定，以獲得優化的結果
%參考《非數值并行算法--模擬退火算法》科學出版社

function [e,c]=GCPanneal2(L,s,t,dt,lamda,d,b,n)

e=zeros(1,n);
e=e+1;
c=zeros(1,d+1);
c(1)=n;
f=0;
%設定初始解

s0=0;
while 1
    a=0;
    for k=1:L
        
        u=1+fix(unifrnd(0,n));
        i=e(u);
        while 1
            j=1+fix(unifrnd(0,d+1));
            if j~=1&(c(i)>1|c(j)>0)
                break;
            end
        end
        df=0;
        for v=1:n
            if e(v)==j
                df=df+b(v,u);
            end
            if e(v)==i
                df=df-b(v,u);
            end
        end
        df=lamda*df;
        if c(i)==1
            df=df-1;
        end
        if c(j)==0
            df=df+1;
        end
        if df==0
            df=1;
        end
        
        
        if df<0|(df/t<88&exp(-df/t)>rand(1,1))
            e(u)=j;c(i)=c(i)-1;c(j)=c(j)+1;f=f+df;
            a=1;
        end
    end
    fprintf('圖中各點被著的顏色號\n');
    disp(e);
    fprintf('著以各顏色的頂點個數\n');
    disp(c);
    t=t*dt
    if a==0
        s0=s0+1;
    else 
        s0=0;
    end
    if s0==s
        break;
    end
end