結果
| 問題 | No.1073 無限すごろく |
| ユーザー |
 ganmodokix
|
| 提出日時 | 2020-06-06 09:23:13 |
| 言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
| 結果 |
AC
|
| 実行時間 | 2 ms / 2,000 ms |
| コード長 | 9,904 bytes |
| コンパイル時間 | 5,090 ms |
| コンパイル使用メモリ | 208,668 KB |
| 実行使用メモリ | 4,380 KB |
| 最終ジャッジ日時 | 2023-08-24 19:47:18 |
| 合計ジャッジ時間 | 6,079 ms |
|
ジャッジサーバーID （参考情報） |
judge11 / judge15 |
（要ログイン）
テストケース
テストケース表示| 入力 | 結果 | 実行時間 実行使用メモリ |
|---|---|---|
| testcase_00 | AC | 1 ms
4,376 KB |
| testcase_01 | AC | 1 ms
4,380 KB |
| testcase_02 | AC | 1 ms
4,376 KB |
| testcase_03 | AC | 1 ms
4,376 KB |
| testcase_04 | AC | 1 ms
4,380 KB |
| testcase_05 | AC | 2 ms
4,376 KB |
| testcase_06 | AC | 2 ms
4,380 KB |
| testcase_07 | AC | 1 ms
4,376 KB |
| testcase_08 | AC | 2 ms
4,376 KB |
| testcase_09 | AC | 1 ms
4,376 KB |
| testcase_10 | AC | 1 ms
4,380 KB |
| testcase_11 | AC | 1 ms
4,376 KB |
| testcase_12 | AC | 1 ms
4,380 KB |
| testcase_13 | AC | 2 ms
4,380 KB |
| testcase_14 | AC | 2 ms
4,376 KB |
| testcase_15 | AC | 2 ms
4,380 KB |
| testcase_16 | AC | 2 ms
4,380 KB |
| testcase_17 | AC | 1 ms
4,376 KB |
| testcase_18 | AC | 2 ms
4,380 KB |
| testcase_19 | AC | 1 ms
4,376 KB |
| testcase_20 | AC | 2 ms
4,380 KB |
| testcase_21 | AC | 2 ms
4,380 KB |
| testcase_22 | AC | 1 ms
4,376 KB |
| testcase_23 | AC | 2 ms
4,376 KB |
| testcase_24 | AC | 2 ms
4,380 KB |
| testcase_25 | AC | 2 ms
4,376 KB |
| testcase_26 | AC | 1 ms
4,380 KB |
| testcase_27 | AC | 1 ms
4,380 KB |
| testcase_28 | AC | 2 ms
4,380 KB |
| testcase_29 | AC | 1 ms
4,376 KB |
| testcase_30 | AC | 2 ms
4,376 KB |
| testcase_31 | AC | 1 ms
4,380 KB |
| testcase_32 | AC | 2 ms
4,380 KB |
ソースコード
// May this submission get accepted
#include <bits/stdc++.h>
// エイリアス
using ll = long signed long;
using ull = long unsigned long;
using ld = long double;
using namespace std;
// AtCoder/Codeforces 用 デバッグ検知
#ifdef ONLINE_JUDGE
constexpr bool DEBUG_MODE = false;
#else
constexpr bool DEBUG_MODE = true;
#endif
// エイリアス (補完・コンパイルが重くなる)
// #include <boost/multiprecision/cpp_int.hpp>
// using mll = boost::multiprecision::cpp_int;
// 汎用マクロ
#define ALL_OF(x) (x).begin(), (x).end()
#define REP(i,n) for (long long i=0, i##_len=(n); i<i##_len; i++)
#define RANGE(i,is,ie) for (long long i=(is), i##_end=(ie); i<=i##_end; i++)
#define DSRNG(i,is,ie) for (long long i=(is), i##_end=(ie); i>=i##_end; i--)
#define STEP(i, is, ie, step) for (long long i=(is), i##_end=(ie), i##_step = (step); i<=i##_end; i+=i##_step)
#define UNIQUE(v) do { sort((v).begin(), (v).end()); (v).erase(unique((v).begin(), (v).end()), (v).end()); } while (false)
#define FOREACH(i,q) for (auto &i : q)
template<class T> bool chmax(T &a, const T b) { if (a < b) {a = b; return true;} return false; }
template<class T> bool chmin(T &a, const T b) { if (a > b) {a = b; return true;} return false; }
constexpr int INF = numeric_limits<int>::max();
constexpr long long LINF = numeric_limits<long long>::max();
constexpr long double EPS = 1e-10L;
#define Yes(q) ((q) ? "Yes" : "No")
#define YES(q) ((q) ? "YES" : "NO")
#define Possible(q) ((q) ? "Possible" : "Impossible")
#define POSSIBLE(q) ((q) ? "POSSIBLE" : "IMPOSSIBLE")
#define IIF(q,t,f) ((q) ? (t) : (f))
#define DUMP(q) DUMP_FUNC(q, #q, __FILE__, __LINE__)
template <typename T> void DUMP_PROC(T x) { if (is_integral<T>() || is_floating_point<T>()) cerr << "\e[32m" << x << "\e[m"; else cerr << x; }
template<> void DUMP_PROC<char>(char x) { cerr << "\e[36m\'" << x << "\'\e[m"; }
template<> void DUMP_PROC<string>(string x) { cerr << "\e[33m\"" << x << "\"\e[m"; }
template <typename T, typename U> void DUMP_PROC(pair<T, U> x) { cerr << "{"; DUMP_PROC(x.first); cerr << ", "; DUMP_PROC(x.second); cerr << "}"; }
template <typename ...T, typename U, U... Seq> void DUMP_PROC(tuple<T...> &x, integer_sequence<U, Seq...>) { (void)(int[]){(cerr << ((const char*[]){"", ", "})[!!Seq] << (DUMP_PROC(get<Seq>(x)), ""), 0)...}; }
template <typename ...T> void DUMP_PROC(tuple<T...> x) {cerr << "{"; DUMP_PROC(x, index_sequence_for<T...>()); cerr << "}";}
template <typename T> void DUMP_PROC(vector<T> x) { cerr << "["; for (auto &xi : x) { DUMP_PROC(xi); cerr << (&xi != &*x.rbegin()?", ":""); } cerr << "]"; }
template <typename T> void DUMP_FUNC(T x, const char* name, const char* fn, int ln) { cerr << "\e[32m[DEBUG]\e[m " << name << ": "; DUMP_PROC(x); cerr << " @ " << fn << "(" << ln << ")" << endl; }
// gcc拡張マクロ
#define popcount __builtin_popcount
#define popcountll __builtin_popcountll
// 標準入出力
struct qin { // query input
size_t sz;
qin(size_t _sz = 1) : sz(_sz) {}
template <typename T> operator T () const { T a; cin >> a; return a; }
template <typename T> operator vector<T> () const { vector<T> a(sz); for (size_t i = 0; i < sz; i++) cin >> a[i]; return a; }
template <typename T, typename U> operator pair<T, U> () const { T f; U s; cin >> f >> s; return pair<T, U>(f, s); }
};
qin in1; // input one
template <typename T> void say(const T x, const char* end = "\n") { cout << x << end; }
void say(const ld x, const char* end = "\n") { cout << setprecision(30) << x << end; }
template <typename T> void say(const vector<T> x, const char* sep = " ", const char* end = "\n") { REP(i, x.size()) { cout << x[i] << (i+1 == i_len ? end : sep); } }
template <typename T> void say(const vector<vector<T>> x, const char* sep = " ", const char* end = "\n") { REP(i, x.size()) { say(x[i], sep, end); } }
// モジュール
// [[LIBRARY]]
// 行列演算
// pdivの余りを取りたくない場合はoperator+=/-=/*のコメントアウトを適宜
constexpr ll pdiv = 1000000007LL;
template <typename T = ll>
class matrix {
private:
size_t col_size, row_size;
vector<T> raw;
public:
matrix() : col_size(0), row_size(0), raw() {}
matrix(size_t n) : col_size(1), row_size(n), raw(n, 0) {}
matrix(size_t n, size_t m) : col_size(m), row_size(n), raw(n * m, 0) {}
matrix(const vector<vector<T>> &a) : col_size(a.size() ? a[0].size() : 0), row_size(a.size()), raw(col_size * row_size) {
for (size_t i = 0; i < row_size; i++) {
for (size_t j = 0; j < col_size; j++) {
at(i, j) = a[i][j];
}
}
}
matrix(initializer_list<T> init) : col_size(init.size()), row_size(1), raw(init) {}
matrix(initializer_list<initializer_list<T>> init) : col_size(0), row_size(0), raw() {
const size_t n = init.size();
const size_t m = [&]() {
size_t maxm = 0;
for (const auto& initi : init) {
const size_t mi = initi.size();
if (mi > maxm) maxm = mi;
}
return maxm;
}();
tie(row_size, col_size) = make_pair(n, m);
raw.resize(n * m);
size_t i = 0, j = 0;
for (const auto &initi : init) {
for (const auto &v : initi) {
at(i, j++) = v;
}
i++; j = 0;
}
}
size_t rows() const { return row_size; }
size_t cols() const { return col_size; }
const T& at(const size_t i, const size_t j) const { return raw.at(i * col_size + j); }
const T& at(const size_t i) const { return raw.at(i); }
T& at(const size_t i, const size_t j) { return raw.at(i * col_size + j); }
T& at(const size_t i) { return raw.at(i); }
T& operator() (size_t i, size_t j) { return raw[i * col_size + j]; }
T& operator() (size_t i) { return raw[i]; }
T& operator[] (size_t i) { return raw[i]; }
const T& operator() (size_t i, size_t j) const { return raw[i * col_size + j]; }
const T& operator() (size_t i) const { return raw[i]; }
const T& operator[] (size_t i) const { return raw[i]; }
matrix<T>& operator+= (const matrix<T> &a) {
for (size_t i = 0; i < row_size; i++) {
for (size_t j = 0; j < col_size; j++) {
// at(i, j) += a.at(i, j);
if ((at(i, j) += a.at(i, j)) >= pdiv) at(i, j) -= pdiv;
}
}
return *this;
}
matrix<T>& operator-= (const matrix<T> &a) {
for (size_t i = 0; i < row_size; i++) {
for (size_t j = 0; j < col_size; j++) {
// at(i, j) -= a.at(i, j);
if ((at(i, j) += pdiv - a.at(i, j)) >= pdiv) at(i, j) -= pdiv;
}
}
return *this;
}
matrix<T> operator* (const matrix<T> &a) {
assert(col_size == a.row_size);
const size_t n = row_size, m = col_size, q = a.col_size;
matrix<T> r(n, q);
for (size_t i = 0; i < n; i++) {
for (size_t k = 0; k < m; k++) {
for (size_t j = 0; j < q; j++) {
// r.at(i, j) += at(i, k) * a.at(k, j);
(r.at(i, j) += at(i, k) * a.at(k, j) % pdiv) %= pdiv;
}
}
}
return r;
}
matrix<T> operator+ (const matrix<T> &a) const {
matrix<T> r = *this;
return r += a;
}
matrix<T> operator- (const matrix<T> &a) const {
matrix<T> r = *this;
return r -= a;
}
matrix<T>& operator*= (const matrix<T> &a) {
return *this = move(*this * a);
}
operator T () const {
assert(col_size == 1 && row_size == 1);
return raw[0];
}
static matrix<T> eye(size_t n) {
matrix<T> r(n, n);
for (size_t i = 0; i < n; i++)
r.at(i, i) = 1;
return r;
}
matrix<T> transpose() const {
matrix<T> r(col_size, row_size);
for (size_t i = 0; i < col_size; i++) {
for (size_t j = 0; j < row_size; j++) {
r.at(i, j) = at(j, i);
}
}
return r;
}
matrix<T> pow(long long int n) const {
assert(row_size == col_size);
matrix<T> r = eye(row_size);
matrix<T> a = *this;
while (n) {
if (n & 1) r *= a;
n /= 2;
if (n) a *= a;
}
return r;
}
};
template <typename T>
void DUMPMTR(matrix<T> m) {
cerr << "[DEBUG] DUMPMTR: " << endl;
for (size_t i = 0; i < m.rows(); i++) {
cerr << " ["[!i] << "[";
for (size_t j = 0; j < m.cols(); j++) {
cerr << setw(4) << m(i, j) << " ]"[j+1 == m.cols()];
}
if (i+1 == m.rows()) cerr << ']';
cerr << endl;
}
}
// 剰余演算ライブラリ (小): 剰余と逆元だけ欲しいときに使う軽量ライブラリ
// mod pdiv 上での a^n
ll modpow(const ll a, ll n, const ll pdiv) {
if (a == 0) return 0;
ll r = 1;
((n %= pdiv-1) += pdiv-1) %= pdiv-1; // a^(p-1) === a^0 mod pdiv
for (ll b = (a % pdiv + pdiv) % pdiv; n; n >>= 1, (b *= b) %= pdiv) if (n & 1) (r *= b) %= pdiv;
return r;
}
inline ll modinv(const ll a, const ll pdiv) { return modpow(a, pdiv-2, pdiv); }
// 処理内容
int main() {
ios::sync_with_stdio(false); // stdioを使うときはコメントアウトすること
cin.tie(nullptr); // インタラクティブ問題ではコメントアウトすること
ll n = in1;
ll sixth = modinv(6, pdiv);
matrix<ll> m = {
{sixth, sixth, sixth, sixth, sixth, sixth},
{1LL, 0LL, 0LL, 0LL, 0LL, 0LL},
{0LL, 1LL, 0LL, 0LL, 0LL, 0LL},
{0LL, 0LL, 1LL, 0LL, 0LL, 0LL},
{0LL, 0LL, 0LL, 1LL, 0LL, 0LL},
{0LL, 0LL, 0LL, 0LL, 1LL, 0LL}
};
matrix<ll> x0 = matrix<ll>{1LL, 0LL, 0LL, 0LL, 0LL, 0LL}.transpose();
say((m.pow(n) * x0).at(0));
}