1742 开心的小Q

如果一个数字存在一个约数是完全平方数，那么小Q就认为这个数是有趣的。

小Q喜欢收集有趣的数字，每找到一个有趣的数，小Q就会变得很开心。

小Q发现12是有趣的，18也是有趣的，它们都是36的约数，而在36的约数中，还有3个数是有趣的，它们是4、9、36。

小Q很好奇，在a~b里每个数字各有多少个有趣的约数，由于a和b太大了，所以他只想知道这些个数之和是多少。

例如4有1个有趣的约数，8有2个有趣的约数，9有1个有趣的约数，所以1~10里每个数的有趣约数个数之和是4。

Input

输入数据包括2个数：a, b，中间用空格分隔。（1≤a≤b≤10^9）

Output

输出a~b里每个数字的有趣约数个数之和。

Input示例

1 10

Output示例

4
思路：分块+莫比乌斯;
根据莫比乌斯abs(mul[i]) = 0表示i这个数有平方项因子，abs（mul[i]） = 1表示这个数没有平方项因子。

，那么答案就是S（m）-S(n-1);

然后分块计算，前sqrt（n）直接用公式，然后后面的因为（n/i）在一段区间内是相同的值，那么我们只要算出这个区间的长度，这个容斥下就行了，算出区间内符合是含有平方因子的数。

#include<bits/stdc++.h>
typedef long long LL;
using namespace std;
bool prime[100005];
LL ak[100005];
LL mul[100005];
const int BufferSize=1<<16;
char buffer[BufferSize],*head,*tail;
inline char Getchar() {
    if(head==tail) {
        int l=fread(buffer,1,BufferSize,stdin);
        tail=(head=buffer)+l;
    }
    return *head++;
}
inline int read() {
    int x=0,f=1;char c=Getchar();
    for(;!isdigit(c);c=Getchar()) if(c=='-') f=-1;
    for(;isdigit(c);c=Getchar()) x=x*10+c-'0';
    return x*f;
}
LL slove(LL n);
int akk[100005];
int main(void)
{
    int cn = 0;
    mul[1] = 1;int ap = 0;
    for(int i = 2; i <= 100000; i++)
    {
        if(!prime[i])
        {
            ak[cn++] = i;
            mul[i] = -1;
        }
        for(int j = 0; j < cn&&(LL)ak[j]*i<=100000; j++)
        {
            if(i%ak[j])
            {
                prime[i*ak[j]] = true;
                mul[i*ak[j]] = -mul[i];
            }
            else
            {
                prime[i*ak[j]] = true;
                mul[i*ak[j]] = 0;
                break;
            }
        }
        if(mul[i])akk[ap++] = i;
    }
    LL n,m;
   n = read(),m=read();
    LL c = sqrt(n); LL sum = 0;
    sum = slove(m)-slove(n-1);
    printf("%lld\n",sum);
    return 0;
}
LL slove(LL n)
{   if(n == 0)return 0;
    LL c = sqrt(n);LL sum = 0;
    for(int i = 2;i <= c;i++)
    {
        if(!mul[i])
        {
            sum += n/i;
        }
    }
    for(int i = 1;i <= n/c-1;i++)
    {
        LL a,b;
        a = n/(i+1);
        b = n/i;
        a++;LL v = 0;
        for(int j = 0;akk[j] <= sqrt(b);j++)
        {
            sum -= (mul[akk[j]]*((b)/(akk[j]*akk[j])-(a-1)/(akk[j]*akk[j])))*i;
        }
    }
    return sum;
}