X-Git-Url: https://ruin.nu/git/?a=blobdiff_plain;f=src%2Fgenealgorithms.h;h=6bc133fa96174478746cd828659a476e02ea52ec;hb=cc3e3bd109d72d82826629b96f8e6721d5243cf7;hp=4c519a7e3167a40384893a86dea243f20d1ac6d9;hpb=74e3465418ae60633025efb581e399308ffdf4b1;p=germs.git diff --git a/src/genealgorithms.h b/src/genealgorithms.h index 4c519a7..6bc133f 100644 --- a/src/genealgorithms.h +++ b/src/genealgorithms.h @@ -25,6 +25,26 @@ class GeneOrder; +struct Component{ + Component(int b,int e,int s):begin(b),end(e),sign(s){} + bool operator==(const Component& c){ + return begin == c.begin && end == c.end && sign == c.sign; + } + int begin; + int end; + int sign; +}; + +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; +}; + /** * Returns the length of the longest increasing sequence and the longest * decreasing sequence. @@ -35,5 +55,30 @@ std::pair longestSequences(const GeneOrder& go); */ std::vector > robinsonSchensted(const GeneOrder& go); +/** + * Counts the number of cycles in the gene order. + */ +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. + */ +std::vector findComponents(const GeneOrder& go); + +/** + * Find intervals. + */ +std::vector findIntervals(const GeneOrder& go); + +/** + * Creates a list with the intervals at each point. + */ +std::vector findIntervalsAtPoints(const std::vector& intervals); + #endif