Submission #15064

#TimeUsernameProblemLanguageResultExecution timeMemory
15064hongjun7씽크스몰 (kriii3_TT)C++98
20 / 30
1233 ms43924 KiB
#include <cstdio> #define _USE_MATH_DEFINES #include <math.h> #include <complex> #include <vector> #include <algorithm> using namespace std; #define sz(v) ((int)(v).size()) #define all(v) (v).begin(),(v).end() typedef complex<long double> base; typedef long long ll; void fft(vector <base> &a, bool invert) { int n = sz(a); for (int i = 1, j = 0; i<n; i++) { int bit = n >> 1; for (; j >= bit; bit >>= 1) j -= bit; j += bit; if (i < j) swap(a[i], a[j]); } for (int len = 2; len <= n; len <<= 1) { long double ang = 2 * M_PI / len*(invert ? -1 : 1); base wlen(cos(ang), sin(ang)); for (int i = 0; i<n; i += len) { base w(1); for (int j = 0; j<len / 2; j++) { base u = a[i + j], v = a[i + j + len / 2] * w; a[i + j] = u + v; a[i + j + len / 2] = u - v; w *= wlen; } } } if (invert) { for (int i = 0; i<n; i++) a[i] /= n; } } void multiply(const vector<ll> &a, const vector<ll> &b, vector<ll> &res) { vector <base> fa(all(a)), fb(all(b)); int n = 1; while (n < sz(a) + sz(b)) n <<= 1; fa.resize(n); fb.resize(n); fft(fa, 0); fft(fb, 0); for (int i = 0; i < n; i++) fa[i] *= fb[i]; fft(fa, 1); res.resize(n); for (int i = 0; i < n; i++) res[i] = (ll)(fa[i].real() + (fa[i].real()>0 ? 0.5 : -0.5)); ll r = res[0]; for (int i = 1; i < n; i++) r ^= res[i]; printf("%lld", r); } int n, m, x; int main() { scanf("%d%d", &n, &m); vector <ll> a(n + 1), b(m + 1), res; for (int i = 0; i <= n; i++) scanf("%lld", &a[i]); for (int i = 0; i <= m; i++) scanf("%lld", &b[i]); if (n < m) for (int i = n; i <= m; i++) a.push_back(0); else for (int i = m; i <= n; i++) b.push_back(0); multiply(a, b, res); return 0; }
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...