[JZOJ4392] 幂

Description

已知$x\in [1,N], y\in [1,M]$，求$x^y$有多少种取值。

Constraint

$N, M \le 10^9$

Analysis

若$A^C = B^D, C\neq 1, D\neq 1$，
则一定存在$r$使得$r^a=A, r^b=B$。

考虑到$A=\prod p_i^{a_i}, B=\prod p_i^{b_i}$，
因为$a_i \cdot C = b_i \cdot D = [C,D] \cdot k_i$，
然后$a_i = \frac{D}{(C,D)} \cdot k_i$，
我们提取出$r = \prod p_i^{k_i}$，
这样$r^{D/(C,D)}=A, r^{C/(C,D)}=B$。

那我们考虑这些不能被表示成$a^b, b\neq 1$的$r$，
考虑以$r^i$为底作出的贡献。
因为$r$之间的结果不相交，所以直接加起来就好。

假设$r^{p} \le N,r^{p+1} > N$，
实际上$p$相等的$r$的贡献不需要考虑底数，
即${(r^i)}^j, 1\le i\le p, 1\le j\le M$只用考虑$i\times j$的不同取值。

考虑$x\in [1, pM]$合法的情况是存在$t\in [1,p]$使得$t \mid x$。
把$[1,pM]$分成$p-1$段，每段是$[(i-1)M+1, iM]$，
考虑区间$[l, r]$里面的数合法，要满足$t\in [i, p], t\mid x$。
这样容斥就能计算，但$[i,p]$的规模可能达到$[1,30]$，
但$t_1, t_2$若满足$t_1 \mid t_2$，$t_2$不需要保留，
实际中规模下降到$15$个数。

时间复杂度：$O(2^{15}\log^2 N)$。