bzoj3944 Sum 杜教筛

模板题，注意杜教筛化式子的技巧：

1、构造求和函数

2、利用积性函数约数性质

3、做差算出单项

4、约数倍数转化贡献

5、分块求商

码：

#include<iostream>
#include<cstdio>
#include<map>
#include<cmath>
using namespace std;
#define M 3000000
#define ll long long
int tot,i,j,T,su[M+5];
ll ans,phe[M+5],n[999],N,mhe[M+5],mu[M+5],phi[M+5];
bool he[M+5];
map<int,ll>pma;
map<int,int>mma;
void eular()
{mu[1]=phi[1]=1;for(i=2;i<=N;i++){if(!he[i]){su[++tot]=i;           phi[i]=i-1;mu[i]=-1;}for(j=1;j<=tot&&su[j]*i<=N;j++){int k=su[j]*i;he[k]=1;if(i%su[j]==0){phi[k]=su[j]*phi[i];mu[k]=0;break;}else{phi[k]=phi[i]*phi[su[j]];mu[k]=-mu[i];}  }}for(i=1;i<=N;i++)phe[i]=phe[i-1]+phi[i],mhe[i]=mhe[i-1]+mu[i];
}
ll dfs1(ll o)
{if(o<=N)return phe[o]; if(pma[o])return pma[o];
ll lin=o*(1+o)/2;
ll l;
for(l=2;l<=o;l++)
{
ll oo=o/l;
lin-=((o/oo)-l+1)*dfs1(oo) ;
l=o/oo;
}
pma[o]=lin;
return lin;
}
int dfs2(int o)
{   if(o<=N)return mhe[o]; if(mma[o])return mma[o];
int lin=1;
unsigned int l;
for(l=2;l<=o;l++)
{
int oo=o/l;
lin-=((o/oo)-l+1)*dfs2(oo);
l=o/oo;
}
mma[o]=lin;
return lin;
}
int main()
{scanf("%d",&T);
for(i=1;i<=T;i++)
{
scanf("%lld",&n[i]);
N=max(n[i],N);
}
N=pow(N,2.000/3) ;
eular();
for(i=1;i<=T;i++)
{
printf("%lld ",dfs1(n[i]));
printf("%d\n",dfs2(n[i]));
}
}

bzoj3944 Sum 杜教筛相关推荐

bzoj3944: Sum//杜教筛
bzoj3944: Sum 前几天的blog,没写完,现在补上吧- 题意求莫比乌斯函数和欧拉函数的前缀和.N<=2^31-1.T(<=10)组询问. 题解杜教筛. 然而这篇博客不是写题 ...
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 杜教筛相关推荐

最新文章

热门文章