GCD(莫比乌斯+去重）

题目链接

莫比乌斯反演模板题， 去重即可；

我们可以发现只有在区间重叠部分才会有重复且为cal(e, e, k)/2;（e表示b, d中较小的一个）；

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<string>
#include<vector>
#include<stack>
#include<bitset>
#include<cstdlib>
#include<cmath>
#include<set>
#include<list>
#include<deque>
#include<map>
#include<queue>
using namespace std;
typedef long long ll;
const int maxn = ;

int a,b,c,d,k;
ll pri[maxn],vis[maxn],mu[maxn];

void init()
{

    mu[] = ;
    int cnt  = ;
    for(int i = ; i <= maxn; i++)
    {
        if(vis[i] == ) {
            pri[cnt++] = i;
            mu[i] = -;
        }
        for(int j = ; j < cnt&&i*pri[j] <= maxn; j++)
        {
            ll k = i*pri[j];
            vis[k] = ;
            if(i%pri[j] == ) {mu[k]= ; break;}
            else mu[k] = -mu[i];
        }
    }
    for(int i = ; i <= maxn; i++)
        mu[i] += mu[i-];
}

ll cal(int l, int r, int k)  // 神奇的分块
{
    if(!l||!r||!k) return ;
    if(l > r) swap(l, r);
    l /= k, r /= k;
    int last = ;
    ll ans  = ;
    for(int i = ; i <= l; i = last+)
    {
        last = min( l/(l/i), r/(r/i) );
        ans += ( mu[last] - mu[i-] )*(l/i)*(r/i);
    }
    return ans;
}

int main()
{
    ios::sync_with_stdio(false);
    init();
    int n;
    cin >> n;
    int i = ;
    while(i <= n)
    {
        cin >> a >> b >> c >> d >> k;
        int e  = min(b,d);
        ll ans = cal(b, d, k) - cal(e, e, k)/;
      //  printf("Case %d: %lld\n", i, ans);
       cout << "Case " << i  << ": " << ans << endl;
        i++;
    }

    return ;
}