UVa 12858 Even Distribution

最新推荐文章于 2026-08-18 23:43:40 发布

原创最新推荐文章于 2026-08-18 23:43:40 发布 · 421 阅读

3 ·

本内容遵循CC 4.0 BY-SA版权协议

GEO检测

收录于

图遍历

题目描述

$Endre\texttt{Endre}$ 有很多侄子和侄女。每年，他都会带其中一些孩子去一个群岛旅行，那里有一家船公司在一些岛屿之间提供双向服务。由于 $Endre\texttt{Endre}$ 和孩子们可以飞往或直接从任何岛屿返回，任何行程都可以描述为一个非空的岛屿序列 $i1,i2,…,ini_1, i_2, \ldots, i_n$ ，使得每对连续的岛屿 $i_j$ 和 $i_{j+1}$ 之间有船运服务。行程的第一个和最后一个岛屿可以是同一个，也可以不同，行程中可以多次访问同一个岛屿。

群岛中的每个岛屿都生产一种不同的特色糖果，并且每次会给到达的每个团体固定数量的糖果。 $Endre\texttt{Endre}$ 自己不喜欢糖果，但孩子们几乎会立刻吃完。为了避免争吵，每次团队到达一个岛屿并收到糖果时，他会把糖果均匀地分给所有孩子。

你可能好奇 $Endre\texttt{Endre}$ 如何总能均匀分配他们在每个岛屿收到的糖果。其实答案很简单。每年，旅行社都会提前把行程计划（序列 $i1,i2,…,ini_1, i_2, \ldots, i_n$ ）发给他。由于他想带尽可能多的侄子和侄女，他计算出在遵守糖果均匀分配规则的前提下，他可以带的最大孩子数量 $k$ 。注意，每个行程计划唯一地决定了带孩子的人数。

这种情况已经持续了很多年，每次 $Endre\texttt{Endre}$ 最终带的孩子数量都不同。他想知道他能带多少种不同的孩子数量，也就是有多少个整数 $k$ ，使得存在一个行程计划，他最终会带 $k$ 个孩子。现在 $Endre\texttt{Endre}$ 正忙于准备今年的旅行，你能帮他回答这个问题吗？

输入格式

输入包含多个测试用例；每个测试用例格式如下。第一行包含两个整数 $I$ 和 $S$ （ $\le I, S \le 10^4$ ），分别表示岛屿数量和船运服务数量。岛屿用从 $1$ 到 $I$ 的不同整数标识。第二行包含 $I$ 个整数 $C1,C2,…,CIC_1, C_2, \ldots, C_I$ ，其中 $C_i$ 表示到达岛屿 $i$ 时团体收到的糖果数量（ $\le C_i \le 10^5$ ）。接下来的 $S$ 行，每行描述一个不同的船运服务，包含两个整数 $A$ 和 $B$ （ $\le A < B \le I$ ），表示可以在岛屿 $A$ 和 $B$ 之间双向旅行。没有两对岛屿之间有重复的船运服务。

输出格式

对于输入中的每个测试用例，输出一行，包含一个整数，表示使得存在某个行程计划且 $Endre\texttt{Endre}$ 最终带 $k$ 个孩子的整数 $k$ 的数量。

样例

输入

2 1
1 9
1 2
4 2
1 2 3 4
1 3
1 2
4 3
30 42 105 70
2 4
1 2
2 3

输出

2
4
11

题目分析

本题的核心是：给定一个无向图，每个顶点有一个正整数权值 $C_i$ 。一条行程对应一个顶点序列（游走）， $Endre\texttt{Endre}$ 会取该游走经过的所有顶点的 $C$ 值的最大公约数作为可带孩子的数量 $k$ 。我们需要统计所有可能的游走（非空，长度任意，可重复访问顶点）所对应的最大公约数的不同取值个数。

设某次游走经过的顶点集合为 $V^{'}$ （非空），则对应的 $\gcd\{C_v \mid v \in V'\}$ 。由于游走可以重复访问顶点，实际上只要顶点集合在原图中是连通的（即存在一条游走恰好经过这些顶点），该集合就能被实现。因此，问题转化为：对于每个连通分量，考虑所有非空连通顶点子集，计算其顶点权值的最大公约数，求所有不同值的总数。

直接枚举所有连通子集是不可能的，因为子集数量指数级。我们需要利用权值的取值范围（ $Ci≤105C_i \le 10^5$ ）以及 $GCD\texttt{GCD}$ 的性质来设计高效算法。

解题思路

核心观察

对于一个顶点集合 $V^{'}$ ，其最大公约数 $g$ 必定是 $V^{'}$ 中每个顶点权值的公约数。特别地， $g$ 一定是某个 $C_v$ 的约数。由于 $Ci≤105C_i \le 10^5$ ，每个数的约数个数最多约为 $128$ 个，因此所有可能出现的 $k$ 值数量是有限的。

我们考虑这样一个事实：如果从某个顶点 $u$ 出发，沿任意路径到达顶点 $v$ ，路径上所有顶点权值的 $GCD\texttt{GCD}$ 可以通过逐步取 $GCD\texttt{GCD}$ 得到。因此，从不同起点出发，沿不同路径到达同一个顶点时，可能得到不同的 $GCD\texttt{GCD}$ 值。这些值构成了该顶点在游走中“结尾”时的可能 $GCD\texttt{GCD}$ 值集合。

