Number Theory

因为github上不支持对数学公式的渲染,所以就放到了csdn上,感觉还不错.(持续更新)

Number Theory

(目录坏掉了TAT)

• FFT
• 线性筛
• 莫比乌斯反演:${\sum }_{i=1}^n{\sum }_{j=1}^m[\gcd (i,j)==k]$
• 莫比乌斯反演:${\sum }_{b_1=1}^{a_1}{\sum }_{b_2=1}^{a_2}...{\sum }_{b_n=1}^{a_n}[\gcd (b_1,b_2,...,b_n)\ge 2]$
• 莫比乌斯反演:${\sum }_{k=1}^{\infty }{\sum }_{i=1}^n{\sum }_{j=1}^m[\gcd (i,j)==k^2]$
• 莫比乌斯反演:${\sum }_{i=1}^n{\sum }_{j=1}^m\gcd (i,j)$
• 莫比乌斯反演:${\sum }_{i=1}^n{\sum }_{j=1}^{m}lcm(i,j)$
• 莫比乌斯反演:${\sum }_{i=1}^n{\sum }_{j=1}^{m}d(ij)$
• 莫比乌斯反演:${\sum }_{i=1}^n[i不含有平方因子]$

FFT

推荐博客

FFT有很多形式,但最终都可以归结到这一种:

$$h=f\otimes g:h(x)=\sum _{i=0}^{x}f(i)\times g(x-i)$$

hdu1402 f[i]=s’[i]-‘0’,g[i]=t’[i]-‘0’,h=f⊗g,逆序输出h[i],h[i]含义:第i位乘积和 AC代码

bzoj3527 $f[flen-1-i]=q_i,g[i]=\frac{1}{(n-1-i)\times |n-1-i|},g[n-1]=0.$ AC代码

$$h=f\otimes g,h[k+flen-1]=\sum _{i=0}^{flen-1}f[flen-1-i]g[i+k]$$
uvalive6886 f[0]=1,f[a[i]]=1,h=f⊗f,f[i]的含义是走0或1步到达i的方案数,h[i]的含义是走0~2步到达i的方案数 AC代码

hdu4609 f[a[i]]=c[i]代表有c[i]个方案选出长度为a[i]树枝,h=f⊗f,h[i] :选两个树枝长度为i有多少种方案 AC代码

直接的h[i]=(b,c)的方案数,a是待会要枚举的a[i],b和c是被加的树枝长,需要减去一些东西

1.减去b=c情况,即枚举每根树枝h[a[i]*2]–

2.减去b>c情况,即每一个h[i]/=2

剩下的就是b<c的情况,枚举a[i],假设a>c,总数为
${\sum }_{i=a[i]+1}^{\infty }h[i]$
,需要减去如下情况

1.a=b<c 或 b<a=c,一共有n-1种情况

2.a<b<c,一共有
$C_{n-1-i}^2=\frac{(n-1-i)(n-2-i)}{2}$
种情况

3.b<a<c,一共有i*(n-1-i)种情况

把所有情况数相加就是答案

uvalive4671 求A串中alen-blen+1个长如B串的子串和B串的Hamming距离小于等于k的种数 AC代码

首先需要求出a串中长如B串的子串们与B串的Hamming距离

$$定义f(x)是a串中以x为起始下标长如B串的子串与B串的Hamming距离$$
$$即f(x)=\sum _{i=0}^{blen-1}[A(x+i)!=B(i)]$$
$$设i=blen-1-i,C:reverse(B,B+blen),\frac{それでは}{那么}$$
$$f(x)=\sum _{blen-i=0}^{blen-1}[A(x+blen-1-i)!=B(blen-1-i)]$$
$$f(x)=\sum _{i=0}^{blen-1}[A(x+blen-1-i)!=C(i)],根据FFT形式有:$$
$$f(x)=(A\otimes C)(x+blen-1)$$
$$现在的问题是如何计算A\otimes C,可以拆分为两个卷积的和:$$
$$A\otimes C=(a(i)={=}^{\prime }{a}^{\prime })\otimes (b(i)={=}^{\prime }{b}^{\prime })+(a(i)={=}^{\prime }{b}^{\prime })\otimes (b(i)={=}^{\prime }{a}^{\prime })$$

然后在用字符串Hash求不重复合法子串个数(不是重点,见代码)

