n=100;

D=randi(30,100);   %n¸ö1¤èËíê3én¸öèÎÎñËùDèòaμÄê±¼ä￡¨¶¼Ôú30minòÔÄú￡©ò»ò»¶Ôó|1ØÏμ¾ØÕó

S1=randperm(100);  %Ëæ»úéú3éò»¸ö3õê¼½a

T0=1800;           %3õê¼Î¶è

Tend=1e-3;         %ÖÕÖ1ζè

L=250;             %ÿ¸öζèÏÂμü′ú′Îêy

q=0.99;            %½μÎÂËù¶è

Time = 0;

m = size(D,1);

for i = 1 : m

    Time = Time + D(i, S1(i));

end                               %¼ÆËãS1¶Ôó|μÄ1¤×÷×üê±¼ä￡¬¼Ç¼μ½bestTime,×÷Îaμ±Ç°×îóÅoˉêyÖμ

bestarr=S1;                       %½«S1×÷Îaμ±Ç°×îóŽa

while T0>Tend

    temp=zeros(L,n+1);

    for k=1:L

        S2=NewAnswer(S1);                 %¸ù¾Yμ±Ç°½a￡¬2úéúD½a

        [S1,R]=Metropolis(S1,S2,D,T0);    %¸ù¾Ymatropolis·¨Ôò￡¬½«D½a￡¨S1»òS2￡©¸3¸øS1￡¬×÷ÎaÏÂò»′Îμü′úμÄμ±Ç°×îóŽa

        temp(k,:)=[S1,R];                 %¼Ç¼S1￡¬¼°Æä¶Ôó|μÄ×ü1¤×÷ê±¼ä

    end

    [time0,index]=min(temp(:,end));       %¼Ç¼¸ÃζèÏÂμü′úL′ÎoóμÃμ½μÄ×î¶ìê±¼ä

    if Time>time0

        Time=time0;                       %¸üDÂ×îóÅoˉêyÖμòÔ¼°¶Ôó|μĽa

        bestarr=temp(index,1:n);

    end

    T0=q*T0;                              %¿aê¼ÏÂò»¸öζèμÄμü′ú

end

disp('×î¶ìê±¼ä￡o')

disp(Time)

disp(bestarr)