Create randomLam.m2

function which inputs integers (n,k) and returns a random partition of size k with length at most n
Macaulay2 · Jun 2, 2024 · ab0d03d · ab0d03d
1 parent ec10638
commit ab0d03d
Showing 1 changed file with 14 additions and 0 deletions.
diff --git a/randomLam.m2 b/randomLam.m2
@@ -0,0 +1,14 @@
+randomLam = method();
+randomLam(ZZ, ZZ) := (n,k) -> (
+    L=new MutableList from {};
+    sumsofar = 0;
+    for i from 0 to (n-2) do(
+        x = random(0, k-sumsofar);
+        L#i=x;
+        sumsofar = sumsofar + x;
+        );
+    L#(n-1) = k - sumsofar;
+    L=toList L;
+    L=rsort L;
+    return delete(0,L)
+    )