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

本文共 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://qeaso.baihongyu.com/

你可能感兴趣的文章
《嵌入式 Linux应用程序开发标准教程(第2版)》——2.2 Linux启动过程详解
查看>>
使用Karate轻松实现自动API测试
查看>>
React
查看>>
TopK 问题
查看>>
Linux-下载、安装和使用Mozilla Firefox浏览器及LumaQQ软件
查看>>
我的友情链接
查看>>
树莓派3音频配置
查看>>
SQLi filter evasion cheat sheet (MySQL)
查看>>
我的友情链接
查看>>
CentOS -bash: warning: setlocale: LC_MESSAGES: cannot change locale (en_US.UTF-8)
查看>>
bootstrapTable editable select empty
查看>>
spring boot 安全控制
查看>>
编写一个基本的Android应用程序
查看>>
我的友情链接
查看>>
【深入浅出Node.js系列五】初探Node.js的异步I/O实现
查看>>
我的友情链接
查看>>
exchange 2010 快速安装
查看>>
ZendStudio自定义代码补全,自定义代码段
查看>>
nginx+php-fpm内存不释放
查看>>
命名空间
查看>>