[CF1290C] Prefix Enlightenment:并查集

题意

题目链接

给定一个 n n n 位 01 串 S S S

给定 k k k 个集合 A 1 , A 2 , ⋯   , A k A_1, A_2, \cdots, A_k A1,A2,,Ak,任意三个集合满足 A i ∩ A j ∩ A k = ∅ A_i\cap A_j\cap A_k=\empty AiAjAk=

一次操作被定义为选择一个集合 A i A_i Ai,对集合中的每个元素 x x x 执行 S x ← S x ⊕ 1 S_x\leftarrow S_x\oplus 1 SxSx1。保证能通过若干次操作使得 S S S 变为全 1 串。

定义 m i m_i mi 为使 S S S 的前 i i i 个位置都变成 1 的最小操作次数,我们需要求出 m 1 , m 2 , ⋯   , m n m_1, m_2, \cdots, m_n m1,m2,,mn

数据范围: 1 ≤ n , k ≤ 3 × 1 0 5 1\le n,k\le 3\times 10^5 1n,k3×105

题解

首先注意到题目竟然没有给 ∑ ∣ A i ∣ \sum |A_i| Ai 的数据范围,这让我们感到疑惑这题出错了,于是我们考虑任意三个集合交集为空这个莫名其妙的条件。稍微思考一下即可发现,它的意思是对 S S S 中任意一个位置,它最多只受两个集合的控制。所以显然 ∑ ∣ A i ∣ = O ( n ) \sum |A_i|=O(n) Ai=O(n)

知道了这件事情以后,我们就需要考虑从每一位的角度看待问题了:

先设 x 1 , x 2 , ⋯   , x k x_1, x_2, \cdots, x_k x1,x2,,xk 表示集合 A 1 , A 2 , ⋯   , A k A_1, A_2, \cdots, A_k A1,A2,,Ak 的状态(选,不选,是 01 变量)

考虑 S i S_i Si

  • i i i 不受集合约束,不管(因为保证了答案存在)
  • i i i 只受一个集合 a a a 约束,则 x a = 1 ⊕ S i x_a=1\oplus S_i xa=1Si
  • i i i 受两个集合 a , b a,b a,b 约束,则 x a ⊕ x b = S i x_a\oplus x_b=S_i xaxb=Si,即 S i = 0 S_i=0 Si=0 x a ≠ x b x_a\not= x_b xa=xb S i = 1 S_i=1 Si=1 x a = x b x_a=x_b xa=xb

这个给定等于和不等于的条件很容易让我们想到 关押罪犯 那套模型啊。。就是并查集拆点,我们将每个集合 i i i 的选项分为两个点 x i , 0 , x i , 1 x_{i,0},x_{i,1} xi,0,xi,1 分别表示两种选择,则我们的操作是:强制选点,连接 ( x a , 0 , x b , 0 ) , ( x a , 1 , x b , 1 ) (x_{a,0},x_{b,0}),(x_{a,1},x_{b,1}) (xa,0,xb,0),(xa,1,xb,1) 或是连接 ( x a , 0 , x b , 1 ) , ( x a , 1 , x b , 0 ) (x_{a,0},x_{b,1}),(x_{a,1},x_{b,0}) (xa,0,xb,1),(xa,1,xb,0)

然后接下来我就不太会做了所以我考场没做出来。。我们要维护的答案是一个有些连通块取 min 有些按照强制选点来取然后还要求和的玩意,但事实上我们有一个巧妙的做法(来自 Mohammad_Yasser),我们建一个权为 + ∞ +\infty + 的点,假设我们 A i A_i Ai 强制选上,只需要将 x i , 0 x_{i,0} xi,0 连向这个新点。这样我们就没有了强制选点的操作,接下来就好做了。

时间复杂度 O ( n log ⁡ n ) O(n\log n) O(nlogn) (我没按秩合并,按秩合并是 O ( n α ( n ) ) O(n\alpha(n)) O(nα(n)) 的)

代码

#include 
using namespace std;

typedef long long ll;

const int MAXN = 600005, INF = 1000000000;
int N, K, a[MAXN]; vector<int> V[MAXN]; char S[MAXN];

int fa[MAXN], val[MAXN], ans;
int get(int x) {
    if (fa[x] != x) fa[x] = get(fa[x]);
    return fa[x];
}
void merge(int x, int y) {
    x = get(x), y = get(y);
    if (x != y) fa[x] = y, val[y] += val[x];
}
inline int va(int x) {
    return min(val[get(x)], val[get(x + K)]);
}

int main() {
    scanf("%d%d%s", &N, &K, S + 1); int i, j;
    for (i = 1; i <= K; ++i) {
        int c, x; scanf("%d", &c);
        while (c--) scanf("%d", &x), V[x].push_back(i);
    }
    for (i = 1; i <= K * 2; ++i) fa[i] = i, val[i] = i > K;
    fa[K * 2 + 1] = K * 2 + 1, val[K * 2 + 1] = INF;
    for (i = 1; i <= N; ++i) {
        if (V[i].size() == 1) {
            int x = V[i][0] + (S[i] == '1') * K;
            ans -= va(V[i][0]), merge(x, K * 2 + 1), ans += va(V[i][0]);
        } else if (V[i].size() == 2) {
            int x = V[i][0], y = V[i][1];
            if (S[i] == '0') {
                if (get(x) != get(y + K))
                    ans -= va(x) + va(y), merge(x, y + K), merge(x + K, y), ans += va(x);
            }
            else {
                if (get(x) != get(y))
                    ans -= va(x) + va(y), merge(x, y), merge(x + K, y + K), ans += va(x);
            }
        }
        printf("%d\n", ans);
    }
    return 0;
}

你可能感兴趣的:([CF1290C] Prefix Enlightenment:并查集)