#include <iostream.h>
#include <math.h>
#include <stdlib.h>
void main()
{
int n;
float *X,*Y,*V,*R;
cout<<"请输入结点个数：";
cin>>n;
X=new float[n];
Y=new float[n];
V=new float[n];
R=new float[n];
cout<<"请输入各结点的横坐标：";
for (int i=0;i<n;i++)
cin>>X[i];
cout<<"请输入各结点的纵坐标：";
for (i=0;i<n;i++)
cin>>Y[i];
cout<<"请输入各结点的供货量或需求量：";
for (i=0;i<n;i++)
cin>>V[i];
cout<<"请输入各结点的内/外向运输费率：";
for (i=0;i<n;i++)
cin>>R[i];
float sumx=0.0,sumy=0.0,sumqr=0.0;
float avex,avey;
for (i=0;i<n;i++)
{
sumx+=V[i]*R[i]*X[i];
sumy+=V[i]*R[i]*Y[i];
sumqr+=V[i]*R[i];
}
avex=sumx/sumqr;
avey=sumy/sumqr;
cout<<"初始仓库地址为："<<"X="<<avex<<' '<<"Y="<<avey<<endl;
float *d;
d=new float[i];
float xzsumx=0.0,xzsumy=0.0,xzsumqr=0.0; //xz表示修正
float xzavex,xzavey;
int j=0;
while (1)
{
for (i=0;i<n;i++)
 
d[i]=sqrt((avex-X[i])*(avex-X[i])+(avey-Y[i])*(avey-Y[i]));
for (i=0;i<n;i++)
{
xzsumx+=V[i]*R[i]*X[i]/d[i];
xzsumy+=V[i]*R[i]*Y[i]/d[i];
xzsumqr+=V[i]*R[i]/d[i];
}
xzavex=xzsumx/xzsumqr;
xzavey=xzsumy/xzsumqr;
if (j==60)
break;
else
{
cout<<"第"<<j+1<<"次迭代："<<"tempX="<<xzavex<<' '<<"tempY="<<xzavey<<endl;
avex=xzavex;
avex=xzavey;
xzsumx=0.0;
xzsumy=0.0;
xzsumqr=0.0;
j++;
}
}
cout<<"仓库最佳地址为："<<"X="<<xzavex<<' '<<"Y="<<xzavey<<endl;
}
}
}
cout<<"仓库最佳地址为："<<"X="<<xzavex<<' '<<"Y="<<xzavey<<endl;
}