結果
問題 |
No.444 旨味の相乗効果
|
ユーザー |
|
提出日時 | 2016-12-27 22:09:19 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 4 ms / 2,500 ms |
コード長 | 897 bytes |
コンパイル時間 | 1,590 ms |
コンパイル使用メモリ | 166,980 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-06-24 21:00:26 |
合計ジャッジ時間 | 2,495 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 23 |
ソースコード
#include <bits/stdc++.h> using namespace std; const int mod = 1e9 + 7; int modpow(int a, long long b) { int ret = 1; while (b > 0) { if (b & 1) { ret = 1LL * ret * a % mod; } a = 1LL * a * a % mod; b >>= 1; } return ret; } int main() { int n; long long c; cin >> n >> c; long long foo = 0; static long long a[80]; for (int i = 0; i < n; i++) { cin >> a[i]; foo += mod - modpow(a[i], c); } // doubling // ref: https://yukicoder.me/problems/no/137 static long long dp[160]; dp[0] = 1; while (c > 0) { for (int i = 0; i < 80; i++) { for (int j = 158; j >= 0; j--) { (dp[j + 1] += dp[j] * a[i]) %= mod; } } for (int j = (c & 1); j < 160; j += 2) { dp[j / 2] = dp[j]; } for (int j = 80; j < 160; j++) { dp[j] = 0; } for (int i = 0; i < 80; i++) { a[i] = a[i] * a[i] % mod; } c >>= 1; } cout << (dp[0] + foo) % mod << endl; }