結果

問題 No.2062 Sum of Subset mod 999630629
ユーザー ygussanyygussany
提出日時 2022-08-21 18:53:12
言語 C
(gcc 12.3.0)
結果
RE  
実行時間 -
コード長 1,549 bytes
コンパイル時間 226 ms
コンパイル使用メモリ 32,000 KB
実行使用メモリ 14,576 KB
最終ジャッジ日時 2024-10-13 12:38:32
合計ジャッジ時間 8,105 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 2 ms
5,248 KB
testcase_05 AC 2 ms
5,248 KB
testcase_06 AC 2 ms
5,248 KB
testcase_07 AC 2 ms
5,248 KB
testcase_08 AC 22 ms
5,376 KB
testcase_09 AC 19 ms
5,248 KB
testcase_10 AC 17 ms
5,248 KB
testcase_11 RE -
testcase_12 RE -
testcase_13 WA -
testcase_14 RE -
testcase_15 WA -
testcase_16 RE -
testcase_17 RE -
testcase_18 TLE -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>

const int Mod = 998244353, Modd = 999630629;

long long pow_mod(int n, long long k)
{
	long long N, ans = 1;
	for (N = n; k > 0; k >>= 1, N = N * N % Mod) if (k & 1) ans = ans * N % Mod;
	return ans;
}

void merge_sort(int n, int x[])
{
	static int y[100001] = {};
	if (n <= 1) return;
	merge_sort(n / 2, &(x[0]));
	merge_sort((n + 1) / 2, &(x[n/2]));
	
	int i, p, q;
	for (i = 0, p = 0, q = n / 2; i < n; i++) {
		if (p >= n / 2) y[i] = x[q++];
		else if (q >= n) y[i] = x[p++];
		else y[i] = (x[p] < x[q])? x[p++]: x[q++];
	}
	for (i = 0; i < n; i++) x[i] = y[i];
}

int main()
{
	int i, N, A[100001], sum = 0;
	scanf("%d", &N);
	for (i = 0; i < N; i++) {
		scanf("%d", &(A[i]));
		sum += A[i];
	}
	merge_sort(N, A);
	
	int j, tmp, cur, prev, dp[2][370000] = {}, q[2][370000], head, tail, M = sum - Modd;
	long long ans = sum * pow_mod(2, N - 1) % Mod;
	if (M >= 0) {
		q[0][0] = 0;
		dp[0][0] = 1;
		tail = 1;
		for (i = N - 1, cur = 1, prev = 0; i >= 0; i--, cur ^= 1, prev ^= 1) {
			for (head = tail - 1, tail = 0; head >= 0; head--) {
				j = q[prev][head];
				tmp = dp[prev][j];
				dp[prev][j] = 0;
				
				if (dp[cur][j] == 0) q[cur][tail++] = j;
				dp[cur][j] += tmp;
				if (dp[cur][j] >= Mod) dp[cur][j] -= Mod;
				j += A[i];
				if (j <= M) {
					if (dp[cur][j] == 0) q[cur][tail++] = j;
					dp[cur][j] += tmp;
					if (dp[cur][j] >= Mod) dp[cur][j] -= Mod;
				}
				dp[prev][j] = 0;
			}
		}
		for (j = 0; j <= M; j++) ans += Mod - dp[prev][j];
	}
	printf("%lld\n", ans % Mod);
	fflush(stdout);
	return 0;
}
0