diff --git a/cpp/LinearFactor.h b/cpp/LinearFactor.h
index ff4f7a55e..cd1ecbc49 100644
--- a/cpp/LinearFactor.h
+++ b/cpp/LinearFactor.h
@@ -14,6 +14,7 @@
 #include <set>
 #include <ostream>
 #include <boost/shared_ptr.hpp>
+#include <vector>
 
 #include "Matrix.h"
 #include "Factor.h"
@@ -82,6 +83,15 @@ public:
 		As.insert(make_pair(key3, A3));
 	}
 
+	/** Construct an n-ary factor */
+	CONSTRUCTOR
+	LinearFactor(const std::vector<std::pair<std::string, Matrix> > &terms,
+	    const Vector &b_in) :
+	    b(b_in) {
+	  for(unsigned int i=0; i<terms.size(); i++)
+	    As.insert(terms[i]);
+	}
+
 	/** Construct from Conditional Gaussian */
 	CONSTRUCTOR
 	LinearFactor(const std::string& key, const boost::shared_ptr<