結果
| 問題 |
No.444 旨味の相乗効果
|
| コンテスト | |
| ユーザー |
square1001
|
| 提出日時 | 2018-05-19 21:19:40 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 42 ms / 2,500 ms |
| コード長 | 1,070 bytes |
| コンパイル時間 | 531 ms |
| コンパイル使用メモリ | 65,976 KB |
| 実行使用メモリ | 5,376 KB |
| 最終ジャッジ日時 | 2024-06-28 14:29:17 |
| 合計ジャッジ時間 | 1,791 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 5 |
| other | AC * 23 |
ソースコード
#include <iostream>
using namespace std;
const int mod = 1000000007;
int n, a[85], x[85][85], y[85][85], z[85][85]; long long b;
int main() {
cin >> n >> b; b--;
for (int i = 0; i < n; i++) {
cin >> a[i];
x[i][i] = a[i];
for (int j = 0; j <= i; j++) y[j][i] = a[i];
}
while (b) {
if (b & 1) {
for (int i = 0; i < n; i++) {
for (int j = i; j < n; j++) {
z[i][j] = 0;
for (int k = i; k <= j; k++) {
z[i][j] = (z[i][j] + 1LL * x[i][k] * y[k][j]) % mod;
}
}
}
for (int i = 0; i < n; i++) {
for (int j = i; j < n; j++) {
x[i][j] = z[i][j];
}
}
}
for (int i = 0; i < n; i++) {
for (int j = i; j < n; j++) {
z[i][j] = 0;
for (int k = i; k <= j; k++) {
z[i][j] = (z[i][j] + 1LL * y[i][k] * y[k][j]) % mod;
}
}
}
for (int i = 0; i < n; i++) {
for (int j = i; j < n; j++) {
y[i][j] = 1LL * z[i][j];
}
}
b >>= 1;
}
int ret = 0;
for (int i = 0; i < n; i++) {
for (int j = i + 1; j < n; j++) {
ret = (ret + x[i][j]) % mod;
}
}
cout << ret << '\n';
return 0;
}
square1001