HDU

链接

求[1,n] 和 [1,m]中有多少对数的GCD的素因子个数小于等于p

直接暴力做特定超时，所以我们想办法预处理，对于p大于18（1到5e5的最大素数因子个数）的情况，每一对都满足条件，O（1）得结果。

p<=18时，预处理sum【i】【j】sum[i][j]表示所有能整除i的质因子个数<=j的(x,y)的对数，然后求一个前缀和，ans=∑⌊n/d⌋∗⌊m/d⌋∗∑k|dμ(d/k)

//#pragma comment(linker, "/stack:200000000")
//#pragma GCC optimize("Ofast,no-stack-protector")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
//#pragma GCC optimize("unroll-loops")
#include<bits/stdc++.h>
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define pi acos(-1.0)
#define ll long long
#define mod 1000000007
#define C 0.5772156649
#define ls l,m,rt<<1
#define rs m+1,r,rt<<1|1
#define pil pair<int,ll>
#define pii pair<int,int>
#define ull unsigned long long
#define base 1000000000000000000
#define fio ios::sync_with_stdio(false);cin.tie(0)

using namespace std;

const double g=10.0,eps=1e-12;
const int N=500000+10,maxn=400000+10,inf=0x3f3f3f3f;

int mu[N],prime[N],sum[N][20];
bool mark[N];
int num[N];
void init()
{
    mu[1]=1;
    int cnt=0;
    for(int i=2;i<N;i++)
    {
        if(!mark[i])prime[++cnt]=i,mu[i]=-1,num[i]=1;
        for(int j=1;j<=cnt;j++)
        {
            int t=i*prime[j];
            if(t>N)break;
            mark[t]=1;
            num[t]=num[i]+1;
            if(i%prime[j]==0){mu[t]=0;break;}
            else mu[t]=-mu[i];
        }
    }
    for(int i=1;i<N;i++)
        for(int j=i;j<N;j+=i)
            sum[j][num[i]]+=mu[j/i];
    for(int i=1;i<N;i++)
        for(int j=1;j<19;j++)
            sum[i][j]+=sum[i][j-1];
    for(int i=1;i<N;i++)
        for(int j=0;j<19;j++)
            sum[i][j]+=sum[i-1][j];
}
int main()
{
    init();
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int n,m,p;
        scanf("%d%d%d",&n,&m,&p);
        if(p>18)
        {
            printf("%lld
",(ll)n*m);
            continue;
        }
        ll ans=0;
        for(int i=1,last;i<=min(n,m);i=last+1)
        {
            last=min(n/(n/i),m/(m/i));
            ans+=(ll)(sum[last][p]-sum[i-1][p])*(n/i)*(m/i);
        }
        printf("%lld
",ans);
    }
    return 0;
}
/********************

********************/