题意
求 ∑ni=1μ(i),n<=1010
杜教筛
有一篇文必须得转
author: skywalkert
original article: http://blog.csdn.net/skywalkert/article/details/50500009
last update time : 2017-04-01
推一波式子 (下面的除法都是整除)
由
∑ni=1∑d|iμ(d)=1
可推
∑ni=1∑nij=1μ(j)(等价于枚举a=ij,给a累加μ(j))
分出i=1的情况,
∑ni=1μ(i)+∑ni=2∑nij=1μ(j)=1
也就是
Ans(n)=∑ni=1μ(i)=−∑ni=2Ans(ni)+1
其中
ni
分块。
下面分析这样递归求解的时间复杂度。
设
T(n)
为求解
Ans(n)
的复杂度,
T(n)=n√(分块)+∑n√i=1T(i)+T(ni)(求解子问题)
T(n)=n√+∑n√i=1i√+ni−−√
(只展开一层,考虑最高次项)
又因为
∑n√i=1i√+ni−−√=O(n34)
,所以总的时间复杂度就是
O(n34)
算过的不需要再算,用map存起来。不用会T比较多点。
还可以进行优化,可以先用线筛预处理出前k个的Ans,经推导k取 n23 时有最小复杂度 O(n23) 。 这样实际上不仅减少了理论复杂度,更大大缩减了hashmap的寻址时间。
核心思想是数学变形+交换主体+利用狄利克雷卷积性质
如果能通过狄利克雷卷积构造一个更好计算前缀和的函数,且用于卷积的另一个函数也易计算,则可以简化计算过程。
经试验,不加map或不预处理都会T!
#include <iostream>
#include <cstdio>
#include <cmath>
#include <map>
#define N 10000001ll
typedef long long ll;
using namespace std;
ll a,b,mu[N],ini;
map<ll,ll> h;
bool is[N];
int p[N];
void init(ll n) {
ini=n=min(n,N-1);
mu[1]=1;
for (ll i=2; i<=n; i++) {
if (is[i]==0) p[++p[0]]=i,mu[i]=-1;
for (ll j=1; j<=p[0] && i*p[j]<=n; j++) {
is[i*p[j]]=1;
if (i%p[j]) mu[i*p[j]]=-mu[i];
else {
mu[i*p[j]]=0;
break;
}
}
mu[i]+=mu[i-1];
}
}
ll solve(ll x) {
if (x<=ini) return mu[x];
if (h.count(x)) return h[x];
ll &k=h[x],ret=1,r;
for (ll l=2; l<=x; l=r+1) {
r=x/(x/l);
ret-=(r-l+1)*solve(x/l);
}
return k=ret;
}
int main() {
freopen("1244.in","r",stdin);
cin>>a>>b;
init(pow(b,2/(double)3));
printf("%lld",solve(b)-solve(a-1));
}