結果
問題 | No.695 square1001 and Permutation 4 |
ユーザー |
![]() |
提出日時 | 2018-05-26 17:21:26 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 3,386 ms / 7,000 ms |
コード長 | 1,265 bytes |
コンパイル時間 | 905 ms |
コンパイル使用メモリ | 71,948 KB |
実行使用メモリ | 42,540 KB |
最終ジャッジ日時 | 2024-07-22 19:14:26 |
合計ジャッジ時間 | 20,373 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 12 |
ソースコード
#include <iostream>#include <vector>using namespace std;int N, M, X[8], dp[10000009]; vector<int>J;long long solve(int mod) {int G = N / 2 + 1;for (int i = 0; i <= G; i++) dp[i] = 0;dp[0] = 1;for (int i = 1; i <= N; i++) {dp[i % G] = 0;for (int j = 0; j < J.size(); j++) { dp[i % G] += (dp[(i - J[j] + G) % G]); dp[i % G] %= mod; }}long long ret = dp[N%G];for (int i = 0; i <= G; i++) dp[i] = 0;dp[0] = 1;for (int i = 1; i <= G; i++) {dp[i] = 0;for (int j = 0; j < J.size(); j++) { dp[i] += (dp[(i - J[j] + G) % G]); dp[i] %= mod; }}for (int i = 1; i <= M; i++) {long long G = X[i]; if (X[i] <= N / 2) continue;for (int j = 0; j <= N - G; j++) { ret += 1LL * dp[j] * dp[(N - G) - j]; ret %= mod; }}return ret;}int main() {cin >> N >> M; N--;for (int i = 1; i <= M; i++) { cin >> X[i]; if (X[i] <= N / 2) J.push_back(X[i]); }long long G1 = solve(17);long long G2 = solve(9920467);long long G3 = solve(592951213);long long R = G3;for (int i = 0; i < 9920467; i++) { long long B = 592951213LL * i + R; if (B % 9920467LL == G2) { R = B; break; } }for (int i = 0; i < 17; i++) { long long B = (592951213LL * 9920467LL) * i + R; if (B % 17LL == G1) { R = B; break; } }cout << R << endl;return 0;}