#ifndef ONLINE_JUDGE #define _GLIBCXX_DEBUG #endif #include using namespace std; #define pass (void)0 #define INF (1<<30)-1 #define INFLL (1LL<<60)-1 #define rep(i, n) for (int i = 0; i < (int)(n); i++) #define repr(i, n) for (int i = (int)(n) - 1; i >= 0; i--) #define rep2(i, a, b) for (int i = (int)(a); i < (int)(b); i++) #define repr2(i, a, b) for (int i = (int)(b) - 1; i >= (int)(a); i--) #define all(x) (x).begin(), (x).end() #define rall(x) (x).rbegin(), (x).rend() #define sz(x) ((int)(x).size()) #define YesNo(cond) cout << ((cond) ? "Yes\n" : "No\n") #define YESNO(cond) cout << ((cond) ? "YES\n" : "NO\n") using ll = long long; using pii = pair; using pll = pair; using vi = vector; using vl = vector; using vvi = vector; using vvl = vector; template void print(const T& value) { cout << value << "\n"; } template void print(const vector& vec) { for (auto& v : vec) cout << v << " "; cout << "\n"; } template void input(vector& vec) { for (auto& v : vec) cin >> v; }; template bool chmin(T& a, const T& b) { if (a > b) { a = b; return true; } return false; } template bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } return false; } #include using namespace atcoder; using mint = modint998244353; using vm = vector; struct Eratosthenes { int n; vector isPrime; vector minfactor; Eratosthenes(int n) : n(n) { isPrime.assign(n + 1, true); minfactor.assign(n + 1, -1); if (n >= 0) isPrime[0] = false; if (n >= 1) { isPrime[1] = false; minfactor[1] = 1; } for (int i = 2; i <= n; i++) { if (isPrime[i]) { minfactor[i] = i; for (long long j = 1LL * i * 2; j <= n; j += i) { isPrime[j] = false; if (minfactor[j] == -1) minfactor[j] = i; } } } } vector> factorize(long long x) const { vector> res; while (x > 1) { int p = minfactor[x]; int cnt = 0; while (minfactor[x] == p) { x /= p; cnt++; } res.emplace_back(p, cnt); } return res; } vector divisor(long long x) const { vector ans = {1}; auto pf = factorize(x); for (auto [p, c] : pf) { int L = ans.size(); for (int i = 0; i < L; i++) { long long v = 1; for (int t = 0; t < c; t++) { v *= p; ans.push_back(ans[i] * v); } } } sort(all(ans)); return ans; } }; vector> C; ll N; ll func(ll limit) { ll ans = 0; rep (i, N) { ll left = -1; ll right = sz(C[i]); while (left+1 < right) { ll mid = (left+right)/2; if (C[i][mid].first < limit) { left = mid; } else { right = mid; } } ans += C[i][right].second; } return ans; } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(10); ll K; cin >> N >> K; vl A(N); vl B(N); input(A); input(B); Eratosthenes E(100000); rep (i, N) { ll GCD = gcd(A[i], B[i]); vl div = E.divisor(A[i]); vector stack; stack.push_back({GCD, 0}); for (auto d : div) { if (d <= GCD) continue; ll cnt = (((d-B[i])%d+d)%d); while (!stack.empty() && cnt <= stack[sz(stack)-1].second) { stack.pop_back(); } stack.push_back({d, cnt}); } C.push_back(stack); } rep (i, N) sort(all(C[i])); ll left = 1; ll right = INFLL; for (auto a : A) chmin(right, a); right ++; while (left+1 < right) { ll mid = (left+right)/2; if (func(mid) <= K) { left = mid; } else { right = mid; } } print(left); }