HDU - 4746 Mophues （莫比乌斯反演+分块）

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

题目大意：题目定义如果数C是数P的幸运数的话，当且仅当C的质因子数小于等于P。现在给出n，m和p，问你能选多少对组合（a,b）满足1<=a<=n，1<=b<=m，同时gcd(a,b)是P的幸运数。

题目思路：对于这个题目，我们可以很直观的列出以下两个式子：

我们令表示同时满足是P的幸运数的数量，令表示是的倍数同时满足是P的幸运数的数量。

那么，且

根据莫比乌斯反演可得。

按照通常的套路，接下来直接进行枚举即可得到如下的式子

但这题由于组数较大，我们无法用这种常规的方法来解决这个题。所以我们要进行分块加速。

我们令，那么上面的式子就可以化为

。

我们现在想快速的求出的值，我们可以用一个前缀和来记录表示前 i 项中，质因子数小于等于 j 莫比乌斯系数的前缀和，通过打表可以看出在5e5的数据范围内，一个数最多只有19个质因子，所以只需要前缀和记录 j<=19的情况即可，当P大于19时就是n*m了。

由于这里的除法都是向下取整的整数除法，所以我们知道对于所有的，的值都是相同的，在这里可以借助这个性质进行分段优化，借助上面求出的前缀和快速求解。

具体实现看代码：

#include <bits/stdc++.h>
#define fi first
#define se second
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define lowbit(x) x&-x
#define pb push_back
#define MP make_pair
#define clr(a) memset(a,0,sizeof(a))
#define _INF(a) memset(a,0x3f,sizeof(a))
#define FIN freopen("in.txt","r",stdin)
#define fuck(x) cout<<"["<<x<<"]"<<endl
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int>pii;
//head
const int MX=5e5+7;

bool isprime[MX];
int num[MX],sum[MX][22];
int prime[MX],mu[MX]= {0,1},primecount;
void mubi() {
    memset(isprime,true, sizeof(isprime));
    primecount = 0;
    for (int i = 2; i < MX; ++i) {
        if (isprime[i]) {
            prime[++primecount] = i;
            num[i] = 1;
            mu[i] = -1;
        }
        for (int j = 1; i * prime[j] <MX;  ++j) {
            isprime[i * prime[j]] = false;
            num[i * prime[j]] = num[i] + 1;
            if(i%prime[j]==0) {
                mu[i*prime[j]] = 0;
                break;
            } else mu[i*prime[j]] = -mu[i];
        }
    }
}
void init(){
	mubi();
	for(int i = 1;i < MX; i++){
    	for(int j = i;j < MX; j+=i){
    		sum[j][num[i]] += mu[j/i];
    	}
    }
    for(int i = 1;i < MX; i++){
    	for(int j = 1;j < 20; j++){
    		sum[i][j] += sum[i][j-1];
    	}
    }
    for(int i = 1;i < MX; i++){
    	for(int j = 0;j < 20; j++)
    		sum[i][j] += sum[i-1][j];
    }
}
int _,n,m,p;

int main(){
	//FIN;
	init();
	for(scanf("%d",&_);_;_--){
		scanf("%d%d%d",&n,&m,&p);
		if(n>m) swap(n,m);
		ll ans=0;
		p=min(p,19);
		for(int i=1,j;i<=n;i=j+1){
			j=min(n/(n/i),m/(m/i));
			ans+=(ll)(n/j)*(m/j)*(sum[j][p]-sum[i-1][p]);
		}
		printf("%lld\n",ans);
	}
	return 0;
}