【BZOJ4176】Lucas的数论 莫比乌斯反演
【BZOJ4176】Lucas的数论
Description
去年的Lucas非常喜欢数论题,但是一年以后的Lucas却不那么喜欢了。
在整理以前的试题时,发现了这样一道题目“求Sigma(f(i)),其中1<=i<=N”,其中 表示i的约数个数。他现在长大了,题目也变难了。
求如下表达式的值:
![](http://www.lydsy.com/JudgeOnline/upload/201507/abc.jpg)
其中 表示ij的约数个数。
他发现答案有点大,只需要输出模1000000007的值。
Input
第一行一个整数n。
Output
一行一个整数ans,表示答案模1000000007的值。
Sample Input
2
Sample Output
8
HINT
对于100%的数据n <= 10^9。
题解:前置技能:
然后直接上莫比乌斯反演
用杜教筛处理μ(d),然后喜闻乐见的分块~
#include <cstdio> #include <cstring> #include <iostream> #include <map> #define mod 1000000007 using namespace std; const int m=10000000; typedef long long ll; int n,num; ll ans; int mu[m+10],sm[m+10],pri[m+10]; bool np[m+10]; map<ll,ll> mp; ll getsm(ll x) { if(x<=m) return sm[x]; if(mp[x]) return mp[x]; ll ret=1,i,last; for(i=2;i<=x;i=last+1) { last=x/(x/i); ret=(ret-(last-i+1)*getsm(x/i)%mod+mod)%mod; } mp[x]=ret; return ret; } ll getf(ll x) { ll ret=0,i,last; for(i=1;i<=x;i=last+1) { last=x/(x/i); ret=(ret-(last-i+1)*(x/i)%mod+mod)%mod; } return ret*ret%mod; } int main() { scanf("%d",&n); ll i,j,last; sm[1]=mu[1]=1; for(i=2;i<=m;i++) { if(!np[i]) pri[++num]=i,mu[i]=-1; sm[i]=sm[i-1]+mu[i]; for(j=1;j<=num&&i*pri[j]<=m;j++) { np[i*pri[j]]=1; if(i%pri[j]==0) { mu[i*pri[j]]=0; break; } mu[i*pri[j]]=-mu[i]; } } for(i=1;i<=n;i=last+1) { last=n/(n/i); ans=(ans+(getsm(last)-getsm(i-1)+mod)%mod*getf(n/i)%mod)%mod; } printf("%lld",ans); return 0; }
| 欢迎来原网站坐坐! >原文链接<