你这个循环结果应该有sum(1:length(hb[,"city"]))个(抱歉写这么罗嗦，因为论坛显示美元那个符号可能会不正常)，但是因为随着i增加，每次j循环的结果都会少一个，所以处理起来有点小麻烦。
要把结果形成向量，不如先设个i*i的矩阵，然后再形成向量。
我不改你的代码了，写了个简单的，你参考一下，再琢磨琢磨自己怎么改吧。至于求和就不说了。