#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;
}

貌似是c++的吧！

是c++的

怎一个乱字了得：

#include <iostream>
#include <math.h>
#include <stdlib.h>

void main()
{
    int n,i;
    float *X,*Y,*V,*R;
    std::cout<<"请输入结点个数：";
    std::cin>>n;
    X=new float[n];
    Y=new float[n];
    V=new float[n];
    R=new float[n];
    std::cout<<"请输入各结点的横坐标：";
    for (i=0;i<n;i++)
        std::cin>>X[i];
    std::cout<<"请输入各结点的纵坐标：";
    for (i=0;i<n;i++)
        std::cin>>Y[i];
    std::cout<<"请输入各结点的供货量或需求量：";
    for (i=0;i<n;i++)
        std::cin>>V[i];
    std::cout<<"请输入各结点的内/外向运输费率：";
    for (i=0;i<n;i++)
        std::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;
    std::cout<<"初始仓库地址为："<<"X="<<avex<<' '<<"Y="<<avey<<std::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
            {
                std::cout<<"第"<<j+1<<"次迭代："<<"tempX="<<xzavex<<' '<<"tempY="<<xzavey<<std::endl;
                avex=xzavex;
                avex=xzavey;
                xzsumx=0.0;
                xzsumy=0.0;
                xzsumqr=0.0;
                j++;
            }
    }
    std::cout<<"仓库最佳地址为："<<"X="<<xzavex<<' '<<"Y="<<xzavey<<std::endl;
    std::cout<<"仓库最佳地址为："<<"X="<<xzavex<<' '<<"Y="<<xzavey<<std::endl;
}

编译通过，
建议楼主看C++ primer

[ 本帖最后由 heroinearth 于 2013-10-7 23:53 编辑 ]