diff --git a/Project.toml b/Project.toml
index c2e2a3e..f4a422a 100644
--- a/Project.toml
+++ b/Project.toml
@@ -30,6 +30,7 @@ Aqua = "4c88cf16-eb10-579e-8560-4a9242c79595"
 Documenter = "e30172f5-a6a5-5a46-863b-614d45cd2de4"
 Graphs = "86223c79-3864-5bf0-83f7-82e725a168b6"
 HiGHS = "87dc4568-4c63-4d18-b0c0-bb2238e4078b"
+IterTools = "c8e1da08-722c-5040-9ed9-7db0dc04731e"
 JET = "c3a54625-cd67-489e-a8e7-0a5a0ff4e31b"
 JuMP = "4076af6c-e467-56ae-b986-b466b2749572"
 JuliaFormatter = "98e50ef6-434e-11e9-1051-2b60c6c9e899"
@@ -38,4 +39,4 @@ SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [targets]
-test = ["Aqua", "Documenter", "Graphs", "HiGHS", "JET", "JuMP", "JuliaFormatter", "LinearAlgebra", "Test", "SparseArrays"]
+test = ["Aqua", "Documenter", "Graphs", "HiGHS", "JET", "IterTools", "JuMP", "JuliaFormatter", "LinearAlgebra", "SparseArrays", "Test"]
diff --git a/docs/src/algorithms.md b/docs/src/algorithms.md
index 5fd73ca..71f4815 100644
--- a/docs/src/algorithms.md
+++ b/docs/src/algorithms.md
@@ -57,6 +57,15 @@ GraphsOptim.min_cost_assignment!
 graph_matching
 ```
 
+## Coloring
+
+!!! danger "Work in progress"
+    Come back later!
+
+```@docs
+minimum_coloring
+```
+
 ## Utils
 
 ```@docs
diff --git a/src/GraphsOptim.jl b/src/GraphsOptim.jl
index 62e58d4..4806db1 100644
--- a/src/GraphsOptim.jl
+++ b/src/GraphsOptim.jl
@@ -7,6 +7,7 @@ module GraphsOptim
 
 using Graphs: AbstractGraph, is_directed
 using Graphs: vertices, edges, nv, ne, src, dst, inneighbors, outneighbors
+using Graphs: Coloring, maximal_cliques
 using FillArrays: Zeros, Fill
 using HiGHS: HiGHS
 using JuMP: Model
@@ -21,10 +22,12 @@ using OptimalTransport: sinkhorn
 export min_cost_flow
 export min_cost_assignment
 export FAQ, GOAT, graph_matching
+export minimum_coloring
 
 include("utils.jl")
 include("flow.jl")
 include("assignment.jl")
 include("graph_matching.jl")
+include("coloring.jl")
 
 end
diff --git a/src/coloring.jl b/src/coloring.jl
new file mode 100644
index 0000000..103329b
--- /dev/null
+++ b/src/coloring.jl
@@ -0,0 +1,41 @@
+"""
+    minimum_coloring(
+        g, max_nb_colors;
+        optimizer
+    )
+
+Finds a graph coloring using the smallest possible number of colors, assuming that it will not exceed `max_nb_colors`.
+
+Returns a vector of color indices.
+
+Beware: this is an NP-complete problem and so runtime can in the worst case increase exponentially in the size of the graph.
+"""
+function minimum_coloring(
+    g::AbstractGraph, max_nb_colors::Integer; optimizer=HiGHS.Optimizer
+)
+    model = Model(optimizer)
+
+    @variable(model, x[1:nv(g), 1:max_nb_colors], Bin)
+    @variable(model, y[1:max_nb_colors], Bin)
+
+    @objective(model, Min, sum(y))
+
+    for v in 1:nv(g)
+        @constraint(model, sum(x[v, :]) == 1)
+    end
+    for v in 1:nv(g), c in 1:max_nb_colors
+        @constraint(model, x[v, c] <= y[c])
+    end
+    for c in 1:max_nb_colors
+        for e in edges(g)
+            u, v = src(e), dst(e)
+            @constraint(model, x[u, c] + x[v, c] <= 1)
+        end
+    end
+
+    set_silent(model)
+    optimize!(model)
+    # return a vector of color indices
+    c = [argmax(value.(x[v, :])) for v in 1:nv(g)]
+    return c
+end
diff --git a/test/coloring.jl b/test/coloring.jl
new file mode 100644
index 0000000..735989d
--- /dev/null
+++ b/test/coloring.jl
@@ -0,0 +1,38 @@
+using Graphs
+using GraphsOptim
+using IterTools
+using Test
+
+function queens_graph(n::Integer, m::Integer)
+    g = SimpleGraph(n * m)
+    for (ix, iy, jx, jy) in product(1:n, 1:m, 1:n, 1:m)
+        dx = ix - jx
+        dy = iy - jy
+        if dx != 0 || dy != 0
+            if dx == 0 || dy == 0 || dx == dy || dx == -dy
+                add_edge!(g, (ix - 1) * m + iy, (jx - 1) * m + jy)
+            end
+        end
+    end
+    return g
+end
+
+queens_graph(n::Integer) = queens_graph(n, n)
+
+function kneser_graph(n::Integer, k::Integer)
+    ss = collect(subsets(1:n, k))
+    return SimpleGraph([isdisjoint(a, b) for a in ss, b in ss])
+end
+
+# https://oeis.org/A088202
+
+nb_colors(c) = length(unique(c))
+
+@test minimum_coloring(queens_graph(1), 10) |> nb_colors == 1
+@test minimum_coloring(queens_graph(2), 10) |> nb_colors == 4
+@test minimum_coloring(queens_graph(3), 10) |> nb_colors == 5
+@test minimum_coloring(queens_graph(4), 10) |> nb_colors == 5
+@test minimum_coloring(queens_graph(5), 10) |> nb_colors == 5
+@test minimum_coloring(queens_graph(6), 10) |> nb_colors == 7
+
+@test minimum_coloring(kneser_graph(11, 4), 10) |> nb_colors == 11 - 2 * 4 + 2
diff --git a/test/runtests.jl b/test/runtests.jl
index 5e19bc2..32a9a71 100644
--- a/test/runtests.jl
+++ b/test/runtests.jl
@@ -39,4 +39,8 @@ using Test
     @testset verbose = true "Graph matching" begin
         include("graph_matching.jl")
     end
+
+    @testset verbose = true "Coloring" begin
+        include("coloring.jl")
+    end
 end;