#include using namespace std; #pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #define INT(...) \ int __VA_ARGS__; \ IN(__VA_ARGS__) #define LL(...) \ ll __VA_ARGS__; \ IN(__VA_ARGS__) #define STR(...) \ string __VA_ARGS__; \ IN(__VA_ARGS__) #define CHR(...) \ char __VA_ARGS__; \ IN(__VA_ARGS__) #define DBL(...) \ double __VA_ARGS__; \ IN(__VA_ARGS__) #define ll long long #define cout std::cout #define yes cout<<"Yes"<<"\n" #define no cout<<"No"<<"\n" #define rep(i,a,b) for(int i=a;i=b;i--) #define fore(i,a) for(auto &i:a) #define all(x) (x).begin(),(x).end() #define allr(x) (x).rbegin(),(x).rend() #define SUM(v) accumulate(all(v), 0LL) #define MIN(v) *min_element(all(v)) #define MAX(v) *max_element(all(v)) #define lb(c, x) distance((c).begin(), lower_bound(all(c), (x))) #define ub(c, x) distance((c).begin(), upper_bound(all(c), (x))) #define pii pair #define pll pair #define pb push_back #define eb emplace_back #define ff first #define ss second #define vi vector #define vll vector #define vc vector #define vvi vector> #define vec(type, name, ...) vector name(__VA_ARGS__) #define VEC(type, name, size) \ vector name(size); \ IN(name) int scan() { return getchar(); } void scan(int &a) { cin >> a; } void scan(long long &a) { cin >> a; } void scan(char &a) { cin >> a; } void scan(double &a) { cin >> a; } void scan(string &a) { cin >> a; } template void scan(pair &p) { scan(p.first), scan(p.second); } template void scan(vector &); template void scan(vector &a) { for(auto &i : a) scan(i); } template void scan(T &a) { cin >> a; } void IN() {} template void IN(Head &head, Tail &...tail) { scan(head); IN(tail...); } #define vv(type, name, h, ...) vector> name(h, vector(__VA_ARGS__)) #define VV(type, name, h, w) \ vector> name(h, vector(w)); \ IN(name) #define vvv(type, name, h, w, ...) vector>> name(h, vector>(w, vector(__VA_ARGS__))) #define vvvv(type, name, a, b, c, ...) \ vector>>> name(a, vector>>(b, vector>(c, vector(__VA_ARGS__)))) template using min_priority_queue = priority_queue, greater>; template pair operator-(const pair &x, const pair &y) { return pair(x.ff - y.ff, x.ss - y.ss); } template pair operator+(const pair &x, const pair &y) { return pair(x.ff + y.ff, x.ss + y.ss); } template pair operator&(const pair &l, const pair &r) { return pair(max(l.ff, r.ff), min(l.ss, r.ss)); } template inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } #define UNIQUE(x) sort(all(x)), x.erase(unique(all(x)), x.end()) //座標圧縮 template void zip(vector &x) { vector y = x; UNIQUE(y); for(int i = 0; i < x.size(); ++i) { x[i] = lb(y, x[i]); } } template T ceil(T x, T y) { assert(y >= 1); return (x > 0 ? (x + y - 1) / y : x / y); } template T floor(T x, T y) { assert(y >= 1); return (x > 0 ? x / y : (x + y - 1) / y); } long long POW(long long x, int n) { long long res = 1LL; for(; n; n >>= 1, x *= x) if(n & 1) res *= x; return res; } long long modpow(long long a, long long n, long long mod) { long long res = 1; while (n > 0) { if (n & 1) res = res * a % mod; a = a * a % mod; n >>= 1; } return res; } //return 0<=a&&a> x; return x; } ll lin() { unsigned long long x; cin >> x; return x; } long long sqrtll(long long x) { assert(x >= 0); long long rev = sqrt(x); while(rev * rev > x) --rev; while((rev+1) * (rev+1)<=x) ++rev; return rev; } int logN(long long n){ int ret=1; while((1LL<::max() / 2; const long long INFL = numeric_limits::max() / 2; template void debug(vector a){ rep(i,0,(int)a.size()){ cout< bool isPalin(T s){ rep(i,0,s.size()){ if(s[i]!=s[s.size()-i-1])return false; } return true; } //modint template class modint { using u64 = std::uint_fast64_t; public: u64 a; constexpr modint(const u64 x = 0) noexcept : a(x % Modulus) {} constexpr u64 &val() noexcept { return a; } constexpr const u64 &val() const noexcept { return a; } constexpr modint operator+(const modint rhs) const noexcept { return modint(*this) += rhs; } constexpr modint operator-(const modint rhs) const noexcept { return modint(*this) -= rhs; } constexpr modint operator*(const modint rhs) const noexcept { return modint(*this) *= rhs; } constexpr modint operator/(const modint rhs) const noexcept { return modint(*this) /= rhs; } constexpr modint &operator+=(const modint rhs) noexcept { a += rhs.a; if (a >= Modulus) { a -= Modulus; } return *this; } constexpr modint &operator-=(const modint rhs) noexcept { if (a < rhs.a) { a += Modulus; } a -= rhs.a; return *this; } constexpr modint &operator*=(const modint rhs) noexcept { a = a * rhs.a % Modulus; return *this; } constexpr modint &operator/=(modint rhs) noexcept { u64 exp = Modulus - 2; while (exp) { if (exp % 2) { *this *= rhs; } rhs *= rhs; exp /= 2; } return *this; } }; using mint9=modint<998244353>; using mint1=modint<1000000007>; constexpr pii dx4[4] = {pii{-1, 0}, pii{0, 1}, pii{1, 0}, pii{0, -1}}; constexpr pii dx8[8] = {{1, 0}, {1, 1}, {0, 1}, {-1, 1}, {-1, 0}, {-1, -1}, {0, -1}, {1, -1}}; constexpr pii dx[8] = { {-1, 0}, {0, -1}}; #define el "\n" #define endl "\n" int main(){ INT(n,W); vec(int,v,n);vec(int,w,n); rep(i,0,n)cin>>v[i]>>w[i]; vv(ll,dp1,n+1,W+20100,-INFL); vv(mint9,dp2,n+1,W+20100); dp1[0][10000]=0; dp2[0][10000]=1; rep(i,0,n){ rep(j,0,W+20100){ if(j+w[i]dp1[i+1][j+w[i]]){ dp2[i+1][j+w[i]]=dp2[i][j]; } chmax(dp1[i+1][j+w[i]],dp1[i][j]+v[i]); } if(dp1[i][j]==dp1[i+1][j]){ dp2[i+1][j]+=dp2[i][j]; }else if(dp1[i][j]>dp1[i+1][j]){ dp2[i+1][j]=dp2[i][j]; } chmax(dp1[i+1][j],dp1[i][j]); } } ll val=-INFL; mint9 cnt=0; rep(i,0,W+10001){ if(val==dp1[n][i])cnt+=dp2[n][i]; else if(val