jisuanke26219 剪刀石头布加强版,小串最多赢大串几次 AC代码

$$定义winnum(x,y,z)=(a(i)==x)\otimes (b(i)==y)(z+alen-1),a是小串$$
$$x="SLPKRSLPKR",y="PPRRLLKKSS",根据游戏规则所得$$
则 A n s = max ⁡ i = 0 b l e n − a l e n + 1 { ∑ j = 0 9 w i n n u m ( x [ j ] , y [ j ] , i ) } 则Ans=\max_{i=0}^{blen-alen+1}\{\sum_{j=0}^9winnum(x[j],y[j],i)\} 则Ans=i=0maxblen−alen+1​{j=0∑9​winnum(x[j],y[j],i)}

克制关系见下图:(S:紫剪刀,K:蓝斯波克,L:绿蜥蜴,R:粉石头,P:黄布)

cf954i 最小次数本质是求6点图的联通分量个数,用30次FFT和并查集来确定答案,每一个答案是a的子串和b串对应位置的两个点相连之后,6-联通块数量. AC代码 推荐博客的第6个题,说的很明白

bzoj4827 求的是 min ⁡ C = − m m { min ⁡ k = 0 l e n − 1 { ∑ i = 0 l e n − 1 x i 2 + y i + k 2 − 2 x i y i + k + C 2 + 2 C x i − 2 C y i + k } } \min_{C=-m}^m\{\min_{k=0}^{len-1}\{\sum_{i=0}^{len-1}x_i^2+y_{i+k}^2-2x_iy_{i+k}+C^2+2Cx_i-2Cy_{i+k}\}\} minC=−mm​{mink=0len−1​{∑i=0len−1​xi2​+yi+k2​−2xi​yi+k​+C2+2Cxi​−2Cyi+k​}} AC代码

由于$y_{0+k}到y_{len-1+k}$的数字和$y_0到y_{len-1}$的数字只是顺序上的差别,因此有:

A n s = min ⁡ C = − m m { min ⁡ k = 0 l e n − 1 { ∑ i = 0 l e n − 1 x i 2 + y i 2 − 2 x i y i + k + C 2 + 2 C x i − 2 C y i } } Ans=\min_{C=-m}^m\{\min_{k=0}^{len-1}\{\sum_{i=0}^{len-1}x_i^2+y_i^2-2x_iy_{i+k}+C^2+2Cx_i-2Cy_i\}\} Ans=C=−mminm​{k=0minlen−1​{i=0∑len−1​xi2​+yi2​−2xi​yi+k​+C2+2Cxi​−2Cyi​}}

A n s = ∑ i = 0 l e n − 1 ( x i 2 + y i 2 ) − 2 max ⁡ k = 0 l e n − 1 { ∑ i = 0 l e n − 1 x i y i + k } + min ⁡ C = − m m { l e n × C 2 + 2 C ∑ i = 0 l e n − 1 ( x i − y i ) } Ans=\sum_{i=0}^{len-1}(x_i^2+y_i^2)-2\max_{k=0}^{len-1}\{\sum_{i=0}^{len-1}x_iy_{i+k}\}+\min_{C=-m}^m\{len\times C^2+2C\sum_{i=0}^{len-1}(x_i-y_i)\} Ans=i=0∑len−1​(xi2​+yi2​)−2k=0maxlen−1​{i=0∑len−1​xi​yi+k​}+C=−mminm​{len×C2+2Ci=0∑len−1​(xi​−yi​)}

其中 ${\sum }_{i=0}^{len-1}(x_i^2+y_i^2)$ 和 ${\sum }_{i=0}^{len-1}(x_i-y_i)$ O(n)能算,

关键是 max ⁡ k = 0 l e n − 1 { ∑ i = 0 l e n − 1 x i y i + k } \max_{k=0}^{len-1}\{\sum_{i=0}^{len-1}x_iy_{i+k}\} maxk=0len−1​{∑i=0len−1​xi​yi+k​}如何使用FFT计算

设x:rev(x,x+len),则上式可写成 max ⁡ k = 0 l e n − 1 { ∑ i = 0 l e n − 1 x l e n − 1 − i y i + k } \max_{k=0}^{len-1}\{\sum_{i=0}^{len-1}x_{len-1-i}y_{i+k}\} maxk=0len−1​{∑i=0len−1​xlen−1−i​yi+k​}

