Let x be a one-dimensional n-vector, and let J = np.ones(n,n) be the n-by-n matrix of all ones. Computing y = J @ x does n**2 multiplications and additions. How can you compute y in O(n) arithmetic operations without forming J? (Note that np.ones(n) is the 1-dimensional n-vector of all ones.)