#include using namespace std; using ll = long long; using ld = long double; using pll = pair; using vl = vector; template using vec = vector; template using vv = vec>; template using vvv = vv>; template using minpq = priority_queue, greater>; #define all(a) (a).begin(),(a).end() #define rep(i, n) for (ll i = 0; i < (n); ++i) #define reps(i, l, r) for(ll i = (l); i < (r); ++i) #define rrep(i, l, r) for(ll i = (r)-1; i >= (l); --i) #define sz(x) (ll) (x).size() template bool chmax(T &a, const T& b) { return a < b ? a = b, true : false; } template bool chmin(T &a, const T& b) { return a > b ? a = b, true : false; } struct Edge { ll from, to, cost; Edge (ll from, ll to, ll cost = 1ll) : from(from), to(to), cost(cost) {} }; struct Graph { vector> G; Graph() = default; explicit Graph(ll N) : G(N) {} size_t size() const { return G.size(); } void add(ll from, ll to, ll cost = 1ll, bool direct = 0) { G[from].emplace_back(from, to, cost); if (!direct) G[to].emplace_back(to, from, cost); } vector &operator[](const int &k) { return G[k]; } }; using Edges = vector; const ll mod = 998244353; // 1000000007; struct mint { ll x; mint(ll y = 0) : x(y >= 0 ? y % mod : (mod - (-y) % mod) % mod) {} mint &operator+=(const mint &p) { if ((x += p.x) >= mod) x -= mod; return *this; } mint &operator-=(const mint &p) { if ((x += mod - p.x) >= mod) x -= mod; return *this; } mint &operator*=(const mint &p) { x = (ll)(1ll * x * p.x % mod); return *this; } mint &operator/=(const mint &p) { *this *= p.inv(); return *this; } mint operator-() const { return mint(-x); } mint operator+(const mint &p) const { return mint(*this) += p; } mint operator-(const mint &p) const { return mint(*this) -= p; } mint operator*(const mint &p) const { return mint(*this) *= p; } mint operator/(const mint &p) const { return mint(*this) /= p; } bool operator==(const mint &p) const { return x == p.x; } bool operator!=(const mint &p) const { return x != p.x; } friend ostream &operator<<(ostream &os, const mint &p) { return os << p.x; } friend istream &operator>>(istream &is, mint &a) { ll t; is >> t; a = mint(t); return (is); } mint inv() const { return pow(mod - 2); } mint pow(ll n) const { mint ret(1), mul(x); while (n > 0) { if (n & 1) ret *= mul; mul *= mul; n >>= 1; } return ret; } }; void solve(){ ll N, K; cin >> N >> K; ll MX = 100009; vl a(N); rep(i, N) cin >> a[i]; vl b(N); rep(i, N) cin >> b[i]; ll ok = 1, ng = *min_element(all(a)) + 1; vv dp(MX+1); reps(i, 1, MX+1){ for(ll k = i; k <= MX; k += i){ dp[k].push_back(i); } } while(ng - ok > 1){ ll mid = (ok + ng) / 2; ll cnt = 0; rep(i, N){ ll tmp = 1e18; for(ll x : dp[a[i]]){ if(x < mid) continue; ll y = x - (b[i] % x); if(y == x) y = 0; chmin(tmp, y); } cnt += tmp; } if(cnt <= K) ok = mid; else ng = mid; } cout << ok << endl; } int main(){ cin.tie(nullptr); ios_base::sync_with_stdio(false); cout << fixed << setprecision(20); int t = 1; // cin >> t; while(t--){ solve(); } }