又是一道卡常好题

坑掉我的 \(define \space int \space long \space long\) 感觉出题人并没有获得什么快乐……

Description

link

题意概述:

设 \(d(x)\) 为 \(x\) 的约数个数,求:

\[\sum_{i=1}^n \sum_{j=1}^n d(ij)
\]

多组询问,\(1\leq n,m,T \leq 5 \times10^4\)

Solution

首先给一个引理:

\[d(ij)=\sum_{x|i}\sum_{y|j} [gcd(x,y)=1]
\]

证明:参考\(@Siyuan\)的博客,好像可以用映射法?

\[Begin
\]

我们把式子转化一下:

\[\sum_{i=1}^n \sum_{j=1}^n \sum_{x|i}\sum_{y|j} [gcd(x,y)=1]
\]

(四个\(\sum\)好恐怖?)

然后改变求和的次序,首先枚举因数 \(x\) 和 \(y\) :

(这一步转化可以考虑用\(1\)到\(a\)中\(b\)的倍数个数来理解)

\[\sum_{x=1}^n \sum_{y=1}^n \lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y} \rfloor[gcd(x,y)=1]
\]

换一下,我们把这个 \(x\) 和 \(y\) 换成 \(i\) 和 \(j\)

\[\sum_{i=1}^n \sum_{j=1}^n \lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{j} \rfloor[gcd(i,j)=1]
\]

下一步开始反演,根据套路式,设:

\[f(x)=\sum_{i=1}^n \sum_{j=1}^n \lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{j} \rfloor[gcd(i,j)=1]
\]

\[g(x)=\sum_{x|d} f(d)
\]

所以得到关于 \(g(x)\) 的表达式:

\[g(x)=\sum_{i=1}^n \sum_{j=1}^n \lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{j} \rfloor[x|gcd(i,j)]
\]

把 \(x\) 提出来,同时去掉 \(gcd\)

\[g(x)=\sum_{i=1}^{\frac{n}{x}}\sum_{j=1}^{\frac{m}{x}} \lfloor \frac{n}{ix} \rfloor \lfloor \frac{m}{jx} \rfloor
\]

再根据\(f(x)\)的定义,得到答案为\(f(1)\)

又有:

\[f(x)=\sum_{n|d} \mu(\frac{d}{n}) g(d)
\]

所以有

\[Ans=\sum_{i=1}^n \mu(i)g(i)
\]

我们预处理\(s(x)=\sum^{x}_ {i=1} \lfloor \frac{x}{i}\rfloor\)然后就可以快速求解本题

\[Finished
\]

Code

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define reg register
namespace yspm{
inline int read()
{
int res=0,f=1; char k;
while(!isdigit(k=getchar())) if(k=='-') f=-1;
while(isdigit(k)) res=res*10+k-'0',k=getchar();
return res*f;
}
const int N=5e4+10;
int tot,mu[N],pri[N];ll s[N];
bool fl[N];
inline void prework()
{
mu[1]=1;
for(reg int i=2;i<N;++i)
{
if(!fl[i]) mu[i]=-1,pri[++tot]=i;
for(reg int j=1;j<=tot&&i*pri[j]<N;++j)
{
fl[i*pri[j]]=1;
if(i%pri[j]==0){mu[i*pri[j]]=0; break;}
else mu[i*pri[j]]=-mu[i];
}
} for(reg int i=1;i<N;++i) mu[i]+=mu[i-1];
for(reg int x=1;x<N;++x)
{
ll res=0;
for(reg int l=1,r;l<=x;l=r+1){r=x/(x/l); res+=1ll*(r-l+1)*(x/l);}
s[x]=res;
}
return ;
}
inline void work()
{
int n=read(),m=read(); if(n>m) swap(n,m);
ll ans=0; for(reg int l=1,r;l<=n;l=r+1){r=min(n/(n/l),m/(m/l)); ans+=1ll*(mu[r]-mu[l-1])*s[n/l]*s[m/l];}
return printf("%lld\n",ans),void();
}
signed main()
{
prework(); int T=read(); while(T--) work();
return 0;
}
}
signed main(){return yspm::main();}

