nicer with size_t for cycles, they can't be negative anyway

[germs.git] / src / genealgorithms.h

diff --git a/src/genealgorithms.h b/src/genealgorithms.h
index 4c519a7e3167a40384893a86dea243f20d1ac6d9..b5ed95d01ccb69af8797ec51aa820305e99db376 100644 (file)
--- 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<size_t,size_t> Interval;
+
 /**
  * Returns the length of the longest increasing sequence and the longest
  * decreasing sequence.
@@ -35,5 +46,25 @@ std::pair<int,int> longestSequences(const GeneOrder& go);
  */
 std::vector<std::vector<int> > 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<Component> findComponents(const GeneOrder& go);
+
+/**
+ * Find intervals.
+ */
+std::vector<Interval> findIntervals(const GeneOrder& go);
+
+/**
+ * Creates a list with the intervals at each point.
+ */
+std::vector<Interval> findIntervalsAtPoints(const std::vector<Interval>& intervals);
+
 #endif