即  max ⁡ k = 0 l e n − 1 { x ⊗ y ( l e n − 1 + k ) } \max_{k=0}^{len-1}\{x⊗y(len-1+k)\} maxk=0len−1​{x⊗y(len−1+k)} ,O(nlogn)能算

NAIPC2016E f(i)=s(len-1-i),g(i)=s(i),h=f⊗g,依次输出h[i],$i\in [len,2len-2]$,h[i+len-1]就是距离为i的逆序对个数 AC代码

线性筛

推荐博客

nowcoder27E $f=i{d}^k\times \mu$ AC代码

$f(1)=0,f(p)=p^k-1,f(p^k)=f(p^{k-1})\times p^k$

莫比乌斯反演:${\sum }_{i=1}^n{\sum }_{j=1}^m[\gcd (i,j)==k]$

luogu3455  ∑ i = 1 n ∑ j = 1 m [ gcd ⁡ ( i , j ) = = k ]   ( T ∈ [ 1 , 50 &ThinSpace; 000 ] { n , m ∈ [ 1 , 50 &ThinSpace; 000 ] } ) \sum_{i=1}^n\sum_{j=1}^m[\gcd(i,j)==k] (T\in[1,50\,000]\{n,m\in[1,50\,000]\}) ∑i=1n​∑j=1m​[gcd(i,j)==k] (T∈[1,50000]{n,m∈[1,50000]}) AC代码

hdu1695 ${\sum }_{i=1}^n{\sum }_{j=i+1}^m[\gcd (i,j)==k](a,b\in [1,100000],k\in [0,100000])$ AC代码

bzoj2301  ∑ i = a b ∑ j = c d [ gcd ⁡ ( i , j ) = = k ]   ( n ∈ [ 1 , 50 &ThinSpace; 000 ] { 1 ≤ a ≤ b ≤ 50 &ThinSpace; 000 , 1 ≤ c ≤ d ≤ 50 &ThinSpace; 000 } ) \sum_{i=a}^b\sum_{j=c}^d[\gcd(i,j)==k] (n\in[1,50\,000]\{1\leq a\leq b\leq 50\,000,1\leq c\leq d\leq 50\,000\}) ∑i=ab​∑j=cd​[gcd(i,j)==k] (n∈[1,50000]{1≤a≤b≤50000,1≤c≤d≤50000}) AC代码

$$设f(k,n,m)=\sum _{i=1}^n\sum _{j=1}^m[\gcd (i,j)==k]$$

$$\frac{倍数反演したい}{试图倍数反演}$$
$$设g(k,n,m)=\sum _{i=1}^{\infty }f(ik,n,m)$$
$$g(k,n,m)=\sum _{i=1}^n\sum _{j=1}^m[\gcd (i,j)\%k==0]$$
$$做代换i=ak,j=bk$$
$$\frac{それでは}{那么}g(k,n,m)=\sum _{ak=1}^n\sum _{bk=1}^m[\gcd (ak,bk)\%k==0]$$
$$\frac{そして}{然后}g(k,n,m)=\sum _{a=1}^{n/k}\sum _{b=1}^{m/k}1$$
$$\frac{だから}{所以}g(k,n,m)=\lfloor \frac{n}{k}\rfloor \times \lfloor \frac{m}{k}\rfloor$$
$$倍数反演成功f(k,n,m)=\sum _{i=1}^{\infty }g(ik,n,m)\times \mu (i)$$
$$f(k,n,m)=\sum _{i=1}^{\infty }\lfloor \frac{n}{ik}\rfloor \times \lfloor \frac{m}{ik}\rfloor \times \mu (i)$$
$$\frac{だって}{因为}\lfloor \frac{n}{ik}\rfloor \ge 1才会对总和有贡献,\frac{だから}{所以}i\le \lfloor \frac{n}{k}\rfloor$$
$$f(k,n,m)=\sum _{i=1}^{n/k}\lfloor \frac{n}{ik}\rfloor \times \lfloor \frac{m}{ik}\rfloor \times \mu (i)\Leftarrow 可整除分块O(\sqrt{n})$$
$$上式O(\sqrt{n})可算,需预处理出\mu (i)的前缀和O(n)$$

