BZOJ4804

来源:互联网 发布:农村淘宝合伙人官网 编辑:程序博客网 时间:2024/05/01 16:58

====kφ(k)i=1nj=1n[gcd(i, j)=k]kφ(k)i=1nj=1nkd|gcd(i, j)μ(d)kφ(k)dμ(d)nkd2tk|tφ(k)μ(tk)nt2tf(t)nt2

其中,
f(n)=d|nφ(d)μ(nd)

线性筛:
f(pr)=k=0rφ(pk)μ(prk)=φ(pr)μ(1)+φ(pr1)μ(p)=φ(pr)φ(pr1)

杜教筛:
==i=1nf(i)i=1nj|iφ(j)μ(ij)i=1nφ(i)j=1n/iμ(j)

实测杜教筛远慢于线性筛

0 0
原创粉丝点击