bzoj4659 Lcm（反演）

Description

给出A，B，考虑所有满足l<=a<=A，l<=b<=B，且不存在n>1使得n^2同时整除a和b的有序数对(a，b)，求其lcm(a，b)之和。答案模2^30。

Input

第一行一个整数T表示数据组数。接下来T行每行两个整数A，B表示一组数据。

T ≤ 2000,A,B ≤ 4 × 10^6

Output

对每组数据输出一行一个整数表示答案模2^30的值

Sample Input

5

2 2

4 6

3 4

5 1

23333 33333

Sample Output

7

148

48

15

451085813

分析：
首先我们要把题面转化成数学语言：
重点就是这句话：不存在n>1使得n^2同时整除a和b

同时整除是一个什么概念呢：a和b都有n^2这个因子
实际上就是gcd(a，b)有n^2这个因子
这样题面就可以转化成：

不存在n>1使得n^2整除gcd(a,b)

也就是说gcd(a,b)没有平方因子，

应该能敏感的反应出：mu值不等于0

设p为满足条件的gcd值（mu(p)!=0）

我们能得到式子：

又是一个关于LCM的式子，实际上化简的过程大同小异（最后得到的形式也很相像）

（实际上就是第一维枚举的元素限制条件变了：mu值不等于0）

但是这样的二维枚举并不能满足出题人
我们还需要进一步优化时间：
（还是bzoj2820的奇技淫巧。。。）

注意：p是mu值不等于0的数

如果我们能预处理出
整个式子就可在sqrt(n)的时间内求出

因为p的特殊限制，我们不能直接线筛，也不能分块
所以我们采取最笨的办法：
枚举所有符合条件的p，更新p的倍数
（均摊O(nlogn)）

还需要说一句：

因为这道题的模数比较特别，是2的整数次幂，所以可以直接依靠int的自然溢出（这样会快一点）
如果直接暴力%的话就T掉了

tip

题面和bzoj2694一样，但是数据好像更刁钻。。。

//这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>

using namespace std;

const int p=1<<30;
const int N=4000010;
int n,m,sshu[N],tot=0,mu[N],h[N];
bool no[N];

void make()
{
    mu[1]=1;
    for (int i=2;i<N;i++)
    {
        if (!no[i])
        {
            sshu[++tot]=i;
            mu[i]=-1;
        }
        for (int j=1;j<=tot&&sshu[j]*i<N;j++)
        {
            no[sshu[j]*i]=1;
            if (i%sshu[j]==0)
            {
                mu[i*sshu[j]]=0;
                break;
            }
            mu[i*sshu[j]]=-mu[i];
        }
    }

    for (int i=1;i<N;i++) if (mu[i]!=0)
        for (int t=1;t*i<N;t++)
        {
            int j=t*i;
            h[j]=h[j]+i*t*t*mu[t];
        }
    for (int i=2;i<N;i++) h[i]=h[i]+h[i-1];
}

int sum(int n)
{
    return n*(n+1)/2;
}

int main()
{
    make();
    int T;
    scanf("%d",&T);
    while (T--)
    {
        scanf("%d%d",&n,&m);
        int ans=0;
        int last;
        for (int i=1;i<=min(n,m);i=last+1)
        {
            last=min(n/(n/i),m/(m/i));
            ans=ans+sum(n/i)*sum(m/i)*(h[last]-h[i-1]);
        }
        printf("%d\n",(ans%p+p)%p);
    }
    return 0;
}