擅长:python、mysql、java
<p>与@Divakar非常相似,但不必排序:</p>
<pre><code>seed_repeated = seed[data[:,0]]
dist_to_center = np.sqrt(np.sum((data[:,1:]-seed_repeated[:,1:])**2, axis=1))
diam = np.zeros(len(seed))
np.maximum.at(diam, data[:,0], dist_to_center)
</code></pre>
<p>众所周知,<code>ufunc.at</code>速度很慢,所以看看哪个更快会很有趣。在</p>