NAME
LCS - Longest Common Subsequence
SYNOPSIS
use LCS;
my $lcs = LCS->LCS( [qw(a b)], [qw(a b b)] );
# $lcs now contains an arrayref of matching positions
# same as
$lcs = [
[ 0, 0 ],
[ 1, 2 ]
];
my $all_lcs = LCS->allLCS( [qw(a b)], [qw(a b b)] );
# same as
$all_lcs = [
[
[ 0, 0 ],
[ 1, 1 ]
],
[
[ 0, 0 ],
[ 1, 2 ]
]
];
DESCRIPTION
LCS is an implementation based on the traditional LCS algorithm.
It contains reference implementations working slow but correct.
Also some utility methods are added to reformat the result.
CONSTRUCTOR
new()
Creates a new object which maintains internal storage areas for the
LCS computation. Use one of these per concurrent LCS() call.
METHODS
LCS(\@a,\@b)
Finds a Longest Common Subsequence, taking two arrayrefs as method
arguments. It returns an array reference of corresponding indices,
which are represented by 2-element array refs.
# position 0 1 2
my $a = [qw(a b )];
my $b = [qw(a b b)];
my $lcs = LCS->LCS($a,$b);
# same like
$lcs = [
[ 0, 0 ],
[ 1, 1 ]
];
LLCS(\@a,\@b)
Calculates the length of the Longest Common Subsequence.
my $llcs = LCS->LLCS( [qw(a b)], [qw(a b b)] );
print $llcs,"\n"; # prints 2
# is the same as
$llcs = scalar @{LCS->LCS( [qw(a b)], [qw(a b b)] )};
allLCS(\@a,\@b)
Finds all Longest Common Subsequences. It returns an array reference
of all LCS.
my $all_lcs = LCS->allLCS( [qw(a b)], [qw(a b b)] );
# same as
$all_lcs = [
[
[ 0, 0 ],
[ 1, 1 ]
],
[
[ 0, 0 ],
[ 1, 2 ]
]
];
The purpose is mainly for testing LCS algorithms, as they only return
one of the optimal solutions. If we want to know, that the result is
one of the optimal solutions, we need to test, if the solution is
part of all optimal LCS:
use Test::More;
use Test::Deep;
use LCS;
use LCS::Tiny;
cmp_deeply(
LCS::Tiny->LCS(\@a,\@b),
any(@{LCS->allLCS(\@a,\@b)} ),
"Tiny::LCS $a, $b"
);
lcs2align(\@a,\@b,$LCS)
Returns the two sequences aligned, missing positions are represented
as empty strings.
use Data::Dumper;
use LCS;
print Dumper(
LCS->lcs2align(
[qw(a b)],
[qw(a b b)],
LCS->LCS([qw(a b)],[qw(a b b)])
)
);
# prints
$VAR1 = [
[
'a',
'a'
],
[
'',
'b'
],
[
'b',
'b'
]
];
align(\@a,\@b)
Returns the same as lcs2align() via calling LCS() itself.
sequences2hunks($a, $b)
Transforms two array references of scalars to an array of hunks (two
element arrays).
hunks2sequences($hunks)
Transforms an array of hunks to two arrays of scalars.
use Data::Dumper;
use LCS;
print Dumper(
LCS->hunks2sequences(
LCS->LCS([qw(a b)],[qw(a b b)])
)
);
# prints (reformatted)
$VAR1 = [ 0, 1 ];
$VAR2 = [ 0, 2 ];
align2strings($hunks, $gap_character)
Returns two strings aligned with gap characters. The default gap
character is '_'.
use Data::Dumper;
use LCS;
print Dumper(
LCS->align2strings(
LCS->lcs2align([qw(a b)],[qw(a b b)],LCS->LCS([qw(a b)],[qw(a b b)]))
)
);
$VAR1 = 'a_b';
$VAR2 = 'abb';
fill_strings($string1, $string2, $fill_character)
Returns both strings filling up the shorter with $fill_character to
the same length.
The default $fill_character is '_'.
clcs2lcs($compact_lcs)
Convert compact LCS to LCS.
lcs2clcs($compact_lcs)
Convert LCS to compact LCS.
max($i, $j)
Returns the maximum of two numbers.
EXPORT
None by design.
STABILITY
Until release of version 1.00 the included methods, names of methods
and their interfaces are subject to change.
Beginning with version 1.00 the specification will be stable, i.e. not
changed between major versions.
REFERENCES
Ronald I. Greenberg. Fast and Simple Computation of All Longest Common
Subsequences,
http://arxiv.org/pdf/cs/0211001.pdf
Robert A. Wagner and Michael J. Fischer. The string-to-string
correction problem. Journal of the ACM, 21(1):168-173, 1974.
SOURCE REPOSITORY
http://github.com/wollmers/LCS
AUTHOR
Helmut Wollmersdorfer <
[email protected]>
COPYRIGHT
Copyright 2014- Helmut Wollmersdorfer
LICENSE
This library is free software; you can redistribute it and/or modify it
under the same terms as Perl itself.
SEE ALSO