莫比乌斯反演:${\sum }_{b_1=1}^{a_1}{\sum }_{b_2=1}^{a_2}...{\sum }_{b_n=1}^{a_n}[\gcd (b_1,b_2,...,b_n)\ge 2]$

hdu6053  ∑ b 1 = 1 a 1 ∑ b 2 = 1 a 2 . . . ∑ b n = 1 a n [ gcd ⁡ ( b 1 , b 2 , . . . , b n ) ≥ 2 ]   ( T ∈ [ 1 , 10 ] { n , a i ∈ [ 1 , 100 &ThinSpace; 000 ] } ) \sum_{b_1=1}^{a_1}\sum_{b_2=1}^{a_2}...\sum_{b_n=1}^{a_n}[\gcd(b_1,b_2,...,b_n)\geq 2] (T\in [1,10]\{n,a_i\in[1,100\,000]\}) ∑b1​=1a1​​∑b2​=1a2​​...∑bn​=1an​​[gcd(b1​,b2​,...,bn​)≥2] (T∈[1,10]{n,ai​∈[1,100000]}) AC代码

$$设f(k)=\sum _{b_1=1}^{a_1}\sum _{b_2=1}^{a_2}...\sum _{b_n=1}^{a_n}[\gcd (b_1,b_2,...,b_n)\%k==0]$$
$$根据容斥原理,ans=\sum _{i=2}^{\infty }-\mu (i)\times f(i)$$
$$显然f(k)=\prod _{i=1}^n\lfloor \frac{a_i}{k}\rfloor ,\frac{けれど}{但是}并不好算$$

f ( k ) = ∏ i = 1 m a x { a } / k i g ( i , d ) f(k)=\prod_{i=1}^{max\{a\}/k}i^{g(i,d)} f(k)=i=1∏max{a}/k​ig(i,d)
$$g(x,y)=\sum _{i=1}^n[\lfloor \frac{a_i}{y}\rfloor ==x]$$
$$g(x,y)=\sum _{i=1}^n[xy\le a_i<xy+y]$$
$$\frac{だから}{所以}需要预处理出h(x)=\sum _{i=1}^n[a_i<x]$$
$$\frac{そして}{然后}g(x,y)=h(xy+y)-h(xy)$$
だ か ら 所 以 f ( k ) = ∏ i = 1 m a x { a } / k i h ( i d + d ) − h ( i d ) \frac{だから}{所以}f(k)=\prod_{i=1}^{max\{a\}/k}i^{h(id+d)-h(id)} 所以だから​f(k)=i=1∏max{a}/k​ih(id+d)−h(id)
a n s = ∑ j = 2 ∞ − μ ( j ) × ∏ i = 1 m a x { a } / j i h ( i d + d ) − h ( i d ) ans=\sum_{j=2}^\infty-\mu(j)\times\prod_{i=1}^{max\{a\}/j}i^{h(id+d)-h(id)} ans=j=2∑∞​−μ(j)×i=1∏max{a}/j​ih(id+d)−h(id)
だ つ て 因 为 m a x { a } / j ≥ 1 , だ か ら 所 以 j ≤ m a x { a } \frac{だつて}{因为}max\{a\}/j\geq 1,\frac{だから}{所以}j\leq max\{a\} 因为だつて​max{a}/j≥1,所以だから​j≤max{a}
a n s = ∑ j = 2 m a x { a } − μ ( j ) × ∏ i = 1 m a x { a } / j i h ( i d + d ) − h ( i d ) ans=\sum_{j=2}^{max\{a\}}-\mu(j)\times\prod_{i=1}^{max\{a\}/j}i^{h(id+d)-h(id)} ans=j=2∑max{a}​−μ(j)×i=1∏max{a}/j​ih(id+d)−h(id)
$$上式O(nlo{g}^{2}n)可算,需预处理出\mu (i)和h(i)O(n)$$

莫比乌斯反演:${\sum }_{k=1}^{\infty }{\sum }_{i=1}^n{\sum }_{j=1}^m[\gcd (i,j)==k^2]$

hdu5663 nm- ∑ k = 1 ∞ ∑ i = 1 n ∑ j = 1 m [ gcd ⁡ ( i , j ) = = k 2 ]   ( T ∈ [ 1 , 10 &ThinSpace; 000 ] { n , m ∈ [ 1 , 10 &ThinSpace; 000 &ThinSpace; 000 ] } ) \sum_{k=1}^\infty\sum_{i=1}^n\sum_{j=1}^m[\gcd(i,j)==k^2] (T\in[1,10\,000]\{n,m\in[1,10\,000\,000]\}) ∑k=1∞​∑i=1n​∑j=1m​[gcd(i,j)==k2] (T∈[1,10000]{n,m∈[1,10000000]}) AC代码

