洛谷 P3455 [POI2007]ZAP-Queries

题目链接：ZAP-Queries

题目大意：给你b,d,k四个数，问[1,b]这个区间里面的x和[1,d]这个区间里面的y有多少对gcd(x,y) == k

题目思路：跟hdoj1695类似，这个在[1,b]和[1,d]这个区间里有多少对gcd(x,y) == k可以转化为[1,b/k]和[1,d/k]这个区间里有多少对gcd(x,y) == 1，因为如果gcd(x,y)=1,则gcd(x∗k,y∗k)=k因为如果gcd(x,y)=1,则gcd(x∗k,y∗k)=k，这个很容易得证，所以我们转化为两个区间内gcd(x,y) == 1有多少对，那么我们可以构造两个函数F（t）代表gcd(x,y) == t的倍数有多少对（x,y），f（t）代表gcd(x,y) == t有多少对(x,y)，然后可知$F(n)=\sum_{d|n} f(d)$ 然后根据莫比乌斯第二条公式反演一下可知$f(n)=\sum_{d|n}\mu(\frac{n}{d}) F(d)$ 然后我们知道的是要求f(1)为多少，所以n为1，并且$F(n) = \frac{x}{n}*\frac{y}{n}$,x和y为两个区间的上限，这个很明显，然后我们就可以直接套用公式去反演这个f(1)啦，然后注意！！！！！这题卡常卡的飞起，我们需要先用分块优化一下，因为我们在计算f(1)的时候用了累加和，实际上是不需要的，因为n/d的个数只有根号个，所以，在某一段区间里面n/i是一样的，这时候我们用（n/(n/i)）这个分块优化去算就好了（期间需要预处理mu[i]的前缀和），当n m k分别为10 10 1的时候，下面的j为1，2，3，5，10，因为4和5的值相同，6,7,8,9,10的值相同，这样可以根号计算，然后改一下long long为Int就好了（long long别开多）然后就写完这个傻逼卡常题了

#include <map>
#include <set>
#include <cmath>
#include <vector>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>

using namespace std;
typedef long long ll;
const int maxn = 5e4+10;

int prime[maxn],mob[maxn],vis[maxn],sum[maxn],cnt;

template <class T>
inline bool scan_d(T &ret)
{
    char c;
    int sgn;
    if (c = getchar(), c == EOF)
    {
        return 0; //EOF
    }
    while (c != '-' && (c < '0' || c > '9'))
    {
        c = getchar();
    }
    sgn = (c == '-') ? -1 : 1;
    ret = (c == '-') ? 0 : (c - '0');
    while (c = getchar(), c >= '0' && c <= '9')
    {
        ret = ret * 10 + (c - '0');
    }
    ret *= sgn;
    return 1;
}

void Mobius(){
    memset(prime,0,sizeof(prime));
    memset(mob,0,sizeof(mob));
    memset(vis,0,sizeof(vis));
    mob[1] = 1;
    cnt = 0;
    for(int i = 2;i < maxn; i++){
        if(!vis[i]){
            prime[cnt++] = i;
            mob[i] = -1;
        }
        for(int j = 0;j < cnt&&i*prime[j] < maxn;j++){
            vis[i*prime[j]] = 1;
            if(i%prime[j]) mob[i*prime[j]] = -mob[i];
            else{
                mob[i*prime[j]] = 0;
                break;
            }
        }
    }
    sum[0] = 0;
    for(int i = 1;i <= maxn;i++) sum[i] = sum[i-1]+mob[i];
}

ll cal(int n,int m){
    if(n > m) swap(n,m);
    ll ans = 0,j;
    for(int i = 1;i <= n;i = j+1){
        j = min(n/(n/i),m/(m/i));
        ans += (ll)(sum[j]-sum[i-1])*(n/i)*(m/i);
    }
    return ans;
}

int main(){
    Mobius();
    ll T,b,d,k;
    //scanf("%lld",&T);
    scan_d(T);
    for(ll Case = 1;Case <= T;Case++){
        //scanf("%lld%lld%lld",&b,&d,&k);
        scan_d(b);
        scan_d(d);
        scan_d(k);
        b /= k,d /= k;
        printf("%lld\n",cal(b,d));
    }
    return 0;
}
/*
2
4 5 2
6 4 3
*/