蒟蒻YZH的博客

Orz w_yqts Orz xudyh Orz popoqqq 求∑nx=1∑ny=1d(ij)\sum_{x=1}^n\sum_{y=1}^nd(ij) n<=109n<=10^9首先证明d(nm)=∑i|n∑j|m[gcd(i,j)=1]d(nm)=\sum_{i|n}\sum_{j|m}[gcd(i,j)=1] 若有n=n′∗pk1,m=m′∗pk2n=n'*p^{k1},m=m'

给定n,m，求模10^9+7的值。 听说是杜教筛裸题… 可能是我太鶸了.感觉好难啊.. Orz w_yqts Orz xudyh (i,j)=gcd(i,j)(i,j)=gcd(i,j) s(i,j)=∑i=1mφ(ni)s(i,j)=\sum_{i=1}^mφ(ni) ans=∑i=1ns(i,m)ans=\sum_{i=1}^n s(i,m) 对于n%p=0,有φ(n*p)=φ(

Orz w_yqts Orz xudyh 杜教筛模板题 给定积性函数f(x) 函数s(x)为f(x)的前缀和 f(x)满足Σf(d)=g(x),d为x的约数,g(x)的前缀和可以在较短时间内算出. 求s(n) 从枚举约数到枚举倍数,将离散变为连续… 上面的j枚举的就是倍数,i是当前倍数可以支持的约数,不难发现,可以支持的约数是连续的. 预处理前n^(2/3),后面直接用map记