#include #include using namespace std; using ll = long long; using pii = pair; using pll = pair; using vi = vector; using vl = vector; #define rep3(i, a, b, c) for (ll i = (a); i < (b); i += (c)) #define rep2(i, a, b) rep3(i, a, b, 1) #define rep1(i, n) rep2(i, 0, n) #define rep0(n) rep1(aaaaa, n) #define ov4(a, b, c, d, name, ...) name #define rep(...) ov4(__VA_ARGS__, rep3, rep2, rep1, rep0)(__VA_ARGS__) #define per(i, a, b) for (ll i = (a) - 1; i >= (b); i--) #define fore(e, v) for (auto&& e : v) #define all(a) begin(a), end(a) #define sz(a) (int)(size(a)) #define lb(v, x) (lower_bound(all(v), x) - begin(v)) #define eb emplace_back template bool chmin(T& a, const S& b) { return a > b ? a = b, 1 : 0; } template bool chmax(T& a, const S& b) { return a < b ? a = b, 1 : 0; } const int INF = 1e9 + 100; const ll INFL = 3e18 + 100; #define i128 __int128_t struct _ { _() { cin.tie(0)->sync_with_stdio(0), cout.tie(0); } } __; constexpr int MOD=998244353; struct mod_fac{ vector fac,inv,invf; mod_fac (int n): fac(n),inv(n),invf(n){ fac[0]=fac[1]=1; inv[0]=inv[1]=1; invf[0]=invf[1]=1; for(int i=2;i matmul(vector a,vector b){ int n=sz(a); vector ret(n,vl(n)); rep(i,n){ rep(j,n){ rep(k,n){ ch(ret[i][j],a[i][k]*b[k][j]); } } } return ret; } vector matpow(vector a,ll N){ int n=sz(a); vector ret(n,vl(n)); rep(i,n)ret[i][i]=1; while(N>0){ if(N%2){ ret=matmul(ret,a); } a=matmul(a,a); N/=2; } return ret; } vl vecmul(vl v,vector m){ int n=sz(v); vl ret(n); rep(i,n)rep(j,n)ch(ret[j],v[i]*m[i][j]); return ret; } int main(){ ll K,L,R; cin>>K>>L>>R; vector mat(K+4,vl(K+4)); mod_fac F(K+10); mat[0][0]=K; mat[K+1][0]=1; mat[K+2][0]=1; mat[K+2][K+2]=K; mat[0][K+3]=1; mat[K+3][K+3]=1; rep(i,K+1){ rep(j,i+1){ mat[j+1][i+1]=F.comb(i,j); } } // fore(i,mat){ // fore(j,i)cerr<