【JZOJ 杂题选讲】51nod1222

题目

1<=a,b<=10¹⁰,1<=m<=10⁶

思路

设f(x)表示lcm=x的方案数
考虑f(x)等于什么
设x=∏ai^pi，则f(x)=∏((pi+1)^m-pi^m)
意思是lcm=x的话m个数中的质因子ai系数最高且至少有一个是pi
容斥一下求积
①f(x)是积性函数
可以分开求每个质因子的贡献
②f§可以表示成多项式
f§=2^m-1
③f(p^c)可快速计算
f(p^c)=(c+1)^m-c^m
min25即可

代码

#include<bits/stdc++.h>
using namespace std;

const int N=316227;

typedef long long ll;

int tot,p[N],mu[N+5];

ll l,r,ans;

bool bz[N+5];

char c;

void init()
{
    mu[1]=1;
    for (int i=2;i<=N;i++)
    {
        if (!bz[i])
        {
            mu[i]=-1; p[tot++]=i;
        }
        for (int j=0;j<tot;j++)
        {
            int I=i*p[j];
            if (I>N) break;
            bz[I]=1;
            if (i%p[j]==0)
            {
                mu[I]=0; break;
            }
            mu[I]=-mu[i];
        }
    }
}

ll calc(ll n)
{
    if (!n) return 0;
    ans=0;
    for (int k=1;(ll)k*k<=n;k++) if (mu[k]!=0)
    {
        ll m=n/((ll)k*k),s=0;
        for (int i=1;(ll)i*i*i<=m;i++)
        {
            for (int j=i+1;(ll)i*j*j<=m;j++)
            {
                s+=(m/((ll)i*j)-j)*6+3;
            }
            s+=(m/((ll)i*i)-i)*3+1;
        }
        if (mu[k]==1) ans+=s;else ans-=s;
    }
    return (ans+n)/2;
}

int main()
{
    init();
    scanf("%lld%lld",&l,&r);
    printf("%lld\n",calc(r)-calc(l-1));
    return 0;
}