X-Git-Url: https://ruin.nu/git/?a=blobdiff_plain;f=src%2Fgenealgorithms.h;h=a91f67ddcae2345497c24841ff93dc6d6d5ebf82;hb=4462cded18dd80075c3622e5a87ca6f5f3c84f62;hp=b0e75b9f65790e184eabe46b6bf9ec9dba8d1002;hpb=65f449adad91a229757c0317c27ad9fb87a4d222;p=germs.git diff --git a/src/genealgorithms.h b/src/genealgorithms.h index b0e75b9..a91f67d 100644 --- a/src/genealgorithms.h +++ b/src/genealgorithms.h @@ -24,14 +24,21 @@ #include class GeneOrder; +struct Component; -struct Component{ - Component(int b,int e,int s):begin(b),end(e),sign(s){} - int begin; - int end; - int sign; +/** + * Simple interval struct, holds the indexes + information if the interval is + * oriented or not. + */ +struct Interval{ + Interval(size_t f,size_t s,bool o = false):first(f),second(s),oriented(o){} + bool operator==(const Interval& i){ + return first == i.first && second == i.second && oriented == i.oriented; + } + size_t first; + size_t second; + bool oriented; }; -typedef std::pair Interval; /** * Returns the length of the longest increasing sequence and the longest @@ -40,13 +47,20 @@ typedef std::pair Interval; std::pair longestSequences(const GeneOrder& go); /** + * Uses the Robinson-Schensted algorithm on the gene order. Gives the longest + * increasing and decreasing subsequences. */ std::vector > robinsonSchensted(const GeneOrder& go); /** * Counts the number of cycles in the gene order. */ -int countCycles(const GeneOrder& go); +size_t countCycles(const GeneOrder& go); + +/** + * Calculates the inversion distance for this gene order + */ +size_t inversionDistance(const GeneOrder& go); /** * Finds the components in the gene order.