#include //#include using namespace std; // using namespace atcoder; // using mint = modint1000000007; // const int mod = 1000000007; // using mint = modint998244353; // const int mod = 998244353; // const int INF = 1e9; // const long long LINF = 1e18; #define rep(i, n) for (int i = 0; i < (n); ++i) #define rep2(i, l, r) for (int i = (l); i < (r); ++i) #define rrep(i, n) for (int i = (n)-1; i >= 0; --i) #define rrep2(i, l, r) for (int i = (r)-1; i >= (l); --i) #define all(x) (x).begin(), (x).end() #define allR(x) (x).rbegin(), (x).rend() #define P pair template inline bool chmax(A &a, const B &b) { if (a < b) { a = b; return true; } return false; } template inline bool chmin(A &a, const B &b) { if (a > b) { a = b; return true; } return false; } #ifndef KWM_T_MATH_SIEVE_HPP #define KWM_T_MATH_SIEVE_HPP #include namespace kwm_t::math { /** * @brief エラトステネスの篩(bool版) * is_prime[i] を返す */ inline std::vector sieve(int n) { std::vector is_prime(n + 1, true); if (n >= 0) is_prime[0] = false; if (n >= 1) is_prime[1] = false; for (int p = 2; p * p <= n; ++p) { if (!is_prime[p]) continue; for (int i = p * p; i <= n; i += p) { is_prime[i] = false; } } return is_prime; } /** * @brief 素数列挙 * sieve を利用して primes を返す */ inline std::vector primes(int n) { auto is_prime = sieve(n); std::vector res; for (int i = 2; i <= n; ++i) { if (is_prime[i]) res.push_back(i); } return res; } } #endif #ifndef KWM_T_MATH_CONVOLUTION_INTERNAL_MULTIPLE_TRANSFORM_HPP #define KWM_T_MATH_CONVOLUTION_INTERNAL_MULTIPLE_TRANSFORM_HPP #include /** * @brief 倍数ゼータ / メビウス変換 * * f[i] = Σ_{j: i|j} f[j] * * 典型用途: * - GCD畳み込み * * 計算量: * O(n log n) * * verified: */ namespace kwm_t::math::convolution::internal { template void multiple_zeta_transform(std::vector& f, const std::vector& is_prime) { int n = (int)f.size(); for (int p = 2; p < n; ++p) { if (!is_prime[p]) continue; for (int i = (n - 1) / p; i >= 1; --i) { f[i] += f[i * p]; } } } template void multiple_mobius_transform(std::vector& f, const std::vector& is_prime) { int n = (int)f.size(); for (int p = 2; p < n; ++p) { if (!is_prime[p]) continue; for (int i = 1; i * p < n; ++i) { f[i] -= f[i * p]; } } } } #endif // KWM_T_MATH_CONVOLUTION_INTERNAL_MULTIPLE_TRANSFORM_HPP int main() { std::ios::sync_with_stdio(false); std::cin.tie(nullptr); int n, w; cin >> n >> w; int b = 18; vectora(1 << b); vectorx(n), y(n); rep(i, n)cin >> x[i]; rep(i, n)cin >> y[i]; rep(i, n)a[x[i]] += y[i]; auto prime = kwm_t::math::sieve(1 << b); kwm_t::math::convolution::internal::multiple_zeta_transform(a, prime); long long ans = 0; rep2(i, w, 1 << b)chmax(ans, a[i]); cout << ans << endl; return 0; }