r2c = np.mod(np.arange(np.power(5, 4)), 10).reshape((5,) * 4)