更进一步，由于图是无向的，如果存在一条从 $u$ 到 $v$ 的路径，其顶点集合的 $GCD\texttt{GCD}$ 为 $g$ ，那么反过来从 $v$ 到 $u$ 也可以得到相同的 $g$ （因为路径相同，顶点集合相同）。因此，我们可以通过“传播”的方式，让每个顶点维护一个集合，表示从任意起点出发，经过任意游走到达该顶点时，路径上顶点权值的 $GCD\texttt{GCD}$ 的所有可能取值。

算法设计

我们为每个顶点 $u$ 维护一个集合 $val[u]\textit{val}[u]$ ，存储所有可能的 $GCD\texttt{GCD}$ 值。初始时， $val[u]={Cu}\textit{val}[u] = \{C_u\}$ ，因为只访问顶点 $u$ 本身的游走给出 $GCD=Cu\texttt{GCD} = C_u$ 。

考虑一条边 $(u, v)$ 。如果从某个起点到达 $u$ 时得到了 $GCD\texttt{GCD}$ 值 $g$ ，那么从该起点继续走到 $v$ 后，新的 $GCD\texttt{GCD}$ 为 $gcd(g, C_v)$ 。因此，对于 $val[u]\textit{val}[u]$ 中的每一个值 $g$ ，我们可以将 $gcd(g, C_v)$ 插入到 $val[v]\textit{val}[v]$ 中。由于图是无向的，反过来从 $v$ 到 $u$ 也需要做同样的传播。

这类似于在图上进行“信息传播”，我们可以使用队列来加速：当一个顶点的集合插入新值时，将该顶点入队，然后处理它所有邻居，尝试用新值去更新邻居的集合。由于每个值最多被插入一次，而每个顶点的集合大小不超过 $C_u$ 的约数个数（最多约 $128$ ），因此总状态数不超过 $\times 128$ 。每条边被检查的次数为其两端顶点的集合大小之和，总操作量为 $\times 128)$ ，完全可接受。

正确性证明

充分性：如果某个值 $g$ 被插入到 $val[u]\textit{val}[u]$ 中，那么存在一个游走从某个起点到达 $u$ ，其经过顶点权值的 $GCD\texttt{GCD}$ 恰好为 $g$ 。该游走本身就是一个合法的行程，因此 $g$ 是可行的。
必要性：对于任意合法的行程，其对应的顶点序列是一条游走。设该游走最终到达的顶点为 $u$ ，则整个游走过程中顶点的 $GCD\texttt{GCD}$ 值必然在传播过程中从起点逐步传播到 $u$ ，最终被插入到 $val[u]\textit{val}[u]$ 中。因此，所有可行值都会被算法捕获。
终止性：每个顶点的集合大小有限（不超过 $C_u$ 的约数个数），且每次插入新值才入队，因此队列处理次数有限，算法必然终止。

复杂度分析

每个顶点的集合大小最多为 $C_u$ 的约数个数，约数个数最大值在 $10^5$ 内约为 $128$ 。
总插入次数 $≤I×128≤1.28×106\le I \times 128 \le 1.28 \times 10^6$ 。
每条边被检查的次数为 $size(val[u])+size(val[v])\textit{size}(\textit{val}[u]) + \textit{size}(\textit{val}[v])$ ，总和 $≤∑udeg⁡(u)×128=2S×128≤2.56×106\le \sum_{u} \deg(u) \times 128 = 2S \times 128 \le 2.56 \times 10^6$ 。
总时间复杂度 $\times D)$ ，其中 $D$ 为约数个数上界，实际运行很快。
空间复杂度 $\times D)$ ，用于存储各顶点的集合。

代码实现

// Even Distribution
// UVa ID: 12858
// Verdict: Accepted
// Submission Date: 2026-06-17
// UVa Run Time: 0.570s
//
// 版权所有（C）2026，邱秋。metaphysis # yeah dot net

#include <bits/stdc++.h>
using namespace std;

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int I, S;
    while (cin >> I >> S) {
        vector<int> C(I + 1);
        for (int i = 1; i <= I; ++i) cin >> C[i];
        vector<vector<int>> adj(I + 1);
        for (int i = 0; i < S; ++i) {
            int a, b;
            cin >> a >> b;
            adj[a].push_back(b);
            adj[b].push_back(a);
        }
        // 每个顶点的set存储所有可能的GCD值
        vector<set<int>> values(I + 1);
        vector<bool> inQueue(I + 1, false);
        queue<int> q;
        for (int i = 1; i <= I; ++i) {
            values[i].insert(C[i]);
            q.push(i);
            inQueue[i] = true;
        }
        while (!q.empty()) {
            int u = q.front(); q.pop();
            inQueue[u] = false;
              for (int v : adj[u]) {
                bool updated = false;
                for (int g : values[u]) {
                    int ng = __gcd(g, C[v]);
                    if (values[v].find(ng) == values[v].end()) {
                        values[v].insert(ng);
                        updated = true;
                    }
                }
                if (updated && !inQueue[v]) {
                    q.push(v);
                    inQueue[v] = true;
                }
            }
        }
        set<int> allValues;
        for (int i = 1; i <= I; ++i)
            for (int g : values[i])
                allValues.insert(g);
        cout << allValues.size() << '\n';
    }
    return 0;
}