【bzoj3561】DZY Loves Math VI

题意：（第n次）略

（依然是）lych大大的题解

（最近越来越懒了）

吐槽：我一开始用T=pd带掉了，结果那一串函数发现不是积性的（似乎）。。于是gg了。后来发现想太多了（居然不是套路）。。

#include <bits/stdc++.h>
#define gc getchar()
#define ll long long
#define N 500009
#define mod 1000000007
using namespace std;
int n,m,mu[N],pri[N],pd[N],cnt,ans,sum[N],Ans,a[N];
int read()
{
	int x=1;
	char ch;
	while (ch=gc,ch<'0'||ch>'9') if (ch=='-') x=-1;
	int s=ch-48;
	while (ch=gc,ch>='0'&&ch<='9') s=s*10+ch-48;
	return s*x;
}
int ksm(int y,int x)
{
	int ans=1;
	for (;x;x>>=1,y=(ll)y*y%mod)
		if (x&1) ans=(ll)ans*y%mod;
	return ans;
}
int main()
{
	n=read(),m=read();
	if (n>m) swap(n,m);
	mu[1]=pd[1]=1;
	for (int i=2;i<N;i++)
	{
		if (!pd[i]) pri[++cnt]=i,mu[i]=-1;
		for (int j=1;j<=cnt&&i*pri[j]<N;j++)
		{
			pd[i*pri[j]]=1;
			if (i%pri[j]==0)
			{
				mu[i*pri[j]]=0;
				break;
			}
			mu[i*pri[j]]=-mu[i];
		}
	}
	for (int i=1;i<=m;i++) a[i]=1;
	for (int d=1;d<=n;d++)
	{
		ans=0;
		for (int p=1;p*d<=m;p++)
		{
			a[p]=(ll)a[p]*p%mod;
			sum[p]=((ll)sum[p-1]+a[p])%mod;
		}
		for (int p=1;p*d<=n;p++)
			if (mu[p]) ans=((ll)ans+(ll)mu[p]*a[p]*a[p]%mod*sum[n/p/d]%mod*sum[m/p/d]%mod+mod)%mod;
		ans=(ll)ans*ksm(d,d)%mod;
		Ans=((ll)Ans+ans)%mod;
	}
	printf("%d\n",Ans);
	return 0;
}