結果

問題 No.2441 行列累乗
コンテスト
ユーザー kw_c
提出日時 2023-09-16 08:07:46
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 4,594 bytes
コンパイル時間 2,701 ms
コンパイル使用メモリ 202,360 KB
最終ジャッジ日時 2025-02-16 22:56:45
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define rrep(i,n) for (int i = n-1; i >= (0); --i)
#define IO ios::sync_with_stdio(0),cin.tie(0);
#define fi first
#define se second
#define all(x) x.begin(), x.end()
#define rall(x) x.rbegin(), x.rend()
#define pb push_back
#define pf push_front
#define eb emplace_back
#define ef emplace_front
#define bpop(x) __builtin_popcount(x)
#define bpopll(x) __builtin_popcountll(x)


template<class T>
bool chmin(T &a, const T &b) { if (a > b) { a = b; return true; } return false; }
template<class T>
bool chmax(T &a, const T &b) { if (a < b) { a = b; return true; } return false; }
void print_imp() {}
template<class First, class... Rest>
void print_imp(First f, Rest... r) { cout<<f; if (sizeof...(r)) cout<<" "; print_imp(r...); }
template<class... Args>
void print(Args... args) { print_imp(args...); cout << endl; }
template<class T>
void print_pair(vector<T> &v) { for (auto [a,b] : v) print(a, b); }
template<class T>
void print_vec(vector<T> &v) { rep(i, v.size()) {cout << v[i]; if (i == v.size()-1) cout << endl; else cout << " ";} }
template<class T>
void print_mat(vector<vector<T>> &m) { for (auto v : m) print_vec(v); }
template<class T>
void print_bit(T s, int n) { rep(i, n) {if (s>>i&1) cout<<1; else cout<<0;} cout<<endl; }
void print_double(double x) { cout << fixed << setprecision(15) << x << endl; }

template<class T>
using min_pq = priority_queue<T, vector<T>, greater<T>>;
template<class T>
using vt = vector<T>;
template<class T>
using vvt = vt<vt<T>>;
using ll = long long;
using pii = pair<int, int>;
using pil = pair<int, ll>;
using pli = pair<ll, int>;
using pll = pair<ll, ll>;
using pdd = pair<double, double>;

const ll MOD = 1e9 + 7;
const int INF = 1001001001;
const ll LINF = 4e18;
const double PI = 3.14159265358979323846;

/*
- comb                          modintのnCk,nPk,nHk
- matrix                        行列累乗
- Mo                            Mo's algorithm
- segtree_min/max               segtree 点更新
- RMQ/RMAQ_min/max              lazy_segtree 区間加算/変更、区間最小値/最大値取得
- RSQ/RSAQ                      lazy_segtree 区間加算/変更、区間和取得
- Eratosthenes                  高速素因数分解可能なエラトステネスのふるい
- Factorize                     素因数分解 計算量: O(sqrt(x)) 
- dydx, dydx8                   dy,dx={右,下,左,上}
- grid_continue                 if (i < 0 || j < 0 || i >= h || j >= w) continue;
- Pascal                        パスカルの三角形
- Graph                         グラフ関連のライブラリ
    - dijkstra                  ダイクストラ法
    - bellman_ford              ベルマンフォード、負のサイクル検出
    - topological_sort          トポロジカルソート
    - find_cycle_directed       有向グラフのサイクル検出
    - bfs01                     01-BFS
    - tree_diameter             木の直径
    - euler_tour                オイラーツアー
    - build_LCA, get_LCA        LCA
    - build_tree_dist           木の根からの距離
    - dfs_temp                  DFSのためのテンプレート

    - init_tree_from_parents    親頂点配列でGraph初期化
    - init_graph_from_grid      グリッド文字列からGraph初期化
*/

template<typename T>
struct Matrix {
    int h, w;
    vector<vector<T>> d;
    Matrix() {}
    Matrix(int h, int w, T val=0) : h(h), w(w), d(h, vector<T>(w,val)) {}
    Matrix& unit() {
        assert(h == w);
        rep(i, h) d[i][i] = 1;
        return *this;
    }
    const vector<T>& operator[](int i) const { return d[i]; }
    vector<T>& operator[](int i) { return d[i]; }
    Matrix operator*(const Matrix& a) const {
        assert(w == a.h);
        Matrix r(h, a.w);
        rep(i, h) rep(k, w) rep(j, a.w) {
            r[i][j] += d[i][k] * a[k][j];
        }
        return r;
    }
    Matrix pow(ll t) const {
        assert(h == w);
        if (!t) return Matrix(h, h).unit();
        if (t == 1) return *this;
        Matrix r = pow(t>>1);
        r = r*r;
        if (t&1) r = r*(*this);
        return r;
    }
};

void solve() {
    int a, b, c, d;
    cin >> a >> b >> c >> d;
    Matrix<int> x(2, 2);
    x[0][0] = a;
    x[0][1] = b;
    x[1][0] = c;
    x[1][1] = d;
    
    auto y = x.pow(3);
    rep(i, 2) {
        rep(j, 2) {
            cout << y[i][j] << " ";
        }
        cout << endl;
    }
}

int main() {
    IO
    /* int t; cin >> t; rep(_, t) solve(); */
    solve();
    return 0;
}
0