結果
問題 | No.1073 無限すごろく |
ユーザー | penta |
提出日時 | 2020-06-06 15:00:30 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 4,745 bytes |
コンパイル時間 | 1,634 ms |
コンパイル使用メモリ | 175,352 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-06-06 02:01:23 |
合計ジャッジ時間 | 2,638 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,376 KB |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | AC | 1 ms
5,376 KB |
testcase_05 | AC | 2 ms
5,376 KB |
testcase_06 | AC | 1 ms
5,376 KB |
testcase_07 | AC | 2 ms
5,376 KB |
testcase_08 | AC | 1 ms
5,376 KB |
testcase_09 | AC | 2 ms
5,376 KB |
testcase_10 | AC | 2 ms
5,376 KB |
testcase_11 | AC | 2 ms
5,376 KB |
testcase_12 | AC | 1 ms
5,376 KB |
testcase_13 | AC | 2 ms
5,376 KB |
testcase_14 | AC | 2 ms
5,376 KB |
testcase_15 | AC | 2 ms
5,376 KB |
testcase_16 | AC | 2 ms
5,376 KB |
testcase_17 | AC | 2 ms
5,376 KB |
testcase_18 | AC | 2 ms
5,376 KB |
testcase_19 | AC | 2 ms
5,376 KB |
testcase_20 | AC | 2 ms
5,376 KB |
testcase_21 | AC | 2 ms
5,376 KB |
testcase_22 | AC | 2 ms
5,376 KB |
testcase_23 | AC | 2 ms
5,376 KB |
testcase_24 | AC | 2 ms
5,376 KB |
testcase_25 | AC | 2 ms
5,376 KB |
testcase_26 | AC | 2 ms
5,376 KB |
testcase_27 | AC | 2 ms
5,376 KB |
testcase_28 | AC | 2 ms
5,376 KB |
testcase_29 | AC | 2 ms
5,376 KB |
testcase_30 | AC | 2 ms
5,376 KB |
testcase_31 | AC | 2 ms
5,376 KB |
testcase_32 | AC | 2 ms
5,376 KB |
ソースコード
#include <bits/stdc++.h> #define rep(i,n) for (int i = 0; i < (n); ++i) #define all(x) (x).begin(),(x).end() using namespace std; using ll = long long; using P = pair<int,int>; template <class T> void chmin(T &a, const T &b) noexcept { if (b < a) a = b; } template <class T> void chmax(T &a, const T &b) noexcept { if (a < b) a = b; } void debug_out() { cout << "\n"; } template <class T, class... Args> void debug_out(const T &x, const Args &... args) { cout << x << " "; debug_out(args...);} #ifdef _DEBUG #define debug(...) debug_out(__VA_ARGS__) #else #define debug(...) #endif const int mod = 1000000007; // const int mod = 998244353; struct mint { ll x; mint(ll x=0):x((x%mod+mod)%mod){} mint operator-() const { return mint(-x);} mint& operator+=(const mint a) { if ((x += a.x) >= mod) x -= mod; return *this; } mint& operator-=(const mint a) { if ((x += mod-a.x) >= mod) x -= mod; return *this; } mint& operator*=(const mint a) { (x *= a.x) %= mod; return *this;} mint operator+(const mint a) const { return mint(*this) += a;} mint operator-(const mint a) const { return mint(*this) -= a;} mint operator*(const mint a) const { return mint(*this) *= a;} mint pow(ll t) const { if (!t) return 1; mint a = pow(t>>1); a *= a; if (t&1) a *= *this; return a; } // for prime mod mint inv() const { return pow(mod-2);} mint& operator/=(const mint a) { return *this *= a.inv();} mint operator/(const mint a) const { return mint(*this) /= a;} bool operator==(const mint rhs) const { return x == rhs.x; } bool operator!=(const mint rhs) const { return x != rhs.x; } bool operator<(const mint &a) const{ return x<a.x;}; }; istream& operator>>(istream& is, mint& a) { return is >> a.x;} ostream& operator<<(ostream& os, const mint& a) { return os << a.x;} template<typename E> struct Matrix { vector<vector<E> > A; Matrix() {} Matrix(int n, int m):A(n,vector<E>(m,0)) {} Matrix(int n):A(n,vector<E>(n,0)) {} int height() const { return A.size();} int width() const { return A[0].size();} vector<E>& operator[](int k){ return A.at(k);} const vector<E>& operator[](int k) const { return A.at(k);} static Matrix I(int n) { Matrix mat(n); for (int i = 0; i < n; ++i) mat[i][i] = 1; return mat; } Matrix T() { int n = height(), m = width(); vector<vector<E> > B = A; for (int i = 0; i < n; ++i) for (int j = 0; j < m; ++j) swap(B[i][j],B[j][i]); return B; } Matrix& operator+=(const Matrix B) { int n = height(), m = width(); assert(n == B.height() && m == B.width()); for (int i = 0; i < n; ++i) for (int j = 0; j < m; ++j) A[i][j] += B[i][j]; return *this; } Matrix& operator-=(const Matrix B) { int n = height(), m = width(); assert(n == B.height() && m == B.width()); for (int i = 0; i < n; ++i) for (int j = 0; j < m; ++j) A[i][j] -= B[i][j]; return *this; } Matrix& operator*=(const Matrix B) { int n = height(), r = width(), m = B.width(); assert(r == B.height()); vector<vector<E> > res(n,vector<E>(m,0)); for(int i = 0; i < n; ++i){ for(int k = 0; k < r; ++k){ for(int j = 0; j < m; ++j){ res[i][j] += A[i][k] * B[k][j]; } } } swap(A,res); return *this; } Matrix operator+(const Matrix B) const { return Matrix(*this) += B;} Matrix operator-(const Matrix B) const { return Matrix(*this) -= B;} Matrix operator*(const Matrix B) const { return Matrix(*this) *= B;} Matrix pow(ll t) const { if (!t) return Matrix::I(height()); Matrix a = pow(t>>1); a *= a; if (t&1) a *= *this; return a; } E determinant() { assert(width() == height()); Matrix B(*this); E res = 1; int n = height(); for(int i = 0; i < n; ++i) { int idx = -1; for(int j = i; j < n; ++j) { if(B[j][i] != 0) idx = j; } if(idx == -1) return 0; if(i != idx) { swap(B[i], B[idx]); res *= -1; } res *= B[i][i]; E vv = B[i][i]; for(int j = 0; j < n; ++j) { B[i][j] /= vv; } for(int j = i+1; j < n; ++j) { E a = B[j][i]; for(int k = 0; k < n; k++) { B[j][k] -= B[i][k] * a; } } } return res; } bool operator==(const Matrix rhs) const { return A == rhs.A; } bool operator!=(const Matrix rhs) const { return A != rhs.A; } }; int main() { std::cin.tie(nullptr); std::ios_base::sync_with_stdio(false); std::cout << std::fixed << std::setprecision(15); ll n; cin >> n; Matrix<mint> a(6,6); rep(i,6) a[0][i] = mint(1)/6; for (int i = 1; i < 6; ++i) a[i][i-1] = 1; auto res = a.pow(n); cout << res[0][0] << endl; return 0; }