前言

在算法中,经常需要用到一种与调和级数有关的方法求解,在分析该方法的复杂度时,我们会经常得到\(O(\frac{n}{1}+\frac{n}{2}+\ldots+\frac{n}{n})\)的复杂度,然后我们都知道这个式子是等价于\(O(n\log n)\)的。在筛素数、字符串连续重复子串等很多算法中都有用到,用处之广,性能之优。今天不妨来证明下这个等价式。

\(O(\frac{n}{1}+\frac{n}{2}+\ldots+\frac{n}{n})\)~\(O(n\log n)\)

分析

要证明\(O(\frac{n}{1}+\frac{n}{2}+\ldots+\frac{n}{n})\)~\(O(n\log n)\),只需证\(O(1+\frac{1}{2}+\ldots+\frac{1}{n})\)~\(O(\ln n)\)——式子1.

为了证明式子1,需要证明4个定理:

1. 确界存在定理

2. 单调有界数列必定收敛

3. 数列\(\{(1+\frac{1}{n})^n\}\)单调增加,\(\{(1+\frac{1}{n})^{n+1}\}\)单调减少,两者收敛于同一极限

4. \(b_n=1+\frac{1}{2}+\ldots+\frac{1}{n}-\ln n\) 收敛

确界存在定理——实数系连续性定理

描述:非空有上界的数集必有上确界;非空有下界的数集必有下确界。

证明:<a href="http://baike.baidu.com/link?url=NojNWL0qcJWW20mxekU1GcfeD1Tp-0-JtF4oyRio7w9Th-ifVdybvf3PaSzmZKXZywk9mCGnlQ9mkPk-NySvn1b_DUqfs-Boez0kGdtwwgN9cgLyW4xYJfHUVGhdiIHQh3hilNBwuRJxARBOsGkqrSpDfuuEwxa56NjmcEgwrlSiddNVOjmmH6WirfboiLNMPo0a06RFoeFRasTvJSUaPa

" target="_blank">百度百科关于确界存在定理

单调有界数列必定收敛

描述:单调递增且有上界数列必定收敛,单调递减且有下界数列必定收敛

证明:

不妨设数列\(\{x_n\}\)单调增加且有上界,根据确界存在定理,由\(\{x_n\}\)构成的数集必有上确界\(\beta\),满足:

(1) \(\forall n\in N^+:x_n\leq \beta;\)

(2) \(\forall \epsilon>0,\exists x_{n_0}:x_{n_0}>\beta-\epsilon。\)

取\(N=n_0\),\(\forall n>N:\beta-\epsilon<x_{n_0}\leq x_n\leq \beta\),因而\(\{x_n-\beta\}<\epsilon\),于是得到\(\lim_{n\rightarrow \infty}x_n=\beta\)

当数列单调递减且有下界时,同理。

数列\(\{(1+\frac{1}{n})^n\}\)单调增加,\(\{(1+\frac{1}{n})^{n+1}\}\)单调减少,两者收敛于同一极限

证明:

记\(x_n=\{(1+\frac{1}{n})^n\}\),\(y_n=\{(1+\frac{1}{n})^{n+1}\}\),利用平均不等式\(\sqrt[n]{a_1a_2\ldots a_n}\leq\frac{a_1+a_2+\ldots+a_n}{n}\)得到

\(x_n=\{(1+\frac{1}{n})^n\}\cdot1\leq [\frac{n(1+\frac{1}{n})+1}{n+1}]^{n+1}=x_{n+1}\)

\(\frac{1}{y_n}=(\frac{n}{n+1})^{n+1}\cdot 1\leq[\frac{(n+1)\frac{n}{n+1}+1}{n+2}]^{n+2}=\frac{1}{y_{n+1}}\)

这表示\(\{x_n\}\)单调增加,而\(\{y_n\}\)单调减少。又由于\(2=x_1\leq x_n<y_n\leq y_1=4\),可知数列\(\{x_n\}\),\(\{y_n\}\)都收敛(单调有界数列必收敛)。

因为\(y_n=x_n(1+\frac{1}{n})\),所以它们具有相同的极限。习惯上用字母\(e\)来表示这一极限,即\(\lim_{n\rightarrow\infty}(1+\frac{1}{n})^n=\lim_{n\rightarrow\infty}(1+\frac{1}{n+1})^n=e\)

\(e=2.718\ 281\ 828\ 459\cdots\)是一个无理数。以\(e\)为底的对数称为自然对数,通常即为\(\ln x(=\log_ex)\)。

\(b_n=1+\frac{1}{2}+\ldots+\frac{1}{n}-\ln n\) 收敛

证明:

由上一定理可知,\((1+\frac{1}{n})^n<e<(1+\frac{1}{n})^{n+1}\),由此得到\(\frac{1}{n+1}<\ln \frac{n+1}{n}<\frac{1}{n}\)。于是有:

\(b_{n+1}-b_n=\frac{1}{n+1}-\ln (n+1)+\ln n=\frac{1}{n+1}-\ln \frac{n+1}{n}<0\)

\(b_n=1+\frac{1}{2}+\frac{1}{3}+\ldots+\frac{1}{n}-\ln n>\ln \frac{2}{1}+\ln \frac{3}{2}+\ln \frac{4}{3}+\ldots+\ln \frac{n+1}{n}-\ln n=\ln (n+1)-\ln n>0\)

这说明数列\(\{b_n\}\)单调减少有下界,从而收敛。(单调有界数列必收敛)

总结

