正题

题目链接:https://www.luogu.org/problemnew/show/P1943


题目大意

定义LocalLocalLocal数为一个数且比它前面的数字都要大。

求一个随机长度为nnn的序列中LoaclLoaclLoacl数的期望数量。


解题思路

ans=Local(a)n!ans=\frac{Local(a)}{n!}ans=n!Local(a)​

现在我们分开求,先考虑如何求出所有序列中LocalLocalLocal数的数量

一个数xxx,有x−1x-1x−1个数比他小,那么若一个数成为LocalLocalLocal数的排列数是
∑k=1xPk−1x−1∗Pn−xn−x=Px−1k−1∗(n−x)!\sum_{k=1}^xP_{k-1}^{x-1}*P_{n-x}^{n-x}=P^{k-1}_{x-1}*(n-x)!k=1∑x​Pk−1x−1​∗Pn−xn−x​=Px−1k−1​∗(n−x)!
也就是答案为
∑x=1n∑k=1x(Px−1−1)∗(n−x)!n!\sum_{x=1}^n\frac{\sum_{k=1}^x(P^{-1}_{x-1})*(n-x)!}{n!}x=1∑n​n!∑k=1x​(Px−1−1​)∗(n−x)!​
打表得(我不会证明啊)
⇒∑x=1n1x\Rightarrow \sum_{x=1}^n\frac{1}{x}⇒x=1∑n​x1​

然后过掉n≤1000000n \leq 1000000n≤1000000的点,然后使出调和级数
∑x=1n1x=log⁡(n+1)+r(n→∞)\sum_{x=1}^n\frac{1}{x}=\log (n+1)+r(n\rightarrow \infty)x=1∑n​x1​=log(n+1)+r(n→∞)
rrr为欧拉函数r≈0.5772156649r\approx 0.5772156649r≈0.5772156649

鬼知道考试的时候谁做的出来\color{white}\text{鬼知道考试的时候谁做的出来}鬼知道考试的时候谁做的出来


codecodecode

// luogu-judger-enable-o2
#include<cstdio>
#include<cmath>
using namespace std;
double n,ans;
int main()
{scanf("%lf",&n);if(n<=1000000)for(int i=1;i<=n;i++)ans+=((double)1/i);else ans=log(n+1)+0.5772156649;printf("%.8lf",ans);
}

P1943-LocalMaxima_NOI导刊2009提高(1)【数论】相关推荐

  1. P1950 长方形_NOI导刊2009提高(2)[单调栈][贡献法]

    P1950 长方形_NOI导刊2009提高(2) 题意:给你\(n\times m\)的矩形,求没有*的子矩形数量.\(1\leq n,m\leq 1000\). 数据比较弱的题目是luoguP119 ...

  2. 洛谷——P1951 收费站_NOI导刊2009提高(2)

    https://www.luogu.org/problem/show?pid=1951 题目描述 在某个遥远的国家里,有n个城市.编号为1,2,3,-,n. 这个国家的政府修建了m条双向的公路.每条公 ...

  3. [洛谷P1951]收费站_NOI导刊2009提高(2)

    题目大意:有一张$n$个点$m$条边的图,每个点有一个权值$w_i$,有边权,询问从$S$到$T$的路径中,边权和小于$s$,且$\max\limits_{路径经过k}\{w_i\}$最小,输出这个最 ...

  4. P1800 software_NOI导刊2010提高(06)

    P1800 software_NOI导刊2010提高(06) 题目描述 一个软件开发公司同时要开发两个软件,并且要同时交付给用户,现在公司为了尽快完成这一任务,将每个软件划分成m个模块,由公司里的技术 ...

  5. P1801 黑匣子_NOI导刊2010提高(06)

    P1801 黑匣子_NOI导刊2010提高(06) 题目描述 Black Box是一种原始的数据库.它可以储存一个整数数组,还有一个特别的变量i.最开始的时候Black Box是空的.而i等于0.这个 ...

  6. P1799 数列_NOI导刊2010提高(06)

    P1799 数列_NOI导刊2010提高(06) f[i][j]表示前i个数删去j个数得到的最大价值. if(i-j==x) f[i][j]=max(f[i][j],f[i-1][j]+1); els ...

  7. 洛谷 P1800 software_NOI导刊2010提高(06)(二分答案+DP检验)

    P1800 software_NOI导刊2010提高(06) 标签 二分答案 难度 普及/提高- 题目描述 一个软件开发公司同时要开发两个软件,并且要同时交付给用户,现在公司为了尽快完成这一任务,将每 ...

  8. 方程的解_NOI导刊2010提高(01) 组合数

    题目描述 佳佳碰到了一个难题,请你来帮忙解决. 对于不定方程a1+a2+-+ak-1+ak=g(x),其中k≥2且k∈N,x是正整数,g(x)=x^x mod 1000(即x^x除以1000的余数), ...

  9. 洛谷 P1795 无穷的序列_NOI导刊2010提高(05)

    P1795 无穷的序列_NOI导刊2010提高(05) 题目描述 有一个无穷序列如下: 110100100010000100000- 请你找出这个无穷序列中指定位置上的数字 输入输出格式 输入格式: ...

最新文章

  1. python 基础 day3
  2. 一篇虚拟试穿的论文介绍
  3. iOS11和iPhone X的适配
  4. 《小印记》源码分享--极光推送服务器篇
  5. linux系统各文件夹的作用,linux系统文件夹的作用 good
  6. 关于将一个数输出为原码、反码和补码的解惑
  7. 010_jQuery获取和设置内容属性
  8. R语言应用实战-基于R的C4.5算法和C5.0算法原理解析及应用案例
  9. centos7 dns配置_Centos7.7 安装FreeIPA (三)
  10. leetcode(1)485——最大连续 1 的个数(C++)
  11. [Mac]Python 安装MySQLdb模块
  12. 简约几何艺术海报PSD分层模板,即使简单也足以控制观众的注意力。
  13. [JNI] 开发前言
  14. 免费的mysql图形化界面_免费好用的MySQL数据库客户端图形界面管理工具
  15. 笔记本显示电源已连接但是未充电的简单解决办法
  16. 美团实习面(45min + 35min)
  17. 7. F1方程式冠军
  18. 生育医疗费用出院结算流程
  19. 圣光机大学计算机专业,25位师生顺利完成圣光机大学实习项目
  20. js promise的用法

热门文章

  1. ios把数据传递到另一个页面_iOS 委托 页面之间传递数值
  2. java递归 treeview_treeview递归
  3. css-6 df15,webpack 样式文件的代码分割(15)
  4. 神州云科DCN存储管理IP_干货!DCN校园网大二层解决方案详解
  5. python3怎么安装gmpy2_python2/3 模块gmpy2在linux下安装
  6. 机器学习之格式变化——reshape(-1,1)
  7. [Swagger2]拓展:其他皮肤
  8. C++实现顺序串(完整代码)
  9. LeetCode 872叶子相似的树-简单
  10. poj3981 字符串替换-字符串的基本操作