莫比乌斯反演的公式_莫比乌斯反演公式-程序员宅基地

技术标签: 算法的数学基础  

莫比乌斯反演

由于莫比乌斯反演的应用非常广泛,内容很多但是结论却并不复杂。
一道经典的莫比乌斯反演题:
求:∑ni=1∑mj=1[gcd(i,j)==d]∑i=1n∑j=1m[gcd(i,j)==d]
也就是说有多少对(i,j)的gcd为d。

 莫比乌斯反演公式

莫比乌斯反演公式

 莫比乌斯函数

莫比乌斯函数

程序模板

void mobius()
{
    int i,j; mbs[1] = 1;
    fo(i,2,N)
        {
            if (!vis[i]) {p[++p[0]] = i; mbs[i] = -1;}
            for (j = 1;j <= p[0] && i * p[j] <= N; j++)
                {
                    vis[i*p[j]] = 1;
                    if (i % p[j] == 0) {mbs[i*p[j]] = 0; break;}
                    mbs[i*p[j]] = - mbs[i];
                }
        }
}

构造

这里写图片描述

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/u011590573/article/details/81416298

智能推荐

Spark Streaming:原理介绍 & 实时流计算 & perfect-程序员宅基地

目录1、Spark Streaming简介1.1 概述1.2 术语定义1.3 Storm与Spark Streming比较2、运行原理2.1 Streaming架构2.2 编程模型2.2.1 如何使用Spark Streaming2.2.2 DStream的输入源2.2.3 DStream的操作2.3 容错、持久化和性能调优2.3.1 容错2....

利用griddata进行插值-程序员宅基地

因为最近在做算法优化,所以对数据统一性有一定要求,在最近的研究中主要用一个简单的最近邻插值对数据集进行降尺度处理。主要运用到的函数时scipy里面的griddata第一步:导入相关库import xarray as xr from scipy.interpolate import griddata # 插值函数import numpy as np第二步:给出插值到的经纬度信息(目标经纬度)mask_tmp = xr.open_dataset('G:/China/temperatu_griddata

熟悉常用的Linux操作 cd ls mkdir rmdir cp mv rm tac cat more head tail touch chown touch tar grep-程序员宅基地

Linux 最基本的命令及其补充的知识点 cd ls mkdir rmdir cp mv rm tac cat more head tail touch chown touch tar grep

杰里之蓝牙模式提示语未播放完被连接提示语打断解决方法【篇】_蓝牙 提示语 格式-程序员宅基地

需要在提示音播完之后在回调函数中调用初始化蓝牙函数既可。_蓝牙 提示语 格式

webView 的获取页面title的_react-native-webview document.title获取不到-程序员宅基地

获取h5页面的携带的title中是很简单的 ,昨天脑子一糊涂就研究了一下午都不知道怎么获取。加载Android加载h5中有两个方法第一个 webview.setWebChromeClient 第二个是webview.setwebviewClient.setWebChromeClient辅助WebView处理JavaScript的对话框,网站图标,网站title,加载进度等;..._react-native-webview document.title获取不到

【Codeforces 733D】 Kostya the Sculptor【贪心】-程序员宅基地

题意:给你n个长方形,让你找出2个或1个长方体,使得他们拼接成的长方体的内接圆半径最大(两个矩形拼接的条件是他们有一个面完全相同)题解:很容易想到一个长方体的内接圆半径是由他的最短的那条边决定的,设长方体的边长从大到小分别为a,b,c,只有用a*b这个面和别的长方体拼接,才有可能使得拼接后的长方体内接圆半径更大于是,根据上述我们读入的时候将长方体边长排序然后sort 以a为第一关键字 以b为第二关...

随便推点

CentOS8 安装图形界面后,怎么转换为命令行界面_centos8切换命令行界面_优小U的博客-程序员宅基地

问题安装了CentOS8图形界面版后,每次系统启动时间太长,而且如果VMWare里面安装了多台虚拟机,图形界面的CentOS耗费的内存太多使用Linux 主机不使用命令行就有点不是内味了~命令行形式启动:提高启动速度、降低内存占用命令行与图形界面转换图形界面->命令行第一步,在图形界面任何界面按ctrl+alt+f3,这时就会切换到命令行界面;第二步,在命令行里输入systemctl set-default multi-user.target,这样下次启动就会默认以命令行启动了。_centos8切换命令行界面

ANSYS 有限元分析 坐标系/工作平面_kwpave-程序员宅基地

目 录一、前沿二、总体及局部坐标系2.1 激活坐标系 CSYS2.2 新建坐标系2.2.1 LOCAL2.2.2 CLOCAL2.2.3 坐标系类命令三、结点坐标系 NROTAT四、单元坐标系五、结果坐标系 RSYS六、工作平面 WPCSYS七、尾声八、参考文献阿阳的 ANSYS 使用指南,本文仅用于个人学习,除此之外,无其他任何用途。因个人能力有限,本文难免有所疏漏/错误,不妥之处还请各位批评指正。Ansys 参数化编程与命令 / APDL一、前沿  模型空间背景色为黑色_kwpave

庞果网在线编程---倒水---扩展欧几里得算法-程序员宅基地

水题一道:#include #include #include #include using namespace std;int gcd(int a,int b){ if(a

WPF 基础知识学习简单总结(二)_wpf各个功能模块的上下文对象-程序员宅基地

行为:封装通用用户界面功能获取行为支持:安装Expression Blend System.Windows.Interativity.dll。定义支持行为的基本类。 Microsoft.Expression.Inteactions.dll。通过添加可选的以核心行为类为基础的动作和触发器类创建行为:创建继承自Behavior基类的类 覆盖OnAttached()和OnDetaching()。 可以访问放置行为的元素(通过AssociatedObject属性), 调用OnAttac..._wpf各个功能模块的上下文对象

React: 傻瓜/展示性组件的简化 (函数型组件及其传参)_无状态组件传参-程序员宅基地

React的组件可以分为“聪明组件&傻瓜组件”或“容器组件&展示组件”。这里展示对后者进行简化的案例。_无状态组件传参

Bootstrap碎语-程序员宅基地

这里记录下某段时间Bootstrap的零散碎片。 1、有关Bootstrap的参考网站: ● 官方:http://getbootstrap.com/● 主题:http://bootswatch.com/● Font-Awsome: http://fortawesome.github.io/Font-Awesome/● 幻灯片:lokeshdhakar.com/projects/ligh..._"$(document).delegate('*[data-toggle=\"lightbox\"]"