結果

問題 No.2669 Generalized Hitting Set
コンテスト
ユーザー suisen
提出日時 2024-01-14 04:57:16
言語 C
(gcc 13.3.0)
結果
MLE  
実行時間 -
コード長 3,196 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,251 ms
コンパイル使用メモリ 44,320 KB
実行使用メモリ 468,748 KB
最終ジャッジ日時 2024-09-28 02:54:54
合計ジャッジ時間 23,099 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other TLE * 1 MLE * 30 -- * 10
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")

#include <stdio.h>

#define MAX_N 24
#define HALF MAX_N / 2

int res[1 << (MAX_N - 1)][HALF];
int f[1 << MAX_N] = {};

int S[300000] = {};
char buf[MAX_N + 1];

int G[300001] = {};

int main() {
    int n, m, k;
    scanf("%d %d %d", &n, &m, &k);

    for (int i = 0; i < m; ++i) {
        scanf("%s", buf);
        for (int j = 0; j < n; ++j) {
            S[i] |= (buf[j] - '0') << j;
        }
    }

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

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

        for (int topbit = 0; topbit < 2; ++topbit) {
            for (int i = 0; i < 1 << (MAX_N - 1); ++i) for (int j = 0; j < HALF; ++j) res[i][j] = 0;

            for (int i = 0; i < m; ++i) if ((S[i] >> (n - 1)) == topbit) ++res[S[i] & ~(1 << (n - 1))][0];

            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);
                        }
                    }
                }
            }
            for (int T = 0; T < 1 << n; ++T) {
                int x = (T >> (n - 1)) == 1 && topbit == 1;
                for (int i = 0; i < k - x; ++i) {
                    f[T] -= res[T & ~(1 << (n - 1))][i];
                }
            }
        }
    } else {
        const int r = n - k + 1;
        /*
        [x,x]
        [x,1]
        */

        for (int topbit = 0; topbit < 2; ++topbit) {
            for (int i = 0; i < 1 << (MAX_N - 1); ++i) for (int j = 0; j < HALF; ++j) res[i][j] = 0;
        
            for (int i = 0; i < m; ++i) if ((S[i] >> (n - 1)) == topbit) ++res[S[i] & ~(1 << (n - 1))][0];

        
            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;
                    }
                }
            }
            for (int T = 0; T < 1 << n; ++T) {
                int x = (T >> (n - 1)) == 0 || topbit == 0;
                for (int i = 0; i < r - x; ++i) {
                    f[T] += res[T & ~(1 << (n - 1))][i];
                }
            }
        }
    }

    for (int i = 0; i <= m; ++i) {
        G[i] = n;
    }
    for (int T = 0; T < 1 << n; ++T) {
        int T_size = __builtin_popcount(T);
        if (T_size < G[f[T]]) {
            G[f[T]] = T_size;
        }
    }

    for (int i = m - 1; i >= 0; --i) {
        if (G[i + 1] < G[i]) G[i] = G[i + 1];
    }

    for (int i = 1; i <= m; ++i) printf("%d\n", G[i]);
}
0