综合数论/因子分解/分块/莫比乌斯反演/树状数组HDU4947

Posted skywalker767

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了综合数论/因子分解/分块/莫比乌斯反演/树状数组HDU4947相关的知识,希望对你有一定的参考价值。

HDU-4947

题目链接:https://acm.hdu.edu.cn/showproblem.php?pid=4947

队友刷数据结构发现的好题(指秃头)。

题目大意是这样的,给你一串长为 L L L的数组,有如下操作:

  1. 1   n   d   v 1 ~n~d~ v 1 n d v,指 g c d ( x , n ) = d gcd(x , n) = d gcd(x,n)=d的下标为 x x x a x a_x ax加上 v v v
  2. 2 2 2 x x x ∑ 1 x a i ∑_1^x a_i 1xai

数据范围:
1 ≤ l q ≤ 5 e 4 1 ≤ n , d , v ≤ 2 e 5 , 1 ≤ x ≤ l 1 \\le l\\\\ q \\le 5e4 \\\\ 1 \\le n,d,v \\le 2e5 , 1 \\le x \\le l 1lq5e41n,d,v2e5,1xl
被坑辽, l l l的范围不是5e4,是2e5,第一次re辽
我们首思考一下暴力做法:对于 1 1 1操作,我们转化一下, g c d ( x , n ) = d gcd(x , n ) = d gcd(x,n)=d 可以转化为 g c d ( x / d , n / d ) = 1 gcd(x / d , n / d) = 1 gcd(x/d,n/d)=1,相当于什么呢?相当于:枚举 1 1 1 l l l中与 n / d n / d n/d互质的数乘 d d d加上 v v v,这样的时间复杂度如何呢?枚举所有与 n / d n / d n/d互质的数的时间复杂度为: O ( n / d ) O(n / d) O(n/d)最坏情况为 O ( n ) O(n) O(n),对于操作 2 2 2,我们可以用树状数组来做,那么这样我们发现,不考虑操作 2 2 2,光是修改,我们的时间都不够用。

这样不可取,那我们考虑优化。
[ g c d ( x , n ) = d ]    ⟺    [ g c d ( x / d , n / d ) = 1 ] [gcd(x , n) = d] \\iff [gcd(x / d , n / d) = 1] \\\\ [gcd(x,n)=d][gcd(x/d,n/d)=1]
对于操作 1 1 1,相当于在 1 ≤ x ≤ l 1 \\le x \\le l 1xl满足的加上 v v v,等价于:
[ g c d ( x / d , n / d ) = 1 ] ∗ v 利 用 莫 比 乌 斯 函 数 的 性 质 去 掉 方 括 号 ∑ i ∣ g c d ( x / d , n / d ) μ ( i ) ∗ v ∑ i ∣ x d , i ∣ n d μ ( i ) ∗ v n / d 已 知 , 我 们 枚 举 n / d 每 一 个 因 子 p , p ∗ d 即 为 x . ƒ ( p d ) = ∑ p d ∣ x μ ( p ) ∗ v [gcd(x / d, n / d) = 1] * v\\\\ 利用莫比乌斯函数的性质去掉方括号\\\\ ∑_i | gcd(x / d , n / d) \\mu(i) *v\\\\ ∑_i | \\fracxd , i | \\fracnd\\mu(i) *v \\\\ n / d已知,我们枚举n / d每一个因子p,p * d即为x.\\\\ ƒ(pd) = ∑_p d| x\\mu(p) *v \\\\ [gcd(x/d,n/d)=1]vigcd(x/d,n/d)μ(i)vidxidnμ(i)vn/dn/dp,pdx.ƒ(pd)=pdxμ(p)v
我们定义,
a i = ∑ d ∣ i f ( d ) a_i = ∑_d|if(d) ai=dif(d)
那么
∑ i = 1 x a i = ∑ i = 1 x ∑ j ∣ i f ( j ) ∑ i = 1 x a i = ∑ i = 1 x f ( i ) [ x i ] ∑_i = 1^xa_i = ∑_i = 1^x ∑_j | i f(j)\\\\ ∑_i = 1^xa_i = ∑_i = 1^xf(i)[\\fracxi]\\\\ i=1xai=i=1xjif(j)i=1xai=i=1xf(i)[ix]
利用分块和树状数组就可以求出。


#include <bits/stdc++.h>

#define fi first
#define se second
#define PI acos(-1)
#define mk make_pair
#define lowbit(x) (x & (-x))
#define sz(x) (int)(x).size()
#define all(a) a.begin() , a.end()
#define rep(i, a, b) for (int i = (a); i <= (b); ++i)
#define pre(i, b, a) for (int i = (b); i >= (a); --i)
#define debug freopen("1.in", "r", stdin), freopen("1.out", "w", stdout);

using namespace std;

typedef long long LL;
typedef pair<int , int > PII;
const int INF = 0x3f3f3f3f;
const int N = 2e5 + 10;

bool st[N];
LL tr[N];
int l , q;
vector<int> fact[N];
int primes[N] , mobius[N] , cnt , f[N];

int get(int n , int a) 
    return n / (n / a);


void add(int x , int  v) 
    for (int i = x;i <= l;i += lowbit(i)) tr[i] += v;


LL query(int x) 
    LL ans = 0;
    for (int i = x; i ; i -= lowbit(i)) ans += tr[i];
    return ans;


void init(int n) 
    mobius[1] = 1;
    for (int i = 以上是关于综合数论/因子分解/分块/莫比乌斯反演/树状数组HDU4947的主要内容,如果未能解决你的问题,请参考以下文章

数论入门——莫比乌斯函数,欧拉函数,狄利克雷卷积,线性筛,莫比乌斯反演,杜教筛

莫比乌斯反演+树状数组+分块求和GCD Array

bzoj 3309 DZY Loves Math —— 莫比乌斯反演+数论分块

数论分块之整除分块

从基础数论函数说起3:莫比乌斯反演

数论18——反演定理(莫比乌斯反演)