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