我有一个数组,例如[3 5 7 11]
。我想找出这些元素的独特倍数,即
我目前很困惑如何以有效的方式解决它,以便当我将它扩展到数组中的大量元素时,我不会爆炸我的内存并杀死我的计算机。
我目前的代码:
prime_factor=[3 5 7 11];
for i=1:length(prime_factor)
pfm=unique([pfm cumprod(circshift(prime_factor,i,2))])
endfor
魔术函数是nchoosek
,nchoosek(v,k)
从输入矢量k
产生v
元素的所有组合。只需沿输出行计算产品即可获得所需结果:
v = [3 5 7 11];
prod(nchoosek(v,2),2)
收益率:
ans =
15
21
33
35
55
77
prod(nchoosek(V,3),2)
收益率:
ans =
105
165
231
385
prod(nchoosek(V,4),2)
收益率:
ans =
1155