Джулии: это массив содержит конкретный вложенный массив

Question

В Юле мы можем проверить, если массив содержит значение, например так:

> 6 in [4,6,5] 
true 

Однако это возвращает ложь, при попытке проверить наличие подрешетки в определенном порядке:

> [4,6] in [4,6,5] 
false 

Что такое правильный синтаксис, чтобы проверить, если конкретный суб-массив существует в массиве?

Answer 1

Для т.е. вектора третьего состояния [4,6] появляется как суб-вектор 4,6,5 следующая функция предлагается:

issubvec(v,big) = 
    any([v == slice(big,i:(i+length(v)-1)) for i=1:(length(big)-length(v)+1)]) 

Для второго условия, то есть, дать логический для каждого элемента в els векторах, появляется в set векторе, следующее предлагается:

function vecin(els,set) 
    res = zeros(Bool,size(els)) 
    res[findin(els,set)]=true 
    res 
end 

с вектором в ОП, эти результаты:

julia> vecin([4,6],[4,6,5]) 
2-element Array{Bool,1}: 
true 
true 

julia> issubvec([4,6],[4,6,5]) 
true 

Answer 2

Это занимает немного кода, чтобы сделать функцию, которая выполняет хорошо, но это намного быстрее, чем версия выше issubvec:

function subset2(x,y) 
    lenx = length(x) 
    first = x[1] 
    if lenx == 1 
     return findnext(y, first, 1) != 0 
    end 
    leny = length(y) 
    lim = length(y) - length(x) + 1 
    cur = 1 
    while (cur = findnext(y, first, cur)) != 0 
     cur > lim && break 
     beg = cur 
     @inbounds for i = 2:lenx 
      y[beg += 1] != x[i] && (beg = 0 ; break) 
     end 
     beg != 0 && return true 
     cur += 1 
    end 
    false 
end 

Примечание: также было бы гораздо полезнее, если функция фактически вернула позицию начала подмассива, если она найдена, или 0, если нет, аналогично функциям findfirst/findnext.

информация Timing (второй использует свою функцию subset2):

0.005273 seconds (65.70 k allocations: 4.073 MB) 
    0.000086 seconds (4 allocations: 160 bytes) 

Answer 3

Я использовал это в последнее время, чтобы найти подпоследовательности в массивы целых чисел. Это не так хорошо или быстро, как @ scott's subset2(x,y) ... но он возвращает индексы.

function findsequence(arr::Array{Int64}, seq::Array{Int64}) 
    indices = Int64[] 
    i = 1 
    n = length(seq) 
    if n == 1 
     while true 
      occurrence = findnext(arr, seq[1], i) 
      if occurrence == 0 
       break 
      else 
       push!(indices, occurrence) 
       i = occurrence +1 
      end 
     end 
    else 
     while true 
      occurrence = Base._searchindex(arr, seq, i) 
      if occurrence == 0 
       break 
      else 
       push!(indices, occurrence) 
       i = occurrence +1 
      end 
     end 
    end 
    return indices 
end 

julia> @time findsequence(rand(1:9, 1000), [2,3]) 
    0.000036 seconds (29 allocations: 8.766 KB) 
    16-element Array{Int64,1}: 
    80 
    118 
    138 
    158 
    234 
    243 
    409 
    470 
    539 
    589 
    619 
    629 
    645 
    666 
    762 
    856