尝试编写解决方案以在10维空间中查找向量B={b1, b2..b10}的完整列表，使得（a）每个b1..b10是0到30之间的整数（包括0和30）; （b）t(B) %*% M %*% B的值介于两个标量值0.98*K和K之间，其中M a为10x10矩阵。
如下所示的蛮力方法耗尽了这个大小的数组的空间。
M<- array(dim=c(30, 10))
d <- do.call(expand.grid,lapply(dim(M),seq))等......
是否有人对（a）我可以使用的求解器，或（b）解决计算空间问题的非求解方法有想法？