結果
問題 | No.1810 RGB Biscuits |
ユーザー |
![]() |
提出日時 | 2024-11-30 18:44:59 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 76 ms / 2,000 ms |
コード長 | 1,813 bytes |
コンパイル時間 | 1,999 ms |
コンパイル使用メモリ | 198,632 KB |
最終ジャッジ日時 | 2025-02-26 09:42:46 |
ジャッジサーバーID (参考情報) |
judge1 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 20 |
ソースコード
#include <bits/stdc++.h>#include <atcoder/modint>#define INF 4611686018427387903llusing namespace std;using namespace atcoder;using ll = long long;using mint = modint1000000007;struct Matrix22 {mint v11, v12, v21, v22;Matrix22(void){this->v11 = 1;this->v12 = 0;this->v21 = 0;this->v22 = 1;}Matrix22(mint v11, mint v12, mint v21, mint v22) {this->v11 = v11;this->v12 = v12;this->v21 = v21;this->v22 = v22;}Matrix22 operator*(const Matrix22 other) const {return Matrix22((this->v11 * other.v11) + (this->v12 * other.v21),(this->v11 * other.v12) + (this->v12 * other.v22),(this->v21 * other.v11) + (this->v22 * other.v21),(this->v21 * other.v12) + (this->v22 * other.v22));}pair<mint, mint> operator*(const pair<mint, mint> other) const {return pair<mint, mint>((this->v11 * other.first) + (this->v12 * other.second),(this->v21 * other.first) + (this->v22 * other.second));}Matrix22 pow(ll n){Matrix22 ret = Matrix22();Matrix22 A = *this;while(n){if(n & 1){ret = A * ret;}A = A * A;n >>= 1;}return ret;}};int main(){ll a, b; cin >> a >> b;ll n; cin >> n;Matrix22 A(a, b, 1, 0);pair<mint, mint> v(1, 1);for(int i = 0; i < n; ++i){ll t; cin >> t;ll d = t / 2, r = t % 2;pair<mint, mint> v1 = A.pow(d) * v;pair<mint, mint> v2 = A * v1;mint ans = 0;if(r == 1){ans = v1.first + v1.second + v2.first;}else{ans = v1.first + v1.second;}cout << ans.val() << endl;}}