7 2 2 10
10 -3
-3  2
-2  1

first line: n (# responses), d (# c values), p (# y values), M (# solutions)

second line: y = yTbot

rest: R = HRbot, d-by-p

% python3 -i solve.py
n = 7 , X = [1, 2, 3, 1]
nm = 18 , m = 2.5714285714285716 , vnn = 0 , v = 0.0 , s = 0.0
[2, 0, 4, 1]
[3, 0, 1, 3]
[1, 1, 5, 0]
[2, 1, 2, 2]
[1, 2, 3, 1]
[2, 2, 0, 3]
[0, 3, 4, 0]
[1, 3, 1, 2]
[0, 4, 2, 1]
[0, 5, 0, 2]
search_count = 21 , solution_count = 10
>>> cTlist
[[2, 0], [3, 0], [1, 1], [2, 1], [1, 2], [2, 2], [0, 3], [1, 3], [0, 4], [0, 5]]
>>> yTbot
Matrix([[10, -3]])
>>> HRbot
Matrix([
[-3, 2],
[-2, 1]])
>>> avg(cTlist)
[1.2, 2.1]
>>> avg(xTlist)
[1.2, 2.1, 2.2, 1.5]
>>>