#include //#include //#include //#include #define rep(i, a) for (int i = (int)0; i < (int)a; ++i) #define rrep(i, a) for (int i = (int)a - 1; i >= 0; --i) #define REP(i, a, b) for (int i = (int)a; i < (int)b; ++i) #define RREP(i, a, b) for (int i = (int)a - 1; i >= b; --i) #define repl(i, a) for (ll i = (ll)0; i < (ll)a; ++i) #define pb push_back #define eb emplace_back #define all(x) x.begin(), x.end() #define rall(x) x.rbegin(), x.rend() #define popcount __builtin_popcount #define popcountll __builtin_popcountll #define fi first #define se second using ll = long long; constexpr ll mod = 1e9 + 7; constexpr ll mod_998244353 = 998244353; constexpr ll INF = 1LL << 60; // #pragma GCC target("avx2") // #pragma GCC optimize("O3") // #pragma GCC optimize("unroll-loops") //using lll=boost::multiprecision::cpp_int; //using Double=boost::multiprecision::number>;//仮数部が1024桁 template inline bool chmin(T &a, T b) { if (a > b) { a = b; return true; } return false; } template inline bool chmax(T &a, T b) { if (a < b) { a = b; return true; } return false; } template T mypow(T x, T n, const T &p = -1) { //x^nをmodで割った余り if (p != -1) { x =(x%p+p)%p; } T ret = 1; while (n > 0) { if (n & 1) { if (p != -1) ret = (ret * x) % p; else ret *= x; } if (p != -1) x = (x * x) % p; else x *= x; n >>= 1; } return ret; } using namespace std; //using namespace atcoder; template struct Matrix{ vector> mat; Matrix(int n,int m,T val=0){ mat.assign(n,vector(m,val)); } inline vector& operator[](int k){//変更するとき return mat.at(k); } inline const vector& operator[](int k) const{//変更しないとき return mat.at(k); } Matrix I(int n){ Matrix v(n,n); for(int i=0;i>work(n,vector(m,0)); for(int i=0;iI(n); while(k>0){ if(k&1){ v*=*this; } *this*=*this; k>>=1; } mat.swap(v.mat); return *this; } Matrix operator+(const Matrix &v) const{ return (Matrix(*this)+=v);//thisポインタを渡しているのでコピーコンストラクタが呼ばれている } Matrix operator-(const Matrix &v) const{ return (Matrix(*this)-=v);//thisポインタを渡しているのでコピーコンストラクタが呼ばれている } Matrix operator*(const Matrix &v) const{ return (Matrix(*this)*=v);//thisポインタを渡しているのでコピーコンストラクタが呼ばれている } Matrix operator^(const Matrix &v) const{ return (Matrix(*this)^=v);//thisポインタを渡しているのでコピーコンストラクタが呼ばれている } }; template struct Modint{ int x; Modint():x(0){} Modint(int64_t y):x((y%mod+mod)%mod){} Modint &operator+=(const Modint &p){ if((x+=p.x)>=mod) x -= mod; return *this; } Modint &operator-=(const Modint &p){ if((x+=mod-p.x)>=mod) x -= mod; return *this; } Modint &operator*=(const Modint &p){ x = (1LL * x * p.x) % mod; return *this; } Modint &operator/=(const Modint &p){ *this *= p.inverse(); return *this; } Modint operator-() const { return Modint(-x); } Modint operator+(const Modint &p) const{ return Modint(*this) += p; } Modint operator-(const Modint &p) const{ return Modint(*this) -= p; } Modint operator*(const Modint &p) const{ return Modint(*this) *= p; } Modint operator/(const Modint &p) const{ return Modint(*this) /= p; } bool operator==(const Modint &p) const { return x == p.x; } bool operator!=(const Modint &p) const{return x != p.x;} Modint inverse() const{//非再帰拡張ユークリッド int a = x, b = mod, u = 1, v = 0; while(b>0){ int t = a / b; swap(a -= t * b, b); swap(u -= t * v, v); } return Modint(u); } Modint pow(int64_t n) const{//繰り返し二乗法 Modint ret(1), mul(x); while(n>0){ if(n&1) ret *= mul; mul *= mul; n >>= 1; } return ret; } friend ostream &operator<<(ostream &os,const Modint &p){ return os << p.x; } }; using modint = Modint; using modint2= Modint; void solve() { ll n; cin>>n; Matrixmat(2,2),a(2,2); a[0][0]=1; mat[0][0]=1,mat[0][1]=1; mat[1][0]=1; mat^=n; mat*=a; modint ans=mat[0][0]; ans*=mat[1][0]; cout<