#include using i32 = std::int32_t; using u32 = std::uint32_t; using i64 = std::int64_t; using u64 = std::uint64_t; using isize = std::ptrdiff_t; using usize = std::size_t; class rep { struct Iter { usize itr; constexpr Iter(const usize pos) noexcept : itr(pos) {} constexpr void operator++() noexcept { ++itr; } constexpr bool operator!=(const Iter& other) const noexcept { return itr != other.itr; } constexpr usize operator*() const noexcept { return itr; } }; const Iter first, last; public: explicit constexpr rep(const usize first, const usize last) noexcept : first(first), last(std::max(first, last)) {} constexpr Iter begin() const noexcept { return first; } constexpr Iter end() const noexcept { return last; } }; template constexpr T totient(T x) { T ret = x; for (T i = 2; i * i <= x; ++i) { if (x % i == 0) { ret /= i; ret *= i - 1; while (x % i == 0) x /= i; } } if (x > 1) { ret /= x; ret *= x - 1; } return ret; } template constexpr T rem_euclid(T value, const T& mod) { return (value %= mod) >= 0 ? value : value + mod; } template * = nullptr> class StaticModint { using Mint = StaticModint; static inline constexpr u32 PHI = totient(MOD); u32 v; public: static constexpr u32 mod() noexcept { return MOD; } template and std::is_integral_v>* = nullptr> static constexpr T normalize(const T x) noexcept { return rem_euclid>(x, MOD); } template and std::is_integral_v>* = nullptr> static constexpr T normalize(const T x) noexcept { return x % MOD; } constexpr StaticModint() noexcept : v(0) {} template constexpr StaticModint(const T x) noexcept : v(normalize(x)) {} template static constexpr Mint raw(const T x) noexcept { Mint ret; ret.v = x; return ret; } constexpr u32 get() const noexcept { return v; } constexpr Mint neg() const noexcept { return raw(v == 0 ? 0 : MOD - v); } constexpr Mint inv() const noexcept { return pow(PHI - 1); } constexpr Mint pow(u64 exp) const noexcept { Mint ret(1), mult(*this); for (; exp > 0; exp >>= 1) { if (exp & 1) ret *= mult; mult *= mult; } return ret; } constexpr Mint operator-() const noexcept { return neg(); } constexpr Mint operator~() const noexcept { return inv(); } constexpr Mint operator+(const Mint& rhs) const noexcept { return Mint(*this) += rhs; } constexpr Mint& operator+=(const Mint& rhs) noexcept { if ((v += rhs.v) >= MOD) v -= MOD; return *this; } constexpr Mint operator-(const Mint& rhs) const noexcept { return Mint(*this) -= rhs; } constexpr Mint& operator-=(const Mint& rhs) noexcept { if (v < rhs.v) v += MOD; v -= rhs.v; return *this; } constexpr Mint operator*(const Mint& rhs) const noexcept { return Mint(*this) *= rhs; } constexpr Mint& operator*=(const Mint& rhs) noexcept { v = (u64)v * rhs.v % MOD; return *this; } constexpr Mint operator/(const Mint& rhs) const noexcept { return Mint(*this) /= rhs; } constexpr Mint& operator/=(const Mint& rhs) noexcept { return *this *= rhs.inv(); } constexpr bool operator==(const Mint& rhs) const noexcept { return v == rhs.v; } constexpr bool operator!=(const Mint& rhs) const noexcept { return v != rhs.v; } friend std::ostream& operator<<(std::ostream& stream, const Mint& rhs) { return stream << rhs.v; } }; using Modint1000000007 = StaticModint<1000000007>; using Modint998244353 = StaticModint<998244353>; template struct NumSemiRing { using Type = T; static constexpr T zero() { return T(0); } static constexpr T one() { return T(1); } static constexpr T sum(const T& x, const T& y) { return x + y; } static constexpr T product(const T& x, const T& y) { return x * y; } }; template class SemiRingMatrix { using T = typename S::Type; using Self = SemiRingMatrix; std::vector> data; public: SemiRingMatrix() = default; explicit SemiRingMatrix(const usize h, const usize w, const T& val = S::zero()) : data(h, std::vector(w, val)) {} SemiRingMatrix(const std::vector>& vec) : data(vec) {} SemiRingMatrix(const std::initializer_list>& list) { data.reserve(list.size()); for (const auto& v : list) { data.emplace_back(v.begin(), v.end()); } } template SemiRingMatrix(const std::array, H>& arr) { data.reserve(H); for (const auto& v : arr) { data.emplace_back(v.begin(), v.end()); } } usize height() const { return data.size(); } usize width() const { return data.empty() ? 0 : data[0].size(); } T& operator()(const usize i, const usize j) { assert(i < height()); assert(j < width()); return data[i][j]; } const T& operator()(const usize i, const usize j) const { assert(i < height()); assert(j < width()); return data[i][j]; } Self operator+(const Self& other) const { return Self(*this) += other; } Self& operator+=(const Self& other) { assert(height() == other.height()); assert(width() == other.width()); for (const usize i : rep(0, height())) { for (const usize j : rep(0, width())) { data[i][j] = S::sum(data[i][j], other.data[i][j]); } } return *this; } Self operator*(const Self& other) const { assert(width() == other.height()); Self ret(height(), other.width(), S::zero()); for (const usize i : rep(0, height())) { for (const usize k : rep(0, width())) { for (const usize j : rep(0, other.width())) { ret.data[i][j] = S::sum(ret.data[i][j], S::product(data[i][k], other.data[k][j])); } } } return ret; } Self operator*(const T& other) const { return Self(*this) *= other; } Self& operator*=(const T& other) { for (const usize i : rep(0, height())) { for (const usize j : rep(0, width())) { data[i][j] = S::product(data[i][j], other); } } } Self pow(u64 exp) const { assert(height() == width()); Self ret(height(), width(), S::zero()), mult(*this); for (const usize i : rep(0, height())) { ret.data[i][i] = S::one(); } for (; exp > 0; exp >>= 1) { if (exp & 1) ret = ret * mult; mult = mult * mult; } return ret; } }; template using Vec = std::vector; using Fp = Modint1000000007; using Num = NumSemiRing; using Matrix = SemiRingMatrix; void main_() { u32 C, N, M; std::cin >> N; C = 3; M = 1; Matrix A = {{1, 0, 0, 0}}; Matrix B = {{Fp(1) / Fp(C), 0, 0, 1}, {Fp(C - 1) / Fp(C), 0, 0, 0}, {0, Fp(1) / Fp(C - 1), Fp(1) / Fp(C), Fp(C - 2) / Fp(C - 1)}, {0, 0, Fp(C - 1) / Fp(C), 0}}; const auto P = (A * B.pow(N))(0, 0); std::cout << Fp(1) - (Fp(1) - P).pow(M) << '\n'; } int main() { std::ios_base::sync_with_stdio(false); std::cin.tie(nullptr); usize cases; std::cin >> cases; while (cases--) main_(); return 0; }