1 files changed, 58 insertions, 0 deletions
diff --git a/STL/random.cpp b/STL/random.cpp
new file mode 100644
index 00000000000..57ef645ae64
--- /dev/null
+++ b/STL/random.cpp
@@ -0,0 +1,58 @@
+/*
+ *
+ * Copyright (c) 1994
+ * Hewlett-Packard Company
+ *
+ * Permission to use, copy, modify, distribute and sell this software
+ * and its documentation for any purpose is hereby granted without fee,
+ * provided that the above copyright notice appear in all copies and
+ * that both that copyright notice and this permission notice appear
+ * in supporting documentation.  Hewlett-Packard Company and Microsoft
+ * Corporation make no representations about the suitability of this
+ * software for any purpose.  It is provided "as is" without express or
+ * implied warranty.
+ *
+ */
+
+#include <stddef.h>
+
+#define __SEED 161803398
+
+class __random_generator {
+protected:
+    unsigned long table[55];
+    size_t index1;
+    size_t index2;
+public:
+    unsigned long operator()(unsigned long limit) {
+	index1 = (index1 + 1) % 55;
+	index2 = (index2 + 1) % 55;
+	table[index1] = table[index1] - table[index2];
+	return table[index1] % limit;
+    }
+    void seed(unsigned long j);
+    __random_generator(unsigned long j) { seed(j); }
+};
+
+void __random_generator::seed(unsigned long j) {
+    unsigned long k = 1;
+    table[54] = j;
+    for (size_t i = 0; i < 54; i++) {
+        size_t ii = 21 * i % 55;
+        table[ii] = k;
+        k = j - k;
+        j = table[ii];
+    }
+    for (int loop = 0; loop < 4; loop++) {
+        for (i = 0; i < 55; i++)
+            table[i] = table[i] - table[(1 + i + 30) % 55];
+    }
+    index1 = 0;
+    index2 = 31;
+}
+
+__random_generator rd(__SEED);
+
+unsigned long __long_random(unsigned long limit) {
+    return rd(limit);
+}