#include using namespace std; template struct modint{ long long v; modint(long long x=0):v((x%=mod)<0?x+mod:x){} constexpr long long val()noexcept{return v;} constexpr modint& operator=(const modint o){v=o.v;return *this;} constexpr modint& operator+=(const modint o)noexcept{if((v+=o.v)>=mod)v-=mod;return *this;} constexpr modint& operator-=(const modint o)noexcept{if((v-=o.v)<0)v+=mod;return *this;} constexpr modint& operator*=(const modint o)noexcept{v=v*o.v%mod;return *this;} constexpr modint& operator/=(modint o)noexcept{return *this*=o.inv();} constexpr modint pow(long long n)noexcept{ modint ret=1,base=*this; while(n>0){if(n&1)ret*=base;base*=base;n/=2;} return ret; } constexpr modint inv()noexcept{return pow(mod-2);} constexpr modint operator+(const modint o)noexcept{return modint(*this)+=o;} constexpr modint operator+()noexcept{return *this;} constexpr modint operator-(const modint o)noexcept{return modint(*this)-=o;} constexpr modint operator-()noexcept{return modint()-*this;} constexpr modint operator*(const modint o)noexcept{return modint(*this)*=o;} constexpr modint operator/(const modint o)noexcept{return modint(*this)/=o;} constexpr modint& operator++()noexcept{return *this+=1;} constexpr modint operator++(int)noexcept{return (*this+=1)-1;} constexpr modint& operator--()noexcept{return *this-=1;} constexpr modint operator--(int)noexcept{return (*this-=1)+1;} constexpr bool operator==(const modint o)noexcept{return v==o.v;} constexpr bool operator!=(const modint o)noexcept{return v!=o.v;} friend ostream& operator<<(ostream& os,const modint o)noexcept{os<; int main(){ long long k,l,r; cin>>k>>l>>r; int m=k+4; vector> mat(m,vector(m,0)); mat[0][0]=k+1; mat[0][1]=-k; mat[0][2]=1; mat[0][3]=1; mat[1][0]=1; mat[2][2]=k; mint fact[k+1]; fact[0]=1; for(int i=1;i<=k;i++) fact[i]=fact[i-1]*i; for(int i=0;i<=k;i++) for(int j=0;j<=k;j++){ if(i> matrix1,vector> matrix2){ vector> ret(m,vector(m,0)); for(int i=0;i> e(m,vector(m,0)); for(int i=0;i0){ if(l&1) e_l=mul(mat_l,e_l); mat_l=mul(mat_l,mat_l); l>>=1; } auto mat_r=mat; auto e_r=e; while(r>0){ if(r&1) e_r=mul(mat_r,e_r); mat_r=mul(mat_r,mat_r); r>>=1; } mint sl=0,sr=0; vector v0(m,0); v0[0]=1; v0[2]=1; v0[m-1]=1; for(int i=0;i