結果
問題 | No.314 ケンケンパ |
ユーザー |
|
提出日時 | 2019-03-07 11:49:35 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
CE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 4,160 bytes |
コンパイル時間 | 730 ms |
コンパイル使用メモリ | 78,452 KB |
最終ジャッジ日時 | 2024-11-14 21:17:31 |
合計ジャッジ時間 | 1,283 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。(リジャッジ後のコンパイルエラーは公開されます)
ただし、clay言語の場合は開発者のデバッグのため、公開されます。
ただし、clay言語の場合は開発者のデバッグのため、公開されます。
コンパイルメッセージ
main.cpp:17:39: error: '::numeric_limits' has not been declared 17 | template<class T> constexpr T INF = ::numeric_limits<T>::max()/32*15+208; | ^~~~~~~~~~~~~~ main.cpp:17:55: error: expected primary-expression before '>' token 17 | template<class T> constexpr T INF = ::numeric_limits<T>::max()/32*15+208; | ^ main.cpp:17:61: error: no matching function for call to 'max()' 17 | template<class T> constexpr T INF = ::numeric_limits<T>::max()/32*15+208; | ~~~~~^~ In file included from /home/linuxbrew/.linuxbrew/Cellar/gcc@12/12.3.0/include/c++/12/string:50, from /home/linuxbrew/.linuxbrew/Cellar/gcc@12/12.3.0/include/c++/12/bits/locale_classes.h:40, from /home/linuxbrew/.linuxbrew/Cellar/gcc@12/12.3.0/include/c++/12/bits/ios_base.h:41, from /home/linuxbrew/.linuxbrew/Cellar/gcc@12/12.3.0/include/c++/12/ios:42, from /home/linuxbrew/.linuxbrew/Cellar/gcc@12/12.3.0/include/c++/12/ostream:38, from /home/linuxbrew/.linuxbrew/Cellar/gcc@12/12.3.0/include/c++/12/iostream:39, from main.cpp:1: /home/linuxbrew/.linuxbrew/Cellar/gcc@12/12.3.0/include/c++/12/bits/stl_algobase.h:254:5: note: candidate: 'template<class _Tp> constexpr const _Tp& std::max(const _Tp&, const _Tp&)' 254 | max(const _Tp& __a, const _Tp& __b) | ^~~ /home/linuxbrew/.linuxbrew/Cellar/gcc@12/12.3.0/include/c++/12/bits/stl_algobase.h:254:5: note: template argument deduction/substitution failed: main.cpp:17:61: note: candidate expects 2 arguments, 0 provided 17 | template<class T> constexpr T INF = ::numeric_limits<T>::max()/32*15+208; | ~~~~~^~ /home/linuxbrew/.linuxbrew/Cellar/gcc@12/12.3.0/include/c++/12/bits/stl_algobase.h:300:5: note: candidate: 'templat
ソースコード
#include <iostream>#include <algorithm>#include <iomanip>#include <map>#include <set>#include <queue>#include <stack>#include <numeric>#include <bitset>static const int MOD = 1000000007;using ll = int64_t;using u32 = uint32_t;using namespace std;template<class T> constexpr T INF = ::numeric_limits<T>::max()/32*15+208;template<class T, T M = 1000000007>struct modint{T val;modint(): val(0){}explicit modint(int v): val(v) {}explicit modint(ll t){val = t%M; if(val < 0) val += M;}modint pow(ll k){modint res(1), x(val);while(k){if(k&1) res *= x;x *= x;k >>= 1;}return res;}modint& operator=(int a){return (*this) = modint(a);}modint inv() {return pow(M-2);}modint& operator+=(modint a){ val += a.val; if(val >= M) val -= M; return *this;}modint& operator-=(modint a){ val += M-a.val; if(val >= M) val -= M; return *this;}modint& operator*=(modint a){ val = 1LL*val*a.val%M; return *this;}modint& operator/=(modint a){ return (*this) *= a.inv();}modint operator+(modint a) const {return modint(val) +=a;}modint operator-(modint a) const {return modint(val) -=a;}modint operator*(modint a) const {return modint(val) *=a;}modint operator/(modint a) const {return modint(val) /=a;}modint operator-(){return val ? M-val : val;}bool operator==(const modint a) const {return val == a.val;}bool operator!=(const modint a) const {return val != a.val;}bool operator<(const modint a) const {return val < a.val;}};#include <assert.h>template<class T>struct matrix {vector<vector<T>> A;matrix() = default;matrix(size_t n, size_t m) : A(n, vector<T>(m)) {}explicit matrix(size_t n) : A(n, vector<T> (n)) {};size_t height() const { return (A.size()); }size_t width() const { return (A[0].size()); }const vector<T> &operator [] (int k) const { return A.at(k); }vector<T> &operator[] (int k) { return A.at(k); }static matrix I(size_t n){matrix mat(n);for (int i = 0; i < n; ++i) mat[i][i] = 1;return mat;}matrix &operator+= (const matrix &B){size_t h = height(), w = width();assert(h == B.height() && w == B.width());for (int i = 0; i < h; ++i) {for (int j = 0; j < w; ++j) {(*this)[i][j] += B[i][j];}}}matrix &operator-= (const matrix &B){size_t h = height(), w = width();assert(h == B.height() && w == B.width());for (int i = 0; i < h; ++i) {for (int j = 0; j < w; ++j) {(*this)[i][j] -= B[i][j];}}}matrix &operator*=(const matrix &B){size_t n = height(), m = B.width(), p = width();assert(p == B.height());matrix C (n, m);for (int i = 0; i < n; ++i) {for (int j = 0; j < m; ++j) {for (int k = 0; k < p; ++k) {C[i][j] = (C[i][j] + (*this)[i][k] * B[k][j]);}}}A.swap(C.A);return (*this);}template <class U>matrix &operator%= (const U &m){for (int i = 0; i < height(); ++i) {for (int j = 0; j < width(); ++j) {(*this)[i][j] %= m;}}}matrix pow(ll n) const {matrix a = (*this), res = I(height());while(n > 0){if(n & 1) res *= a;a *= a;n >>= 1;}return res;}matrix operator+(const matrix &A) const {return matrix(*this) += A;}matrix operator-(const matrix &A) const {return matrix(*this) -= A;}matrix operator*(const matrix &A) const {return matrix(*this) *= A;}template <class U>matrix operator%(const U &m) const {return matrix(*this) %= m;}};int main() {int n;cin >> n;matrix<modint<int>> v(3);v[0][1] = 1, v[1][0] = 1, v[1][2] = 1, v[2][0] = 1;auto v2 = v.pow(n-1);cout << (v2[0][0] + v2[1][0] + v2[2][0]).val << "\n";return 0;}