結果
問題 | No.1844 Divisors Sum Sum |
ユーザー |
|
提出日時 | 2022-02-19 03:23:04 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 791 ms / 3,000 ms |
コード長 | 1,214 bytes |
コンパイル時間 | 2,328 ms |
コンパイル使用メモリ | 196,112 KB |
最終ジャッジ日時 | 2025-01-28 00:57:44 |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 38 |
ソースコード
#include <bits/stdc++.h> using namespace std; const long long MOD = 1000000007; array<array<long long, 3>, 3> matmul(array<array<long long, 3>, 3> A, array<array<long long, 3>, 3> B){ array<array<long long, 3>, 3> ans = {array<long long, 3>{0, 0, 0}, array<long long, 3>{0, 0, 0}, array<long long, 3>{0, 0, 0}}; for (int i = 0; i < 3; i++){ for (int j = 0; j < 3; j++){ for (int k = 0; k < 3; k++){ ans[i][k] += A[i][j] * B[j][k]; ans[i][k] %= MOD; } } } return ans; } array<array<long long, 3>, 3> matexp(array<array<long long, 3>, 3> A, long long b){ array<array<long long, 3>, 3> ans = {array<long long, 3>{1, 0, 0}, array<long long, 3>{0, 1, 0}, array<long long, 3>{0, 0, 1}}; while (b > 0){ if (b % 2 == 1){ ans = matmul(ans, A); } A = matmul(A, A); b /= 2; } return ans; } int main(){ int L; cin >> L; vector<int> P(L), e(L); for (int i = 0; i < L; i++){ cin >> P[i] >> e[i]; } long long ans = 1; for (int i = 0; i < L; i++){ array<array<long long, 3>, 3> A = {array<long long, 3>{P[i], 1, 1}, array<long long, 3>{0, 1, 1}, array<long long, 3>{0, 0, 1}}; A = matexp(A, e[i] + 1); ans *= A[0][2]; ans %= MOD; } cout << ans << endl; }