bzoj3512 DZY Loves Math IV 欧拉函数+杜教筛

litble

于 2019-03-04 17:17:52 发布

阅读量324

点赞数 2

分类专栏：数学文章标签：数学欧拉函数杜教筛

本文链接：https://blog.csdn.net/litble/article/details/88126095

版权

数学专栏收录该内容

74 篇文章 3 订阅

订阅专栏

题目分析

好难不会QAQ。

设 $S(n,m)=\sum_{i=1}^m \phi(im)$

设 $w$ 为 $n$ 的所有质因子的（一次方）乘积， $v=\frac{n}{w}$ ， $d = g c d (w, i)$ ，根据欧拉函数的性质，则有：

$S(n,m)=v\sum_{i=1}^m \phi(iw)=v\sum_{i=1}^m \phi(i) \phi(\frac{w}{d})d$

知 $\sum_{t|d} \phi(t)=d$

则 $\sum_{i=1}^m \phi(i) \phi(\frac{w}{d}) \sum_{t|d} \phi(t)=v \sum_{i=1}^m \phi(i) \phi(\frac{w}{d}) \sum_{t|d} \phi(\frac{d}{t})$

看向 $\phi(\frac{w}{d})\sum_{t|d} \phi(\frac{d}{t})$ 这个玩意，因为 $w$ 的每种质因子都只存在一个，所以 $\frac{w}{d}$ 不含 $d$ 中的任何质因子，所以 $\frac{w}{d}$ 和 $\frac{d}{t}$ 互质，所以这玩意等于 $\sum_{t|w且t|i}\phi(\frac{w}{t})$ 。

所以 $\sum_{i=1}^m \phi(i)\sum_{t|w且t|i}\phi(\frac{w}{t})=v\sum_{t|w}\phi(\frac{w}{t})\sum_{i=1}^{\frac{w}{t}} \phi(it)$

也就是 $S(n,m)=v\sum_{t|m} \phi(\frac{w}{t})S(t,\frac{w}{t})$

当 $n = 1$ 的时候用杜教筛，否则递归处理。

复杂度咋证我也不知道。

代码

#include<bits/stdc++.h>
using namespace std;
#define RI register int
const int mod=1e9+7,N=5000000;
int n,m,ans,tot,pri[N+5],is[N+5],phi[N+5],sumphi[N+5],w[N+5];
typedef pair<int,int> PR;
map<int,int> mp2;
map<PR,int> mp1;

int qm(int x) {return x>=mod?x-mod:x;}
void prework() {
	phi[1]=w[1]=1;
	for(RI i=2;i<=N;++i) {
		if(!is[i]) pri[++tot]=i,phi[i]=i-1,w[i]=i;
		for(RI j=1;j<=tot&&pri[j]*i<=N;++j) {
			is[pri[j]*i]=1;
			if(i%pri[j]==0) {
				phi[i*pri[j]]=1LL*phi[i]*pri[j]%mod;
				w[i*pri[j]]=w[i];break;
			}
			else {
				phi[i*pri[j]]=1LL*phi[i]*phi[pri[j]]%mod;
				w[i*pri[j]]=1LL*w[i]*pri[j]%mod;
			}
		}
	}
	for(RI i=1;i<=N;++i) sumphi[i]=qm(sumphi[i-1]+phi[i]);
}
int du_seive(int n) {
	if(n<=N) return sumphi[n];
	if(mp2.count(n)) return mp2[n];
	int re=(1LL*n*(n+1)/2)%mod;
	for(RI i=2,j;i<=n;i=j+1)
		j=n/(n/i),re=qm(re-1LL*(j-i+1)*du_seive(n/i)%mod+mod);
	mp2[n]=re;return re;
}
int work(int n,int m) {
	if(!n||!m) return 0;
	if(n==1) return du_seive(m);
	if(m==1) return phi[n];
	if(mp1.count((PR){n,m})) return mp1[(PR){n,m}];
	int re=0;
	for(RI i=1;i*i<=w[n];++i) if(w[n]%i==0) {
		re=qm(re+1LL*phi[w[n]/i]*work(i,m/i)%mod);
		if(i!=w[n]/i) re=qm(re+1LL*phi[i]*work(w[n]/i,m/(w[n]/i))%mod);
	}
	re=1LL*re*(n/w[n])%mod,mp1[(PR){n,m}]=re;return re;
}
int main()
{
	scanf("%d%d",&n,&m);
	prework();
	for(RI i=1;i<=n;++i) ans=qm(ans+work(i,m));
	printf("%d\n",ans);
	return 0;
}

litble

关注

2
点赞
踩
1

收藏

觉得还不错? 一键收藏
0
评论
bzoj3512 DZY Loves Math IV 欧拉函数+杜教筛

题目分析好难不会QAQ。设S(n,m)=∑i=1mϕ(im)S(n,m)=\sum_{i=1}^m \phi(im)S(n,m)=∑i=1mϕ(im)设www为nnn的所有质因子的（一次方）乘积，v=nwv=\frac{n}{w}v=wn，d=gcd(w,i)d=gcd(w,i)d=gcd(w,i)，根据欧拉函数的性质，则有：S(n,m)=v∑i=1mϕ(iw)=v∑i=1mϕ(i)ϕ...
复制链接

扫一扫