X-Git-Url: https://ruin.nu/git/?a=blobdiff_plain;f=src%2Fgenealgorithms.h;h=b5ed95d01ccb69af8797ec51aa820305e99db376;hb=148a31960aeb2ba127ef4322a12ae3af28ef2988;hp=4c519a7e3167a40384893a86dea243f20d1ac6d9;hpb=74e3465418ae60633025efb581e399308ffdf4b1;p=germs.git diff --git a/src/genealgorithms.h b/src/genealgorithms.h index 4c519a7..b5ed95d 100644 --- a/src/genealgorithms.h +++ b/src/genealgorithms.h @@ -25,6 +25,17 @@ 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; +}; +typedef std::pair Interval; + /** * Returns the length of the longest increasing sequence and the longest * decreasing sequence. @@ -35,5 +46,25 @@ 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); + +/** + * 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