#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using ll = long long; constexpr int INF = 1001001001; constexpr int mod = 1000000007; // constexpr int mod = 998244353; template inline bool chmax(T& x, T y){ if(x < y){ x = y; return true; } return false; } template inline bool chmin(T& x, T y){ if(x > y){ x = y; return true; } return false; } const double pi = acos(-1); struct Complex { double real, imag; Complex(double r = 0, double i = 0) : real(r), imag(i) {} inline Complex& operator=(const Complex& c){ real = c.real; imag = c.imag; return *this; } }; inline Complex operator+(const Complex& x, const Complex& y){ return {x.real + y.real, x.imag + y.imag}; } inline Complex operator-(const Complex& x, const Complex& y){ return {x.real - y.real, x.imag - y.real}; } inline Complex operator*(const Complex& x, const Complex& y){ return {x.real * y.real - x.imag * y.imag, x.real * y.imag + x.imag * y.real}; } inline Complex operator*(const Complex& x, double a){ return {x.real * a, x.imag * a}; } inline Complex operator/(const Complex& x, double a){ return {x.real / a, x.imag / a}; } vector dft(vector& a, int n, int sgn = 1){ if(n == 1) return a; int m = n / 2; vector f(m), g(m); for(int i = 0; i < m; ++i){ f[i] = a[2 * i]; g[i] = a[2 * i + 1]; } // divide and conquer f = dft(f, m, sgn); g = dft(g, m, sgn); // zeta : nth root of 1 Complex zeta(cos(2.0 * pi / n), sin(2.0 * pi / n) * sgn); Complex x = 1; for(int i = 0; i < m; ++i){ a[i] = f[i] + x * g[i]; x = x * zeta; } for(int i = m; i < n; ++i){ a[i] = f[i - m] + x * g[i - m]; x = x * zeta; } return a; } vector inv_dft(vector& a, int n){ a = dft(a, n, -1); for(int i = 0; i < n; ++i) a[i] = a[i] / n; return a; } vector multiply(vector& a, vector& b){ int sz = a.size() + b.size() - 1; int n = 1; while(n < sz) n <<= 1; a.resize(n); b.resize(n); a = dft(a, n); b = dft(b, n); vector f(n); for(int i = 0; i < n; ++i){ f[i] = a[i] * b[i]; } f = inv_dft(f, n); vector res(n); for(int i = 0; i < n; ++i){ res[i] = llround(f[i].real); } return res; } int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); int L, M, N, A, B, Q; cin >> L >> M >> N; vector a(N), b(N); for(int i = 0; i < L; ++i){ cin >> A; a[A - 1] = Complex(1, 0); } for(int i = 0; i < M; ++i){ cin >> B; b[N - B] = Complex(1, 0); } auto c = multiply(a, b); cin >> Q; for(int i = 0; i < Q; ++i){ cout << c[N - 1 + i] << '\n'; } return 0; }