LGOJ3327 【SDOI2015】约数个数和的更多相关文章

  1. BZOJ_3994_[SDOI2015]约数个数和_莫比乌斯反演

    BZOJ_3994_[SDOI2015]约数个数和_莫比乌斯反演 Description  设d(x)为x的约数个数,给定N.M,求   Input 输入文件包含多组测试数据. 第一行,一个整数T,表 ...

  2. P3327/bzoj3994 [SDOI2015]约数个数和(莫比乌斯反演)

    P3327 [SDOI2015]约数个数和 神犇题解(转) 无话可补 #include<iostream> #include<cstdio> #include<cstri ...

  3. 【BZOJ 3994】3994: [SDOI2015]约数个数和(莫比乌斯反演)

    3994: [SDOI2015]约数个数和 Description  设d(x)为x的约数个数,给定N.M,求   Input 输入文件包含多组测试数据. 第一行,一个整数T,表示测试数据的组数. 接 ...

  4. 洛谷 [SDOI2015]约数个数和 解题报告

    [SDOI2015]约数个数和 题目描述 设\(d(x)\)为\(x\)的约数个数,给定\(N,M\),求$ \sum\limits^N_{i=1}\sum\limits^M_{j=1}d(ij)$ ...

  5. BZOJ 3994: [SDOI2015]约数个数和

    3994: [SDOI2015]约数个数和 Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 898  Solved: 619[Submit][Statu ...

  6. 【BZOJ3994】[SDOI2015]约数个数和 莫比乌斯反演

    [BZOJ3994][SDOI2015]约数个数和 Description  设d(x)为x的约数个数,给定N.M,求   Input 输入文件包含多组测试数据. 第一行,一个整数T,表示测试数据的组 ...

  7. 洛谷P3327 - [SDOI2015]约数个数和

    Portal Description 共\(T(T\leq5\times10^4)\)组数据.给出\(n,m(n,m\leq5\times10^4)\),求\[\sum_{i=1}^n\sum_{j= ...

  8. P3327 [SDOI2015]约数个数和 莫比乌斯反演

    P3327 [SDOI2015]约数个数和 莫比乌斯反演 链接 luogu 思路 第一个式子我也不会,luogu有个证明,自己感悟吧. \[d(ij)=\sum\limits_{x|i}\sum\li ...

  9. [BZOI 3994] [SDOI2015]约数个数和(莫比乌斯反演+数论分块)

    [BZOI 3994] [SDOI2015]约数个数和 题面 设d(x)为x的约数个数,给定N.M,求\(\sum _{i=1}^n \sum_{i=1}^m d(i \times j)\) T组询问 ...

  10. 【BZOJ】3994: [SDOI2015]约数个数和

    题意: \(T(1 \le T \le 50000)\)次询问,每次给出\(n, m(1 \le n, m \le 50000)\),求\(\sum_{i=1}^{n} \sum_{j=1}^{m} ...

随机推荐

  1. jdk 的安装教程

    1. 配置的位置 鼠标右击计算机----属性-----高级系统设置---环境变量----在系统变量中配置(推荐) 2.增加环境变量 (1)新建 -----  变量名:   JAVA_HOME 变量值 ...

  2. outlook 2013邮件在服务器保留副本

    用outlook2013来收邮件确实是比较方便,但是它收邮件默认设置是:当outlook2013将在线邮箱的邮件下载至本机计算机之后,它就会删除在线邮箱中的邮件.不知道是不是以前邮箱容量比较小,所以要 ...

  3. nodejs(15)express开启cors跨域

    express开启cors跨域 package.json "dependencies": { "body-parser": "^1.18.3" ...

  4. VUE- 解决插件报错问题

    VUE- 解决插件报错问题 解决办法:先删除node_modules整个文件夹,然后cnpm cache clean,然后cnpm install. cnpm  cache clean时可能会提示 这 ...

  5. CF round #622 (div2)

    CF Round 622 div2 A.简单模拟 B.数学 题意: 某人A参加一个比赛,共n人参加,有两轮,给定这两轮的名次x,y,总排名记为两轮排名和x+y,此值越小名次越前,并且对于与A同分者而言 ...

  6. POJ 1979:Red and Black

    Red and Black Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 26058   Accepted: 14139 D ...

  7. Vue专题-路由系统

    一切分离都是为了更好的结合,本文详细介绍了前后端分离架构之后,前端如果实现路由控制,即Vue路由系统. Vue路由系统 VueRouter实现原理 VueRouter的实现原理是根据监控锚点值的改变, ...

  8. 使用java(jdbc)向mysql中添加数据时出现“unknown column……”错误

    错误情况如题,出现这个错误的原因是这样的: 在数据库中,插入一个字符串数据的时候是需要用单引号引起来的. 而下面的代码,注意看: sta.executeUpdate("INSERT INTO ...

  9. JXCPC 试题册

    JXCPC 试题册 Input file: standard input Output file: standard output Time limit: 1s Memory limit: 256 m ...

  10. SeetaFaceEngine系列2:Face Alignment编译和使用

    前面一篇写了编译人脸检测部分,现在就介绍下人脸配准部分,SeetaFace的Face Alignment通过人脸的五个关键点来配准人脸,也就是双眼.鼻尖.两个嘴角. 这部分的编译也和上一篇一样,步骤如 ...