#line 1 "combined.cpp" #pragma region Macros #include using namespace std; template inline bool chmax(T &a, T b) { if(a < b) { a = b; return 1; } return 0; } template inline bool chmin(T &a, T b) { if(a > b) { a = b; return 1; } return 0; } #ifdef DEBUG template ostream &operator<<(ostream &os, const pair &p) { os << '(' << p.first << ',' << p.second << ')'; return os; } template ostream &operator<<(ostream &os, const vector &v) { os << '{'; for(int i = 0; i < (int)v.size(); i++) { if(i) { os << ','; } os << v[i]; } os << '}'; return os; } void debugg() { cerr << endl; } template void debugg(const T &x, const Args &... args) { cerr << " " << x; debugg(args...); } #define debug(...) \ cerr << __LINE__ << " [" << #__VA_ARGS__ << "]: ", debugg(__VA_ARGS__) #define dump(x) cerr << __LINE__ << " " << #x << " = " << (x) << endl #else #define debug(...) (void(0)) #define dump(x) (void(0)) #endif struct Setup { Setup() { cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(15); } } __Setup; using ll = long long; #define OVERLOAD3(_1, _2, _3, name, ...) name #define ALL(v) (v).begin(), (v).end() #define RALL(v) (v).rbegin(), (v).rend() #define REP1(i, n) for(int i = 0; i < int(n); i++) #define REP2(i, a, b) for(int i = (a); i < int(b); i++) #define REP(...) OVERLOAD3(__VA_ARGS__, REP2, REP1)(__VA_ARGS__) #define UNIQUE(v) sort(ALL(v)), (v).erase(unique(ALL(v)), (v).end()) #define REVERSE(v) reverse(ALL(v)) #define SZ(v) ((int)(v).size()) const int INF = 1 << 30; const ll LLINF = 1LL << 60; constexpr int MOD = 1000000007; constexpr int MOD2 = 998244353; const int dx[4] = {1, 0, -1, 0}; const int dy[4] = {0, 1, 0, -1}; void Case(int i) { cout << "Case #" << i << ": "; } int popcount(int x) { return __builtin_popcount(x); } ll popcount(ll x) { return __builtin_popcountll(x); } #pragma endregion Macros #line 2 "/Users/siro53/kyo-pro/compro_library/math/convolution/convolution.hpp" #line 7 "/Users/siro53/kyo-pro/compro_library/math/convolution/convolution.hpp" namespace fft { template struct Complex { D x, y; Complex(): x(0), y(0) {}; Complex(D x, D y) : x(x), y(y) {} Complex &operator+=(const Complex &c) { x += c.x; y += c.y; return (*this); } Complex &operator-=(const Complex &c) { x -= c.x; y -= c.y; return (*this); } Complex &operator*=(const Complex &c) { D nx = x * c.x - y * c.y; D ny = x * c.y + y * c.x; x = nx, y = ny; return (*this); } Complex &operator/=(const Complex& c) { // (a + bi) / (c + di) // (a + bi)(c - di) / (c^2 + d^2) // (ac + bd) + (bc - ad)i / (c^2 + d^2) D nx = (x * c.x + y * c.y) / (c.x * c.x + c.y * c.y); D ny = (y * c.x - x * c.y) / (c.x * c.x + c.y * c.y); x = nx, y = ny; return (*this); } Complex operator-() const { return Complex(-x, -y); } Complex operator+(const Complex &c) const { return Complex(*this) += c; } Complex operator-(const Complex &c) const { return Complex(*this) -= c; } Complex operator*(const Complex &c) const { return Complex(*this) *= c; } Complex operator/(const Complex &c) const { return Complex(*this) /= c; } }; template constexpr const D PI = std::acos(D(-1)); template inline Complex omega(int k, int n) { return Complex(std::cos(D(k) * 2 * PI / n), std::sin(D(k) * 2 * PI / n)); } inline int revbit(int mask, int bitlen) { int res = 0; while(bitlen--) { res = (res << 1) | (mask & 1); mask >>= 1; } return res; } template void fft(std::vector>& a, int bitlen) { int n = (int)a.size(); int len = n; while(len > 1) { for(int i = 0; i < n; i += len) { int t = len >> 1; for(int j = 0; j < t; j++) { int p = i + j; auto l = a[p]; auto r = a[p + t]; a[p] = l + r; a[p + t] = (l - r) * omega(j, len); } } len >>= 1; } for(int i = 0; i < n; i++) { int j = revbit(i, bitlen); if(i < j) std::swap(a[i], a[j]); } } template void ifft(std::vector>& a, int bitlen) { int n = (int)a.size(); for(int i = 0; i < n; i++) { int j = revbit(i, bitlen); if(i < j) std::swap(a[i], a[j]); } int len = 2; while(len <= n) { for(int i = 0; i < n; i += len) { int t = len >> 1; for(int j = 0; j < t; j++) { int p = i + j; auto l = a[p]; auto r = a[p + t] * omega(-j, len); a[p] = l + r; a[p + t] = l - r; } } len <<= 1; } for(int i = 0; i < n; i++) a[i] /= Complex(n, 0); } template std::vector convolution(const std::vector& a, const std::vector& b) { int m = (int)a.size() + (int)b.size() - 1; int n = 1, bitlen = 0; while(n < m) { n <<= 1; bitlen++; } std::vector> A(n), B(n); for(int i = 0; i < (int)a.size(); i++) A[i] = Complex(a[i], 0); for(int i = 0; i < (int)b.size(); i++) B[i] = Complex(b[i], 0); fft(A, bitlen); fft(B, bitlen); for(int i = 0; i < n; i++) A[i] *= B[i]; ifft(A, bitlen); std::vector res(m); for(int i = 0; i < m; i++) res[i] = A[i].x; return res; } }; // namespace fft #line 78 "combined.cpp" int main() { int L, M, N; cin >> L >> M >> N; vector fa(N+1, 0), fb(N+1, 0); REP(i, L) { int a; cin >> a; fa[a] += 1; } REP(i, M) { int b; cin >> b; fb[N - b] += 1; } auto c = fft::convolution(fa, fb); int Q; cin >> Q; REP(v, Q) cout << (ll)c[N + v] << '\n'; }