Используя также m_ind, n_ind = w.T, просто сделайте M += np.bincount(m_ind, weights=N[n_ind], minlength=len(M))
m_ind, n_ind = w.T
M += np.bincount(m_ind, weights=N[n_ind], minlength=len(M))