欢迎大家来看我写的A+B呀

##A+B

感谢大李老板友情赞助防ak,赛前打包票两个人尝试没有人能过嘻嘻嘻

A = ∑ x 1 = 1 n ⋅ ⋅ ⋅ ∑ x m = 1 n x 1 3 ∑ d ∣ g c d ( x 1 , x 2 , ⋅ ⋅ ⋅ , x m ) μ ( d ) A= \sum_{x_1=1}^{n} \cdot \cdot \cdot \sum_{x_m=1}^{n} x_1^3 \sum_{d|gcd(x_1,x_2,\cdot\cdot\cdot,x_m)}{\mu(d)} A=x1=1nxm=1nx13dgcd(x1,x2,,xm)μ(d)

= ∑ d = 1 n μ ( d ) ∑ x 1 = 1 , d ∣ x 1 n ⋅ ⋅ ⋅ ∑ x m = 1 , d ∣ x m n ∗ x 1 3 = \sum_{d=1} ^ {n} \mu (d) \sum_{x_1=1,d|x_1} ^ {n} \cdot \cdot \cdot \sum_{x_m=1,d|x_m} ^ {n}*x_1 ^ 3 =d=1nμ(d)x1=1,dx1nxm=1,dxmnx13

= ∑ d = 1 n μ ( d ) ∗ ⌊ n d ⌋ m − 1 ∗ ∑ x 1 n d x 1 3 ∗ d 3 =\sum_{d=1}^n {\mu(d)}*\lfloor \frac{n}{d} \rfloor^{m-1} *\sum_{x_1}^{\frac{n}{d}}x_1^3*d^3 =d=1nμ(d)dnm1x1dnx13d3

= ∑ d = 1 n d 3 μ ( d ) ⌊ n d ⌋ m − 1 ∗ f ( ⌊ n d ⌋ ) =\sum_{d=1}^{n}d^3\mu(d)\lfloor \frac{n}{d}\rfloor^{m-1}*f(\lfloor \frac{n}{d}\rfloor) =d=1nd3μ(d)dnm1f(dn)

f ( n ) = n 2 ( n + 1 ) 2 4 f(n)=\frac{n^2(n+1)^2}{4} f(n)=4n2(n+1)2其中 ∑ d 3 μ ( d ) \sum d^3\mu(d) d3μ(d)

用 线 性 筛 + 前 缀 和 预 处 理 o ( n ) 用线性筛 +前缀和预处理o(n) 线+o(n)当然直接暴力好像也不是不可以

⌊ n d ⌋ m − 1 f ( ⌊ n d ⌋ ) \lfloor \frac{n}{d}\rfloor^{m-1}f(\lfloor \frac{n}{d}\rfloor) dnm1f(dn)

$可以用数论分块+快速幂 $

μ ( d ) 是 莫 比 乌 斯 函 数 \mu(d)是莫比乌斯函数 μ(d)

复 杂 度 : o ( n + q ∗ n ∗ log ⁡ m ) 复杂度:o(n+q*\sqrt{n}*\log{m}) o(n+qn logm)

B = ∑ i = 1 m ∑ x 1 = 1 n ⋅ ⋅ ⋅ ∑ x m = 1 n ( x i 2 + x i ∗ ∑ i ≠ j , j = 1 n x j ) ∗ ∑ d ∣ g c d ( x 1 , x 2 , ⋅ ⋅ ⋅ , x m ) μ ( d ) B=\sum_{i=1}^{m}\sum_{x_1=1}^{n} \cdot\cdot\cdot \sum_{x_m=1}^{n}{(x_i^2+x_i*\sum_{i\neq j,j=1}^{n}x_j) }*\sum_{d|gcd(x_1,x_2,\cdot\cdot\cdot,x_m)}\mu(d) B=i=1mx1=1nxm=1n(xi2+xii=j,j=1nxj)dgcd(x1,x2,,xm)μ(d)

