diff --git a/sum-with-multiplicity b/sum-with-multiplicity
new file mode 100644
index 0000000..d41813d
--- /dev/null
+++ b/sum-with-multiplicity
@@ -0,0 +1,25 @@
+# Time:  O(n^2), n is the number of disctinct A[i]
+# Space: O(n)
+
+import collections
+import itertools
+
+
+class Solution(object):
+    def threeSumMulti(self, A, target):
+        """
+        :type A: List[int]
+        :type target: int
+        :rtype: int
+        """
+        count = collections.Counter(A)
+        result = 0
+        for i, j in itertools.combinations_with_replacement(count, 2):
+            k = target - i - j
+            if i == j == k:
+                result += count[i] * (count[i]-1) * (count[i]-2) // 6
+            elif i == j != k:
+                result += count[i] * (count[i]-1) // 2 * count[k]
+            elif max(i, j) < k:
+                result += count[i] * count[j] * count[k]
+        return result % (10**9 + 7)