結果

問題 No.2669 Generalized Hitting Set
ユーザー suisensuisen
提出日時 2024-01-14 02:30:04
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
MLE  
実行時間 -
コード長 5,489 bytes
コンパイル時間 1,047 ms
コンパイル使用メモリ 78,172 KB
実行使用メモリ 528,556 KB
最終ジャッジ日時 2024-09-28 02:39:13
合計ジャッジ時間 55,841 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 2 ms
6,944 KB
testcase_03 AC 18 ms
6,944 KB
testcase_04 AC 16 ms
6,940 KB
testcase_05 AC 13 ms
6,940 KB
testcase_06 AC 22 ms
6,940 KB
testcase_07 AC 5 ms
6,940 KB
testcase_08 AC 6 ms
6,944 KB
testcase_09 AC 16 ms
6,944 KB
testcase_10 AC 9 ms
6,944 KB
testcase_11 AC 5 ms
6,944 KB
testcase_12 AC 19 ms
6,940 KB
testcase_13 MLE -
testcase_14 RE -
testcase_15 AC 4 ms
7,684 KB
testcase_16 AC 48 ms
23,912 KB
testcase_17 RE -
testcase_18 AC 2 ms
6,940 KB
testcase_19 AC 33 ms
13,556 KB
testcase_20 AC 2 ms
6,944 KB
testcase_21 AC 30 ms
23,932 KB
testcase_22 RE -
testcase_23 AC 43 ms
10,588 KB
testcase_24 AC 22 ms
6,940 KB
testcase_25 RE -
testcase_26 MLE -
testcase_27 RE -
testcase_28 MLE -
testcase_29 MLE -
testcase_30 RE -
testcase_31 RE -
testcase_32 RE -
testcase_33 AC 47 ms
14,544 KB
testcase_34 AC 29 ms
8,184 KB
testcase_35 RE -
testcase_36 AC 2 ms
6,940 KB
testcase_37 AC 71 ms
39,136 KB
testcase_38 AC 5 ms
9,464 KB
testcase_39 MLE -
testcase_40 MLE -
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <cassert>
#include <iostream>
#include <string>
#include <vector>

constexpr int MAX_N = 24;
constexpr int HALF = MAX_N / 2;

constexpr int MAX_SIZE = HALF;

int res[1 << (MAX_N - 1)][MAX_SIZE]{};
int f[1 << MAX_N]{};
int cnt[1 << MAX_N]{};

