数学规划学习之配送中心选址问题

学习了数学规划，现在把做的题目进行总结——配送选址问题

某省共有92个城市，城市位置、标号，公路交通网数据见附件1。
某企业在该省标号前20位的城市建立了直销中心，各直销中心负责所在城市的销售，销售量见附件1。该企业欲在该省设立一个配送中心负责给直销中心配送产品，配送中心建设成本为30万元。每吨公里运费2元，每吨产品的销售利润为300元。
试建立数学模型分析研究下面的问题：
（1）为了降低运输成本，配送中心应选在哪个城市？
（2）请为该企业制定一个成本最小的5年产品配送计划：应设立几个配送中心、各设在何处？
（3）如果该企业考虑重新为20个直销中心选址，请给出最佳的5年产品销售、配送计划。
（4）假定没有直销中心城市的客户按就近的原则购买产品，请重新考虑问题（3）。

（1）对于第一个问题，首先考虑题目，降低运输成本，在配送中心费用相同的情况下，考虑最小的公里数，因此，把此题归为最短路径问题。

首先，对于最短路径问题，最常用的就是floyd算法和dijkstra算法，在学习的过程中，已经对这两种算法有所理解，在mydijkstra算法中需要知道邻接矩阵然后要明确自己需要从哪个值到另一个值的最短路径。
对于dijkstra有：

function [mydistance,mypath]=mydijkstra(a,sb,db);
% 输入：a—邻接矩阵，a（i，j)是指i到j之间的距离，可以是有向的
% sb—起点的标号, db—终点的标号
% 输出：mydistance—最短路的距离, mypath—最短路的路径
n=size(a,1); %n取a的行数，即顶点个数
visited(1:n) = 0;visited(sb)=1;%记录顶点是否得到标号
distance(1:n) = inf; distance(sb) = 0; %起点到各顶点距离的初始化
 u=sb;  %u为最新的P标号顶点
parent(1:n) = 0; %前驱顶点的初始化
for i = 1: n-1
     id=find(visited==0); %查找未标号的顶点
     for v = id           
         if  a(u, v) + distance(u) < distance(v)
             distance(v) = distance(u) + a(u, v);  %修改标号值 
             parent(v) = u;                                    
         end            
     end
     temp=distance;
     temp(visited==1)=inf;  %已标号点的距离换成无穷
     [t, u] = min(temp);  %找标号值最小的顶点 
     visited(u) = 1;       %标记已经标号的顶点
 end
mypath = [];
if parent(db) ~= 0   %如果存在路!
    t = db; mypath = [db];
    while t ~= sb
        p = parent(t);
        mypath = [p mypath];
        t = p;      
    end
end
mydistance = distance(db);

而对于floyd可以知道，在floyd算法中并不需要要从哪到哪的要求，只需要算出邻接矩阵即可。
对于floyd算法有：

function [dist,mypath,path]=floyd2(a);
% 输入：a—邻接矩阵，元素(aij)是顶点i到j之间的直达距离，可以是有向的

% 输出：dist—最短路的距离；% mypath—最短路的路径
n=size(a,1); path=zeros(n);
for k=1:n
    for i=1:n
        for j=1:n
            if a(i,j)>a(i,k)+a(k,j)
                a(i,j)=a(i,k)+a(k,j);
                path(i,j)=k;
            end
        end
    end
end
dist=a;

for g=1:n
   for h=1:n
       gg=g;
        while path(gg,h)~=0&&path(path(gg,h),h)~=0 
              path(g,h)=path(path(gg,h),h);
              gg=path(gg,h);      
       end
   end
end
for i=1:length(a)
    for j=1:length(a)
parent=path(i,:); %从起点i到终点j的最短路上各顶点的前驱顶点
parent(parent==0)=i; %path中的分量为0，表示该顶点的前驱是起点
mypath{i,j}=[j]; t=j;
while t~=i
        p=parent(t); mypath{i,j}=[p,mypath{i,j}];
        t=p;
end
    end
end

在对于第一个问题中，可以知道应该算出每个点与点之间的距离，但在计算距离的时候，应该要筛选出sheet2中，有点点并没有相互连接的地方，在matlab编程中，可以使用find等函数进行筛选。
然后因为是练习数学规划，在此，下面是lingo程序的第一个问题的解答。

此时的程序为：

model:
  sets:
  nodesr/n1..n92/:xu;!xu为每个城市的需求量,问题1、2中只需考虑前20个城市;
  nodesc/m1..m92/:y,cost;!y为城市是否为配送中心的0-1变量;
  links(nodesr,nodesc): dist,x;!x(i,j)为判定第j个城市作为配送中心是否为第i个直销中心配送货物的0-1变量;
endsets

  !目标函数;
  min=z;