= ∑ i = 1 m ∑ x 1 = 1 n ⋅ ⋅ ⋅ ∑ x m = 1 n ( x i 2 ) ∗ ∑ d ∣ g c d ( x 1 , x 2 , ⋅ ⋅ ⋅ , x m ) μ ( d ) + ∑ i = 1 m ∑ x 1 = 1 n ⋅ ⋅ ⋅ ∑ x m = 1 n x i ∗ ∑ i ≠ j , j = 1 n x j ∗ ∑ d ∣ g c d ( x 1 , x 2 , ⋅ ⋅ ⋅ , x m ) μ ( d ) =\sum_{i=1}^{m}\sum_{x_1=1}^{n} \cdot\cdot\cdot\sum_{x_m=1}^{n}(x_i^2)*\sum_{d|gcd(x_1,x_2,\cdot\cdot\cdot,x_m)}\mu(d)+\sum_{i=1}^{m}\sum_{x_1=1}^{n} \cdot\cdot\cdot \sum_{x_m=1}^{n}x_i*\sum_{i\neq j,j=1}^{n}x_j*\sum_{d|gcd(x_1,x_2,\cdot\cdot\cdot,x_m)}\mu(d) =i=1mx1=1nxm=1n(xi2)dgcd(x1,x2,,xm)μ(d)+i=1mx1=1nxm=1nxii=j,j=1nxjdgcd(x1,x2,,xm)μ(d)

= P + Q =P+Q =P+Q

由于对称性,知 ∀ i ∈ [ 1 , m ] , 有 \forall i\in [1,m],有 i[1,m],

∑ x 1 = 1 n ⋅ ⋅ ⋅ ∑ x m = 1 n ( x i 2 ) ∑ d ∣ g c d ( x 1 , x 2 , ⋅ ⋅ ⋅ , x m ) μ ( d ) \sum_{x_1=1} ^ {n} \cdot \cdot \cdot \sum_{x_m=1} ^ {n}(x_i ^ 2) \sum_{ d |gcd ( x _ 1,x _ 2,\cdot\cdot\cdot,x_m)}\mu(d) x1=1nxm=1n(xi2)dgcd(x1,x2,,xm)μ(d)

= ∑ x 1 = 1 n ⋅ ⋅ ⋅ ∑ x m = 1 n ( x 1 2 ) ∑ d ∣ g c d ( x 1 , x 2 , ⋅ ⋅ ⋅ , x m ) μ ( d ) =\sum_{x_1=1}^{n}\cdot \cdot \cdot\sum_{x_m=1}^{n}(x_1^2)\sum_{d|gcd(x_1,x_2,\cdot\cdot\cdot,x_m)}\mu(d) =x1=1nxm=1n(x12)dgcd(x1,x2,,xm)μ(d)

P = m ∗ ∑ x 1 = 1 n ⋅ ⋅ ⋅ ∑ x m = 1 n ( x 1 2 ) ∑ d ∣ g c d ( x 1 , x 2 , ⋅ ⋅ ⋅ , x m ) μ ( d ) P=m*\sum_{x_1=1}^{n}\cdot \cdot \cdot\sum_{x_m=1}^{n}(x_1^2)\sum_{d|gcd(x_1,x_2,\cdot\cdot\cdot,x_m)}\mu(d) P=mx1=1nxm=1n(x12)dgcd(x1,x2,,xm)μ(d)

= m ∗ ∑ d = 1 n μ ( d ) ∗ ∑ x 1 = 1 , d ∣ x 1 n ⋅ ⋅ ⋅ ∑ x m = 1 , d ∣ x m n ( x 1 2 ) =m*\sum_{d=1}^{n}\mu(d)*\sum_{x_1=1,d|x_1}^{n}\cdot \cdot \cdot\sum_{x_m=1,d|x_m}^{n}(x_1^2) =md=1nμ(d)x1=1,dx1nxm=1,dxmn(x12)

= m ∗ ∑ d = 1 n μ ( d ) ⌊ n d ⌋ m − 1 ∑ x 1 = 1 n d x 1 2 ∗ d 2 =m*\sum_{d=1}^{n}\mu(d)\lfloor \frac{n}{d}\rfloor^{m-1}\sum_{x_1=1}^{\frac{n}{d}}x_1^2*d^2 =md=1nμ(d)dnm1x1=1dnx12d2

= m ∗ ∑ d = 1 n μ ( d ) d 2 ⌊ n d ⌋ m − 1 ∗ g ( ⌊ n d ⌋ ) =m*\sum_{d=1}^{n}\mu(d)d^2\lfloor \frac{n}{d}\rfloor^{m-1}*g(\lfloor \frac{n}{d}\rfloor) =md=1nμ(d)d2dnm1g(dn)

其 中 , g ( n ) = n ( n + 1 ) ( 2 n + 1 ) 6 其中,g(n)=\frac{n(n+1)(2n+1)}{6} g(n)=6n(n+1)(2n+1)

预 处 理 ∑ d 2 μ ( d ) 预处理 \sum {d^2\mu(d)} d2μ(d)

数 论 分 块 + 快 速 幂 数论分块+快速幂 +

