51nod1222 最小公倍数计数

最新推荐文章于 2018-12-09 12:04:45 发布

sdfzyhx

最新推荐文章于 2018-12-09 12:04:45 发布

阅读量418

点赞数

分类专栏：数学其他oj

本文链接：https://blog.csdn.net/sdfzyhx/article/details/73409737

版权

数学同时被 2 个专栏收录

229 篇文章 0 订阅

订阅专栏

其他oj

66 篇文章 0 订阅

订阅专栏

= = \sum i = 1 n \sum j = 1 n [lcm (i, j) \leq n] \sum d = 1 n \sum i = 1 ⌊ n d ⌋ \sum j = 1 ⌊ n i d ⌋ [g c d (i, j) = 1] \sum x = 1 n μ (x) \sum i \sum j \sum k [i j k \leq n x 2]

$\begin{align} & \sum_{i=1}^n\sum_{j=1}^n[\text{lcm}(i,j)\le n] \\ = & \sum_{d=1}^n\sum_{i=1}^{\lfloor\frac n d\rfloor} \sum_{j=1}^{\lfloor \frac n{id}\rfloor} [gcd(i,j)=1] \\ = & \sum_{x=1}^n\mu(x) \sum_i \sum_j \sum_k[ijk\le \frac n{x^2}] \\ \end{align}$
直接枚举

x $x$ ，这一步只需要枚举到

n‾‾√ $\sqrt n$ ，然后需要计算

f(n)=∑i∑j∑k[ijk≤n] $f(n)=\sum_i\sum_j\sum_k [ijk\le n]$ 。不妨设

i≤j≤k $i\le j\le k$ ，那么

i $i$ 只用枚举到

n‾‾√3 $\sqrt[3] n$ ，

j $j$ 只用枚举到

ni‾‾√ $\sqrt \frac{n}{i}$ 。
下面来推导一下复杂度。
首先看

f(n) $f(n)$ 的复杂度，

O(∫n√30nx‾‾√dx)=O(n23) $O(\int _0^{\sqrt[3]n}\sqrt{\frac n x} dx)=O(n^{\frac 2 3})$ 。
整个问题的复杂度

O(∫n√1(nx2)23dx)=O(n23⋅3(1−1n√6)) $O(\int_1^{\sqrt n}(\frac n{x^2})^{\frac 2 3} dx )=O(n^{\frac 2 3}\cdot 3(1-\frac1{\sqrt[6]n}))$
，右边一部分是常数，因此复杂度

O(n23) $O(n^{\frac 2 3})$ 。

#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
#define LL long long
const int maxn=400010,mod=10000007,lim=10000010;
int vis[maxn],mu[maxn],prm[maxn],
fir[mod+10],ne[lim],
tot;
LL a,b,mx;
/*struct set
{
    int fir[mod+10],ne[lim],tot;
    LL val[lim],ans[lim];
    LL find(LL n)
    {
        for (int i=fir[n%mod];i;i=ne[i])
            if (val[i]==n) return ans[i];
        return 0;
    }
    void ins(LL n,LL x)
    {
        if (tot>lim-5) return;
        int y=n%mod;
        ne[++tot]=fir[y];
        fir[y]=tot;
        val[tot]=n;
        ans[tot]=x;
    }
}h2,h3;
LL sum2(LL n)
{
    LL ret=h2.find(n);
    if (ret) return ret;
    for (LL i=1,j;i<=n;i=j+1)
    {
        j=n/(n/i);
        ret+=(j-i+1)*(n/i);
    }
    h2.ins(n,ret);
    return ret;
}
LL sum3(LL n)
{
    LL ret=h3.find(n);
    if (ret) return ret;
    ret=0;
    for (LL i=1,j;i<=n;i=j+1)
    {
        j=n/(n/i);
        ret+=(j-i+1)*sum2(n/i);
    }
    h3.ins(n,ret);
    return ret;
}*/
LL sum3(LL n)
{
    LL ret=0,x;
    for (LL i=1;i*i*i<=n;i++)
    {
        ret++;
        x=n/(i*i);
        if (x>i) ret+=(x-i)*3;
        for (LL j=i+1;i*j*j<=n;j++)
        {
            x=n/(i*j);
            ret+=(x-j)*6+3;
        }
    }
    return ret;
}
LL solve(LL n)
{
    LL ret=0;
    for (int i=1;(LL)i*i<=n;i++)
        if (mu[i]) ret+=mu[i]*sum3(n/i/i);
    return ret;
}
int main()
{
    //freopen("a.in","r",stdin);
    scanf("%lld%lld",&a,&b);
    mx=sqrt(b+0.5);
    mu[1]=1;
    for (int i=2;i<=mx;i++)
    {
        if (!vis[i])
        {
            prm[++tot]=i;
            mu[i]=-1;
        }
        for (int j=1;j<=tot&&(LL)i*prm[j]<=mx;j++)
        {
            vis[i*prm[j]]=1;
            if (i%prm[j]) mu[i*prm[j]]=-mu[i];
            else
            {
                mu[i*prm[j]]=0;
                break;
            }
        }
    }
    printf("%lld\n",(solve(b)-solve(a-1)+b-a+1)/2);
}