HDU-4746 Mophues

链接：http://acm.hdu.edu.cn/showproblem.php?pid=4746

莫比乌斯反演搞，然后要改变枚举顺序，就可以发现前半部分是整除分块能搞的，后半部分可以预处理（因为数据范围5e5质因数最多的也不会超过19个），然后就可以根号回答每个询问了

具体安利sw大佬博客（Orz）：https://www.luogu.org/blog/spaly/hdu4746-mophues

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=5e5+10;
int zys[N],miu[N],prime[N/3],pnum;
bool vis[N];ll dp[20][N];
void init()
{
	int nw;
	zys[1]=0;miu[1]=1;
	for(int i=2;i<N;i++)
	{
		if(!vis[i])prime[++pnum]=i,miu[i]=-1;
		for(int j=1;j<=pnum&&prime[j]*i<N;j++)
		{
			vis[i*prime[j]]=1;
			if(i%prime[j]==0)break;
			miu[i*prime[j]]=-miu[i];
		}
	}
	for(int i=2;i<N;i++)
	{
		nw=i;
		for(int j=2;j*j<=nw;j++)
			while(nw%j==0)zys[i]++,nw/=j;
		if(nw>1)zys[i]++;
	}
	for(int i=1;i<N;i++)
	{
		nw=zys[i];
		for(int j=1;j*i<N;j++)
		{
			dp[nw][i*j]+=miu[j];
		}
	}
	for(int i=1;i<N;i++)
	{
		for(int j=1;j<=19;j++)
			dp[j][i]+=dp[j-1][i];
		for(int j=0;j<=19;j++)
			dp[j][i]+=dp[j][i-1];
	}
}
int main()
{
	int Q,n,m,p,mn;ll ans;
	init();
	scanf("%d",&Q);
	while(Q--)
	{
		scanf("%d%d%d",&n,&m,&p);
		if(p>19)printf("%lld\n",1LL*n*m);
		else
		{
			mn=min(n,m);ans=0;
			for(int l=1,r;l<=mn;l=r+1)
			{
				r=min(n/(n/l),m/(m/l));
				ans+=(dp[p][r]-dp[p][l-1])*(n/l)*(m/l);
			}
			printf("%lld\n",ans);
		}
	}
}