// includes #include using namespace std; // macros #define pb emplace_back #define mk make_pair #define FOR(i, a, b) for(int i=(a);i<(b);++i) #define rep(i, n) FOR(i, 0, n) #define rrep(i, n) for(int i=((int)(n)-1);i>=0;i--) #define irep(itr, st) for(auto itr = (st).begin(); itr != (st).end(); ++itr) #define irrep(itr, st) for(auto itr = (st).rbegin(); itr != (st).rend(); ++itr) #define all(x) (x).begin(),(x).end() #define sz(x) ((int)(x).size()) #define UNIQUE(v) v.erase(unique(v.begin(), v.end()), v.end()) #define bit(n) (1LL<<(n)) // functions template bool chmax(T &a, const T &b){if(a < b){a = b; return 1;} return 0;} template bool chmin(T &a, const T &b){if(a > b){a = b; return 1;} return 0;} template istream &operator>>(istream &is, vector &vec){for(auto &v: vec)is >> v; return is;} template ostream &operator<<(ostream &os, const vector& vec){for(int i = 0; i < vec.size(); i++){ os << vec[i]; if(i + 1 != vec.size())os << " ";} return os;} template ostream &operator<<(ostream &os, const set& st){for(auto itr = st.begin(); itr != st.end(); ++itr){ os << *itr; auto titr = itr; if(++titr != st.end())os << " ";} return os;} template ostream &operator<<(ostream &os, const unordered_set& st){for(auto itr = st.begin(); itr != st.end(); ++itr){ os << *itr; auto titr = itr; if(++titr != st.end())os << " ";} return os;} template ostream &operator<<(ostream &os, const multiset& st){for(auto itr = st.begin(); itr != st.end(); ++itr){ os << *itr; auto titr = itr; if(++titr != st.end())os << " ";} return os;} template ostream &operator<<(ostream &os, const unordered_multiset& st){for(auto itr = st.begin(); itr != st.end(); ++itr){ os << *itr; auto titr = itr; if(++titr != st.end())os << " ";} return os;} template ostream &operator<<(ostream &os, const pair &p){os << p.first << " " << p.second; return os;} template ostream &operator<<(ostream &os, const map &mp){for(auto itr = mp.begin(); itr != mp.end(); ++itr){ os << itr->first << ":" << itr->second; auto titr = itr; if(++titr != mp.end())os << " "; } return os;} template ostream &operator<<(ostream &os, const unordered_map &mp){for(auto itr = mp.begin(); itr != mp.end(); ++itr){ os << itr->first << ":" << itr->second; auto titr = itr; if(++titr != mp.end())os << " "; } return os;} // types using ll = long long int; using P = pair; // constants const int inf = 1e9; const ll linf = 1LL << 50; const double EPS = 1e-10; const int mod = 1000000007; const int dx[4] = {-1, 0, 1, 0}; const int dy[4] = {0, -1, 0, 1}; // io struct fast_io{ fast_io(){ios_base::sync_with_stdio(false); cin.tie(0); cout << fixed << setprecision(20);} } fast_io_; template T extgcd(T a, T b, T &x, T &y){ T d = a; if(b != 0){ d = extgcd(b, a % b, y, x); y -= (a / b) * x; }else{ x = 1, y = 0; } return d; } template T modinv(T a, T m){ long long x = 0, y = 0; extgcd(a, m, x, y); x %= m; if(x < 0)x += m; return x; } template long long garner(vector b, vector m, T MOD){ m.emplace_back(MOD); vector coef(m.size(), 1); vector consts(m.size(), 0); for(int i = 0; i < b.size(); i++){ long long t = ((b[i] - consts[i]) % m[i]) * modinv(coef[i], m[i]) % m[i]; if(t < 0)t += m[i]; for(int j = i + 1; j < m.size(); j++){ consts[j] = (consts[j] + t * coef[j] % m[j]) % m[j]; coef[j] = coef[j] * m[i] % m[j]; } } return consts.back(); } template T gcd(T a, T b) { if(a > b)return gcd(b, a); if(a == 0)return b; return gcd(b % a, a); } template long long arbitrary_garner(vector &b, vector &m, T MOD){ for(int i = 0; i < b.size(); i++){ for(int j = 0; j < i; j++){ T g = gcd(m[i], m[j]); if((b[i] - b[j]) % g != 0)return -1; m[i] /= g, m[j] /= g; T gi = gcd(m[i], g), gj = g / gi; do{ g = gcd(gi, gj); gi *= g, gj /= g; }while(g != 1); m[i] *= gi, m[j] *= gj; b[i] %= m[i], b[j] %= m[j]; } } return garner(b, m, MOD); } int main(int argc, char const* argv[]) { int n, k; cin >> n >> k; vector v(n); rep(i, n)v[i] = i; rep(i, k){ int x, y; cin >> x >> y; x--; y--; swap(v[x], v[y]); } vector idx(n, -1); vector> ids; int curr = 0; rep(i, n){ if(idx[i] >= 0)continue; int j = i; do{ idx[j] = curr; j = v[j]; }while(j != i); curr++; } rep(i, curr){ vector vec; rep(j, n){ if(idx[j] == i){ vec.pb(v[j]); } } ids.pb(vec); } int q; cin >> q; rep(i_, q){ vector s(n); cin >> s; rep(i, n)s[i]--; vector> qi; rep(j, curr){ vector vec; rep(k, n){ if(idx[k] == j){ vec.pb(s[k]); } } qi.pb(vec); } vector b, m; bool ng = false; rep(i, sz(ids)){ //cerr << ids[i] << endl; //cerr << qi[i] << endl; //cerr << "-------------" << endl; int B = -1; rep(j, sz(ids[i])){ bool flg = true; rep(k, sz(ids[i])){ if(ids[i][k] != qi[i][(k+j)%sz(ids[i])])flg = false; } if(flg){ B = j; break; } } if(B >= 0){ b.pb((B+1+sz(ids[i])) % sz(ids[i])); m.pb(sz(ids[i])); }else{ ng = true; break; } } if(!ng){ ll res = arbitrary_garner(b, m, mod); if(res == 0){ res = m[0]; FOR(i, 1, sz(m)){ res = res / gcd(res, m[i]) * m[i]; } } cout << res << endl; }else{ cout << -1 << endl; } } return 0; }