$$ans=\sum _{k=1}^{\infty }\sum _{i=1}^n\sum _{j=1}^m[\gcd (i,j)==k^2]$$
$$ans=\sum _{k=1}^{\infty }\sum _{d=1}^{\infty }[d==k^2]\sum _{i=1}^n\sum _{j=1}^m[\gcd (i,j)==k]$$
$$ans=\sum _{k=1}^{\infty }\sum _{d=1}^{\infty }[d==k^2]\times \sum _{i=1}^{n/d}\lfloor \frac{n}{id}\rfloor \times \lfloor \frac{m}{id}\rfloor \times \mu (i)$$
$$设f(x)=[x是平方数],\frac{それでは}{那么}有$$
$$ans=\sum _{d=1}^{\infty }f(d)\times \sum _{i=1}^{n/d}\lfloor \frac{n}{id}\rfloor \times \lfloor \frac{m}{id}\rfloor \times \mu (i)$$
$$ans=\sum _{d=1}^{\infty }\sum _{i=1}^{n/d}f(d)\times \lfloor \frac{n}{id}\rfloor \times \lfloor \frac{m}{id}\rfloor \times \mu (i)$$
$$设T=id,\frac{そして}{然后}改变求和次序$$
$$ans=\sum _{T=1}^{\infty }\lfloor \frac{n}{T}\rfloor \times \lfloor \frac{m}{T}\rfloor \sum _{i\mid T}f(\frac{T}{i})\times \mu (i)\Leftarrow 可整除分块O(\sqrt{n})$$
$$上式O(\sqrt{n})可算$$

$$设g(T)=\sum _{i\mid T}f(\frac{T}{i})\times \mu (i),g(T)可以暴力筛出来,复杂度是O(n)$$

莫比乌斯反演:${\sum }_{i=1}^n{\sum }_{j=1}^m\gcd (i,j)$

bzoj2005 ${\sum }_{i=1}^n{\sum }_{j=1}^{m}2\times \gcd (i,j)-1(n,m\in [1,100000])$ AC代码

$$ans=\sum _{i=1}^n\sum _{j=1}^m\gcd (i,j)$$
$$ans=\sum _{d=1}^{n}d\sum _{i=1}^n\sum _{j=1}^m[\gcd (i,j)==d]$$
$$ans=\sum _{d=1}^{n}d\sum _{i=1}^{n/d}\lfloor \frac{n}{id}\rfloor \times \lfloor \frac{m}{id}\rfloor \times \mu (i)$$
$$设f(k,n,m)=\sum _{i=1}^n\sum _{j=1}^m[\gcd (i,j)==k]=\sum _{i=1}^{n/k}\lfloor \frac{n}{ik}\rfloor \times \lfloor \frac{m}{ik}\rfloor \times \mu (i)$$
$$ans=\sum _{d=1}^{n}d\times f(1,\lfloor \frac{n}{d}\rfloor ,\lfloor \frac{n}{d}\rfloor )$$
$$O(n)的做法就出来了,不过通过调整求和次序可以获得O(\sqrt{n})的做法$$
$$ans=\sum _{d=1}^{n}d\sum _{i=1}^{n/d}\lfloor \frac{n}{id}\rfloor \times \lfloor \frac{m}{id}\rfloor \times \mu (i)$$
$$设T=id,\frac{そして}{然后}交换求和次序为\sum _{T=1}^n\sum _{i\mid T}$$
$$ans=\sum _{T=1}^n\lfloor \frac{n}{T}\rfloor \times \lfloor \frac{m}{T}\rfloor \times \sum _{i\mid T}\mu (i)\times \lfloor \frac{T}{i}\rfloor$$
$$上式O(\sqrt{n})可算,需要预处理出\sum _{i\mid T}\mu (i)\times \lfloor \frac{T}{i}\rfloor ,可以用线性筛O(n)时间求出$$

$$设f(T)=\mu (i)\times \lfloor \frac{T}{i}\rfloor ,那么对于线性筛有如下式子$$
$$f(1)=1,f(p)=p-1,f(p^k)=p\times f(p^{k-1})$$

