[BZOJ3512]DZY loves Math IV-数论

【题目地址】

简易题意

给定$n,m$，求下面式子在$\mathrm{m}\mathrm{o}\mathrm{d}10^9+7$意义下的值：

$$\sum _{i=1}^n\sum _{j=1}^m\phi (ij)$$

和这道题很像，但是这里$n⩽10^5,m⩽10^9$，数据范围完全不同了，而且只有单次询问，所以我们考虑$n$很小，所以从后半部分下手。

我们令$S(n,m)=\sum _{i=1}^m\phi (in)$，那么最后的式子就是$\sum _{i=1}^{n}S(i,m)$，所以考虑如何快速求取$S$函数。

现在只看这个式子：

$$\sum _{i=1}^m\phi (in)$$

考虑如何把$\phi (in)$拆开，根据求取$\phi$的式子：

$$\begin{gathered} n=\prod _{p_i|n}{p}_i^{c_i} \\ \phi (n)=n\prod _{p_i|n}\frac{p_i-1}{p_i} \end{gathered}$$

我们提取每种质因子的一个出来，那么就变成了：

$$\begin{gathered} a=\prod _{p_i|n}{p}_i \\ b=\prod _{p_i|n}{p}_i^{c_i-1} \\ n=a\times b \\ \phi (n)=\phi (a)\times b \end{gathered}$$

因为在例如$gcd(i,p)̸=1$时，$\phi (i\times p)=\phi (i)\times p$，而每个质因子之间可以分开独立考虑，所以上面式子成立。

而$\phi (i\times p)=\phi (i)\times p$的证明看这里：

摘取自友链博客【IN-there-by hdxrie】

所以我们把$n$变成$a,b$后，原式变成如下：

$$\begin{array}{rl}S(n,m)=& b\sum _{i=1}^m\phi (ai)\\ =& b\sum _{i=1}^m\phi (\frac{a}{gcd(i,a)})\phi (i)gcd(i,a)\end{array}$$

由于$a$中全部是一次质因子的乘积，所以：

$$\begin{array}{rl}\phi (ai)=& ai\prod _{p_k|ai}\frac{p_k-1}{p_k}\\ =& ai\prod _{p_k|i}\frac{p_k-1}{p_k}\prod _{p_k|a,p_k\nmid i}\frac{p_k-1}{p_k}\\ =& i\prod _{p_k|i}\frac{p_k-1}{p_k}\frac{a}{gcd(i,a)}\prod _{p_k|\frac{a}{gcd(i,a)}}\frac{p_k-1}{p_k}\times gcd(i,a)\\ =& \phi (i)\phi (\frac{a}{gcd(i,a)})gcd(i,a)\end{array}$$

这里，我们考虑将其拆开，一部分是属于$i$的质因子，而只属于$a$的也就是除去公共质因子的部分，相当于属于$\frac{a}{gcd(i,a)}$的，但是最后还少了个$gcd(i,a)$，再乘上即可。而且这里由于$a$的所有质因子都只出现一次，所以除以了$gcd(i,a)$后肯定与$i$互质，不会有重复的$\frac{p_k-1}{p_k}$。

然后我们根据$\phi \ast 1=id$，将后面的$gcd(i,a)={\sum }_{d|gcd(i,a)}\phi (d)$，所以：

$$\begin{array}{rl}S(n,m)=& b\sum _{i=1}^m\phi (\frac{a}{gcd(i,a)})\phi (i)\sum _{d|gcd(i,a)}\phi (d)\end{array}$$

由于$a$中每个质因子只出现一次，所以我们可以将$\frac{a}{gcd(i,a)}$乘到${\sum }_{d|gcd(i,a)}\phi (d)$中的每一项，此时，由于$\frac{a}{gcd(i,a)}$肯定与每个$d$互质，根据积性函数的定义：

$$\frac{a}{gcd(i,a)}\sum _{d|gcd(i,a)}\phi (d)=\sum _{d|gcd(i,a)}\phi (\frac{ad}{gcd(i,a)})$$

其实就相当于等于：

$$\sum _{d|gcd(i,a)}\phi (\frac{a}{\frac{gcd(i,a)}{d}})=\sum _{d|gcd(i,a)}\phi (\frac{a}{d})$$

因为枚举$\frac{gcd(i,a)}{d}$，$d$是一样的。

所以我们继续化简：

