Разделить массив на блоки, где каждый блок является самой длинной восходящей (увеличивающейся) последовательностью

Question

Я пытаюсь взять вектор и разделить на разделы. код требует от меня создания отдельного вектора, когда возрастающий порядок останавливается следующим образом: [3, 5, 9, 21, 27, 15, 10, 13] будет разделен на три блока [3 5 9 21 27], [15] и [10, 13].

Я не должен использовать петлю, поэтому было бы полезно, если кто-то сможет мне помочь. Благодарю.

Answer 1

Функция diff вычисляет разницу между смежными элементами массива. Элемент, который является отрицательным в результате diff, указывает на то место, где восходящий порядок «сломан».

Итак:

v = [3, 5, 9, 21, 27, 15, 10, 13]; 
dv = diff(v); 
nb = sum(dv <= 0); % number of blocks 
b = cell(1, nb); % container of blocks 

ix = find(dv <= 0); % start indices of blocks 
ix0 = 1; 
for k = 1:nb 
     ix1 = ix(k); 
     b{k} = a(ix0:ix1); 
     ix0 = ix1 +1 ; 
end; 

Мое мнение таково, что трудно избежать с помощью петель, однако, приведенный выше код является гораздо более «мягким» относительно вычислительной нагрузки по сравнению с грубым силовым методом нахождения блоков элемента по элементу. не

Answer 2

Нет петли, одна линия, я надеюсь, что все в порядке;)

a = [3, 5, 9, 21, 27, 15, 10, 13]; 

output = accumarray(cumsum([0; diff(a(:))] < 0)+1, a, [], @(x) {x}) 

некоторые пояснения:

%// find out where the vector decreases: 
da = [0; diff(a(:))] 
%// and create a mask 
mask = da < 0 
%// sum the mask up cumulatively and add +1 
subs = cumsum(mask) + 1 
%// use accumarray to gather everything 
output = accumarray(subs,a,[],@(x) {x}) 

---

Если есть последующее такое же количество, как и здесь:

a = [3, 5, 9, 21, 27, 27, 15, 10, 13]; 

решение выше считает seco й 27 к первой группе, если вы хотите быть отдельная группа, изменить маску:

mask = da <= 0 

Answer 3

Это можно сделать очень легко с mat2cell:

x = find(diff(a)<0); 
result = mat2cell(a, 1, [x(1) diff(x) numel(a)-x(end)]); 

Если вы хотите один -liner:

result = mat2cell(a, 1, diff(find(diff([inf a -inf])<0)));