莫比乌斯反演:${\sum }_{i=1}^n{\sum }_{j=1}^{m}lcm(i,j)$

bzoj2154 ${\sum }_{i=1}^n{\sum }_{j=1}^{m}lcm(i,j)(n,m\in [1,10000000])$ AC代码

$$ans=\sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j)$$
$$ans=\sum _{i=1}^n\sum _{j=1}^m\frac{ij}{\gcd (i,j)}$$
$$ans=\sum _{d=1}^n\sum _{i=1}^n\sum _{j=1}^m\frac{ij}{d}[\gcd (i,j)==d]$$
$$做代换i=id,j=jd$$
$$ans=\sum _{d=1}^n\sum _{id=1}^n\sum _{jd=1}^m\frac{idjd}{d}[\gcd (id,jd)==d]$$
$$ans=\sum _{d=1}^{n}d\sum _{i=1}^{n/d}\sum _{j=1}^{m/d}ij[\gcd (i,j)==1]$$
$$设f(k,n,m)=\sum _{i=1}^n\sum _{j=1}^{m}ij\times [\gcd (i,j)==k]$$

$$\frac{倍数反演したい}{试图倍数反演}$$
$$设g(k,n,m)=\sum _{i=1}^{\infty }f(ik,n,m)$$
$$g(k,n,m)=\sum _{i=1}^n\sum _{j=1}^{m}ij\times [\gcd (i,j)\%k==0]$$
$$\frac{それでは}{那么}g(k,n,m)可以\frac{かんたん}{轻易}算出$$
$$//做代换i=ik,j=jk$$
$$\frac{それでは}{那么}g(k,n,m)=\sum _{ik=1}^n\sum _{jk=1}^{m}ikjk\times [\gcd (ik,jk)\%k==0]$$
$$\frac{そして}{然后}g(k,n,m)=k^2\times \sum _{i=1}^{n/k}\sum _{j=1}^{m/k}ij\times [\gcd (i,j)\%1==0]$$
$$\frac{だつて}{因为}[\gcd (i,j)\%1==0]永真$$
$$\frac{だから}{所以}g(k,n,m)=k^2\times \sum _{i=1}^{n/k}\sum _{j=1}^{m/k}ij$$

$$设S(x)=1+2+...+x$$
$$\frac{それでは}{那么}S(x)=\frac{x^2+x}{2},O(1)算出$$
$$\frac{だから}{所以}g(k,n,m)=k^2\times S(\lfloor \frac{n}{k}\rfloor )S(\lfloor \frac{m}{k}\rfloor )$$
$$倍数反演成功,f(k,n,m)=\sum _{d=1}^{\infty }\mu (d)\times g(dk)$$
$$f(k,n,m)=\sum _{d=1}^{\infty }\mu (d)\times (dk{)}^2\times S(\lfloor \frac{n}{dk}\rfloor )S(\lfloor \frac{m}{dk}\rfloor )$$
$$f(k,n,m)=k^2\sum _{d=1}^{\infty }\mu (d)\times d^2\times S(\lfloor \frac{n}{dk}\rfloor )S(\lfloor \frac{m}{dk}\rfloor )$$
$$\frac{だつて}{因为}\lfloor \frac{n}{dk}\rfloor \ge 1,\frac{だから}{所以}d\le \lfloor \frac{n}{k}\rfloor$$
$$f(k,n,m)=k^2\sum _{d=1}^{n/k}\mu (d)\times d^2\times S(\lfloor \frac{n}{dk}\rfloor )S(\lfloor \frac{m}{dk}\rfloor )$$
$$上式O(\sqrt{n})可算,需要预处理\mu (d)\times d^2的前缀和,可以线性求出来O(n)$$
$$ans=\sum _{d=1}^{n}d\times f(1,\lfloor \frac{n}{d}\rfloor ,\lfloor \frac{m}{d}\rfloor )$$
$$上式O(\sqrt{n}\times \sqrt{n})可算,\frac{だから}{所以}\sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j)需要O(n)的时间$$

莫比乌斯反演:${\sum }_{i=1}^n{\sum }_{j=1}^{m}d(ij)$

