From e1bd03f421cd8dc67d0ff7a0ad9d938be5cc78e2 Mon Sep 17 00:00:00 2001
From: Maxime Vorwerk <maxime@4-werk.com>
Date: Fri, 25 Aug 2023 15:43:41 +0200
Subject: [PATCH] finished powersubsets

---
 julia/Powerset.jl | 12 ++++++++++--
 1 file changed, 10 insertions(+), 2 deletions(-)

diff --git a/julia/Powerset.jl b/julia/Powerset.jl
index 132df21..490ba48 100644
--- a/julia/Powerset.jl
+++ b/julia/Powerset.jl
@@ -31,8 +31,16 @@ end
 function getPowerSubSet(setSize::Integer, subSetSize::Integer)
     PowerSet = getPowerSet(setSize)
     setSizes = dropdim(sum(Powerset, dims = 2), dims = 2)
-
-    return 
+    powerSubSetSize = binomial(setSize, subSetSize)
+    powerSubSets = Matrix{Int64}(undef, powerSubSetSize, setSize)
+    j = 1
+    for i ∈ 1:setSize
+        if setSizes[i] == subSetSize
+            powerSubSets[j, :] = PowerSet[i, :]
+            j += 1
+        end
+    end
+    return powerSubSets
 end
 
 end