P1943-LocalMaxima_NOI导刊2009提高(1)【数论】
正题
题目链接: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∑xPk−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∑nn!∑k=1x(Px−1−1)∗(n−x)!
打表得(我不会证明啊)
⇒∑x=1n1x\Rightarrow \sum_{x=1}^n\frac{1}{x}⇒x=1∑nx1
然后过掉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∑nx1=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)【数论】相关推荐
- P1950 长方形_NOI导刊2009提高(2)[单调栈][贡献法]
P1950 长方形_NOI导刊2009提高(2) 题意:给你\(n\times m\)的矩形,求没有*的子矩形数量.\(1\leq n,m\leq 1000\). 数据比较弱的题目是luoguP119 ...
- 洛谷——P1951 收费站_NOI导刊2009提高(2)
https://www.luogu.org/problem/show?pid=1951 题目描述 在某个遥远的国家里,有n个城市.编号为1,2,3,-,n. 这个国家的政府修建了m条双向的公路.每条公 ...
- [洛谷P1951]收费站_NOI导刊2009提高(2)
题目大意:有一张$n$个点$m$条边的图,每个点有一个权值$w_i$,有边权,询问从$S$到$T$的路径中,边权和小于$s$,且$\max\limits_{路径经过k}\{w_i\}$最小,输出这个最 ...
- P1800 software_NOI导刊2010提高(06)
P1800 software_NOI导刊2010提高(06) 题目描述 一个软件开发公司同时要开发两个软件,并且要同时交付给用户,现在公司为了尽快完成这一任务,将每个软件划分成m个模块,由公司里的技术 ...
- P1801 黑匣子_NOI导刊2010提高(06)
P1801 黑匣子_NOI导刊2010提高(06) 题目描述 Black Box是一种原始的数据库.它可以储存一个整数数组,还有一个特别的变量i.最开始的时候Black Box是空的.而i等于0.这个 ...
- 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 ...
- 洛谷 P1800 software_NOI导刊2010提高(06)(二分答案+DP检验)
P1800 software_NOI导刊2010提高(06) 标签 二分答案 难度 普及/提高- 题目描述 一个软件开发公司同时要开发两个软件,并且要同时交付给用户,现在公司为了尽快完成这一任务,将每 ...
- 方程的解_NOI导刊2010提高(01) 组合数
题目描述 佳佳碰到了一个难题,请你来帮忙解决. 对于不定方程a1+a2+-+ak-1+ak=g(x),其中k≥2且k∈N,x是正整数,g(x)=x^x mod 1000(即x^x除以1000的余数), ...
- 洛谷 P1795 无穷的序列_NOI导刊2010提高(05)
P1795 无穷的序列_NOI导刊2010提高(05) 题目描述 有一个无穷序列如下: 110100100010000100000- 请你找出这个无穷序列中指定位置上的数字 输入输出格式 输入格式: ...
最新文章
- python 基础 day3
- 一篇虚拟试穿的论文介绍
- iOS11和iPhone X的适配
- 《小印记》源码分享--极光推送服务器篇
- linux系统各文件夹的作用,linux系统文件夹的作用 good
- 关于将一个数输出为原码、反码和补码的解惑
- 010_jQuery获取和设置内容属性
- R语言应用实战-基于R的C4.5算法和C5.0算法原理解析及应用案例
- centos7 dns配置_Centos7.7 安装FreeIPA (三)
- leetcode(1)485——最大连续 1 的个数(C++)
- [Mac]Python 安装MySQLdb模块
- 简约几何艺术海报PSD分层模板,即使简单也足以控制观众的注意力。
- [JNI] 开发前言
- 免费的mysql图形化界面_免费好用的MySQL数据库客户端图形界面管理工具
- 笔记本显示电源已连接但是未充电的简单解决办法
- 美团实习面(45min + 35min)
- 7. F1方程式冠军
- 生育医疗费用出院结算流程
- 圣光机大学计算机专业,25位师生顺利完成圣光机大学实习项目
- js promise的用法
热门文章
- ios把数据传递到另一个页面_iOS 委托 页面之间传递数值
- java递归 treeview_treeview递归
- css-6 df15,webpack 样式文件的代码分割(15)
- 神州云科DCN存储管理IP_干货!DCN校园网大二层解决方案详解
- python3怎么安装gmpy2_python2/3 模块gmpy2在linux下安装
- 机器学习之格式变化——reshape(-1,1)
- [Swagger2]拓展:其他皮肤
- C++实现顺序串(完整代码)
- LeetCode 872叶子相似的树-简单
- poj3981 字符串替换-字符串的基本操作