BZOJ 2820: YY的GCD
2820: YY的GCD
Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 1705 Solved: 908
[Submit][Status][Discuss]
Description
神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对kAc这种
傻×必然不会了,于是向你来请教……多组输入
Input
第一行一个整数T 表述数据组数接下来T行,每行两个正整数,表示N, M
Output
T行,每行一个整数表示第i组数据的结果
Sample Input
2
10 10
100 100
10 10
100 100
Sample Output
30
2791
2791
HINT
T = 10000
N, M <= 10000000
Source
分析:
我们很容易就可以得出一下这个复杂度爆炸的式子:
\sum_ {prime(p)}^{n} \sum_ {d=1}^{n/p} μ(d) \left\lfloor \frac{n}{pd} \right\rfloor \left\lfloor \frac{m}{pd} \right\rfloor
但是有10000个case,这个复杂度明显爆炸...我们设k=pd
可以得到一下式子:
\sum_ {k=1}^{n} \left\lfloor \frac{n}{k} \right\rfloor \left\lfloor \frac{m}{k} \right\rfloor \sum_ {prime(p) p\mid k} μ(\frac {k}{p})
然后对于后面的μ我们可以预处理前缀和,这样单个case复杂度大概可以降到sqrt(n)+sqrt(m)...
代码:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 | #include<algorithm> #include<iostream> #include<cstring> #include<cstdio> //by NeighThorn #define int long long using namespace std; const int maxn=10000000+5; int n,m,cas,cnt,f[maxn],mu[maxn],pri[maxn],vis[maxn]; inline void prework( void ){ mu[1]=1; for ( int i=2;i<=10000000;i++){ if (!vis[i]) pri[++cnt]=i,mu[i]=-1; for ( int j=1;j<=cnt&&i*pri[j]<=10000000;j++){ vis[i*pri[j]]=1; if (i%pri[j]==0){ mu[i*pri[j]]=0; break ; } mu[i*pri[j]]=-mu[i]; } } for ( int i=1;i<=cnt;i++) for ( int j=1;pri[i]*j<=10000000;j++) f[pri[i]*j]+=mu[j]; for ( int i=1;i<=10000000;i++) f[i]+=f[i-1]; } signed main( void ){ scanf ( "%lld" ,&cas);prework(); while (cas--){ scanf ( "%lld%lld" ,&n,&m); if (n>m) swap(n,m); int ans=0,r; for ( int i=1;i<=n;i=r+1){ r=min(n/(n/i),m/(m/i)); ans+=(f[r]-f[i-1])*(n/i)*(m/i); } printf ( "%lld\n" ,ans); } return 0; } |
By NeighThorn
分类:
数学
标签:
Dirichlet卷积
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 你所不知道的 C/C++ 宏知识
· 聊一聊 操作系统蓝屏 c0000102 的故障分析
· SQL Server 内存占用高分析
· .NET Core GC计划阶段(plan_phase)底层原理浅谈
· .NET开发智能桌面机器人:用.NET IoT库编写驱动控制两个屏幕
· 我干了两个月的大项目,开源了!
· 推荐一款非常好用的在线 SSH 管理工具
· 千万级的大表,如何做性能调优?
· 聊一聊 操作系统蓝屏 c0000102 的故障分析
· .NET周刊【1月第1期 2025-01-05】