[BZOJ3512] DZY Loves Math IV

传送门
昨天花了好久好像是看懂了，那今天早上尝试自己推一遍柿子 顺便水了一篇博客

$\mathbf{D}\mathbf{e}\mathbf{s}\mathbf{c}\mathbf{r}\mathbf{i}\mathbf{p}\mathbf{t}\mathbf{i}\mathbf{o}\mathbf{n}$

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

$1\le n\le 10^5,1\le m\le 10^9$，模 $10^9+7$

$\mathbf{S}\mathbf{o}\mathbf{l}\mathbf{u}\mathbf{t}\mathbf{i}\mathbf{o}\mathbf{n}$

先给个结论，当 $|\mu (n)|=1$ 时，有

$$\phi (ni)=\phi (i)\sum _{d|(n,i)}\phi (\frac{n}{d})$$

简单证明一下（想了一节晚自习才想通呢QAQ）
（gcd的表示有点乱，不要在意）

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

推的时候要注意 $n$ 的性质，还要熟悉 $\phi$，不然就推不粗来QAQ

然后令 $S(n,m)={\sum }_{i=1}^m\phi (ni)$，令 $P$ 为 $n$ 所有质因子的乘积，$Q=\frac{n}{P}$，然后可得

$$S(n,m)=Q\cdot \sum _{i=1}^m\phi (i)\sum _{d|(P,i)}\phi (\frac{P}{d})$$

先枚举因数，得到

$$S(n,m)=Q\cdot \sum _{d|P}\phi (\frac{P}{d})\sum _{i=1}^{\lfloor \frac{m}{d}\rfloor }\phi (di)$$
$$=Q\cdot \sum _{d|P}\phi (\frac{P}{d})S(d,\lfloor \frac{m}{d}\rfloor )$$

然后就阔以记忆化啦？
当 $n=1$ 的时候就是 $\phi$ 的前缀和，杜教筛板子啦？
复杂度玄学？不会证啦。
代码？还没写出来啦。
写出来再贴吧。

---

写出来啦，十分暴力QAQ

//#pragma GCC optimize(2)
//#pragma GCC optimize(3)
#include<bits/stdc++.h>
#include<tr1/unordered_map>
#define LL long long
#define re register
#define fr(i,x,y) for(int i=(x);i<=(y);i++)
#define rf(i,x,y) for(int i=(x);i>=(y);i--)
#define frl(i,x,y) for(int i=(x);i<(y);i++)
using namespace std;
using namespace tr1;
const int N=2000002;
const int INF=2147483647;
const int p=1e9+7;
int n,m;

void read(int &x){ scanf("%d",&x); }

void Add(int &x,int y){
	x+=y;
	while(x<0) x+=p;
	while(x>=p) x-=p;
}

int phi[N];
int pri[N/10],b[N],L;
void init(){
	phi[1]=1;
	frl(i,2,N){
		if (!b[i]) pri[++L]=i,phi[i]=i-1;
		for(int j=1;j<=L&&i*pri[j]<N;j++){
			b[i*pri[j]]=1;
			if (i%pri[j]==0){
				phi[i*pri[j]]=phi[i]*pri[j];
				break;
			}
			phi[i*pri[j]]=phi[i]*(pri[j]-1);
		}
	}
	frl(i,2,N) Add(phi[i],phi[i-1]);
}

map<int,int> s,mp[N];
int S(int n,int m){
	//cout<<n<<' '<<m<<endl;
	if (n==0||m==0) return 0;
	if (n==1){
		if (m<N) return phi[m];
		if (s.count(m)) return s[m];
		int ans=(1LL*m*(m+1)/2)%p;
		for(re int L=2,r=2;L<=m;L=r+1)
		 r=m/(m/L),Add(ans,-1LL*S(n,m/L)*(r-L+1)%p);
		return s[m]=ans;
	}
	if (mp[n].count(m)) return mp[n][m];
	int ans=0;
	for(re int i=1;i*i<=n;i++)
	 if (n%i==0){
	 	Add(ans,1LL*(phi[n/i]-phi[n/i-1])*S(i,m/i)%p);
	 	if (n/i!=i) Add(ans,1LL*(phi[i]-phi[i-1])*S(n/i,m/(n/i))%p);
	 }
	return mp[n][m]=ans;
}

int main(){
	//freopen("1.in","r",stdin);
	init();
	read(n);read(m);
	int ans=0;
	fr(xx,1,n){
		int s=1,x=xx;
		for(int i=2;i*i<=x;i++)
		 if (x%i==0){
		 	s*=i;
		 	while(x%i==0) x/=i;
		 }
		s*=x;
		Add(ans,1LL*S(s,m)*(xx/s)%p);
		//cout<<S(s,m)*(xx/s)<<endl;
	}
	cout<<ans<<endl;
	//cout<<S(s,m)<<endl;
	//printf("%lld\n",1LL*S(s,m)*(n/s)%p);
	return 0;
}