#include using namespace std; #define rep2(i, m, n) for(int i=int(m); i struct ModInt { int x; ModInt(int _x = 0) : x( (0 <= _x && _x < MOD) ? _x : (_x%MOD+MOD)%MOD ) {} constexpr ModInt operator-() const noexcept { return ModInt(MOD-x); } constexpr ModInt operator+(const ModInt y) const noexcept { return ModInt(*this) += y; } constexpr ModInt operator-(const ModInt y) const noexcept { return ModInt(*this) -= y; } constexpr ModInt operator*(const ModInt y) const noexcept { return ModInt(*this) *= y; } constexpr ModInt &operator+=(const ModInt y) noexcept { x += y.x; if (x >= MOD) x -= MOD; return *this; } constexpr ModInt &operator-=(const ModInt y) noexcept { x -= y.x; if (x < 0) x += MOD; return *this; } constexpr ModInt &operator*=(const ModInt y) noexcept { x = ll(x) * y.x % MOD; return *this; } }; using mint = ModInt; using Vm = vector; using VVm = vector; istream &operator>>(istream &is, mint &a) { return is >> a.x; } ostream &operator<<(ostream &os, const mint &a) { return os << a.x; } int main() { int k, n, m; cin >> k >> n >> m; Vm a(k+n), c(k+1); c[0] = -1; rep(i, k) cin >> a[i]; rep(i, k) cin >> c[i+1]; rep2(i, k, n) rep2(j, 1, k+1) a[i] += c[j]*a[i-j]; VVm d(n, Vm(k)); rep(i, n) { d[i][0] = c[0]*a[i]; rep2(j, 1, k) d[i][j] += d[i][j-1] + c[j]*a[i-j]; } Vm y(n); rep(_, m) { int l, r; cin >> l >> r; rep2(i, l, min(l+k, n)) y[i] += d[i-l][i-l]; rep2(i, r, min(r+k, n)) y[i] -= d[i-l][i-r]; } Vm x(n); rep(i, n) { x[i] = - y[i]; rep2(j, 1, min(i, k)+1) x[i] += c[j]*x[i-j]; cout << x[i] << '\n'; } return 0; }