51 nod 1222 最小公倍数计数

最新推荐文章于 2019-06-10 18:59:49 发布

原创最新推荐文章于 2019-06-10 18:59:49 发布 · 891 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

题解

数学题

莫比乌斯反演

本文介绍了一种算法，用于计算指定区间内所有可能的不同二元组数量，这些二元组的最小公倍数位于该区间。通过数学变换和枚举策略实现了高效的计算。

Description

定义F(n)表示最小公倍数为n的二元组的数量。
即：如果存在两个数（二元组）X，Y(X <= Y)，它们的最小公倍数为N，则F(n)的计数加1。
例如：F(6) = 5，因为[2,3] [1,6] [2,6] [3,6] [6,6]的最小公倍数等于6。

给出一个区间[a,b]，求最小公倍数在这个区间的不同二元组的数量。
例如：a = 4，b = 6。符合条件的二元组包括：
[1,4] [2,4] [4,4] [1,5] [5,5] [2,3] [1,6] [2,6] [3,6] [6,6]，共10组不同的组合。

Input

输入数据包括2个数：a, b，中间用空格分隔（1 <= a <= b <= 10^11)。

Output

输出最小公倍数在这个区间的不同二元组的数量。

Input示例

4 6

Output示例

Solution

忽略x<=y这个条件

a n s = \sum i = 1 n \sum j = 1 n [i j g c d ( i , j ) < = n]

$ans=\sum_{i=1}^n\sum_{j=1}^n [{ij\over gcd(i,j)}<=n]$
双sigma+gcd？直接考虑反演
为方便考虑，下取整忽略不写，即除法默认下取整
套路得

a n s = \sum d = 1 n \sum i = 1 n \sum j = 1 n [g c d (i, j) = 1] [i j d < = n]

$ans=\sum_{d=1}^n\sum_{i=1}^n\sum_{j=1}^n [gcd(i,j)=1][{ijd}<=n]$
设f[d]表示gcd(i,j)=d满足上面的条件的个数

f [d] = \sum i = 1 n d \sum j = 1 n d [i j d 2 < = n]

$f[d]=\sum_{i=1}^{n\over d}\sum_{j=1}^{n\over d}[ijd^2<=n]$
代回

a n s = \sum k = 1 n \sum d = 1 n μ (d) \sum i = 1 n d \sum j = 1 n d [k i j d 2 < = n]

$ans=\sum_{k=1}^n\sum_{d=1}^n\mu(d)\sum_{i=1}^{n\over d}\sum_{j=1}^{n\over d} [kijd^2<=n]$

a n s = \sum d = 1 n μ (d) \sum k = 1 n \sum i = 1 n d \sum j = 1 n d [k i j < = n d 2]

$ans=\sum_{d=1}^n\mu(d)\sum_{k=1}^n\sum_{i=1}^{n\over d}\sum_{j=1}^{n\over d} [kij<={n\over d^2}]$
即求

k∗i∗j<n $k*i*j<n$
先假设

k<i<j $k<i<j$ 那么结果*6
再假设

k=i<j和k<i=j $k=i<j 和 k<i=j$ 结果*3
最后再加上X<=Y这个限制即ans=(ans+n)/2

Code

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define N 501000
#define fo(i,a,b) for(ll i=a;i<=b;i++)
#define ll long long
using namespace std;
ll mu[N],p[N],bz[N],s[N];
ll calc(ll n)
{
    ll ans=0;
    for(ll d=1;(ll)d*d<=n;d++)
    if(mu[d]!=0)
    {
        ll t=(n/d/d),an=0;
        for(ll a=1;(ll)a*a*a<=t;a++)
        {
            for(ll b=a+1;(ll)a*b*b<=t;b++) 
            an+=((t/a/b)-b)*6+3;
            an++;
            an+=((t/a/a)-a)*3;
        }
        ans+=an*mu[d];
    }
    return (ans+n)/2;
}
int main()
{
    mu[1]=1;
    fo(i,2,N-1)
    {
        if (bz[i]==0) p[++p[0]]=i,mu[i]=-1;
        fo(j,1,p[0])
        {
            int x=i*p[j];
            if(x>N-1) break;
            bz[x]=1;
            if(i%p[j]==0) {mu[x]=0;break;}
            mu[x]=-mu[i];
        }
    }
    ll n,m;
    scanf("%lld %lld",&n,&m);
    ll ans=calc(m);
    ans-=calc(n-1);
    printf("%lld",ans);
}