使用vertcat的Matlab错误内存不足

Question

我正在运行以下代码。 它适用于较小的NumberOfVariables，但由于内存不足，不适用于8或更大的NumberOfVariables。 我真的只需要第一个不包含零的AllAnswers，但如果可能的话，需要所有AllAnswers。

NumberOfVariables = 9;
k=NumberOfVariables^2-NumberOfVariables+1;
integers = 0:k-1;
numbers = 1:k-1;
tic
s = combnk(integers,NumberOfVariables);
AllAnswers = [];
for i = 1:size(s,1)
G=combnk(s(i,:),2);
G = [(G(:,1)'-G(:,2)') (G(:,2)'-G(:,1)')];
G = sort(mod(G,k));
if (isequal(G,numbers))
AllAnswers = [AllAnswers;s(i,:)];
end
end
toc

Answer 1

 s = combnk(integers,NumberOfVariables);

是所有Number of Variables大小的integers组合的列表。 所以这是一个清单

ķ！/（（NumberOfVariables！（K-NumberOfVariables）！）

或者，因为k = NumberOfVariables ^ 2-NumberOfVariables

（NumberOfVariables ^ 2-NumberOfVariables）！/（（NumberOfVariables！（NumberOfVariables ^ 2-2 * NumberOfVariables）！）

使用x表示NumberOfVariables，以保持可读性：

      (x²-x)!        x²!            x²!            x²!     
#s= ----------  > ---------- > ------------ = ------------ 
     x!(x²-2x)!   x!(x²-2x)!   x!(x²-2x+1)!    x!((x-1)²!)

你明白了。 如果您实际上是通过为#s元素分配内存来接近它的话，那么这不是您的朋友。 实际上，看着这个，对于x> = 4，它的增长速度将比e ^ x快。

对于从10开始的NumberOfVariables，它变得非常大。

算一算！

https://www.wolframalpha.com/share/clip?f=d41d8cd98f00b204e9800998ecf8427euo0stpao37