#include using namespace std; #include #include #include #include #include template inline bool chmax(T& a, T b) { if (a <= b) { a = b; return 1; } return 0; } template inline bool chmin(T& a, T b) { if (a >= b) { a = b; return 1; } return 0; } #define rep(i,n) for (int i = 0; i < (n); ++i) typedef long long ll; using P=pair; const int INF=1001001001; //const int mod=998244353; // 負の数にも対応した mod // 例えば -17 を 5 で割った余りは本当は 3 (-17 ≡ 3 (mod. 5)) // しかし単に -17 % 5 では -2 になってしまう inline long long mod(long long a, long long m) { return (a % m + m) % m; } // 拡張 Euclid の互除法 // ap + bq = gcd(a, b) となる (p, q) を求め、d = gcd(a, b) をリターンします long long extGcd(long long a, long long b, long long &p, long long &q) { if (b == 0) { p = 1; q = 0; return a; } long long d = extGcd(b, a%b, q, p); q -= a/b * p; return d; } // 中国剰余定理 // リターン値を (r, m) とすると解は x ≡ r (mod. m) // 解なしの場合は (0, -1) をリターン pair ChineseRem(long long b1, long long m1, long long b2, long long m2) { long long p, q; long long d = extGcd(m1, m2, p, q); // p is inv of m1/d (mod. m2/d) if ((b2 - b1) % d != 0) return make_pair(0, -1); long long m = m1 * (m2/d); // lcm of (m1, m2) long long tmp = (b2 - b1) / d * p % (m2/d); long long r = mod(b1 + m1 * tmp, m); return make_pair(r, m); } void solve(){ int n,m; cin>>n>>m; vectora(n),b(m); rep(i,n){cin>>a[i];} rep(i,m){cin>>b[i];} ll mn=1e18; rep(i,n){ rep(j,m){ if(a[i]!=b[j]){continue;} auto t=ChineseRem(i,n,j,m); if(t.first==0&&t.second==-1){continue;} chmin(mn,t.first+1); } } if(mn==1e18){mn=-1;} cout<