#include #include using namespace std; namespace my{ using ml=atcoder::modint998244353; auto&operator>>(istream&i,ml&x){int t;i>>t;x=t;return i;} auto&operator<<(ostream&o,const ml&x){return o<<(int)x.val();} #define LL(...) ll __VA_ARGS__;lin(__VA_ARGS__) #define VL(n,...) vec__VA_ARGS__;setsize({n},__VA_ARGS__);lin(__VA_ARGS__) #define FO(n) for(ll ij=n;ij-->0;) #define FOR(i,...) for(auto[i,i##stop,i##step]=range(0,__VA_ARGS__);isync_with_stdio(0);cout<r{0,0,1};ll I=0;((r[I++]=a),...);if(!s&&I==1)swap(r[0],r[1]);r[0]-=s;if(s)r[2]*=-1;return r;} constexpr char newline=10; constexpr char space=32; template>auto&sort(auto&a,F f={}){ranges::sort(a,f);return a;} templateostream&operator<<(ostream&o,const unordered_map&m){fe(m,e)o<auto pack_kth(const auto&...a){return get(make_tuple(a...));} templateauto pack_slice(const auto&...a){return[&](index_sequence){return array{get(forward_as_tuple(a...))...};}(make_index_sequence{});} templateconcept vectorial=is_base_of_v,V>; templatestruct vec_attr{using core_type=T;static constexpr int rank=0;}; templatestruct vec_attr{using core_type=typename vec_attr::core_type;static constexpr int rank=vec_attr::rank+1;}; templateusing core_t=vec_attr::core_type; templateistream&operator>>(istream&i,vector&v){fe(v,e)i>>e;return i;} templateostream&operator<<(ostream&o,const vector&v){fe(v,e)o<?newline:space);return o;} templatestruct vec; templatestruct tensor_helper{using type=vec::type>;}; templatestruct tensor_helper<0,T>{using type=T;}; templateusing tensor=typename tensor_helper::type; templatestruct vec:vector{ using vector::vector; vec(const vector&v){vector::operator=(v);} templaterequires(sizeof...(A)>=3)vec(A...a){const ll n=sizeof...(a)-1;auto t=pack_slice(a...);ll s[n];fo(i,n)s[i]=t[i];*this=make_vec(s,pack_kth(a...));} templatestatic auto make_vec(const ll(&s)[n],T x){if constexpr(i==n-1)return vec(s[i],x);else{auto X=make_vec(s,x);return vec(s[i],X);}} vec&operator^=(const vec&u){this->insert(this->end(),u.begin(),u.end());return*this;} vec operator^(const vec&u)const{return vec{*this}^=u;} vec&operator+=(const vec&u){vec&v=*this;fo(i,v.size())v[i]+=u[i];return v;} vec&operator-=(const vec&u){vec&v=*this;fo(i,v.size())v[i]-=u[i];return v;} vec operator+(const vec&u)const{return vec{*this}+=u;} vec operator-(const vec&u)const{return vec{*this}-=u;} vec&operator++(){fe(*this,e)++e;return*this;} vec&operator--(){fe(*this,e)--e;return*this;} vec operator-()const{vec v=*this;fe(v,e)e=-e;return v;} vec&operator%=(auto M){vec&v=*this;fe(v,e)e%=M;return v;} vec operator%(auto M)const{return vec{*this}%=M;} ll size()const{return vector::size();} templateauto iota()const{vecr(size());fo(i,size())r[i]=i;return r;} void emplace_front(const V&x={}){this->emplace(this->begin(),x);} template>auto sort(F f={})const{vec v=*this;ranges::sort(v,f);return v;} template>auto sort_index(F f={})const{const vec&v=*this;return iota().sort([&](ll i,ll j){return!f(v[i],v[j])&&!f(v[j],v[i])?i>auto rank(F f={})const{auto o=sort_index(f);vecr(size());fo(i,size())r[o[i]]=i;return r;} auto pre_enumerate()const requires(!vectorial){ unordered_mappre; vecres(size(),-1); fo(i,size()){ auto e=(*this)[i]; if(pre.contains(e))res[i]=pre[e]; pre[e]=i; } return res; } auto rectangle_sum(ll a,ll b,ll c,ll d)const requires(vec_attr>::rank==2){ static_assert(is_same_v>>); assert(0<=a&&a<=b&&b<=size()); assert(0<=c&&c<=d&&d<=(*this)[0].size()); core_tr{}; if(b&&d)r+=(*this)[b-1][d-1]; if(a&&d)r-=(*this)[a-1][d-1]; if(b&&c)r-=(*this)[b-1][c-1]; if(a&&c)r+=(*this)[a-1][c-1]; return r; } }; templaterequires(sizeof...(A)>=2)vec(A...a)->vec(declval>()))>>>; vec(ll)->vec; templatevoid setsize(const ll(&l)[n],A&...a){((a=vec::make_vec(l,core_t{})),...);} void lin(auto&...a){(cin>>...>>a);} templatevoid pp(const auto&...a){ll n=sizeof...(a);((cout<0,c)),...);cout<=0)dp[i][j]+=su.rectangle_sum(pre[i],i,j-a[s[i]],j-a[s[i]]+1); if(j-b[s[i]]>=0)dp[i][j]+=su.rectangle_sum(pre[i],i,j-b[s[i]],j-b[s[i]]+1); su[i][j]=dp[i][j]+(i?su[i-1][j]:0)+(j?su[i][j-1]:0)-(i&&j?su[i-1][j-1]:0); } } pp(su[N][K]); }}