已证明\(O(1+\frac{1}{2}+\ldots+\frac{1}{n})\)~\(O(\ln n)\),因此可知\(O(\frac{n}{1}+\frac{n}{2}+\ldots+\frac{n}{n})\)~\(O(n\log n)\)

以后可能会附上用此公式的算法题目\(\ldots\)(待续)

证明 O(n/1+n/2+…+n/n)=O(nlogn)的更多相关文章

  1. PAT复杂度_最大子列和问题、最大子列和变种

    01-复杂度1. 最大子列和问题 给定K个整数组成的序列{ N1, N2, ..., NK },“连续子列”被定义为{ Ni, Ni+1, ..., Nj },其中 1 <= i <= j ...

  2. k-d tree 学习笔记

    以下是一些奇怪的链接有兴趣的可以看看: https://blog.sengxian.com/algorithms/k-dimensional-tree http://zgjkt.blog.uoj.ac ...

  3. BZOJ2733 [HNOI2012]永无乡

    直接平衡树启发式合并就好了...貌似是个很高端的东西.. 貌似可以证明splay的启发式合并是均摊$O(nlogn)$的...而其他平衡树都不行,所以其他的复杂度都是$O(nlog^2n)的$的 所以 ...

  4. 评估Divide and Conquer算法时间复杂度的几种策略

    算法导论的第四章对于divide-conquer进行了阐述, 感觉这本书特别在,实际给出的例子并不多,更多其实是一些偏向数学性质的分析, 最重要的是告诉你该类算法分析的一般性策略. 估计 首先是估计算 ...

  5. 算法(第四版)C# 习题题解——2.4

    写在前面 整个项目都托管在了 Github 上:https://github.com/ikesnowy/Algorithms-4th-Edition-in-Csharp 查找更方便的版本见:https ...

  6. Luogu P3384 【【模板】树链剖分】

    转载请注明出处,部分内容引自banananana大神的博客 ~~别说你不知道什么是树~~╮(─▽─)╭(帮你百度一下) 先来回顾两个问题:1,将树从x到y结点最短路径上所有节点的值都加上z 这也是个模 ...

  7. 几道splay

    hdu 1890 题意:每次将第i位到第i小数字所在的位置之间的位置翻转,每次输出第i小数字所在的位置 分析: 简单的splay处理区间翻转问题 有三点需要注意: 1.区间是1~n+2 2.此题里的查 ...

  8. Educational Codeforces Round 75 (Rated for Div. 2)

    知识普及: Educational使用拓展ACM赛制,没有现场hack,比赛后有12h的全网hack时间. rank按通过题数排名,若通过题数相等则按罚时排名. (罚时计算方式:第一次通过每题的时间之 ...

  9. 【学术篇】CF833B TheBakery 分治dp+主席树

    题目の传送门~ 题目大意: 将\(n\)个蛋糕分成恰好\(k\)份, 求每份中包含的蛋糕的种类数之和的最大值. 这题有两种做法. 第一种是线段树优化dp, 我还没有考虑. 另一种就是分治+主席树. 然 ...

随机推荐

  1. Serv-U 的升级及数据备份和迁移【转】

    Serv-U 配置备份   在serv-u7.x及以上版本安装目录下,有一个文件Serv-U.Archive是serv-u的配置文件,有一个users文件夹是Serv-U的域和用户的信息,那么我们只需 ...

  2. C#连接MySQL 操作步骤

    1.工具安装: 安装 MySQL For Windows,这个不多说,上官网下载: 安装mysql-connector-net,这个是MySQL数据库.NET开发驱动,因为C#是.NET架构的,所以需 ...

  3. 【自用】bat ftp下载前一天备份

    @echo off rem 指定FTP用户名 set ftpUser=app rem 指定FTP密码 set ftpPass=app rem 指定FTP服务器地址 set ftpIP=192.168. ...

  4. 03 Editor plugins and IDEs 编辑器插件和 ide

    Editor plugins and IDEs  编辑器插件和 ide Introduction  介绍 Options 选项   Introduction 介绍 This document list ...

  5. redis入门——redis常用命令

    http://blog.csdn.net/wclxyn/article/details/8449082 https://jingyan.baidu.com/article/90bc8fc87ce8e2 ...

  6. java基础80 JSP基础知识点(网页知识)

    1.jsp的引入 Servlet的作用:用java语言开发动态资源技术!!!    Jsp的作用:用java语言(+html语言)开发的动态资源技术!!!     Jsp就是servlet 问题:为什 ...

  7. ZooKeeper常见问题

    转载自原文:zookeeper(二)常见问题汇总 一.为什么zookeeper要部署基数台服务器? 所谓的zookeeper容错是指,当宕掉几个zookeeper服务器之后,剩下的个数必须大于宕掉的个 ...

  8. 20165203 2017-2018-2 《Java程序设计》课程总结

    20165203 2017-2018-2 <Java程序设计>课程总结 一.每周作业及实验报告链接汇总 我期望的师生关系(预备作业一):浅谈一下对师生关系的看法和对自己未来学习和生活的期望 ...

  9. ubuntu16.04 安装 caffe cuda 相关流程

    不多说了,经历了很多莫名其妙的错误最后终于安装好了,直接放安装脚本: #!/bin/bash #安装时要注意有些库可能安装失败以及安装caffe有和protobuf相关错误时可能需要重新对protob ...

  10. java解析Xml格式的字符串

    最近在工作中,需要调别的接口,接口返回的是一个字符串,而且内容是xml格式的,结果在解析json的时候报错,最终修改了接口的返回方式,以Map返回, 才得以接收到这个xml的字符串,然后通过dom4j ...