What is the Mathematica equivalent of the following Python code with the vectors' broadcast addition?
import numpy as np
a = np.random.rand(5000, 1, 5);
b = np.random.rand(1, 500, 5);
result = a + b #shape: (5000, 500, 5)
import numpy as np
a = np.random.rand(5000, 1, 5);
b = np.random.rand(1, 500, 5);
result = a + b #shape: (5000, 500, 5)