bzoj3994  ∑ i = 1 n ∑ j = 1 m d ( i j )   ( T ∈ [ 1 , 50 &ThinSpace; 000 ] { n , m ∈ [ 1 , 50 &ThinSpace; 000 ] } ) \sum_{i=1}^n\sum_{j=1}^md(ij) (T\in[1,50\,000]\{n,m\in[1,50\,000]\}) ∑i=1n​∑j=1m​d(ij) (T∈[1,50000]{n,m∈[1,50000]}) AC代码

$$ans=\sum _{i=1}^n\sum _{j=1}^{m}d(ij)$$
$$ans=\sum _{i=1}^n\sum _{j=1}^m\sum _{x|i}\sum _{y|j}[\gcd (x,y)==1]$$

$$稍微枚举几项即可知:\sum _{i=1}^n\sum _{x|i}any=\sum _{x=1}^n\lfloor \frac{n}{x}\rfloor any$$
$$ans=\sum _{x=1}^n\sum _{y=1}^m\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor [\gcd (x,y)==1]$$
$$设f(k,n,m)=\sum _{x=1}^n\sum _{y=1}^m\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor [\gcd (x,y)==k]$$
$$\frac{倍数反演したい}{试图倍数反演}$$
$$设g(k,n,m)=\sum _{i=1}^{\infty }f(ik,n,m)$$
$$\frac{それでは}{那么}g(k,n,m)=\sum _{x=1}^n\sum _{y=1}^m\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor [\gcd (x,y)\%k==0]$$
$$//做代换x=xk,y=yk$$
$$g(k,n,m)=\sum _{xk=1}^n\sum _{yk=1}^m\lfloor \frac{n}{xk}\rfloor \lfloor \frac{m}{yk}\rfloor [\gcd (xk,yk)\%k==0]$$
$$g(k,n,m)=\sum _{x=1}^{n/k}\sum _{y=1}^{m/k}\lfloor \frac{n}{xk}\rfloor \lfloor \frac{m}{yk}\rfloor$$

$$设S(n)=\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor$$
$$S(n)等价于\sum _{i=1}^{n}d(i)//d(i):i的因数个数$$
$$d(i)是一个积性函数,所以可以用线性筛出来$$
$$需要一个辅助数组a(i)=i的最小因子的次幂$$
$$S(1)=1,a(1)=0,S(p)=2,a(p)=1,S(p^k)=\frac{S(p^{k-1})}{1+a(p^{k-1})}\times (2+a(p^{k-1})),a(p^k)=a(p^{k-1})+1$$
$$g(k,n,m)=S(\lfloor \frac{n}{k}\rfloor )\times S(\lfloor \frac{m}{k}\rfloor )$$
$$倍数反演成功,f(k,n,m)=\sum _{i=1}^{\infty }g(ik,n,m)\times \mu (i)$$
$$f(k,n,m)\sum _{i=1}^{\infty }\mu (i)\times S(\lfloor \frac{n}{ik}\rfloor )\times S(\lfloor \frac{m}{ik}\rfloor )$$
$$\frac{だつて}{因为}\lfloor \frac{n}{ik}\rfloor \ge 1,\frac{だから}{所以}i\le \lfloor \frac{n}{k}\rfloor$$
$$f(k,n,m)\sum _{i=1}^{n/k}\mu (i)\times S(\lfloor \frac{n}{ik}\rfloor )\times S(\lfloor \frac{m}{ik}\rfloor )$$
$$上式可O(\sqrt{n})计算,\frac{けれど}{但是}需要预处理出\mu (i)的前缀和以及S(n)的每一项,需要O(n)的时间$$

莫比乌斯反演:${\sum }_{i=1}^n[i不含有平方因子]$

bzoj2440  求 第 k 个 不 含 平 方 因 子 的 数 ( 1 不 算 )   ( T ∈ [ 1 , 50 ] { k ∈ [ 1 , 1 e 9 ] } ) 求第k个不含平方因子的数(1不算) (T\in[1,50]\{k\in[1,1e9]\}) 求第k个不含平方因子的数(1不算) (T∈[1,50]{k∈[1,1e9]}) AC代码

$$直接容斥,f(n)=\sum _{i=1}^n[i不含有平方因子]=\sum _{i=1}^{\sqrt{n}}\mu (i)\times \frac{n}{i^2}$$
$$上式O(\sqrt{n})可算,对于bzoj2440来说可以二分找到答案$$