Problem - G - Codeforces
問題
整数 $x$ が与えられる. $x$ は素数列 $(p_i)$ の積で表される.
$x$ の約数の集合 $S$ であって,どの相異なる 2 要素も約数の関係にないものの考える. $|S|$ の最大値を求めよ.
制約
- $1 \leq n \leq 2 \cdot 10^5$
- $2 \leq p_i \leq 3 \cdot 10^6$
考察
まずどのような $S$ が最適かを考えるのだが,これがかなり難しい.結論を言うと, $a$ を素因数分解したときの指数部の和を $deg(a)$ としたとき, $deg(d) = \lfloor n / 2 \rfloor$ なる $x$ の約数 $d$ を全部持ってくるのが最善.
まず同じ $deg$ のもの同士は互いに約数の関係にないことは明らかだろう.その上で Hasse図 を眺めると, $deg$ が真ん中辺りになるほど点の数が増えていくことから最適っぽいことが分かる. ここ にちゃんとした証明があるらしいので気になる人はどうぞ.
後は $deg$ の分布が求まれば良い.これは,素因数 $p^q$ が増えるときに $0, 1, \cdots, q$ 全ての幅で一気に遷移するようなナップザック DP で求まるが,これは愚直にやると $O(n^2)$ かかる.
一方これを母関数で考えると,この遷移は $1 + x + \cdots + x^q$ を掛けることに等しい.よってこの掛け算に NTT を適用することで計算量が落ちる.
適当な順番で掛け算を行うと計算量は悪いままだが,次数の小さい方からマージすると $O(n \log^2 n)$ になる.
実装例
Submission #68987208 - Codeforces
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
#include <iostream>
#include <vector>
#include <queue>
#include <map>
template <class T>
using MaxHeap = std::priority_queue<T>;
template <class T>
using MinHeap = std::priority_queue<T, std::vector<T>, std::greater<T>>;
template <int MOD>
struct ModInt { ... };
constexpr int MOD = 998244353;
using mint = ModInt<MOD>;
template <int MOD, int Root>
struct NumberTheoreticalTransform { ... };
using mints = std::vector<mint>;
using NTT = NumberTheoreticalTransform<MOD, 3>;
const NTT Ntt;
// priority queue用の比較関数
struct Compare {
bool operator()(const mints& a, const mints& b) {
return a.size() > b.size();
}
};
void solve() {
int n;
std::cin >> n;
std::map<int, int> cnt;
for (int i = 0; i < n; ++i) {
int x;
std::cin >> x;
if (!cnt.count(x)) cnt[x] = 0;
++cnt[x];
}
// 多項式をpriority queueに突っ込む
std::priority_queue<mints, std::vector<mints>, Compare> que;
for (auto p : cnt) {
int q = p.second;
auto v = std::vector<mint>(q + 1, 1);
que.emplace(v);
}
// サイズの小さい方からマージ
while (que.size() > 1) {
auto f = que.top();
que.pop();
auto g = que.top();
que.pop();
mints h = Ntt.ntt(f, g);
que.push(h);
}
auto f = que.top();
std::cout << f[n / 2] << std::endl;
}
|