Hmm, mogę wymyślić dwa możliwe algorytmy: liniowy skan sekwencji A lub zbudowanie słownika z ciągłym wyszukiwaniem indeksów.
Jeśli testujesz wiele potencjalnych podsekwencji B względem pojedynczej większej sekwencji A , sugerowałbym użycie wariantu ze słownikiem.
Skan liniowy
Opis
Utrzymujemy kursor dla sekwencji A . Następnie iterację wszystkich elementów w podciąg B . Dla każdego elementu przesuwamy kursor do przodu w A, aż znajdziemy pasujący element. Jeśli nie znaleziono pasującego elementu, B nie jest podsekwencją.
To zawsze działa w O (rozmiar sekw . ) .
Pseudo kod
Tryb rozkazujący:
def subsequence? seq, subseq:
i = 0
for item in subseq:
i++ while i < seq.size and item != seq[i]
return false if i == seq.size
return true
Funkcjonalny styl:
let rec subsequence? = function
| _ [] -> true
| [] _ -> false
| cursor::seq item::subseq ->
if cursor = item
then subsequence? seq subseq
else subsequence? seq item::subseq
Przykładowa implementacja (Perl):
use strict; use warnings; use signatures; use Test::More;
sub is_subsequence_i ($seq, $subseq) {
my $i = 0;
for my $item (@$subseq) {
$i++ while $i < @$seq and $item != $seq->[$i];
return 0 if $i == @$seq;
}
return 1;
}
sub is_subsequence_f ($seq, $subseq) {
return 1 if @$subseq == 0;
return 0 if @$seq == 0;
my ($cursor, @seq) = @$seq;
my ($item, @subseq) = @$subseq;
return is_subsequence_f(\@seq, $cursor == $item ? \@subseq : $subseq);
}
my $A = [1, 2, 3, 4];
my $B = [1, 3];
my $C = [1, 3, 4];
my $D = [3, 1];
my $E = [1, 2, 5];
for my $is_subsequence (\&is_subsequence_i, \&is_subsequence_f) {
ok $is_subsequence->($A, $B), 'B in A';
ok $is_subsequence->($A, $C), 'C in A';
ok ! $is_subsequence->($A, $D), 'D not in A';
ok ! $is_subsequence->($A, $E), 'E not in A';
ok $is_subsequence->([1, 2, 3, 4, 3, 5, 6], [2, 3, 6]), 'multiple nums';
}
done_testing;
Wyszukiwanie w słowniku
Opis
Mapujemy elementy sekwencji A na ich indeksy. Następnie wyszukujemy odpowiednie indeksy dla każdego elementu w B , pomijamy te, które są zbyt małe, i wybieramy najmniejszy możliwy indeks jako dolny limit. Gdy nie znaleziono żadnych wskaźników, B nie jest podsekwencją.
Działa w coś takiego jak O (rozmiar podsekcji · k) , gdzie k opisuje liczbę zduplikowanych liczb seq. Plus narzut O (rozmiar sek.)
Zaletą tego rozwiązania jest to, że decyzja negatywna może być podjęta znacznie szybciej (aż do stałego czasu), po zapłaceniu kosztów ogólnych związanych z budowaniem tabeli odnośników.
Pseudo kod:
Tryb rozkazujący:
# preparing the lookup table
dict = {}
for i, x in seq:
if exists dict[x]:
dict[x].append(i)
else:
dict[x] = [i]
def subsequence? subseq:
min_index = -1
for x in subseq:
if indices = dict[x]:
suitable_indices = indices.filter(_ > min_index)
return false if suitable_indices.empty?
min_index = suitable_indices[0]
else:
return false
return true
Funkcjonalny styl:
let subsequence? subseq =
let rec subseq-loop = function
| [] _ -> true
| x::subseq min-index ->
match (map (filter (_ > min-index)) data[x])
| None -> false
| Some([]) -> false
| Some(new-min::_) -> subseq-loop subseq new-min
in
subseq-loop subseq -1
Przykładowa implementacja (Perl):
use strict; use warnings; use signatures; use Test::More;
sub build_dict ($seq) {
my %dict;
while (my ($i, $x) = each @$seq) {
push @{ $dict{$x} }, $i;
}
return \%dict;
}
sub is_subsequence_i ($seq, $subseq) {
my $min_index = -1;
my $dict = build_dict($seq);
for my $x (@$subseq) {
my $indices = $dict->{$x} or return 0;
($min_index) = grep { $_ > $min_index } @$indices or return 0;
}
return 1;
}
sub is_subsequence_f ($seq, $subseq) {
my $dict = build_dict($seq);
use feature 'current_sub';
return sub ($subseq, $min_index) {
return 1 if @$subseq == 0;
my ($x, @subseq) = @$subseq;
my ($new_min) = grep { $_ > $min_index } @{ $dict->{$x} // [] } or return 0;
__SUB__->(\@subseq, $new_min);
}->($subseq, -1);
}
my $A = [1, 2, 3, 4];
my $B = [1, 3];
my $C = [1, 3, 4];
my $D = [3, 1];
my $E = [1, 2, 5];
for my $is_subsequence (\&is_subsequence_i, \&is_subsequence_f) {
ok $is_subsequence->($A, $B), 'B in A';
ok $is_subsequence->($A, $C), 'C in A';
ok ! $is_subsequence->($A, $D), 'D not in A';
ok ! $is_subsequence->($A, $E), 'E not in A';
ok $is_subsequence->([1, 2, 3, 4, 3, 5, 6], [2, 3, 6]), 'multiple nums';
}
done_testing;
Wariant wyszukiwania w słowniku: Kodowanie jako skończona maszyna stanowa
Opis
Możemy jeszcze bardziej zmniejszyć złożoność algorytmiczną do O (rozmiar podsekw.), Jeśli wymienimy więcej pamięci. Zamiast mapować elementy na ich indeksy, tworzymy wykres, w którym każdy węzeł reprezentuje element w jego indeksie. Krawędzie pokazują możliwe przejścia, np. Sekwencja a, b, abędzie miała krawędzie a@1 → b@2, a@1 → a@3, b@2 → a@3. Ten wykres jest równoważny skończonej maszynie stanów.
Podczas wyszukiwania utrzymujemy kursor, który początkowo jest pierwszym węzłem drzewa. Następnie spacer krawędzi dla każdego elementu w podmenu B . Jeśli taka krawędź nie istnieje, B nie jest listą podrzędną. Jeśli po wszystkich elementach kursor zawiera prawidłowy węzeł, wówczas B jest listą podrzędną.
Pseudo kod
Tryb rozkazujący:
# preparing the graph
graph = {}
for x in seq.reverse:
next_graph = graph.clone
next_graph[x] = graph
graph = next_graph
def subseq? subseq:
cursor = graph
for x in subseq:
cursor = graph[x]
return false if graph == null
return true
Funkcjonalny styl:
let subseq? subseq =
let rec subseq-loop = function
| [] _ -> true
| x::subseq graph -> match (graph[x])
| None -> false
| Some(next-graph) -> subseq-loop subseq next-graph
in
subseq-loop subseq graph
Przykładowa implementacja (Perl):
use strict; use warnings; use signatures; use Test::More;
sub build_graph ($seq) {
my $graph = {};
for (reverse @$seq) {
$graph = { %$graph, $_ => $graph };
}
return $graph;
}
sub is_subsequence_i ($seq, $subseq) {
my $cursor = build_graph($seq);
for my $x (@$subseq) {
$cursor = $cursor->{$x} or return 0;
}
return 1;
}
sub is_subsequence_f ($seq, $subseq) {
my $graph = build_graph($seq);
use feature 'current_sub';
return sub ($subseq, $graph) {
return 1 if @$subseq == 0;
my ($x, @subseq) = @$subseq;
my $next_graph = $graph->{$x} or return 0;
__SUB__->(\@subseq, $next_graph);
}->($subseq, $graph);
}
my $A = [1, 2, 3, 4];
my $B = [1, 3];
my $C = [1, 3, 4];
my $D = [3, 1];
my $E = [1, 2, 5];
for my $is_subsequence (\&is_subsequence_i, \&is_subsequence_f) {
ok $is_subsequence->($A, $B), 'B in A';
ok $is_subsequence->($A, $C), 'C in A';
ok ! $is_subsequence->($A, $D), 'D not in A';
ok ! $is_subsequence->($A, $E), 'E not in A';
ok $is_subsequence->([1, 2, 3, 4, 3, 5, 6], [2, 3, 6]), 'multiple nums';
}
done_testing;