結果

問題 No.695 square1001 and Permutation 4
ユーザー e869120
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
}
0