プロファイラだと、　　XY[r] += b[i]*b[j];がボドルネック。