#include #include #include #include #include #include #include #include #include #include #include #include #define debug_value(x) cerr << "line" << __LINE__ << ":<" << __func__ << ">:" << #x << "=" << x << endl; #define debug(x) cerr << "line" << __LINE__ << ":<" << __func__ << ">:" << x << endl; 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; } using namespace std; typedef long long ll; const ll MOD = 998244353; class ModInt{ public: ll v; ModInt(ll _v = 0){ if(_v < 0) _v = (v%MOD)+MOD; if(_v >= MOD) _v %= MOD; v = _v; } ModInt operator+(ll n){ return ModInt((v+n)%MOD); } ModInt operator-(ll n){ return ModInt((v-n+MOD)%MOD); } ModInt operator*(ll n){ if(n >= MOD) n %= MOD; return ModInt((v*n)%MOD); } ModInt operator/(ll n){ return ModInt((ModInt(n).inv()*v).v%MOD); } void operator+=(ll n){ v = (v+n)%MOD; } void operator-=(ll n){ v = (v-n+MOD)%MOD; } void operator*=(ll n){ v = (v*n+MOD)%MOD; } ModInt operator+(ModInt n){ return ModInt((v+n.v)%MOD); } ModInt operator-(ModInt n){ return ModInt((v-n.v+MOD)%MOD); } ModInt operator*(ModInt n){ return ModInt((v*n.v)%MOD); } ModInt operator/(ModInt n){ return ModInt((n.inv()*v).v%MOD); } void operator+=(ModInt n){ v = (v+n.v)%MOD; } void operator-=(ModInt n){ v = (v-n.v+MOD)%MOD; } void operator*=(ModInt n){ v = (v*n.v)%MOD; } void operator=(ModInt n){ v = n.v; } bool operator==(ModInt n){ return v == n.v; } bool operator!=(ModInt n){ return v != n.v; } void operator=(ll n){ v = n%MOD; } ModInt inv(){ if(v == 1) return ModInt(1); else return ModInt(MOD-ModInt(MOD%v).inv().v*(MOD/v)%MOD); } }; ostream& operator<<(ostream& os, const ModInt& m){ os << m.v; return os; } istream & operator >> (istream &in, ModInt &m){ in >> m.v; return in; } using mint = ModInt; mint dp[2005][2005]; mint dp_cumsum[2005][2005]; // i刻みで持っておく mint cumsum[2005][2005]; // i刻みで持っておく mint sect_sum(int i, int l, int r){ // dp[i][l] + ... + dp[i][r]; assert((r-l)%i == 0); if(l-i < 0) return dp_cumsum[i][r]; return dp_cumsum[i][r]-dp_cumsum[i][l-i]; } mint sect_cumsum(int i, int l, int r){ if(l-i < 0) return cumsum[i][r]; return cumsum[i][r]-cumsum[i][l-i]; } int main(){ ios::sync_with_stdio(false); cin.tie(0); cout << setprecision(10) << fixed; int n, L, R; cin >> n >> L >> R; vector a(n+1); for(int i = 0; i < n; i++) cin >> a[i+1]; mint ans = 0; for(int i = 1; i <= n; i++) cumsum[i][0] = 1; for(int j = 1; j <= R; j++){ mint sum = 0; for(int i = 1; i <= n; i++){ if(j < i) continue; int l = j-min(a[i], j/i)*i; int r = j-i; dp[i][j] += sect_cumsum(i, l, r)-sect_sum(i, l, r); // cout << i << ' ' << j << ' ' << l << ' ' << r << ' ' << sect_sum(l, r) << '-' << sect_sum(i, l, r) << endl; dp_cumsum[i][j] = dp_cumsum[i][j-i]+dp[i][j]; sum += dp[i][j]; if(j >= L) { ans += dp[i][j]; } } for(int i = 1; i <= n; i++){ if(i > j) cumsum[i][j] = sum; else cumsum[i][j] = cumsum[i][j-i] + sum; } } cout << ans << endl; }