diff --git a/packages/random/src/index.ts b/packages/random/src/index.ts
index d9752cb2e5..140ce47e88 100644
--- a/packages/random/src/index.ts
+++ b/packages/random/src/index.ts
@@ -12,6 +12,7 @@ export * from "./xsadd";
 export * from "./coin";
 export * from "./random-bytes";
 export * from "./random-id";
+export * from "./unique-indices";
 export * from "./uuid";
 export * from "./weighted-random";
 
diff --git a/packages/random/src/unique-indices.ts b/packages/random/src/unique-indices.ts
new file mode 100644
index 0000000000..59572586d3
--- /dev/null
+++ b/packages/random/src/unique-indices.ts
@@ -0,0 +1,61 @@
+import { assert, Fn0 } from "@thi.ng/api";
+import type { IRandom } from "./api";
+import { SYSTEM } from "./system";
+
+/**
+ * Attempts to draw `k` unique values from given zero-arg function `fn`
+ * (presumably a PRNG of sorts) and adds them to `existing` array of unique
+ * samples (or creates a new one). Returns the array. Gives up after
+ * `maxTrials`.
+ *
+ * @param k
+ * @param fn
+ * @param existing
+ * @param maxTrials
+ */
+export const uniqueValuesFrom = (
+    k: number,
+    fn: Fn0<number>,
+    existing: number[] = [],
+    maxTrials = 100
+) => {
+    let n = 0;
+    while (n < k) {
+        let i: number;
+        let trials = maxTrials;
+        do {
+            i = fn();
+        } while (existing.includes(i) && --trials > 0);
+        if (trials <= 0) break;
+        existing.push(i);
+        n++;
+    }
+    return existing;
+};
+
+/**
+ * Similar to (and based o) {@link uniqueValuesFrom}. Attempts to add `k` unique
+ * integer indices in the `[0, max)` interval to the (optional) array of
+ * pre-`existing` indices (which will never be picked again and new indices will
+ * be added to). Returns updated array.
+ *
+ * @remarks
+ * Candidates are drawn from the provided `rnd` {@link IRandom} (default:
+ * {@link SYSTEM}) and only `maxTrials` are attempted before giving up.
+ *
+ * @param k
+ * @param max
+ * @param existing
+ * @param maxTrials
+ * @param rnd
+ */
+export const uniqueIndices = (
+    k: number,
+    max: number,
+    existing?: number[],
+    maxTrials?: number,
+    rnd: IRandom = SYSTEM
+) => {
+    assert(k >= 0 && k <= max, `k must be in [0, ${max}] interval`);
+    return uniqueValuesFrom(k, () => rnd.int() % max, existing, maxTrials);
+};