#include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using ll = long long; class range {private: struct I{int x;int operator*(){return x;}bool operator!=(I& lhs){return x ostream& operator<<(ostream& os, const pair& p){ return os << "{" << p.first << ", " << p.second << "}"; } template ostream& operator<<(ostream& os, const vector& obj) { os << "{"; for (const auto& e : obj) os << e << ", "; return os << "}"; } template ostream& operator<<(ostream& os, const set& obj) { os << "set{"; for (const auto& e : obj) os << e << ", "; return os << "}"; } template ostream& operator<<(ostream& os, const map& obj) { os << "map{"; for (const auto& e : obj) os << e << ", "; return os << "}"; } template void take(vector& vec, int n) { vec.resize(n); for (int i = 0; i < n; ++i) cin >> vec[i]; } #ifdef ONLINE_JUDGE #define dump(expr) ; #else #define dump(expr) { cerr << "\033[33m#L" << __LINE__ << ": " << expr << "\033[39m" << endl; } #endif struct UnionFindWithWeights { vector data, weights; UnionFindWithWeights(int n) : data(n, -1), weights(n, 1) {} int root(int x) { return data[x] < 0 ? x : data[x] = root(data[x]); } bool find(int x, int y) { return root(x) == root(y); } void uni(int x, int y) { x = root(x); y = root(y); if (x != y) { int new_weight = weights[x] + weights[y]; if (data[x] < data[y]) swap(x, y); data[x] += data[y]; data[y] = x; weights[x] = new_weight; } } int weight(int x) { return weights[root(x)]; } }; namespace solver { template struct In2 { T1 a; T2 b; friend std::istream& operator>>(std::istream& is, In2& obj) { T1 t1; T2 t2; is >> t1 >> t2; obj = {t1, t2}; return is; } }; template struct In3 { T1 a; T2 b; T3 c; friend std::istream& operator>>(std::istream& is, In3& obj) { T1 t1; T2 t2; T3 t3; is >> t1 >> t2 >> t3; obj = {t1, t2, t3}; return is; } }; int n; vector vd, vw; void read() { cin >> n; take(vd, n); take(vw, n); } using RetType = string; RetType run() { UnionFindWithWeights uf(n); vector selfloop; for (int i : range(n)) { int j = (i + vd[i]) % n; int k = ((i - vd[i]) % n + n) % n; if (j == k) selfloop.push_back(j); else uf.uni(j, k); } vector parity(n); vector ignore(n); for (int i : range(n)) { parity[uf.root(i)] ^= vw[i]; } for (int i : selfloop) { ignore[uf.root(i)] = 1; } for (int i : range(n)) { int r = uf.root(i); if (ignore[uf.root(i)]) continue; if (parity[r] != (uf.weight(r) & 1)) return "No"; } return "Yes"; } } // namespace template void run(F f) { if constexpr (std::is_same_v) f(); else cout << f() << endl; } int main(int argc, char** argv) { cerr << fixed << setprecision(12); cout << fixed << setprecision(12); int testcase = 1; if (argc > 1) testcase = atoi(argv[1]); while (testcase--) { solver::read(); } run(solver::run); }