bzoj2693 jzptab

最新推荐文章于 2019-07-28 10:39:00 发布

Euryale_

最新推荐文章于 2019-07-28 10:39:00 发布

阅读量500

点赞数

分类专栏：数论文章标签：莫比乌斯反演分块

本文链接：https://blog.csdn.net/u010336344/article/details/67631485

版权

数论专栏收录该内容

25 篇文章 0 订阅

订阅专栏

题目链接：bzoj2693
题目大意：
多组数据。给定n、m，求 $\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i*j)$ 。
T <= 10000 N, M<=10000000

题解：
莫比乌斯反演
这个n,m的范围是一样的，但是是多组数据。可以说是bzoj2154的加强版。
那么我们根据上次的式子继续化，上次化到的式子是

\sum i = 1 n \sum j = 1 m l c m (i * j) = \sum d = 1 m i n (n, m) d \sum t = 1 m i n (⌊ n d ⌋, ⌊ m d ⌋) t 2 * μ (t) * s u m (⌊ n d t ⌋, ⌊ m d t ⌋)

$\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i*j)=\sum_{d=1}^{min(n,m)}d\sum_{t=1}^{min(\lfloor\frac{n}{d}\rfloor,\lfloor\frac{m}{d}\rfloor)}t^2*\mu(t)*sum({\lfloor\frac{n}{dt}\rfloor},{\lfloor\frac{m}{dt}\rfloor})$
改变求和指标。

\sum i = 1 n \sum j = 1 m l c m (i * j) = \sum D = 1 m i n (n, m) s u m (⌊ n D ⌋, ⌊ m D ⌋) * D \sum t | D t * μ (t)

$\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i*j)=\sum_{D=1}^{min(n,m)}sum(\lfloor\frac{n}{D}\rfloor,\lfloor\frac{m}{D}\rfloor)*D\sum_{t|D}t*\mu(t)$
设

f(x)=∑i|xi∗μ(i) $f(x)=\sum_{i|x}i*\mu(i)$ ，

g(i)=i∗f(i) $g(i)=i*f(i)$
那么式子就变成了

\sum i = 1 n \sum j = 1 m l c m (i * j) = \sum D = 1 m i n (n, m) s u m (⌊ n D ⌋, ⌊ m D ⌋) * g (D)

$\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i*j)=\sum_{D=1}^{min(n,m)}sum(\lfloor\frac{n}{D}\rfloor,\lfloor\frac{m}{D}\rfloor)*g(D)$
对于前面的看到

nD $\frac{n}{D}$ 这种东西肯定套路分块啦。然后再看后面的

g(i) $g(i)$
观察

f $f$ ，可以发现它是个积性函数。
即对于互质的两个数就有

f[i∗prime[j]]=f[i]∗f[prime[j]] $f[i*prime[j]]=f[i]*f[prime[j]]$ ，于是我们就会想把它线性筛出来。但是对于不互质的两个数怎么办呢？注意到

f $f$ 里面有个

μ $\mu$ ，就是指数有大于1的话

μ $\mu$ 是0。所以如果那两个数不互质，其乘积一定有指数大于等于2的，所以那部分是不作贡献的，即

f[i∗prime[j]]=f[i] $f[i*prime[j]]=f[i]$ 。
然后再for一遍就能处理出

g <script type="math/tex" id="MathJax-Element-60">g</script>了。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
#define N 10000100

const LL ny=50000005;
const LL mod=100000009;
int mymin(int x,int y){return (x<y)?x:y;}
bool ispri[N];int cnt,pri[N/5];LL D[N],f[N];
void pre(int lim)
{
    cnt=0;f[1]=1;
    for (int i=2;i<=lim;i++)
    {
        if (!ispri[i]) {pri[++cnt]=i;f[i]=(1-i+mod)%mod;}
        for (int j=1;j<=cnt && i*pri[j]<=lim;j++)
        {
            int k=i*pri[j];
            ispri[k]=true;
            if (i%pri[j]==0)
            {
                f[k]=f[i];
                break;
            }f[k]=(f[i]*f[pri[j]])%mod;
        }
    }
    D[1]=f[1];for (int i=2;i<=lim;i++) D[i]=(f[i]*i%mod+D[i-1])%mod;
}
LL sum(int n,int m)
{
    LL s1=(LL)(n)*(LL)(n+1)%mod*ny%mod;
    LL s2=(LL)(m)*(LL)(m+1)%mod*ny%mod;
    return s1*s2%mod;
}
int main()
{
    //freopen("a.in","r",stdin);
    //freopen("a.out","w",stdout);
    int T,n,m,i,lim,r;LL ans;
    pre(10000000);scanf("%d",&T);
    while (T--)
    {
        scanf("%d%d",&n,&m);
        lim=mymin(n,m);ans=0;
        for (i=1;i<=lim;i=r+1)
        {
            r=mymin(n/(n/i),m/(m/i));
            ans+=(D[r]-D[i-1]+mod)%mod*sum(n/i,m/i)%mod;
            ans%=mod;
        }while (ans<0) ans+=mod;
        printf("%lld\n",ans);
    }
    return 0;
}