点此看题面

大致题意: 求\(\sum_{x=1}^N\sum_{y=1}^M[gcd(x,y)==d]\)。

一道类似的题目

推荐先去做一下这道题:【洛谷2257】YY的GCD,来初步了解一下莫比乌斯反演

再来看这题,就非常简单了。

一些定义

按照上面提到的那题的思路,首先,我们可以定义\(f(d)\)和\(F(d)\)如下:

\[f(d)=\sum_{i=1}^N\sum_{j=1}^M[gcd(i,j)==d]
\]

\[F(d)=\sum_{i=1}^N\sum_{j=1}^M[d|gcd(i,j)]
\]

通过定义,不难发现:

\[F(n)=\sum_{n|d}f(d)=\lfloor\frac Nn\rfloor\lfloor\frac Mn\rfloor
\]

由于莫比乌斯反演的某些性质,我们又可以得到:

\[f(n)=\sum_{n|d}\mu(\lfloor\frac dn\rfloor)F(d)
\]

公式化简

我们应该不难想到:

\[answer=f(d)
\]

貌似比 YY的GCD 简单了许多。

然后就是一波化简。

先来一波莫比乌斯反演

\[answer=\sum_{p|d}\mu(\lfloor\frac dp\rfloor)F(d)
\]

但是,这样有点难以处理。

于是,我们改成枚举\(\lfloor\frac dp\rfloor\),于是原式就变成了这样:

\[answer=\sum_{p=1}^{min(\lfloor\frac Nd\rfloor,\lfloor\frac Md\rfloor)}\mu(p)\lfloor\frac N{d·p}\rfloor\frac M{d·p}\rfloor
\]

这样就很容易求解了。

求解答案

不难想到,我们可以用除法分块

不难发现,在一定范围内\(\lfloor\frac Ni\rfloor\)的值是保持不变的(\(\lfloor\frac Mi\rfloor\)同理),则\(\lfloor\frac NG\rfloor\lfloor\frac MG\rfloor\)其实最多只有\(\sqrt N+\sqrt M\),而对于\(\lfloor\frac NG\rfloor\lfloor\frac MG\rfloor\)相同的值,我们可以一起求解,于是就能想到用\(sum_i\)来表示\(\sum_{i=1}^i \mu(i)\),这样就能快速求解了。

代码

#include<bits/stdc++.h>
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)<(y)?(x):(y))
#define uint unsigned int
#define LL long long
#define ull unsigned long long
#define swap(x,y) (x^=y,y^=x,x^=y)
#define abs(x) ((x)<0?-(x):(x))
#define INF 1e9
#define Inc(x,y) ((x+=(y))>=MOD&&(x-=MOD))
#define ten(x) (((x)<<3)+((x)<<1))
#define N 50000
using namespace std;
int n,m,k;
class FIO
{
private:
#define Fsize 100000
#define tc() (FinNow==FinEnd&&(FinEnd=(FinNow=Fin)+fread(Fin,1,Fsize,stdin),FinNow==FinEnd)?EOF:*FinNow++)
#define pc(ch) (FoutSize<Fsize?Fout[FoutSize++]=ch:(fwrite(Fout,1,FoutSize,stdout),Fout[(FoutSize=0)++]=ch))
int f,FoutSize,OutputTop;char ch,Fin[Fsize],*FinNow,*FinEnd,Fout[Fsize],OutputStack[Fsize];
public:
FIO() {FinNow=FinEnd=Fin;}
inline void read(int &x) {x=0,f=1;while(!isdigit(ch=tc())) f=ch^'-'?1:-1;while(x=ten(x)+(ch&15),isdigit(ch=tc()));x*=f;}
inline void read_char(char &x) {while(isspace(x=tc()));}
inline void read_string(string &x) {x="";while(isspace(ch=tc()));while(x+=ch,!isspace(ch=tc())) if(!~ch) return;}
inline void write(LL x) {if(!x) return (void)pc('0');if(x<0) pc('-'),x=-x;while(x) OutputStack[++OutputTop]=x%10+48,x/=10;while(OutputTop) pc(OutputStack[OutputTop]),--OutputTop;}
inline void write_char(char x) {pc(x);}
inline void write_string(string x) {register int i,len=x.length();for(i=0;i<len;++i) pc(x[i]);}
inline void end() {fwrite(Fout,1,FoutSize,stdout);}
}F;
class Class_Mobius//莫比乌斯反演
{
private:
int Prime_cnt,mu[N+5],Prime[N+5];bool IsNotPrime[N+5];
public:
LL sum[N+5];
Class_Mobius()//预处理
{
register int i,j;
for(mu[1]=1,i=2;i<=N;++i)//求出莫比乌斯函数
{
if(!IsNotPrime[i]) Prime[++Prime_cnt]=i,mu[i]=-1;
for(j=1;j<=Prime_cnt&&i*Prime[j]<=N;++j)
if(IsNotPrime[i*Prime[j]]=true,i%Prime[j]) mu[i*Prime[j]]=-mu[i];else break;
}
for(i=1;i<=N;++i) sum[i]=sum[i-1]+mu[i];//求前缀和
}
}Mobius;
int main()
{
register int i,nxt,lim,T;register LL ans;F.read(T);
while(T--)
{
F.read(n),F.read(m),F.read(k),lim=min(n,m)/k;
for(ans=0,i=1;i<=lim;i=nxt+1) nxt=min(n/(n/i),m/(m/i)),ans+=1LL*(n/(i*k))*(m/(i*k))*(Mobius.sum[nxt]-Mobius.sum[i-1]);//除法分块
F.write(ans),F.write_char('\n');//输出答案
}
return F.end(),0;
}

