とりあえず、素直に数式を落とし込むとこんな感じでしょうか。
Python
1import numpy as np
2
3A=np.reshape(np.arange(1,7), (2,3))
4B=np.reshape(np.arange(1,7), (3,2))
5C=np.array([[2,3],[1,2]])
6x=np.array([-1,1,3])
7
8M,N = A.shape[0], B.shape[1]
9
10ret = np.zeros((M,N))
11for i in range(M):
12 for j in range(N):
13 ret[i,j] = sum(A[i] * x * B[:,j]) - C[j,i]
14
15print(ret)
16# => [[ 48. 59.]
17# [ 98. 118.]]
18