From 9819b394e494bdfaaf55893a3f90b71ddb52227e Mon Sep 17 00:00:00 2001
From: Noobie99 <noob999noob999@gmail.com>
Date: Tue, 30 Jan 2024 22:01:40 +0100
Subject: add multiply function for fft

---
 math/transforms/multiply.cpp | 12 ++++++++++++
 1 file changed, 12 insertions(+)
 create mode 100644 math/transforms/multiply.cpp

diff --git a/math/transforms/multiply.cpp b/math/transforms/multiply.cpp
new file mode 100644
index 0000000..59c394e
--- /dev/null
+++ b/math/transforms/multiply.cpp
@@ -0,0 +1,12 @@
+vector<ll> mul(vector<ll>& a, vector<ll>& b) {
+	int n = 1 << (__lg(sz(a) + sz(b)) + 1);
+	vector<cplx> a2(all(a)), b2(all(b));
+	a2.resize(n), b2.resize(n);
+	fft(a2), fft(b2);
+	for (int i=0; i<n; i++) a2[i] *= b2[i];
+	fft(a2, true);
+
+	vector<ll> ans(n);
+	for (int i=0; i<n; i++) ans[i] = llround(a2[i].real());
+	return ans;
+}
-- 
cgit v1.2.3