o ( n + q n l o g m ) o(n+q\sqrt n log{m}) o(n+qn logm)

$ 同理Q$

Q = m ( m − 1 ) ∑ x 1 = 1 n ⋅ ⋅ ⋅ ∑ x m = 1 n x 1 x 2 ∑ d ∣ g c d ( x 1 , x 2 , ⋅ ⋅ ⋅ , x m ) μ ( d ) Q=m(m-1)\sum_{x_1=1}^{n} \cdot\cdot\cdot \sum_{x_m=1}^{n}x_1x_2\sum_{d|gcd(x_1,x_2,\cdot\cdot\cdot,x_m)}\mu(d) Q=m(m1)x1=1nxm=1nx1x2dgcd(x1,x2,,xm)μ(d)

不 理 解 的 可 以 手 推 ( x + y + z ) 2 不理解的可以手推(x+y+z)^2 (x+y+z)2

∑ x = 1 n ∑ y = 1 n ∑ z = 1 n ( x + y + z ) 2 \sum_{x=1}^{n}\sum_{y=1}^{n}\sum_{z=1}^{n}(x+y+z)^2 x=1ny=1nz=1n(x+y+z)2

= 3 ∑ x = 1 n ∑ y = 1 n ∑ z = 1 n x 2 + 6 ∗ ∑ x = 1 n ∑ y = 1 n ∑ z = 1 n x y =3\sum_{x=1}^{n}\sum_{y=1}^{n}\sum_{z=1}^{n}x^2+6*\sum_{x=1}^{n}\sum_{y=1}^{n}\sum_{z=1}^{n}xy =3x=1ny=1nz=1nx2+6x=1ny=1nz=1nxy

或者借助完全图来理解变量之间的对称

Q = m ( m − 1 ) ∑ x 1 = 1 n ⋅ ⋅ ⋅ ∑ x m = 1 n x 1 x 2 ∑ d ∣ g c d ( x 1 , x 2 , ⋅ ⋅ ⋅ , x m ) μ ( d ) Q=m(m-1)\sum_{x_1=1}^{n} \cdot\cdot\cdot \sum_{x_m=1}^{n}x_1x_2\sum_{d|gcd(x_1,x_2,\cdot\cdot\cdot,x_m)}\mu(d) Q=m(m1)x1=1nxm=1nx1x2dgcd(x1,x2,,xm)μ(d)

= m ( m − 1 ) ∑ d = 1 n μ ( d ) ⌊ n d ⌋ m − 2 ∗ ∑ x 1 = 1 n d ∑ x 2 = 1 n d d 2 x 1 x 2 =m(m-1)\sum_{d=1}^n\mu(d)\lfloor \frac{n}{d}\rfloor^{m-2}*\sum_{x_1=1}^{\frac{n}{d}}\sum_{x_2=1}^{\frac{n}{d}}d^2x_1x_2 =m(m1)d=1nμ(d)dnm2x1=1dnx2=1dnd2x1x2

= m ( m − 1 ) ∑ d = 1 n μ ( d ) d 2 ⌊ n d ⌋ m − 2 ∗ ( ∑ x 1 = 1 n d x 1 ) 2 =m(m-1)\sum_{d=1}^n\mu(d)d^2\lfloor \frac{n}{d}\rfloor^{m-2}*(\sum_{x_1=1}^{\frac{n}{d}}x_1)^2 =m(m1)d=1nμ(d)d2dnm2(x1=1dnx1)2

= m ( m − 1 ) ∑ d = 1 n μ ( d ) d 2 ⌊ n d ⌋ m − 2 ∗ ( ⌊ n d ⌋ ( 1 + ⌊ n d ⌋ ) 2 ) 2 =m(m-1)\sum_{d=1}^n\mu(d)d^2\lfloor \frac{n}{d}\rfloor^{m-2}*(\frac{\lfloor\frac{n}{d}\rfloor(1+\lfloor\frac{n}{d}\rfloor) }{2})^2 =m(m1)d=1nμ(d)d2dnm2(2dn(1+dn))2

= m ( m − 1 ) ∑ d = 1 n μ ( d ) d 2 ⌊ n d ⌋ m ( ( 1 + ⌊ n d ⌋ ) 2 ) 2 =m(m-1)\sum_{d=1}^n\mu(d)d^2\lfloor \frac{n}{d}\rfloor^{m}(\frac{(1+\lfloor\frac{n}{d}\rfloor) }{2})^2 =m(m1)d=1nμ(d)d2dnm(2(1+dn))2

数 论 分 块 + 快 速 幂 数论分块+快速幂 +

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值