牛客练习赛77 E小G的GLS图 割点+建图优化

题意

有n个点，每个点都有权值，两个点之间有一条边当且仅当gcd(a, b) != 1，每次可以去掉一个点当且仅当连通块个数增加，然后移除改点和改点相连的所有边，问最多可以移除几个点。

分析

去掉一个点连通块增加这很明显就是割点的定义，但本题如果用n^2的复杂度去暴力建图肯定是不行的，因此要考虑怎么去优化建图。

如果两个点可以相连说明他们的gcd不为1，什么时候gcd可以不为1呢，就是他们有公共质因子的时候。这样我们就可以建几个虚点代表质数，然后向虚点去连边，最后求割点。

这时有几个细节需要注意

1. 质数肯定是割点，因此不能算进去
2. 如果只有一个点连某条质数，那那个点必是割点，就出问题了，比如只有一个点6，6连2,6连3，这样6就是割点了，但6其实不是割点，因此这种情况不需要建边

Code

#include <bits/stdc++.h>
using namespace std;
#define fi first
#define se second
typedef long long ll;
typedef pair<int, int> PII;
const int N = 2e5 + 100, M = 5e5 + 5, INF = 0x3f3f3f3f;
const int MOD = 1e9 + 7;

int dfn[N], low[N], idx, cut[N], a[N];
vector<int> g[N];
void tarjan(int u, int root, int fa) {
    dfn[u] = low[u] = ++idx;
    int child = 0;
    for (auto v : g[u]) {
        if (!dfn[v]) {
            tarjan(v, root, u);
            low[u] = min(low[u], low[v]);
            if (low[v] >= dfn[u] && u != root) cut[u] = 1;
            if (u == root) child++;
        }
        else if(v != fa)
            low[u] = min(low[u], dfn[v]);
    }
    if (child >= 2 && u == root) cut[u] = 1;
}
vector<int> tmp[N];
void solve() {
    int n; cin >> n;
    for (int i = 1; i <= n; i++) cin >> a[i];
    map<int, int> mp;
    for (int i = 1; i <= n; i++) {
        int x = a[i];
        for (int j = 2; j * j <= x; j++) {
            if (x % j == 0) tmp[i].push_back(j), mp[j]++;
            while (x % j == 0) x /= j;
        }
        if (x != 1) tmp[i].push_back(x), mp[x]++;
    }
    int cnt = n;
    map<int, int> id;
    for (int i = 1; i <= n; i++) {
        for (auto it : tmp[i]) {
            if (mp[it] == 1) continue;
            if (!id[it]) id[it] = ++cnt;
            g[i].push_back(id[it]);
            g[id[it]].push_back(i);
        }
    }
    int ans = 0;
    for (int i = 1; i <= n; i++) if (!dfn[i]) tarjan(i, i, 0);
    for (int i = 1; i <= n; i++) if (cut[i]) ans += cut[i];
    cout << ans << endl;
}

signed main() {
    ios_base::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
#ifdef ACM_LOCAL
    freopen("input", "r", stdin);
    freopen("output", "w", stdout);
#endif
    solve();
}