int main() {
    std::ios::sync_with_stdio(false);
    std::cin.tie(nullptr);

    int n, m, k;
    std::cin >> n >> m >> k;

    for (int i = 0; i < m; ++i) {
        std::string S;
        std::cin >> S;
        int val = 0;
        for (int j = 0; j < n; ++j) val |= ((S[j] - '0') << j);
        ++cnt[val];
    }

    if (n <= HALF) {
        for (int i = 0; i < 1 << n; ++i) {
            res[i][0] = cnt[i];
        }
        /*
        [1,1]
        [1,x]
        */
        for (int B = 1; B < 1 << n; B <<= 1) {
            for (int L = 0; L < 1 << n; L += 2 * B) {
                for (int p = 0; p < B; ++p) {
                    int* u = res[L + p], * v = res[L + B + p];
                    for (int t = k; t--;) {
                        int tmp = u[t];
                        u[t] = tmp + v[t];
                        v[t] = tmp + (t ? v[t - 1] : 0);
                    }
                }
            }
        }

        std::vector<int> g(m + 1, n);
        for (int T = 0; T < 1 << n; ++T) {
            int T_size = __builtin_popcount(T);
            int f_T = m;
            for (int i = 0; i < k; ++i) f_T -= res[T][i];
            g[f_T] = std::min(g[f_T], T_size);
        }

        for (int i = m - 1; i >= 0; --i) g[i] = std::min(g[i], g[i + 1]);

        for (int i = 1; i <= m; ++i) std::cout << g[i] << '\n';
        return 0;
    }

    if (k <= HALF) {
        /*
        [1,1]
        [1,x]
        */

        for (int i = 0; i < 1 << n; ++i) {
            f[i] = m;
        }

        for (int topbit : { 0, 1 }) {
            for (int i = 0; i < 1 << (n - 1); ++i) {
                for (int j = 0; j < HALF; ++j) res[i][j] = 0;
                res[i][0] = cnt[i | (topbit << (n - 1))];
            }
            for (int B = 1; B < 1 << (n - 1); B <<= 1) {
                for (int L = 0; L < 1 << (n - 1); L += 2 * B) {
                    for (int p = 0; p < B; ++p) {
                        int* u = res[L + p], * v = res[L + B + p];
                        for (int t = k; t--;) {
                            int tmp = u[t];
                            u[t] = tmp + v[t];
                            v[t] = tmp + (t ? v[t - 1] : 0);
                        }
                    }
                }
            }
            if (topbit == 0) {
                for (int T = 0; T < 1 << n; ++T) {
                    for (int i = 0; i < k; ++i) {
                        f[T] -= res[T & ~(1 << (n - 1))][i];
                    }
                }
            } else {
                for (int T = 0; T < 1 << n; ++T) {
                    if (T >> (n - 1)) {
                        for (int i = 0; i < k - 1; ++i) {
                            f[T] -= res[T & ~(1 << (n - 1))][i];
                        }
                    } else {
                        for (int i = 0; i < k; ++i) {
                            f[T] -= res[T][i];
                        }
                    }
                }
            }
        }

        std::vector<int> g(m + 1, n);
        for (int T = 0; T < 1 << n; ++T) {
            int T_size = __builtin_popcount(T);
            g[f[T]] = std::min(g[f[T]], T_size);
        }

        for (int i = m - 1; i >= 0; --i) g[i] = std::min(g[i], g[i + 1]);

        for (int i = 1; i <= m; ++i) std::cout << g[i] << '\n';
    } else {
        assert(false);

        const int r = n - k + 1;
        assert(r <= MAX_SIZE);
        /*
        [x,x]
        [x,1]
        */

        for (int topbit : { 0, 1 }) {
            for (int i = 0; i < 1 << (n - 1); ++i) {
                for (int j = 0; j < HALF; ++j) res[i][j] = 0;
                res[i][0] = cnt[i | (topbit << (n - 1))];
            }
            for (int B = 1; B < 1 << (n - 1); B <<= 1) {
                for (int L = 0; L < 1 << (n - 1); L += 2 * B) {
                    for (int p = 0; p < B; ++p) {
                        int* u = res[L + p], * v = res[L + B + p];
                        for (int t = r; t-- > 1;) {
                            u[t] = u[t - 1] + v[t - 1];
                            v[t] += u[t - 1];
                        }
                        u[0] = 0;
                    }
                }
            }
            if (topbit == 0) {
                for (int T = 0; T < 1 << n; ++T) {
                    for (int i = 0; i < r - 1; ++i) {
                        f[T] += res[T & ~(1 << (n - 1))][i];
                    }
                }
            } else {
                for (int T = 0; T < 1 << n; ++T) {
                    if (T >> (n - 1)) {
                        for (int i = 0; i < r - 1; ++i) {
                            f[T] += res[T & ~(1 << (n - 1))][i];
                        }
                    } else {
                        for (int i = 0; i < r; ++i) {
                            f[T] += res[T][i];
                        }
                    }
                }
            }
        }

        std::vector<int> g(m + 1, n);
        for (int T = 0; T < 1 << n; ++T) {
            int T_size = __builtin_popcount(T);
            g[f[T]] = std::min(g[f[T]], T_size);
        }

        for (int i = m - 1; i >= 0; --i) g[i] = std::min(g[i], g[i + 1]);

        for (int i = 1; i <= m; ++i) std::cout << g[i] << '\n';
    }
}
0