$$\begin{array}{rl}S(n,m)=& b\sum _{i=1}^m\phi (i)\sum _{d|gcd(i,a)}\phi (\frac{a}{d})\end{array}$$

我们交换枚举顺序即可得到：

$$\begin{array}{rl}S(n,m)=& b\sum _{d|a}\phi (\frac{a}{d})\sum _{d|i}^m\phi (i)\\ =& b\sum _{d|a}\phi (\frac{a}{d})\sum _{i=1}^{\lfloor \frac{m}{d}\rfloor }\phi (di)\end{array}$$

然后可以发现后面部分就是一个子问题了，所以递归求就好啦。

$$\begin{array}{rl}S(n,m)=& b\sum _{d|a}\phi (\frac{a}{d})S(d,\lfloor \frac{m}{d}\rfloor )\end{array}$$

复杂度和杜教筛类似。

代码~~

#include<vector>
#include<map>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
#define lowbit(a) ((a)&(-(a)))
using namespace std;
const int M=1e6+10,N=1e5+10,S=1<<20;
const ll Mod=1e9+7;
ll prime[N],phi[M],F[M],last[M],cnt;
bool vis[M];
int bit[S];
void init(){
	phi[1]=1;
	bit[1]=0;
	for(int i=1;i<=19;i++)bit[1<<i]=i;
	for(ll i=2;i<M;i++){
		if(!vis[i]){
			prime[++cnt]=i;
			phi[i]=i-1;last[i]=1;\\last记录上次的质因子
		}
		for(ll j=1,v;j<=cnt&&i*prime[j]<M;j++){
			v=i*prime[j];
			vis[v]=1;last[v]=i;
			if(!(i%prime[j])){
				phi[v]=phi[i]*prime[j];
				break;
			}
			phi[v]=phi[i]*phi[prime[j]];
		}
	}
	for(int i=1;i<M;i++)F[i]=(F[i-1]+phi[i])%Mod;
}
map <ll,ll> mp_phi;
map <ll,ll> mp_cal[N];
ll calc_phi(ll n){
	if(n<M) return F[n];//杜教筛求phi的前缀和
	if(mp_phi.count(n)) return mp_phi[n];
	ll ans=((n*(n+1ll))>>1)%Mod;
	for(ll i=2,j;i<=n;i=j+1){
		j=(n/(n/i));
		ans=(ans-((j-i+1)%Mod)*calc_phi(n/i)%Mod)%Mod;
	}
	if(ans<0)ans+=Mod;
	return mp_phi[n]=ans;
}
map <ll,ll> mp_p;
ll Getp(ll a){
	if(a<M) return phi[a];
	if(mp_p.count(a)) return mp_p[a];
    ll t=a,tmp=a;
    for(ll j=2;j*j<=a;j++){
        if(t%j==0) {
            a=a/j*(j-1);
            while(t%j==0) t/=j;
        }
    }if(t>1) a=a/t*(t-1);
    return mp_p[tmp]=a;	
}//求单独一个phi
ll calc(ll n,ll m){
	if(n==1) return calc_phi(m);//n=1时就是phi的前缀和
	if(m==1) return Getp(n);//m=1就是一个phi
	if(!m||!n) return 0;//为0就返回0
	if(mp_cal[n].count(m)) return mp_cal[n][m];//记忆化
	vector <ll> div;div.clear();
	ll ans=0,x=1,y,w=1,tc=n;
	while(n>1){
		y=n/last[n];w*=y;n=last[n];
		div.push_back(y);
		while(!(n%y))n/=y,x*=y;
	}
	//w为a，x为b
	n=tc;
	int sze=div.size(),up=1<<sze;
	for(int i=0;i<up;i++){
		ll d=1ll;
		for(int j=i;j>0;j-=lowbit(j))d*=div[bit[lowbit(j)]];
		ans=(ans+(Getp(w/d)*calc(d,m/d))%Mod)%Mod;
	}//2进制枚举约数的质因子组成
	ans=(ans*x)%Mod;//最后乘以b
	return mp_cal[n][m]=ans;
}
int n,m;ll ans;
int main(){
	init();
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)ans=(ans+calc(i,m))%Mod;//累加
	printf("%lld\n",ans);
	return 0;
}

另一种实现方式，较为简单，原理差不多：
挂个链接Orz%%%-IN

---

End