#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include // #include "utiltime.hpp" using namespace std; typedef long long int ll; typedef vector vi; typedef vector> vvi; typedef pair P; typedef pair Pll; typedef vector vll; typedef vector> vvll; typedef complex cdouble; const int INFL = (int)1e9; const ll INFLL = (ll)1e18; const double INFD = numeric_limits::infinity(); const double PI = 3.14159265358979323846; #define Loop(i, n) for(int i = 0; i < (int)n; i++) #define Loopll(i, n) for(ll i = 0; i < (ll)n; i++) #define Loop1(i, n) for(int i = 1; i <= (int)n; i++) #define Loopll1(i, n) for(ll i = 1; i <= (ll)n; i++) #define Loopr(i, n) for(int i = (int)n - 1; i >= 0; i--) #define Looprll(i, n) for(ll i = (ll)n - 1; i >= 0; i--) #define Loopr1(i, n) for(int i = (int)n; i >= 1; i--) #define Looprll1(i, n) for(ll i = (ll)n; i >= 1; i--) #define Loopitr(itr, container) for(auto itr = container.begin(); itr != container.end(); itr++) #define printv(vector) Loop(i, vector.size()) { cout << vector[i] << " "; } cout << endl; #define printmx(matrix) Loop(i, matrix.size()) { Loop(j, matrix[i].size()) { cout << matrix[i][j] << " "; } cout << endl; } #define quickio() ios::sync_with_stdio(false); cin.tie(0); #define bitmanip(m,val) static_cast>(val) ll rndf(double x) { return (ll)(x + (x >= 0 ? 0.5 : -0.5)); } ll floorsqrt(double x) { ll m = (ll)sqrt(x); return m + (m * m <= (ll)(x) ? 0 : -1); } ll ceilsqrt(double x) { ll m = (ll)sqrt(x); return m + ((ll)x <= m * m ? 0 : 1); } ll rnddiv(ll a, ll b) { return (a / b + (a % b * 2 >= b ? 1 : 0)); } ll ceildiv(ll a, ll b) { return (a / b + (a % b == 0 ? 0 : 1)); } ll gcd(ll m, ll n) { if (n == 0) return m; else return gcd(n, m % n); } /*******************************************************/ // a.size() should be 2^m vector FFT(vector a) { int n = int(a.size()); vector f(n); if (n == 1) f = a; else { vector a0(n / 2), a1(n / 2); Loop(i, n / 2) { a0[i] = a[i * 2]; a1[i] = a[i * 2 + 1]; } vector q = FFT(a0); vector r = FFT(a1); Loop(i, n / 2) { cdouble x = exp(cdouble({ 0, 2 * PI * i / n })); f[i] = q[i] + x * r[i]; f[i + n / 2] = q[i] - x * r[i]; } } return f; } // f.size() should be 2^m vector IFFT(vector f) { int n = int(f.size()); vector a(n); if (n == 1) a = f; else { vector q(n / 2), r(n / 2); Loop(i, n / 2) { cdouble x = exp(cdouble({ 0, 2 * PI * i / n })); q[i] = (f[i] + f[i + n / 2]) / 2.0; r[i] = (f[i] - f[i + n / 2]) / (x * 2.0); } vector a0 = IFFT(q); vector a1 = IFFT(r); Loop(i, n / 2) { a[i * 2] = a0[i]; a[i * 2 + 1] = a1[i]; } } return a; } int legal_size_of(int n) { int ret = 1 << (int)log2(n); if (ret < n) ret <<= 1; return ret; } int main() { quickio(); int L, M, N; cin >> L >> M >> N; int n = legal_size_of(N * 2 - 1); vector a(n, 0), b(n, 0); Loop(i, L) { int abuf; cin >> abuf; abuf--; a[abuf] = 1; } Loop(i, M) { int bbuf; cin >> bbuf; bbuf--; if (bbuf == 0) b[0] = 1; else b[n - bbuf] = 1; } vector g = FFT(a); vector h = FFT(b); vector gh(n); Loop(i, gh.size()) { gh[i] = g[i] * h[i]; } vector c = IFFT(gh); int q; cin >> q; Loop(i, q) { cout << rndf(c[i].real()) << endl; } }