diff --git a/src/day01.rs b/src/day01.rs
index ef5de76..0815a9b 100644
--- a/src/day01.rs
+++ b/src/day01.rs
@@ -4,24 +4,18 @@ use crate::utils;
 use anyhow::{anyhow, Result};
 
 fn day01_part1(data: &[i64]) -> Result<i64> {
-    for i in 0..data.len() - 2 {
-        for j in i + 1..data.len() - 1 {
-            if data[i] + data[j] == 2020 {
-                return Ok(data[i] * data[j]);
-            }
+    for c in utils::combinations(2, data.len()) {
+        if data[c[0]] + data[c[1]] == 2020 {
+            return Ok(data[c[0]] * data[c[1]]);
         }
     }
     Err(anyhow!("no result found"))
 }
 
 fn day01_part2(data: &[i64]) -> Result<i64> {
-    for i in 0..data.len() - 3 {
-        for j in i + 1..data.len() - 2 {
-            for k in j + 1..data.len() - 1 {
-                if data[i] + data[j] + data[k] == 2020 {
-                    return Ok(data[i] * data[j] * data[k]);
-                }
-            }
+    for c in utils::combinations(3, data.len()) {
+        if data[c[0]] + data[c[1]] + data[c[2]] == 2020 {
+            return Ok(data[c[0]] * data[c[1]] * data[c[2]]);
         }
     }
     Err(anyhow!("no result found"))
diff --git a/src/utils.rs b/src/utils.rs
index 5780853..3c20968 100644
--- a/src/utils.rs
+++ b/src/utils.rs
@@ -13,3 +13,65 @@ pub fn read_i64s(filename: &str) -> Result<Vec<i64>> {
 pub fn read_file(filename: &str) -> Result<String> {
     Ok(fs::read_to_string(filename)?)
 }
+
+// Lexicographic combinations.
+//
+// This algorithm generates all t-combinations of n numbers {0, 1, ..., n-1},
+// given n >= t >= 0. Additional variables c[t] and c[t+1] are used as
+// sentinels.
+//
+// Source: Donald E. Knuth, The Art of Computer Programming, VOLUME 4A,
+//         Combinatorial Algorithms, Part 1, Section 7.2.1.3, Algorithm L,
+//         Page 358, First printing, January 2011.
+pub struct CombinationsState {
+    t: usize,
+    c: Vec<usize>,
+    done: bool,
+}
+
+impl Iterator for CombinationsState {
+    type Item = Vec<usize>;
+
+    fn next(&mut self) -> Option<Vec<usize>> {
+        // L2. Visit.
+        let ret = if self.done {
+            None
+        } else {
+            Some((&self.c[0..self.c.len() - 2]).to_vec())
+        };
+
+        // L3. Find j.
+        let mut j = 0;
+        while (self.c[j] + 1) == self.c[j + 1] {
+            self.c[j] = j;
+            j = j + 1;
+        }
+
+        // L4. Done?
+        if j >= self.t {
+            self.done = true;
+        };
+
+        // L5. Increate c[j].
+        self.c[j] = self.c[j] + 1;
+
+        // return Some(visited) or None from L2 above
+        ret
+    }
+}
+
+pub fn combinations(t: usize, n: usize) -> CombinationsState {
+    let mut c = CombinationsState {
+        t: t,
+        c: vec![0; t + 2],
+        done: false,
+    };
+
+    // L1. Initialize
+    for i in 0..t {
+        c.c[i] = i;
+    }
+    c.c[t + 0] = n;
+    c.c[t + 1] = 0;
+    return c;
+}