#pragma GCC target("avx") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include // #include // #include // #include // using namespace __gnu_pbds; // #include // namespace multiprecisioninteger = boost::multiprecision; // using cint=multiprecisioninteger::cpp_int; using namespace std; using ll=long long; #define double long double using datas=pair; using ddatas=pair; using tdata=pair; using vec=vector; using mat=vector; using pvec=vector; using pmat=vector; // using llset=tree,rb_tree_tag,tree_order_statistics_node_update>; #define For(i,a,b) for(i=a;i<(ll)b;++i) #define bFor(i,b,a) for(i=b,--i;i>=(ll)a;--i) #define rep(i,N) For(i,0,N) #define rep1(i,N) For(i,1,N) #define brep(i,N) bFor(i,N,0) #define brep1(i,N) bFor(i,N,1) #define all(v) (v).begin(),(v).end() #define allr(v) (v).rbegin(),(v).rend() #define vsort(v) sort(all(v)) #define vrsort(v) sort(allr(v)) #define endl "\n" #define eb emplace_back #define print(v) cout< inline bool chmax(T& a,T b){bool x=a inline bool chmin(T& a,T b){bool x=a>b;if(x)a=b;return x;} void startupcpp(){ cin.tie(0); ios::sync_with_stdio(false); cout<0){ if(n&1)res=res*a%mod; a=a*a%mod; n>>=1; } return res; } ll gcd(ll a,ll b){if(!b)return abs(a);return (a%b==0)?abs(b):gcd(b,a%b);} ll lcm(ll a,ll b){return a/gcd(a,b)*b;} ll countdigits(ll n){ ll ans=0; while(n){n/=10;ans++;} return ans; } ll sumdigits(ll n){ ll ans=0; while(n){ans+=n%10;n/=10;} return ans; } void topcoder(vector& v){ string s; while(1){ cin>>s; int i=s[0]=='{',x=0; while(s[i]>='0'&&s[i]<='9'){ x=x*10+s[i]-'0'; ++i; } v.eb(x); if(s[i]=='}')break; } } void topcoder(string& s){ string t; cin>>t; int i=1; while(t[i]!='"'){ s+=t[i++]; } } class matrix{ mat a; ll H,W; public: matrix(mat& g):a(g){ H=g.size(); W=g[0].size(); } matrix(ll i,ll j):a(i,vec(j,0)){H=i;W=j;} matrix(ll n):a(n,vec(n,0)){H=W=n;} inline vec& operator [](int k){ return a.at(k); } // matrix operator =(matrix b){ // this->a.swap(b->a); // this->H=b.H; // this->W=b.W; // return (*this); // } matrix operator +=(matrix b){ ll i,j; rep(i,this->H)rep(j,this->W)(*this)[i][j]+=b[i][j]; return (*this); } matrix operator -=(matrix b){ ll i,j; rep(i,this->H)rep(j,this->W)(*this)[i][j]-=b[i][j]; return (*this); } matrix operator *=(matrix b){ ll i,j,k; assert(this->W==b.H); matrix c(this->H,b.W); rep(i,this->H)rep(j,b.W){ c[i][j]=0; rep(k,this->W)(c[i][j]+=(*this)[i][k]*b[k][j]%mod)%=mod; } (*this)=c; return (*this); } matrix operator ^=(ll K){ assert(this->H==this->W); matrix c(this->H); ll i; rep(i,this->H)c[i][i]=1; if(K&1)c*=(*this); while(K){ K>>=1; (*this)*=(*this); if(K&1)c*=(*this); } this->a.swap(c.a); return (*this); } matrix operator +(matrix c){ return matrix(*this)+=c; } matrix operator -(matrix c){ return matrix(*this)-=c; } matrix operator *(matrix c){ return matrix(*this)*=c; } matrix operator ^(ll K){ return matrix(*this)^=K; } void out(){ for(auto x:a)output(x); } }; int main(){ // startupcpp(); // int codeforces;cin>>codeforces;while(codeforces--){ ll i,j,N,X,ans=0; cin>>N; matrix calc(6),v(6,1); rep(i,6)calc[0][i]=166666668; rep(i,5)calc[i+1][i]=1; v[0][0]=1; calc^=N; calc*=v; print(calc[0][0]); }