我有一个例程,它将整数列表作为向量返回。
这些整数来自连续数字组;例如,它可能如下所示:
vector = 6 7 8 12 13 14 15 26 27 28 29 30 55 56
请注意,上面有四个'运行'数字(6-8,12-15,26-30和55-56)。我想做的是转发最长'运行'数字到一个新的向量。在这种情况下,这将是26-30运行,所以我想生产:
newVector = 26 27 28 29 30
这个计算必须在不同的向量上执行很多次,所以我能做得越多越好!任何智慧都会被感激地接受。
答案 0 :(得分:36)
你可以试试这个:
v = [ 6 7 8 12 13 14 15 26 27 28 29 30 55 56];
x = [0 cumsum(diff(v)~=1)];
v(x==mode(x))
这导致
ans =
26 27 28 29 30
答案 1 :(得分:2)
这是让球滚动的解决方案。 。 。
vector = [6 7 8 12 13 14 15 26 27 28 29 30 55 56]
d = [diff(vector) 0]
maxSequence = 0;
maxSequenceIdx = 0;
lastIdx = 1;
while lastIdx~=find(d~=1, 1, 'last')
idx = find(d~=1, 1);
if idx-lastIdx > maxSequence
maxSequence = idx-lastIdx;
maxSequenceIdx = lastIdx;
end
d(idx) = 1;
lastIdx=idx;
end
output = vector(1+maxSequenceIdx:maxSequenceIdx+maxSequence)
在此示例中,diff命令用于查找连续的数字。当数字是连续的时,差值为1.然后使用while循环查找最长的1组,并存储该连续组的索引。但是,我相信这可以进一步优化。
答案 2 :(得分:1)
没有循环 - 应该更快
temp = find ( ([(vector(2:end) - vector(1:end-1))==1 0])==0);
[len,ind]=max(temp(2:end)-temp(1:end-1));
vec_out = vector(temp(ind)+1:temp(ind)+len)
答案 3 :(得分:1)
使用diff
无循环:
vector = [6 7 8 12 13 14 15 26 27 28 29 30 55 56];
seqGroups = [1 find([1 diff(vector)]~=1) numel(vector)+1]; % beginning of group
[~, groupIdx] = max( diff(seqGroups)); % bigger group index
output = vector( seqGroups(groupIdx):seqGroups(groupIdx+1)-1)
输出向量是
ans =
26 27 28 29 30