#include using namespace std; #define ll long long #define ull unsigned long long #define ld long double using LL = long long; using ULL = unsigned long long; using VI = vector; using VVI = vector; using VVVI = vector; using VL = vector; using VVL = vector; using VVVL = vector; using VB = vector; using VVB = vector; using VVVB = vector; using VD = vector; using VVD = vector; using VVVD = vector; using VC = vector; using VS = vector; using VVC = vector; using PII = pair; using PLL = pair; using PDD = pair; using PIL = pair; using MII = map; using MLL = map; using SI = set; using SL = set; using MSI = multiset; using MSL = multiset; template using MAXPQ = priority_queue; template using MINPQ = priority_queue< T, vector, greater >; const ll MOD = 1000000007; const ll MOD2 = 998244353; const ll INF = 1LL << 60; #define PI 3.14159265358979323846 #define FOR(i, a, b) for(int i = (a); i < (b); ++i) #define REP(i, n) FOR(i, 0, n) #define EACH(e, v) for(auto &e : v) #define RITR(it, v) for(auto it = (v).rbegin(); it != (v).rend(); ++it) #define ALL(v) v.begin(),v.end() vector x8={1,1,1,0,0,-1,-1,-1},y8={1,0,-1,1,-1,1,0,-1}; int dx4[4]={1,-1,0,0}, dy4[4]={0,0,1,-1}; /* memo -uf,RMQ(segtree),BIT,BIT2,SegTree,SegTreeLazy -isprime,Eratosthenes,gcdlcm,factorize,divisors,modpow,moddiv nCr(+modnCr,inverse,extend_euclid.powmod),tobaseB,tobase10 -dijkstra,Floyd,bellmanford,sccd,topological,treediamiter -compress1,compress2,rotate90 -co,ci,fo1,fo2,fo3,fo4 -bitsearch,binaryserach -bfs -SegTreedec,SegTreeLazydec */ int main(){ cin.tie(0); ios_base::sync_with_stdio(0); ll N,K; cin >> N >> K; VL a(N); for(ll i = 0; i < N; i++) cin >> a[i]; VL b(N); for(ll i = 0; i < N; i++) cin >> b[i]; //cout << -1 << endl; VVL div(N); ll M = INF; for(ll i = 0; i < N; i++){ for(ll j = 1; j*j <= a[i]; j++){ if(a[i] % j == 0){ ll d = a[i]/j; div[i].push_back(j); if(j!=d) div[i].push_back(d); } } sort(ALL(div[i])); M = min(M,div[i].back()); } //cout << -1 << endl; auto ok = [&](ll X){ ll res = 0; for(ll i = 0; i < N; i++){ ll need = 1e10+1; for(auto v : div[i]){ if(v >= X){ ll p = b[i] % v; need = min((v-p)%v,need); } } res += need; } return res <= K; }; ll l = 0, r = M+1; while(r-l>1){ ll m = (l+r)/2; if(ok(m)) l = m; else r = m; } cout << l << '\n'; }