题意
给定 n , m n, m n,m,求 ∏ X ∈ [ 1 , m ] n l c m ( X 1 , ⋯ , X n ) gcd ( X 1 , ⋯ , X n ) \prod_{X\in [1,m]^n}\mathrm{lcm}(X_1,\cdots,X_n)^{\gcd(X_1,\cdots,X_n)} X∈[1,m]n∏lcm(X1,⋯,Xn)gcd(X1,⋯,Xn)
X ∈ [ 1 , m ] n X\in [1,m]^n X∈[1,m]n表示 X X X取遍所有长度为 n n n的序列,满足 X i X_i Xi为 1 1 1到 m m m之间的整数。
n ≤ 1 0 9 , m ≤ 1 0 8 n\le 10^9,m\le 10^8 n≤109,m≤108
分析
先考虑这个问题的弱化版:求
S ( m ) = ∏ X ∈ [ 1 , m ] n l c m ( X 1 , ⋯ , X n ) S(m)=\prod_{X\in [1,m]^n}\mathrm{lcm}(X_1,\cdots,X_n) S(m)=X∈[1,m]n∏lcm(X1,⋯,Xn)
第一种算法是
S ( m ) = ∏ p ∈ P , p k ≤ m p m n − ( m − ⌊ m p k ⌋ ) n S(m)=\prod_{p\in \mathbb{P},p^k\le m}p^{m^n-(m-\lfloor\frac{m}{p^k}\rfloor)^n} S(m)=p∈P,pk≤m∏pmn−(m−⌊pkm⌋)n
可以预处理
f ( n ) = { p n = p k , p ∈ P 1 o t h e r w i s e f(n)=\begin{cases} p & n = p^k,p\in \mathbb{P}\\ 1 & otherwise \end{cases} f(n)={ p1n=pk,p∈Potherwise
的前缀积,然后通过分块 O ( m log P ) O(\sqrt m\log P) O(mlogP)求出答案。
第二种算法是通过min-max反演得到
l c m ( X 1 , ⋯ , X n ) = ∏ k = 1 n ∏ 1 ≤ y 1 < ⋯ < y k ≤ n gcd ( X y 1 , ⋯ , X y k ) ( − 1 ) k − 1 \mathrm{lcm}(X_1,\cdots,X_n)=\prod_{k=1}^n\prod_{1\le y_1<\cdots <y_k\le n}\gcd(X_{y_1},\cdots,X_{y_k})^{(-1)^{k-1}} lcm(X1,⋯,Xn)=k=1∏n1≤y1<⋯<yk≤n∏gcd(Xy1,⋯,Xyk)(−1)k−1
如果令
g ( p ) = ∑ k = 1 n ∑ 1 ≤ y 1 < ⋯ < y k ≤ n , 1