#ifdef ONLINE_JUDGE #pragma GCC optimize("Ofast") #pragma GCC optimize("unroll-loops") #endif #include #include using namespace std; using namespace atcoder; typedef long long ll; typedef unsigned long long ull; typedef long double ld; #define rep(i, n) for (ll i = 0; i < (ll)(n); i++) #define rrep(i,start,end) for (ll i = start;i >= (ll)(end);i--) #define repn(i,end) for(ll i = 0; i <= (ll)(end); i++) #define reps(i,start,end) for(ll i = start; i < (ll)(end); i++) #define repsn(i,start,end) for(ll i = start; i <= (ll)(end); i++) #define each(p,a) for(auto &p:a) typedef vector vll; typedef vector> vpll; typedef vector>> vvpll; typedef vector> vvll; typedef vector>> vvvll; typedef set sll; typedef map mpll; typedef pair pll; typedef tuple tpl3; typedef tuple tpl4; typedef tuple tpl5; typedef tuple tpl6; #define LL(...) ll __VA_ARGS__; input(__VA_ARGS__) #define LD(...) ld __VA_ARGS__; input(__VA_ARGS__) #define Str(...) string __VA_ARGS__; input(__VA_ARGS__) #define Ch(...) char __VA_ARGS__; input(__VA_ARGS__) #define all(a) (a).begin(),(a).end() #define UNIQUE(v) v.erase( unique(v.begin(), v.end()), v.end() ); #define sz(x) (ll)x.size() // << std::fixed << std::setprecision(10) const ll INF = 1LL << 60; const ld EPS = 1e-9; ll lceil(ll a,ll b){if(a%b==0){return a/b;}if(a>=0){return (a/b)+1;}else{return -((-a)/b);}} ll lfloor(ll a,ll b){if(a%b==0){return a/b;}if(a>=0){return (a/b);}else{return -((-a)/b)-1;}} inline ll positive_mod(ll a,ll m){return (a % m + m)%m;} inline ll popcnt(ull a){ return __builtin_popcountll(a);} //0indexed inline ll topbit(ll a){assert(a != 0);return 63 - __builtin_clzll(a);} inline ll smlbit(ll a){assert(a != 0);return __builtin_ctzll(a);} template bool chmin(T& a, T b){if(a > b){a = b;return true;}return false;} template bool chmax(T& a, T b){if(a < b){a = b;return true;}return false;} template std::istream &operator>>(std::istream&is,std::vector&v){for(T &in:v){is>>in;}return is;} template std::ostream &operator<<(std::ostream&os,const std::vector&v){for(auto it=std::begin(v);it!=std::end(v);){os<<*it<<((++it)!=std::end(v)?" ":"");}return os;} templatestd::ostream &operator<< (std::ostream &os, std::pair p){os << "{" << p.first << "," << p.second << "}";return os;} templatevoid input(T&... a){(cin >> ... >> a);} void print(){cout << endl;} templatevoid print(const T& a, const Ts&... b){cout << a;((cout << ' ' << b), ...);cout << endl;} template void pspace(const T& a){ cout << a << ' ';} void perr(){cerr << endl;} templatevoid perr(const T& a, const Ts&... b){cerr << a;((cerr << ' ' << b), ...);cerr << endl;} void yes(bool i = true){ return print(i?"yes":"no"); } void Yes(bool i = true){ return print(i?"Yes":"No"); } void YES(bool i = true){ return print(i?"YES":"NO"); } template vector &operator++(vector &v) {for(auto &e : v) e++;return v;} template vector operator++(vector &v, signed) {auto res = v;for(auto &e : v) e++;return res;} template vector &operator--(vector &v) {for(auto &e : v) e--;return v;} template vector operator--(vector &v, signed) {auto res = v;for(auto &e : v) e--;return res;} //grid探索用 vector _ta = {0,0,1,-1,1,1,-1,-1}; vector _yo = {1,-1,0,0,1,-1,1,-1}; bool isin(ll now_i,ll now_j,ll h,ll w){return (0<=now_i && now_i < h && 0 <= now_j && now_j < w);} ll lpow(ll x,ll n){ll ans = 1;while(n >0){if(n & 1)ans *= x;x *= x;n >>= 1;}return ans;} ll Modlpow(ll x,ll n,ll m){ll ans = 1;ll a = x%m;while(n >0){if(n & 1){ans *= a;ans%= m;}a *= a;a %= m;n >>= 1;}return ans;} const ll MOD9 = 998244353LL; const ll MOD10 = 1000000007LL; //約数列挙 vector divisor(ll n){ vector ret; for(ll i = 1;i * i <= n;i++){ if(n % i == 0){ ret.push_back(i); if(i * i != n)ret.push_back(n/i); } } //昇順 sort(all(ret)); return ret; } int main(){ ios::sync_with_stdio(false);cin.tie(nullptr); LL(n,k); vll a(n),b(n);cin >> a >> b; // ll ans = INF; vvll rets(n); vvll mins(n); rep(i,n){ vll ret = divisor(a[i]); rets[i] = ret; mins[i].resize(sz(ret)+1,INF); rep(j,sz(ret)){ ll p = ret[j]; if(p > b[i]){ mins[i][j] = p-b[i]; }else if(p == b[i]){ mins[i][j] = 0; }else{ // b[i] > p ll bi = b[i] % p; mins[i][j] = (p-bi)%p; } // if(p > b[i]){ // if(k >= p-b[i]){ // mins[i][j] = p-b[i]; // } // }else{ // ll num = lceil(b[i]-p,a[i]); // ll nai = a[i] * num + p; // if(k >= nai-b[i]){ // mins[i][j] = nai-b[i]; // } // } } rrep(j,sz(ret)-1,0){ chmin(mins[i][j],mins[i][j+1]); } } ll ng = lpow(10,5)+10;//条件を満たさない値 ll ok = 0;//条件を満たす値 auto f = [&](ll x)->bool { ll need = 0; rep(i,n){ ll id = lower_bound(all(rets[i]),x) - rets[i].begin(); need += mins[i][id]; if(need > k)return false; } return need <= k; }; while(abs(ok-ng) > 1){ ll mid = (ok + ng)/2; if(f(mid)){ ok = mid; }else{ ng = mid; } } cout << ok << endl; }