timebegin=cputime;
%初始化定義全局變量
global customer;
global depot;
global demand;
global bestx;
global bestl;
global bestSC;
global best_iter;
global UD;
global c;
global cm;
global h;
global A;
global hs;
global FC;
global W;
global Q;
global Max;
global ET;
global LT;
global S;
  global La;
global Ea;
global sudu;


La=10;
Ea=5;
sudu=20;
c=20;
cm=1;
h=1;
A=15;
hs=2;
Q=150;Max=1000;


depot=4;
customer=15;

demand=[12.8  8.5 245
18.4  3.4  345
15.4  16.6 360
15.5  11.6 451
10.6  7.6 523
12.5  2.1 396
13.8  5.2 589
6.7  16.9 492
14.8  2.6 450
1.8  8.7 360
17.1  11.0 389
0.2  2.8  461
11.9  19.8 574
6.4  5.6 326
9.6  14.8 465
18.9  15.2 0
8.6  8.4  0
7.4  1.0  0
13.2  15.1 0];

ELT=[4.7  10.5
1.5  6
4.7  10.2
3.7   8.9
6.7  12.3
0.6  5.7
2.6  6.8
2.5  8.1
4.1  10.1
3.4   8.1
0  6.0
2.1  6.3
6.8  12.0
6.0  10.4
5.4  9.6
5.8  11.7
0  12
0  12
0  12
0  12];
    ELT=ELT';
    ET=ELT(1,:);
    LT=ELT(2,:);
    S=zeros(depot+customer,0);
    for j=1:customer
        S(j)=0.05;
    end
UD=[5 4 8 9 6 7 6 5 8 7 5 6 8 7 9 ]; 
  


FC=[250 430 150 240];
W=[2 3 2 4];
x0=zeros(customer,2);
l0=zeros(1,customer);
%計算需求點和配送中心之間的距離
for j=1:depot
    for k=1:customer
        Dd(k,j)=sqrt((demand(k,1)-demand(j+customer,1))^2+(demand(k,2)-demand(j+customer,2))^2);
    end
end
for i=1:customer
   x0(i,2)=i;
   l0(i)=1;
   [dis,n]=min(Dd(i,:));
   x0(i,1)=n;
end


SCx0=SC(x0,l0);
%'1'
bestx=x0;
bestl=l0;
best_iter=1;
bestSC(best_iter)ans =
 
[ 24399.590296639882581075653433800, 23712.790087948884320212528109550, 22618.995555551806319272145628929, 22297.884686486595455789938569069, 22128.230722577744018053635954857, 21719.815654735539283137768507004, 21587.065041518435464240610599518, 21368.914124218230426777154207230, 21117.032004308723116992041468620, 21018.582059241343813482671976089, 20766.699939331840141676366329193, 20419.747079433192993747070431709, 20348.460673141806182684376835823, 19772.468454482132074190303683281, 19286.050510775770817417651414871, 18794.168390866263507632538676262, 18449.798897460110310930758714676, 18377.660738007714826380833983421, 18034.118208943356876261532306671, 18006.447376724674541037529706955, 17680.285367364267585799098014832, 17664.925701112199021736159920692, 17021.758863451457727933302521706, 17021.758863451457727933302521706, 16944.393699856522289337590336800, 16927.155252213880885392427444458, 16868.241564150095655350014567375, 16847.086757511886389693245291710, 16847.086757511886389693245291710, 16843.062456037012452725321054459, 16831.754463902710995171219110489, 16797.151757167124742409214377403, 16752.492792313943937188014388084, 16749.833599823607073631137609482, 16749.833599823607073631137609482, 16745.091792038674611831083893776, 16745.091792038674611831083893776, 16715.257806492874806281179189682, 16696.194348822500614915043115616, 16696.194348822500614915043115616, 16664.656891076978354249149560928, 16664.656891076978354249149560928, 16641.236242129889433272182941437, 16641.236242129889433272182941437, 16636.152091801388451131060719490, 16627.046982419127743924036622047, 16627.046982419127743924036622047, 16619.502968642977066338062286377, 16606.851766400883207097649574280, 16602.575100419573573162779211998, 16587.822092048936610808596014977]
 =SCx0;


    one=x0(:);
  tabu=zeros(7,length(one));
  tabu=updatetabu(x0,tabu);
for count1=1:5
    [x1,l1,tabu]=location_improve(x0,l0,tabu);

    [x2,l2,tabu]=routing_improve(x1,l1,tabu);

end

timeend=cputime;
time=timeend-timebegin
%tabu
bestx
bestl
vpa(bestSC)