bzoj3944: Sum//杜教筛

bzoj3944: Sum
前几天的blog，没写完，现在补上吧…

题意

求莫比乌斯函数和欧拉函数的前缀和。N<=2^31-1。T(<=10)组询问。

题解

杜教筛。
然而这篇博客不是写题解的。这主要是因为
1.网上有很多详细题解。
2.不想打公式。
3.表达能力拙计。

注意事项

今天真是涨姿势…居然还有这种挂题方法…
首先CE的时候显示 File size limit exceeded。问了一圈，不解其意…
改了几个地方，把N调小了一点，莫名过了编译，但是RE了？
查错，发现l会变成负的，于是从

long long l l=(long long)n*(n+1)>>1;

改成了

long long l=(long long)n*(long long)(n+1)>>1;

诶怎么还是负的？想了想发现自己弱智……
再交，又CE，还是File size limit exceeded（FSLE??），~~二分了几发~~试了一下内存，发现好像是MLE①的另一种表达方式……真是十分奇怪的体验

UPD (2018.3.28)：昨天又CE了一道题后…发现File size limit exceeded是因为数组用花括号初始化会占内存 QAQ

①：bz上MLE的表达方法还有TLE…大概

代码

#include<cstdio>
#include<map>
#define N 5125000
using namespace std;
typedef long long ll;
typedef pair<ll,int> pr;
pr b;
ll f[N]={0,1};
int T,a,p[N],g[N]={0,1},tot;
bool pri[N];
map<int,pr>m;
map<int,pr>::iterator s;
pr G(int n)
{if(n<N)return make_pair(f[n],g[n]);if((s=m.find(n))!=m.end())return s->second;ll l=n*((ll)n+1)>>1,t;int r=1;for(ll i=2;i<n;i=t)t=n/(n/i)+1,b=G(n/i),l-=(t-i)*b.first,r-=(t-i)*b.second;return m[n]=make_pair(l,r);
}
int main()
{for(int i=2;i<N;i++){if(!pri[i])p[tot++]=i,g[i]=-1,f[i]=i-1;for(int j=0;i*p[j]<N;j++){pri[a=i*p[j]]=1,f[a]=f[i]*p[j];if(!(i%p[j]))break;g[i*p[j]]=-g[i],f[a]-=f[i];}g[i]+=g[i-1],f[i]+=f[i-1];}scanf("%d",&T);while(T--)scanf("%d",&a),b=G(a),printf("%lld %d\n",b.first,b.second);
}

bzoj3944: Sum//杜教筛相关推荐

bzoj3944 Sum 杜教筛
模板题,注意杜教筛化式子的技巧: 1.构造求和函数 2.利用积性函数约数性质 3.做差算出单项 4.约数倍数转化贡献 5.分块求商码: #include<iostream> #inclu ...
BZOJ3944 Sum 杜教筛
看完题一副不可做的样子默默点开了题解发现是杜教筛就花了半天学习了一下先说下杜教筛可以在优于线性的复杂度内求出积性函数的前缀和下求 ∑ni=1f(i) \sum_{i=1}^nf(i) 令 ...
bzoj 3944: Sum 杜教筛
本来以为这种东西只能O(N)线性筛,但是大千世界,无(sang)奇(xin)不(bing)有(kuang),确实存在更快的算法. 省选的时候rzz讲这种东西在国内OI称为杜教筛,用来求数论函数的前缀和 ...
2019-ACM-ICPC-南京区网络赛-E. K Sum（莫比乌斯反演 + 杜教筛）
K Sum 推式子 Fn(k)=∑l1=1n∑l2=1n⋯∑lk=1n(gcd(l1,l2,-,lk))2=∑d=1nd2∑l1=1nd∑l2=1nd⋯∑lk=1nd(gcd(l1,l2,-,lk)= ...
杜教筛模板（P4213 【模板】杜教筛（Sum））
P4213 [模板]杜教筛(Sum) 套路推式子求s(n)=∑i=1nf(i)∑i=1n(f∗g)(i)=∑i=1n∑d∣if(d)g(id)=∑d=1n∑i=1⌊nd⌋f(i)g(d)=∑d=1n ...
P4213-[模板]杜教筛(Sum)
正题题目链接:https://www.luogu.com.cn/problem/P4213 题目大意给出nnn,求∑i=1nφ(i)\sum_{i=1}^n\varphi(i)i=1∑nφ(i) ...
luoguP4213 【模板】杜教筛（Sum）杜教筛
链接 luogu 思路为了做hdu来学杜教筛. 杜教筛模板题. 卡常数,我加了register居然跑到不到800ms. 太深了. 代码 // luogu-judger-enable-o2 #incl ...
【bzoj4176】Lucas的数论莫比乌斯反演+杜教筛
题目描述去年的Lucas非常喜欢数论题,但是一年以后的Lucas却不那么喜欢了. 在整理以前的试题时,发现了这样一道题目"求Sigma(f(i)),其中1<=i<=N" ...
数论函数 - 莫比乌斯函数与莫比乌斯反演 - 基础杜教筛
原文链接http://www.cnblogs.com/zhouzhendong/p/8627380.html 省选后发现我数学好差.于是先从数论开始学习. 如果发现本文有任何错误,欢迎留言指正. 本文 ...

bzoj3944: Sum//杜教筛

题意

题解

注意事项

代码

bzoj3944: Sum//杜教筛相关推荐

最新文章

热门文章