z=@sum(nodesr(I)|I#LE#20:@sum(nodesc(J): dist(I,J)*x(I,J)*xu(I)))*2;

  !约束条件;
@for(nodesr(I)|I#LE#20:@sum(nodesc(J): x(I,J))=1);!每个直销中心的需求都得到满足;
@sum(nodesc(J):y(J))=1;    !只设置1个直销中心;
@for(nodesc(J):@for(nodesr(I)|I#LE#20:x(I,J)<=y(J)));!反映x和y的对应关系,其实这儿可以写成'=';
@for(nodesc(J):x(J,J)=y(J)); !j恰好同时为直销中心和配送中心;
@for(links:@bin(x));
@for(nodesc:@bin(y));
data:
dist=  0.0000000e+000  1.8987490e+001  3.8838842e+001  4.5352167e+001  9.3742889e+001  9.5375178e+001  1.1500349e+002  9.0226247e+001  9.2254378e+001  1.4649568e+002  1.9087929e+002  2.2236153e+002  2.2001753e+002  1.6028474e+002  1.4249318e+002  9.2868122e+001  3.5912053e+001  2.5645722e+001  1.7583464e+001  5.2631990e+001  1.9293439e+002  2.1096215e+002  2.2501753e+002  2.2893203e+002  2.1090427e+002  1.8187929e+002  1.8931232e+002  1.9001160e+002  1.9515806e+002  1.2083445e+002  1.1281154e+002  1.0360174e+002  9.8502720e+001  9.7279316e+001  8.8011738e+001  9.0823632e+001  9.5922652e+001  5.8809349e+001  5.5809349e+001  3.8131680e+001  4.4412053e+001  2.6063196e+001  1.8000938e+001  2.8474323e+001  8.1303534e+001  8.0925709e+001  1.0830311e+002  1.1850115e+002  8.8742889e+001  8.6067321e+001  8.2259434e+001  8.0655768e+001  8.2034685e+001  5.9231176e+001  5.1497831e+001  7.6413127e+001  6.4033709e+001  7.1533709e+001  7.9343959e+001  6.2744611e+001  9.7457721e+001  4.8852167e+001  3.5044403e+001  2.5989018e+001  2.3599092e+001  2.0436814e+001  1.6194173e+001  1.2071068e+001  5.0000000e+000  1.0385165e+001  1.1403124e+001  1.6403124e+001  1.0296111e+001  6.2649820e+000  9.3005376e+000  1.2836072e+001  1.6403124e+001  6.4031242e+000  1.3111328e+001  1.7583464e+001  3.2353926e+001  3.6439214e+001  3.1030887e+001  4.0879744e+001  4.8159854e+001  5.6237542e+001  5.1942002e+001  4.7910873e+001  4.3879744e+001  4.5171338e+001  4.9914755e+001  6.9939739e+001
  1.8987490e+001  0.0000000e+000  2.1116536e+001  5.6850675e+001  7.8337113e+001  9.8420771e+001  9.7281187e+001  7.2503941e+001  7.4532072e+001  1.2877337e+002  1.7315698e+002  2.0463922e+002  2.0103004e+002  1.4129725e+002  1.2477088e+002  7.3880632e+001  2.5911116e+001  4.3847598e+001  3.6570954e+001  7.0833867e+001  1.7394690e+002  1.9197466e+002  2.0603004e+002  2.1120972e+002  1.9318197e+002  1.6415698e+002  1.7159002e+002  1.7228930e+002  1.7743576e+002  1.0311214e+002  9.5089235e+001  8.5879433e+001  8.0780413e+001  7.9557010e+001  7.0289431e+001  7.1836142e+001  7.6935162e+001  3.9821859e+001  3.6821859e+001  1.9144190e+001  3.4411116e+001  1.6062258e+001  8.0000000e+000  9.4868330e+000  6.3581228e+001  6.3203403e+001  9.2897333e+001  1.0309537e+002  7.3337113e+001  8.3281962e+001  7.9474076e+001  7.5172913e+001  6.6628909e+001  4.3825401e+001  3.3775525e+001  7.9415554e+001  7.5532217e+001  8.3032217e+001  8.2389552e+001  7.4243119e+001  1.0895623e+002  6.0350675e+001  4.6542911e+001  3.7487526e+001  3.1656574e+001  2.8494297e+001  2.4251656e+001  2.1058558e+001  1.3987490e+001  8.6023253e+000  2.0390614e+001  1.6062258e+001  2.4124515e+001  2.5252472e+001  2.5586250e+001  2.9121784e+001  3.3593920e+001  2.5390614e+001  3.2098818e+001  3.6570954e+001  5.0555802e+001  5.4641090e+001  4.9232763e+001  5.9081621e+001  6.6361731e+001  7.4439418e+001  7.0143879e+001  6.6112750e+001  6.2081621e+001  6.3373215e+001  6.8116631e+001  8.0727959e+001
  3.8838842e+001  2.1116536e+001  0.0000000e+000  4.0433852e+001  5.7220577e+001  7.7304235e+001  7.6164651e+001  5.1387404e+001  5.3415536e+001  1.0765683e+002  1.5204044e+002  1.8352268e+002  1.8740507e+002  1.2767227e+002  1.0365434e+002  6.0255658e+001  4.7027652e+001  5.8949096e+001  4.1942567e+001  8.5935365e+001  1.6032193e+002  1.7834968e+002  1.9240507e+002  1.9009319e+002  1.7206543e+002  1.4304044e+002  1.5047348e+002  1.5117276e+002  1.5631922e+002  8.1995603e+001  7.3972698e+001  6.4762897e+001  5.9663877e+001  5.8440474e+001  4.9172895e+001  5.4172895e+001  5.9271915e+001  6.0938396e+001  5.7938396e+001  4.0260726e+001  5.5527652e+001  3.7178794e+001  2.9116536e+001  1.1629703e+001  4.2464691e+001  4.2086867e+001  7.1780797e+001  8.1978836e+001  5.2220577e+001  6.2165426e+001  5.8357539e+001  5.4056377e+001  4.5512373e+001  2.2708865e+001  1.2658989e+001  5.8299017e+001  5.9115393e+001  6.6615393e+001  6.1273015e+001  5.7826296e+001  9.2539406e+001  4.3933852e+001  3.0126088e+001  2.1070703e+001  1.5239751e+001  1.8402028e+001  2.2644669e+001  2.6767775e+001  3.3838842e+001  2.9718862e+001  4.0241967e+001  3.7178794e+001  4.5241052e+001  4.5103824e+001  3.1157107e+001  2.7621573e+001  3.2093709e+001  4.2093709e+001  4.6414702e+001  5.0886838e+001  6.5657300e+001  6.9742588e+001  6.4334261e+001  7.4183119e+001  8.1463229e+001  8.9540916e+001  8.5245377e+001  8.1214248e+001  7.7183119e+001  7.8474713e+001  8.3218129e+001  1.0184450e+002
  4.5352167e+001  5.6850675e+001  4.0433852e+001  0.0000000e+000  4.9200435e+001  5.0023011e+001  7.6566903e+001  8.3272828e+001  8.9866681e+001  1.4410798e+002  1.8849159e+002  2.1997383e+002  2.0981794e+002  1.5008515e+002  1.1475067e+002  8.2668530e+001  7.4705246e+001  6.3843619e+001  4.6837089e+001  6.7988800e+001  1.8273480e+002  2.0076256e+002  2.1481794e+002  2.2654433e+002  2.0851657e+002  1.7949159e+002  1.8692462e+002  1.6226909e+002  1.5535338e+002  8.1029762e+001  9.9673358e+001  8.7968658e+001  9.1549301e+001  9.4891619e+001  8.5624041e+001  8.0624041e+001  8.5723060e+001  4.8609758e+001  4.5609758e+001  6.3287427e+001  8.3205246e+001  6.4856389e+001  5.6794131e+001  4.7363843e+001  7.9972291e+001  7.3972291e+001  6.3760655e+001  7.3958694e+001  5.0555808e+001  4.0715154e+001  3.6907267e+001  3.5303601e+001  4.3847605e+001  3.4494537e+001  4.4544413e+001  3.1060960e+001  1.8681542e+001  2.6181542e+001  3.3991791e+001  1.7392444e+001  5.2105554e+001  3.5000000e+000  1.0307764e+001  1.9363149e+001  2.5194101e+001  2.8356379e+001  3.2599019e+001  3.6722125e+001  4.3793193e+001  4.9178358e+001  5.0196317e+001  5.5196317e+001  5.5648278e+001  5.1617149e+001  3.6051630e+001  3.2516096e+001  3.6988232e+001  4.6988232e+001  5.1309225e+001  5.5781361e+001  7.0551823e+001  7.4637111e+001  6.9228784e+001  7.0796774e+001  6.3516664e+001  7.1594352e+001  8.1859032e+001  7.7827903e+001  7.3796774e+001  7.7332308e+001  8.0869284e+001  1.0089427e+002
  9.3742889e+001  7.8337113e+001  5.7220577e+001  4.9200435e+001  0.0000000e+000  2.9426289e+001  2.7366468e+001  3.5356854e+001  4.6954268e+001  1.0041606e+002  1.4479967e+002  1.7628191e+002  1.8655062e+002  1.2969629e+002  6.5550234e+001  6.2279671e+001  1.0424823e+002  1.1223434e+002  9.5227811e+001  1.1703955e+002  1.6234594e+002  1.7749524e+002  1.9155062e+002  1.8285241e+002  1.6482466e+002  1.3579967e+002  1.4323271e+002  1.1306865e+002  1.0615294e+002  3.1829327e+001  5.0472922e+001  3.8768223e+001  4.3633327e+001  5.1199700e+001  5.1196908e+001  5.6196908e+001  5.7780281e+001  9.4211191e+001  9.1211191e+001  9.7481303e+001  1.1274823e+002  9.4399371e+001  8.6337113e+001  6.8850280e+001  5.0657392e+001  4.4657392e+001  1.4560220e+001  2.4758259e+001  5.0000000e+000  8.4852814e+000  1.2293168e+001  1.6594331e+001  1.1708204e+001  3.4511712e+001  4.4561588e+001  2.0836971e+001  3.0518894e+001  2.3018894e+001  1.5208644e+001  3.8658304e+001  5.3758259e+001  5.2550748e+001  5.8698486e+001  6.7753871e+001  7.2460328e+001  7.5622605e+001  7.9865246e+001  8.3988352e+001  9.1059419e+001  8.6939439e+001  9.7462544e+001  9.4399371e+001  1.0246163e+002  1.0000787e+002  8.4442351e+001  8.0906817e+001  8.5378953e+001  9.5378953e+001  9.9699947e+001  1.0417208e+002  1.1894254e+002  1.2302783e+002  1.1761951e+002  1.1984752e+002  1.1256741e+002  1.2064510e+002  1.3090978e+002  1.2687865e+002  1.2284752e+002  1.2638306e+002  1.2992003e+002  1.4994502e+002
  9.5375178e+001  9.8420771e+001  7.7304235e+001  5.0023011e+001  2.9426289e+001  0.0000000e+000  2.7672317e+001  3.5662703e+001  4.7260117e+001  1.0072191e+002  1.4510552e+002  1.7658776e+002  1.8685647e+002  1.3000213e+002  6.5856083e+001  6.2585520e+001  1.2433189e+002  1.1386663e+002  9.6860100e+001  1.1786212e+002  1.6265179e+002  1.7780109e+002  1.9185647e+002  1.8315826e+002  1.6513050e+002  1.3610552e+002  1.4353855e+002  1.1337450e+002  1.0645879e+002  3.2135176e+001  5.0778771e+001  3.9074071e+001  4.3939176e+001  5.1505549e+001  5.1502757e+001  5.6502757e+001  5.8086130e+001  9.4517040e+001  9.1517040e+001  1.0919471e+002  1.3283189e+002  1.1448303e+002  1.0642077e+002  8.8933938e+001  5.0963241e+001  4.4963241e+001  1.4866069e+001  2.5064108e+001  3.3194889e+001  2.2754582e+001  1.8946695e+001  2.3247858e+001  3.1791862e+001  5.4595370e+001  6.4645246e+001  2.7490499e+001  3.1341469e+001  2.3841469e+001  1.6031220e+001  3.9480880e+001  5.4064108e+001  5.3373324e+001  6.0330775e+001  6.9386160e+001  7.5217112e+001  7.8379390e+001  8.2622030e+001  8.6745136e+001  9.3816204e+001  9.9201369e+001  1.0021933e+002  1.0521933e+002  1.0567129e+002  1.0164016e+002  8.6074640e+001  8.2539107e+001  8.7011243e+001  9.7011243e+001  1.0133224e+002  1.0580437e+002  1.2057483e+002  1.2466012e+002  1.1925179e+002  1.2067010e+002  1.1338999e+002  1.2146768e+002  1.3173236e+002  1.2770123e+002  1.2367010e+002  1.2720563e+002  1.3074261e+002  1.5076759e+002
  1.1500349e+002  9.7281187e+001  7.6164651e+001  7.6566903e+001  2.7366468e+001  2.7672317e+001  0.0000000e+000  2.4777246e+001  2.9092085e+001  7.3283506e+001  1.1766712e+002  1.4914936e+002  1.5941807e+002  1.0901221e+002  3.8183766e+001  4.1595595e+001  1.1508414e+002  1.3511375e+002  1.1810722e+002  1.4440602e+002  1.4166187e+002  1.5036268e+002  1.6441807e+002  1.5571986e+002  1.3769210e+002  1.0866712e+002  1.1610015e+002  8.5702184e+001  8.0154569e+001  5.8309519e+000  2.3106454e+001  1.1401754e+001  1.6500774e+001  2.4067147e+001  3.3334725e+001  3.5512832e+001  3.0413813e+001  7.3527115e+001  7.0527115e+001  8.8204784e+001  1.1360516e+002  1.1334344e+002  1.0528119e+002  8.7794354e+001  4.0042929e+001  3.4077784e+001  1.2806248e+001  1.2902020e+001  3.2366468e+001  3.5851750e+001  3.9659636e+001  4.3960799e+001  3.9074672e+001  6.1878181e+001  6.3505662e+001  4.8203440e+001  5.7885362e+001  5.0385362e+001  4.2575112e+001  6.6024772e+001  4.1902020e+001  7.9917216e+001  8.6064954e+001  9.5120339e+001  9.1404402e+001  9.4566679e+001  9.8809320e+001  1.0293243e+002  1.1000349e+002  1.0588351e+002  1.1640662e+002  1.1334344e+002  1.2140570e+002  1.2126848e+002  1.0732176e+002  1.0378622e+002  1.0825836e+002  1.1825836e+002  1.2257935e+002  1.2705149e+002  1.4182195e+002  1.4590724e+002  1.4049891e+002  1.4721399e+002  1.3993388e+002  1.4801157e+002  1.5827625e+002  1.5424512e+002  1.5021399e+002  1.5374952e+002  1.5728650e+002  1.5992201e+002
  9.0226247e+001  7.2503941e+001  5.1387404e+001  8.3272828e+001  3.5356854e+001  3.5662703e+001  2.4777246e+001  0.0000000e+000  1.1597414e+001  6.5059205e+001  1.0944282e+002  1.4092506e+002  1.5119377e+002  9.4339432e+001  5.4761836e+001  2.6922817e+001  9.8415056e+001  1.1033650e+002  9.3329971e+001  1.3732277e+002  1.2698909e+002  1.4213838e+002  1.5619377e+002  1.4749556e+002  1.2946780e+002  1.0044282e+002  1.0787585e+002  1.0228025e+002  1.0493182e+002  3.0608198e+001  2.5080192e+001  1.3375492e+001  8.2764727e+000  1.5842846e+001  1.5840054e+001  2.0840054e+001  2.5939074e+001  5.8854337e+001  5.5854337e+001  7.3532007e+001  9.8932386e+001  8.8566199e+001  8.0503941e+001  6.3017108e+001  1.5300538e+001  9.3005376e+000  2.0796634e+001  3.0994673e+001  4.0356854e+001  4.3842135e+001  4.7650022e+001  5.1951185e+001  4.7065058e+001  4.8778291e+001  3.8728416e+001  5.6193825e+001  6.5875748e+001  5.8375748e+001  5.0565498e+001  7.4015158e+001  5.9994673e+001  8.6772828e+001  7.2965064e+001  7.2458107e+001  6.6627155e+001  6.9789433e+001  7.4032073e+001  7.8155179e+001  8.5226247e+001  8.1106266e+001  9.1629371e+001  8.8566199e+001  9.6628456e+001  9.6491229e+001  8.2544511e+001  7.9008977e+001  8.3481113e+001  9.3481113e+001  9.7802107e+001  1.0227424e+002  1.1704470e+002  1.2112999e+002  1.1572167e+002  1.2557052e+002  1.3285063e+002  1.4092832e+002  1.3663278e+002  1.3260165e+002  1.2857052e+002  1.2986212e+002  1.3460553e+002  1.4524923e+002
  9.2254378e+001  7.4532072e+001  5.3415536e+001  8.9866681e+001  4.6954268e+001  4.7260117e+001  2.9092085e+001  1.1597414e+001  0.0000000e+000  5.4241297e+001  9.8624909e+001  1.3010715e+002  1.4037586e+002  8.2742018e+001  5.0238807e+001  1.5325403e+001  8.8813953e+001  1.1236463e+002  9.5358102e+001  1.3935090e+002  1.1539167e+002  1.3132047e+002  1.4537586e+002  1.3667765e+002  1.1864989e+002  8.9624909e+001  9.7057943e+001  9.7757224e+001  1.0724406e+002  3.4923036e+001  2.0557162e+001  1.7690330e+001  1.2591311e+001  5.0249378e+000  4.2426407e+000  9.2426407e+000  1.4341660e+001  4.7256923e+001  4.4256923e+001  6.1934593e+001  8.7334972e+001  9.0594330e+001  8.2532072e+001  6.5045239e+001  1.0950845e+001  1.6950845e+001  3.2394048e+001  4.1994104e+001  5.1954268e+001  5.5439549e+001  5.9247436e+001  6.3548598e+001  5.8662472e+001  5.6428598e+001  4.6378723e+001  6.7791239e+001  7.7473161e+001  6.9973161e+001  6.2162912e+001  8.5612571e+001  7.0994104e+001  9.3366681e+001  8.0615371e+001  7.4486238e+001  6.8655286e+001  7.1817564e+001  7.6060205e+001  8.0183310e+001  8.7254378e+001  8.3134397e+001  9.3657503e+001  9.0594330e+001  9.8656588e+001  9.8519360e+001  8.4572643e+001  8.1037109e+001  8.5509245e+001  9.5509245e+001  9.9830238e+001  1.0430237e+002  1.1907284e+002  1.2315812e+002  1.1774980e+002  1.2759865e+002  1.3487876e+002  1.4295645e+002  1.3866091e+002  1.3462978e+002  1.3059865e+002  1.3189025e+002  1.3663366e+002  1.3365182e+002
  1.4649568e+002  1.2877337e+002  1.0765683e+002  1.4410798e+002  1.0041606e+002  1.0072191e+002  7.3283506e+001  6.5059205e+001  5.4241297e+001  0.0000000e+000  4.4383612e+001  7.5865852e+001  8.6134563e+001  1.2775659e+002  9.4430228e+001  6.9566700e+001  1.4305525e+002  1.6660593e+002  1.4959940e+002  1.9359220e+002  9.5106934e+001  7.7079177e+001  9.1134563e+001  8.2436353e+001  6.4408596e+001  3.5383612e+001  4.2816646e+001  1.4194865e+002  1.5143548e+002  7.9114458e+001  6.4748584e+001  6.1881752e+001  5.6782732e+001  4.9216359e+001  5.8483938e+001  6.3483938e+001  6.8582957e+001  1.0149822e+002  9.8498220e+001  1.1617589e+002  1.4157627e+002  1.4483563e+002  1.3677337e+002  1.1928654e+002  6.5192141e+001  7.1192141e+001  8.5855839e+001  8.6185526e+001  1.0541606e+002  1.0890134e+002  1.1270923e+002  1.1701039e+002  1.1212426e+002  1.1066990e+002  1.0062002e+002  1.2125303e+002  1.3093495e+002  1.2343495e+002  1.1562470e+002  1.3907436e+002  1.1518553e+002  1.4760798e+002  1.3485667e+002  1.2872754e+002  1.2289658e+002  1.2605886e+002  1.3030150e+002  1.3442461e+002  1.4149568e+002  1.3737569e+002  1.4789880e+002  1.4483563e+002  1.5289788e+002  1.5276066e+002  1.3881394e+002  1.3527841e+002  1.3975054e+002  1.4975054e+002  1.5407154e+002  1.5854367e+002  1.7331413e+002  1.7739942e+002  1.7199109e+002  1.8183995e+002  1.8912006e+002  1.9719775e+002  1.9290221e+002  1.8887108e+002  1.8483995e+002  1.8613155e+002  1.9087496e+002  1.8789311e+002
  1.9087929e+002  1.7315698e+002  1.5204044e+002  1.8849159e+002  1.4479967e+002  1.4510552e+002  1.1766712e+002  1.0944282e+002  9.8624909e+001  4.4383612e+001  0.0000000e+000  3.7913528e+001  4.1750951e+001  8.3372977e+001  1.3881384e+002  1.1395031e+002  1.8743886e+002  2.1098954e+002  1.9398301e+002  2.3797581e+002  5.0723322e+001  3.2695565e+001  4.6750951e+001  3.8052741e+001  2.0024984e+001  9.0000000e+000  1.6433034e+001  1.8633226e+002  1.9581909e+002  1.2349807e+002  1.0913220e+002  1.0626536e+002  1.0116634e+002  9.3599971e+001  1.0286755e+002  1.0786755e+002  1.1296657e+002  1.4588183e+002  1.4288183e+002  1.6055950e+002  1.8595988e+002  1.8921924e+002  1.8115698e+002  1.6367015e+002  1.0957575e+002  1.1557575e+002  1.3023945e+002  1.3056914e+002  1.4979967e+002  1.5328495e+002  1.5709284e+002  1.6139400e+002  1.5650787e+002  1.5505351e+002  1.4500363e+002  1.6563664e+002  1.7531856e+002  1.6781856e+002  1.6000831e+002  1.8345797e+002  1.5956914e+002  1.9199159e+002  1.7924028e+002  1.7311115e+002  1.6728020e+002  1.7044247e+002  1.7468511e+002  1.7880822e+002  1.8587929e+002  1.8175931e+002  1.9228241e+002  1.8921924e+002  1.9728150e+002  1.9714427e+002  1.8319755e+002  1.7966202e+002  1.8413415e+002  1.9413415e+002  1.9845515e+002  2.0292728e+002  2.1769774e+002  2.2178303e+002  2.1637471e+002  2.2622356e+002  2.3350367e+002  2.4158136e+002  2.3728582e+002  2.3325469e+002  2.2922356e+002  2.3051516e+002  2.3525857e+002  2.3227672e+002
  2.2236153e+002  2.0463922e+002  1.8352268e+002  2.1997383e+002  1.7628191e+002  1.7658776e+002  1.4914936e+002  1.4092506e+002  1.3010715e+002  7.5865852e+001  3.7913528e+001  0.0000000e+000  5.9770021e+001  1.1950282e+002  1.7029608e+002  1.4543255e+002  2.1892110e+002  2.4247178e+002  2.2546525e+002  2.6945805e+002  8.6853163e+001  6.8825406e+001  6.4770021e+001  3.5916300e+001  1.7888544e+001  4.0482240e+001  3.3049206e+001  2.1781450e+002  2.2730133e+002  1.5498031e+002  1.4061444e+002  1.3774760e+002  1.3264858e+002  1.2508221e+002  1.3434979e+002  1.3934979e+002  1.4444881e+002  1.7736407e+002  1.7436407e+002  1.9204174e+002  2.1744212e+002  2.2070148e+002  2.1263922e+002  1.9515239e+002  1.4105799e+002  1.4705799e+002  1.6172169e+002  1.6205138e+002  1.8128191e+002  1.8476719e+002  1.8857508e+002  1.9287624e+002  1.8799011e+002  1.8653575e+002  1.7648587e+002  1.9711888e+002  2.0680080e+002  1.9930080e+002  1.9149055e+002  2.1494021e+002  1.9105138e+002  2.2347383e+002  2.1072252e+002  2.0459339e+002  1.9876244e+002  2.0192471e+002  2.0616735e+002  2.1029046e+002  2.1736153e+002  2.1324155e+002  2.2376465e+002  2.2070148e+002  2.2876374e+002  2.2862651e+002  2.1467979e+002  2.1114426e+002  2.1561639e+002  2.2561639e+002  2.2993739e+002  2.3440952e+002  2.4917998e+002  2.5326527e+002  2.4785695e+002  2.5770580e+002  2.6498591e+002  2.7306360e+002  2.6876806e+002  2.6473693e+002  2.6070580e+002  2.6199740e+002  2.6674081e+002  2.6375896e+002
  2.2001753e+002  2.0103004e+002  1.8740507e+002  2.0981794e+002  1.8655062e+002  1.8685647e+002  1.5941807e+002  1.5119377e+002  1.4037586e+002  8.6134563e+001  4.1750951e+001  5.9770021e+001  0.0000000e+000  5.9732797e+001  1.8056479e+002  1.2714941e+002  2.0876521e+002  2.4487764e+002  2.2934764e+002  2.7186391e+002  2.7083142e+001  9.0553851e+000  5.0000000e+000  2.3853721e+001  4.1881477e+001  5.0750951e+001  5.8183985e+001  2.2808321e+002  2.3757004e+002  1.6524902e+002  1.5088315e+002  1.4801631e+002  1.4291729e+002  1.3535092e+002  1.3823217e+002  1.3323217e+002  1.3833119e+002  1.6120818e+002  1.6420818e+002  1.8188585e+002  2.0128623e+002  2.1709230e+002  2.0903004e+002  1.9903477e+002  1.4494038e+002  1.5094038e+002  1.7199040e+002  1.7232009e+002  1.9155062e+002  1.9503590e+002  1.9884379e+002  2.0314495e+002  1.9825883e+002  1.9041813e+002  1.8036826e+002  2.0738759e+002  2.1706952e+002  2.0956952e+002  2.0175927e+002  2.2520893e+002  2.0132009e+002  2.1331794e+002  2.1460491e+002  2.0847577e+002  2.0264482e+002  2.0580710e+002  2.1004974e+002  2.1417284e+002  2.1501753e+002  2.0963237e+002  2.2142066e+002  2.1709230e+002  2.2515456e+002  2.2628252e+002  2.1856218e+002  2.1502664e+002  2.1949878e+002  2.2642066e+002  2.3312886e+002  2.3760100e+002  2.4898959e+002  2.5401452e+002  2.5026281e+002  2.6011167e+002  2.6739178e+002  2.7546946e+002  2.6898064e+002  2.6714279e+002  2.6311167e+002  2.6274665e+002  2.6749006e+002  2.4760308e+002
  1.6028474e+002  1.4129725e+002  1.2767227e+002  1.5008515e+002  1.2969629e+002  1.3000213e+002  1.0901221e+002  9.4339432e+001  8.2742018e+001  1.2775659e+002  8.3372977e+001  1.1950282e+002  5.9732797e+001  0.0000000e+000  1.3298082e+002  6.7416615e+001  1.4903242e+002  1.8514485e+002  1.6961484e+002  2.1213111e+002  3.2649655e+001  5.0677412e+001  6.4732797e+001  8.3586518e+001  1.0161427e+002  9.2372977e+001  9.9806012e+001  1.8049924e+002  1.8916678e+002  1.1484316e+002  1.0329918e+002  1.0043235e+002  9.5333329e+001  8.7766956e+001  7.8499378e+001  7.3499378e+001  7.8598397e+001  1.0147539e+002  1.0447539e+002  1.2215306e+002  1.4155344e+002  1.5735950e+002  1.4929725e+002  1.3930198e+002  8.5207582e+001  9.1207582e+001  1.1513607e+002  1.2191423e+002  1.3469629e+002  1.3818157e+002  1.4198945e+002  1.4629062e+002  1.4140449e+002  1.3068534e+002  1.2063546e+002  1.5053326e+002  1.6021518e+002  1.5271518e+002  1.4490493e+002  1.6747759e+002  1.5091423e+002  1.5358515e+002  1.5487211e+002  1.4874298e+002  1.4291202e+002  1.4607430e+002  1.5031694e+002  1.5444005e+002  1.5528474e+002  1.4989957e+002  1.6168786e+002  1.5735950e+002  1.6542176e+002  1.6654972e+002  1.5882938e+002  1.5529385e+002  1.5976598e+002  1.6668786e+002  1.7339607e+002  1.7786820e+002  1.8925679e+002  1.9428173e+002  1.9053001e+002  2.0037887e+002  2.0765898e+002  2.1573667e+002  2.0924784e+002  2.0741000e+002  2.0337887e+002  2.0301385e+002  2.0775727e+002  1.8787028e+002
  1.4249318e+002  1.2477088e+002  1.0365434e+002  1.1475067e+002  6.5550234e+001  6.5856083e+001  3.8183766e+001  5.4761836e+001  5.0238807e+001  9.4430228e+001  1.3881384e+002  1.7029608e+002  1.8056479e+002  1.3298082e+002  0.0000000e+000  6.5564210e+001  1.3905276e+002  1.6260344e+002  1.4559691e+002  1.8258978e+002  1.6563048e+002  1.7150941e+002  1.8556479e+002  1.7686658e+002  1.5883882e+002  1.2981384e+002  1.3724687e+002  4.7518417e+001  5.7005250e+001  4.4014718e+001  2.9681644e+001  4.1386344e+001  4.6485364e+001  4.5213869e+001  5.4481447e+001  5.9481447e+001  6.4580467e+001  9.7495730e+001  9.4495730e+001  1.1217340e+002  1.3757378e+002  1.4083314e+002  1.3277088e+002  1.1528405e+002  6.1189651e+001  6.4062374e+001  5.0990015e+001  5.1085786e+001  7.0550234e+001  7.4035516e+001  7.7843402e+001  8.2144565e+001  7.7258438e+001  1.0006195e+002  9.3490252e+001  8.6387206e+001  9.6069128e+001  8.8569128e+001  8.0758878e+001  1.0420854e+002  8.0085786e+001  1.1810098e+002  1.2424872e+002  1.2472504e+002  1.1889409e+002  1.2205637e+002  1.2629901e+002  1.3042212e+002  1.3749318e+002  1.3337320e+002  1.4389631e+002  1.4083314e+002  1.4889539e+002  1.4875817e+002  1.3481145e+002  1.3127592e+002  1.3574805e+002  1.4574805e+002  1.5006904e+002  1.5454118e+002  1.6931164e+002  1.7339693e+002  1.6798860e+002  1.7783746e+002  1.7811765e+002  1.8619533e+002  1.8889972e+002  1.8486859e+002  1.8083746e+002  1.8212905e+002  1.8687247e+002  1.8389062e+002
  9.2868122e+001  7.3880632e+001  6.0255658e+001  8.2668530e+001  6.2279671e+001  6.2585520e+001  4.1595595e+001  2.6922817e+001  1.5325403e+001  6.9566700e+001  1.1395031e+002  1.4543255e+002  1.2714941e+002  6.7416615e+001  6.5564210e+001  0.0000000e+000  8.1615802e+001  1.1772823e+002  1.0219822e+002  1.4471450e+002  1.0006627e+002  1.1809403e+002  1.3214941e+002  1.5100313e+002  1.3397530e+002  1.0495031e+002  1.1238335e+002  1.1308263e+002  1.2175016e+002  4.7426547e+001  3.5882566e+001  3.3015734e+001  2.7916714e+001  2.0350341e+001  1.1082763e+001  6.0827625e+000  1.1181782e+001  3.4058773e+001  3.7058773e+001  5.4736442e+001  7.4136822e+001  8.9942890e+001  8.1880632e+001  7.1885361e+001  1.7790966e+001  2.3790966e+001  4.7719451e+001  5.4497614e+001  6.7279671e+001  7.0764952e+001  7.4572839e+001  7.8874001e+001  7.3987875e+001  6.3268720e+001  5.3218844e+001  8.3116642e+001  9.2798564e+001  8.5298564e+001  7.7488315e+001  1.0006097e+002  8.3497614e+001  8.6168530e+001  8.7455493e+001  8.1326360e+001  7.5495408e+001  7.8657686e+001  8.2900327e+001  8.7023432e+001  8.7868122e+001  8.2482957e+001  9.4271246e+001  8.9942890e+001  9.8005147e+001  9.9133104e+001  9.1412764e+001  8.7877230e+001  9.2349366e+001  9.9271246e+001  1.0597945e+002  1.1045159e+002  1.2184017e+002  1.2686511e+002  1.2311340e+002  1.3296225e+002  1.4024236e+002  1.4832005e+002  1.4183122e+002  1.3999338e+002  1.3596225e+002  1.3559724e+002  1.4034065e+002  1.2045367e+002
  3.5912053e+001  2.5911116e+001  4.7027652e+001  7.4705246e+001  1.0424823e+002  1.2433189e+002  1.1508414e+002  9.8415056e+001  8.8813953e+001  1.4305525e+002  1.8743886e+002  2.1892110e+002  2.0876521e+002  1.4903242e+002  1.3905276e+002  8.1615802e+001  0.0000000e+000  4.6932575e+001  5.3495518e+001  6.7003800e+001  1.8168207e+002  1.9970983e+002  2.1376521e+002  2.2549160e+002  2.0746385e+002  1.7843886e+002  1.8587190e+002  1.8657118e+002  1.9523871e+002  1.2091510e+002  1.0937112e+002  1.0650428e+002  1.0140526e+002  9.3838891e+001  8.4571312e+001  7.9571312e+001  8.4670332e+001  4.7557030e+001  4.4557030e+001  2.6879360e+001  8.5000000e+000  9.8488578e+000  1.7911116e+001  3.5397949e+001  8.9492343e+001  8.9114519e+001  1.1880845e+002  1.2798616e+002  9.9248229e+001  1.0919308e+002  1.0538519e+002  1.0108403e+002  9.2540025e+001  6.9736516e+001  5.9686641e+001  1.0532667e+002  9.3386788e+001  1.0088679e+002  1.0830067e+002  9.2097690e+001  1.2681080e+002  7.8205246e+001  6.4397482e+001  5.5342097e+001  4.9511145e+001  4.6348868e+001  4.2106227e+001  3.7983121e+001  3.0912053e+001  2.55