博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
莫比乌斯反演初步与实际应用
阅读量:5922 次
发布时间:2019-06-19

本文共 1326 字,大约阅读时间需要 4 分钟。

clipboard.png

clipboard.png

clipboard.png

clipboard.png

#include 
#include
#include
#include
using namespace std ;typedef long long LL ;const int MAXN = 10000010 ;long long T, N, M, V[MAXN], P[MAXN], Mu[MAXN], Tot ;long long S[MAXN], Ans, G[MAXN] ;inline long long Read() { long long X = 0, F = 1 ; char ch = getchar() ; while (ch > '9' || ch < '0') F = (ch == '-' ? - 1 : 1), ch = getchar() ; while (ch >= '0' && ch <= '9') X=(X<<1)+(X<<3)+(ch^48), ch = getchar() ; return X * F ;}inline void MU() { memset(V, 0, sizeof(V)) ; Mu[1] = 1 ; Tot = 0 ; for (int i = 2 ; i < MAXN ; i ++) { if (! V[i]) P[Tot ++] = i, Mu[i] = - 1, G[i] = 1 ; for (int j = 0 ; j < Tot && i * P[j] < MAXN ; j ++) { V[i * P[j]] = 1 ; if (i % P[j]) Mu[i * P[j]] = - Mu[i], G[i * P[j]] = Mu[i] - G[i] ; else { Mu[i * P[j]] = 0 ; G[i * P[j]] = Mu[i] ; break ; } } } for (int i = 1 ; i < MAXN ; i ++) S[i] = S[i - 1] + G[i] ;}int main() { MU() ; T = Read() ; while (T --) { Ans = 0 ; N = Read(), M = Read() ; for (int i = 1, j ; i <= min(N, M) ; i = j + 1) { j = min(N / (N / i), M / (M / i)) ; Ans += (N / i) * (M / i) * (S[j] - S[i - 1]) ; } printf("%lld\n", Ans) ; }}

转载地址:http://swsvx.baihongyu.com/

你可能感兴趣的文章
C++——this指针
查看>>
关于MVVM和RAC的一些拓展阅读
查看>>
day7
查看>>
python3 爬取百合网的女人们和男人们
查看>>
模式对象管理
查看>>
Java学习方法浅谈
查看>>
Linux 命令 "cp" 代码实现简介
查看>>
开启我的编程之路
查看>>
POJ-3186_Treats for the Cows
查看>>
bzoj1010: [HNOI2008]玩具装箱toy
查看>>
四叶草——计划会议
查看>>
UVa 1339,紫书P73,词频
查看>>
JavaScript-事件
查看>>
MarkDown换行
查看>>
【温故而知新】HTTP 报文
查看>>
部分关键字段重复的记录(保留最后记录)
查看>>
svn问题终极解决办法
查看>>
快速排序——Java
查看>>
Banana PI D1(Lamobo-D1)
查看>>
音频自动增益 与 静音检测 算法 附完整C代码【转】
查看>>