题目:
http://acm.hdu.edu.cn/showproblem.php?pid=6833
Given you
n
n
n,
x
x
x,
k
k
k , find the value of the following formula:
∑
a
1
=
1
n
.
.
.
∑
a
x
=
1
n
(
∏
j
=
1
x
a
j
k
)
f
(
g
c
d
(
a
1
,
a
2
.
.
.
a
x
)
)
⋅
g
c
d
(
a
1
,
a
2
.
.
.
a
x
)
\sum_{a_1=1}^{n}...\sum_{a_x=1}^{n}\left(\prod_{j=1}^{x}a_j^k\right)f(gcd(a_1,a_2...a_x))\cdot gcd(a_1,a_2...a_x)
a1=1∑n...ax=1∑n(j=1∏xajk)f(gcd(a1,a2...ax))⋅gcd(a1,a2...ax)
The function
f
(
x
)
f(x)
f(x) is defined as follows:
If there exists an ingeter
k
(
k
>
1
)
k (k>1)
k(k>1) , and
k
2
k^2
k2 is a divisor of
x
x
x,
then
f
(
x
)
=
0
f(x)=0
f(x)=0, else
f
(
x
)
=
1
f(x)=1
f(x)=1.
输入:
The first line contains three integers
t
,
k
,
x
(
1
≤
t
≤
1
0
4
,
1
≤
k
≤
1
0
9
,
1
≤
x
≤
1
0
9
)
t,k,x (1≤t≤10^4,1≤k≤10^9,1≤x≤10^9)
t,k,x(1≤t≤104,1≤k≤109,1≤x≤109)
Then t t t test cases follow. Each test case contains an integer n ( 1 ≤ n ≤ 2 × 1 0 5 ) n (1≤n≤2×10^5) n(1≤n≤2×105)
思路:
∑
a
1
=
1
n
.
.
.
∑
a
x
=
1
n
(
∏
j
=
1
x
a
j
k
)
f
(
g
c
d
(
a
1
,
a
2
.
.
.
a
x
)
)
⋅
g
c
d
(
a
1
,
a
2
.
.
.
a
x
)
=
∑
e
=
1
n
∑
a
1
=
1
n
.
.
.
∑
a
x
=
1
n
[
g
c
d
(
a
1
,
a
2
.
.
.
a
x
)
=
e
]
(
∏
j
=
1
x
a
j
k
)
f
(
e
)
⋅
e
=
∑
e
=
1
n
∑
a
1
=
1
⌊
n
e
⌋
.
.
.
∑
a
x
=
1
⌊
n
e
⌋
[
g
c
d
(
a
1
,
a
2
.
.
.
a
x
)
=
1
]
(
∏
j
=
1
x
(
a
j
⋅
e
)
k
)
f
(
e
)
⋅
e
=
∑
e
=
1
n
∑
a
1
=
1
⌊
n
e
⌋
.
.
.
∑
a
x
=
1
⌊
n
e
⌋
∑
d
∣
(
a
1
,
a
2
.
.
.
a
x
)
μ
(
d
)
(
∏
j
=
1
x
(
a
j
⋅
e
)
k
)
f
(
e
)
⋅
e
=
∑
e
=
1
n
∑
d
=
1
⌊
n
e
⌋
∑
a
1
=
1
⌊
⌊
n
e
⌋
d
⌋
.
.
.
∑
a
x
=
1
⌊
⌊
n
e
⌋
d
⌋
μ
(
d
)
(
∏
j
=
1
x
(
a
j
⋅
e
⋅
d
)
k
)
f
(
e
)
⋅
e
=
∑
e
=
1
n
∑
d
=
1
⌊
n
e
⌋
μ
(
d
)
⋅
f
(
e
)
⋅
e
k
x
+
1
⋅
d
k
x
∑
a
1
=
1
⌊
⌊
n
e
⌋
d
⌋
.
.
.
(
∏
j
=
1
x
−
1
a
j
k
)
∑
a
x
=
1
⌊
⌊
n
e
⌋
d
⌋
a
x
k
令
S
U
M
(
x
)
=
∑
i
=
1
x
i
k
.
.
.
=
∑
e
=
1
n
∑
d
=
1
⌊
n
e
⌋
μ
(
d
)
⋅
f
(
e
)
⋅
e
k
x
+
1
⋅
d
k
x
⋅
S
U
M
(
⌊
⌊
n
e
⌋
d
⌋
)
x
=
∑
e
=
1
n
f
(
e
)
⋅
e
k
x
+
1
∑
d
=
1
⌊
n
e
⌋
μ
(
d
)
⋅
d
k
x
⋅
S
U
M
(
⌊
⌊
n
e
⌋
d
⌋
)
x
令
T
=
e
d
=
∑
T
=
1
n
S
U
M
(
⌊
n
T
⌋
)
x
⋅
T
k
x
∑
e
∣
T
f
(
e
)
⋅
e
⋅
μ
(
T
e
)
令
S
u
m
(
T
)
=
T
k
x
∑
e
∣
T
f
(
e
)
⋅
e
⋅
μ
(
T
e
)
=
∑
T
=
1
n
S
U
M
(
⌊
n
T
⌋
)
x
⋅
S
u
m
(
T
)
\begin{aligned} &\sum_{a_1=1}^{n}...\sum_{a_x=1}^{n}\left(\prod_{j=1}^{x}a_j^k\right)f(gcd(a_1,a_2...a_x))\cdot gcd(a_1,a_2...a_x)\\ =&\sum_{e=1}^{n}\sum_{a_1=1}^{n}...\sum_{a_x=1}^{n}[gcd(a_1,a_2...a_x)=e]\left(\prod_{j=1}^{x}a_j^k\right)f(e)\cdot e\\ =&\sum_{e=1}^{n}\sum_{a_1=1}^{\lfloor\frac{n}{e}\rfloor}...\sum_{a_x=1}^{\lfloor\frac{n}{e}\rfloor}[gcd(a_1,a_2...a_x)=1]\left(\prod_{j=1}^{x}(a_j\cdot e)^k\right)f(e)\cdot e\\ =&\sum_{e=1}^{n}\sum_{a_1=1}^{\lfloor\frac{n}{e}\rfloor}...\sum_{a_x=1}^{\lfloor\frac{n}{e}\rfloor}\sum_{d|(a_1,a_2...a_x)}\mu(d)\left(\prod_{j=1}^{x}(a_j\cdot e)^k\right)f(e)\cdot e\\ =&\sum_{e=1}^{n}\sum_{d=1}^{\lfloor\frac{n}{e}\rfloor}\sum_{a_1=1}^{\lfloor\frac{\lfloor\frac{n}{e}\rfloor}{d}\rfloor}...\sum_{a_x=1}^{\lfloor\frac{\lfloor\frac{n}{e}\rfloor}{d}\rfloor}\mu(d)\left(\prod_{j=1}^{x}(a_j\cdot e\cdot d)^k\right)f(e)\cdot e\\ =&\sum_{e=1}^{n}\sum_{d=1}^{\lfloor\frac{n}{e}\rfloor}\mu(d)\cdot f(e)\cdot e^{kx+1}\cdot d^{kx}\sum_{a_1=1}^{\lfloor\frac{\lfloor\frac{n}{e}\rfloor}{d}\rfloor}...\left(\prod_{j=1}^{x-1}a_j^k\right)\sum_{a_x=1}^{\lfloor\frac{\lfloor\frac{n}{e}\rfloor}{d}\rfloor}a_x^k\quad令SUM(x)=\sum_{i=1}^{x}i^k\\ ...\\ =&\sum_{e=1}^{n}\sum_{d=1}^{\lfloor\frac{n}{e}\rfloor}\mu(d)\cdot f(e)\cdot e^{kx+1}\cdot d^{kx}\cdot SUM(\lfloor\frac{\lfloor\frac{n}{e}\rfloor}{d}\rfloor)^x\\ =&\sum_{e=1}^{n}f(e)\cdot e^{kx+1}\sum_{d=1}^{\lfloor\frac{n}{e}\rfloor}\mu(d)\cdot d^{kx}\cdot SUM(\lfloor\frac{\lfloor\frac{n}{e}\rfloor}{d}\rfloor)^x\quad 令T=ed\\ =&\sum_{T=1}^{n}SUM(\lfloor\frac{n}{T}\rfloor)^x\cdot T^{kx}\sum_{e|T}f(e)\cdot e\cdot \mu(\frac{T}{e})\quad 令Sum(T)=T^{kx}\sum_{e|T}f(e)\cdot e\cdot \mu(\frac{T}{e})\\ =&\sum_{T=1}^{n}SUM(\lfloor\frac{n}{T}\rfloor)^x\cdot Sum(T)\\ \end{aligned}
=====...====a1=1∑n...ax=1∑n(j=1∏xajk)f(gcd(a1,a2...ax))⋅gcd(a1,a2...ax)e=1∑na1=1∑n...ax=1∑n[gcd(a1,a2...ax)=e](j=1∏xajk)f(e)⋅ee=1∑na1=1∑⌊en⌋...ax=1∑⌊en⌋[gcd(a1,a2...ax)=1](j=1∏x(aj⋅e)k)f(e)⋅ee=1∑na1=1∑⌊en⌋...ax=1∑⌊en⌋d∣(a1,a2...ax)∑μ(d)(j=1∏x(aj⋅e)k)f(e)⋅ee=1∑nd=1∑⌊en⌋a1=1∑⌊d⌊en⌋⌋...ax=1∑⌊d⌊en⌋⌋μ(d)(j=1∏x(aj⋅e⋅d)k)f(e)⋅ee=1∑nd=1∑⌊en⌋μ(d)⋅f(e)⋅ekx+1⋅dkxa1=1∑⌊d⌊en⌋⌋...(j=1∏x−1ajk)ax=1∑⌊d⌊en⌋⌋axk令SUM(x)=i=1∑xike=1∑nd=1∑⌊en⌋μ(d)⋅f(e)⋅ekx+1⋅dkx⋅SUM(⌊d⌊en⌋⌋)xe=1∑nf(e)⋅ekx+1d=1∑⌊en⌋μ(d)⋅dkx⋅SUM(⌊d⌊en⌋⌋)x令T=edT=1∑nSUM(⌊Tn⌋)x⋅Tkxe∣T∑f(e)⋅e⋅μ(eT)令Sum(T)=Tkxe∣T∑f(e)⋅e⋅μ(eT)T=1∑nSUM(⌊Tn⌋)x⋅Sum(T)
注意:
多组数据只有
n
n
n。
预处理
S
u
m
(
i
)
,
1
≤
i
≤
n
Sum(i),1\le i\le n
Sum(i),1≤i≤n,然后分块