题目链接: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;
}