【BZOJ1101】[POI2007] Zap(莫比乌斯反演)的更多相关文章

  1. BZOJ1101: [POI2007]Zap(莫比乌斯反演)

    1101: [POI2007]Zap Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 2951  Solved: 1293[Submit][Status ...

  2. Bzoj1101: [POI2007]Zap 莫比乌斯反演+整除分块

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1101 莫比乌斯反演 1101: [POI2007]Zap 设 \(f(i)\) 表示 \(( ...

  3. 1101: [POI2007]Zap(莫比乌斯反演)

    1101: [POI2007]Zap Time Limit: 10 Sec Memory Limit: 162 MB Description FGD正在破解一段密码,他需要回答很多类似的问题:对于给定 ...

  4. 【BZOJ1101】Zap [莫比乌斯反演]

    Zap Time Limit: 10 Sec  Memory Limit: 162 MB[Submit][Status][Discuss] Description 对于给定的整数a,b和d,有多少正整 ...

  5. BZOJ 1101: [POI2007]Zap( 莫比乌斯反演 )

    求 answer = ∑ [gcd(x, y) = d] (1 <= x <= a, 1 <= y <= b) . 令a' = a / d, b' = b / d, 化简一下得 ...

  6. BZOJ1101 POI2007 Zap 【莫比乌斯反演】

    BZOJ1101 POI2007 Zap Description FGD正在破解一段密码,他需要回答很多类似的问题:对于给定的整数a,b和d,有多少正整数对x,y,满足x<=a,y<=b, ...

  7. 【题解】Zap(莫比乌斯反演)

    [题解]Zap(莫比乌斯反演) 裸题... 直接化吧 [P3455 POI2007]ZAP-Queries 所有除法默认向下取整 \[ \Sigma_{i=1}^x\Sigma_{j=1}^y[(i, ...

  8. [BZOJ1101][POI2007]Zap

    [BZOJ1101][POI2007]Zap 试题描述 FGD正在破解一段密码,他需要回答很多类似的问题:对于给定的整数a,b和d,有多少正整数对x,y,满足x<=a,y<=b,并且gcd ...

  9. 【莫比乌斯反演】BZOJ1101 [POI2007]zap

    Description 回答T组询问,有多少组gcd(x,y)=d,x<=a, y<=b.T, a, b<=4e5. Solution 显然对于gcd=d的,应该把a/d b/d,然 ...

  10. P3455 [POI2007]ZAP-Queries(莫比乌斯反演)

    题目 P3455 [POI2007]ZAP-Queries 解析 莫比乌斯反演. 给定\(n\),\(m\),\(d\),求\[\sum_{i=1}^{n}\sum_{j=1}^{m}[gcd(i,j ...

随机推荐

  1. 解读人:李琼,Metabolic profiling of tumors, sera and skeletal muscles from an orthotopic murine model of gastric cancer associated-cachexia(胃癌相关恶病质的原位小鼠模型中肿瘤,血清和骨骼肌的代谢谱分析)

    发表时间:(2019年4月) IF:3.950 单位: 厦门大学 厦门理工大学 物种:小鼠的肿瘤,血清和骨骼肌 技术:核磁共振代谢组学分析 一. 概述:(用精炼的语言描述文章的整体思路及结果) 本研究 ...

  2. 《权限系列shiro+cas》---封装公共验证模块

    操作系统:Windows8.1 显卡:Nivida GTX965M 开发工具:Visual Studio 2017 Introduction 我们现在可以将任意属性传递给每个顶点的顶点着色器使用.但是 ...

  3. LightOJ 1336 - Sigma Function

    原题链接 基础数论中很经典的一道题 题意 给出了σ(n)的计算公式,让你找出整数1-n中有多少对应σ(n)的值是偶数. 思路 观察σ(n)的公式发现,每一个乘项都是 (piei+1 - 1) / (p ...

  4. 关于php命名空间的理解

    以phpmailer这个类库为例,composer自动加载好该类库,有用的文件都放在src这个目录下 这些文件的命名空间都是这个:namespace PHPMailer\PHPMailer; 如果我们 ...

  5. git教程3-分支

    https://git-scm.com/book/zh/v1/Git-%E5%88%86%E6%94%AF-%E5%88%86%E6%94%AF%E7%9A%84%E7%AE%A1%E7%90%86 ...

  6. day04 ---Linux安装Python3

    如何linux上安装python3 1.下载源代码,方式有2个, 1.在windows上下载,下载完成后,通过lrzsz工具,或者xftp工具,传输到linux服务器中 2.在linux中直接下载 c ...

  7. js 数组,字符串,json互相转换(在select实现多个输入的时候与后台交互常使用)

    数组转字符串 var arr = [1,2,3,4,'巴德','merge']; var str = arr.join(','); console.log(str); // 1,2,3,4,巴德,me ...

  8. Luogu P3223 [HNOI2012]排队 组合

    本来做了一道  P4901 排队 后来发现自己做错题了...到也都是数学qwq 这题最恶心的就是两只(雾)老师. 那我们分类讨论: 1.两个老师之间是男生: $ A(n,n)*A(n+1,2)*A(n ...

  9. 一些有关PyCharm使用总结

    目前在这里,你能看见 license server Python版本配置 添加另外版本的Python 设置字体大小 关于编码 关于模版 安装好之后,第一个问题就是 license server 问题, ...

  10. FTP服务器FileZilla Server配置及使用方法

    FileZilla Server下载安装完成后,安装过程不写说明了,网上一抓一大把,直接从配置开始记录. 1.创建服务器 ²  Password:栏位中输入本服务器Filezilla服务的密码, ²  ...