Problem Description

百小度最近在逛博客,然后发现了一个有趣的问题。

如下图所示,是一个12 
位数014326951987 
, 它的数字先逐渐变大, 然后变小,再变大,接着变小,又变大,最后变小。我们就称,其共包含6 
个单调区间。

现在问题来了:一个n位数平均包含多少个单调区间?单调区间的平均长度又是多少?

因为我们考虑到这样的数样本太大,有10 n  
这么多,所以百小度决定缩小样本,假定任意两位相邻数字不能相同,而且这个n 
位数允许以0 
开头。现在我已经将样本大小已经被缩小到10∗9 n−1  
,接下来把这个问题交给你,请你开启大脑挖掘机,挖挖答案在哪里。

Input

第一行为T 
,表示输入数据组数。

下面T行,每行包含一个正整数n 
,n 
为不大于100000 
的正整数。

Output

对第i组数据,输出

Case #i:

然后输出两个实数,用空格隔开,分别为平均单调区间数和单调区间平均长度,结果保留六位小数。

Sample Input
2
2
12
Sample Output
Case #1:
1.000000 2.000000
Case #2:
8.037037 2.368664

原题是IMO候选题:

考虑这么一个  位数  ,如图所示,它的数字先逐渐变大,然后开始变小,再变大,再变小,再变大,再变小。我们就说,它一共包含了  个单调区间。我们的问题就是:一个 n 位数平均有多少个单调区间?为了避免歧义,我们假设任意两位相邻的数字都不相同,因而像  这样的数我们就不考虑了。另外,大家可能已经注意到了,我们允许这个 n 位数以数字  开头。因而,更精确地说,我们的问题是:相邻数字都不相同的、允许以  开头的所有 n 位数当中,平均有多少个单调区间?

    这个题目来自  年 IMO 候选题。

    让我们把所有这种 n 位数的个数记作 N 。那么 N 等于多少?这个 n 位数的第一位有  种选择,今后的每一位都只有  种选择(因为要跟前一位不一样),因而 n 位数一共有 N =  · 9n- 个。接下来,我们要求的就是,所有 n 位数当中的所有单调区间一共有多少个。我们换一种方法来累计这些单调区间:先算所有从第一位开始的单调区间,再算所有从第二位开始的单调区间,等等,最后算所有从第 n 位开始的单调区间。如果用 ri 来表示所有从第 i 位开始的单调区间的数目,那么我们要求的平均单调区间数就是 (r1 + r2 + … + rn) / N ,也就是 r1 / N + r2 / N + … + rn / N 。注意到其中的每一项 ri / N 其实就是从 N 个合法的 n 位数中任取一个后,存在以第 i 位数打头的单调区间的概率。因此,我们只需要求出这 n 个概率值,加起来便是我们想要的答案了。
显然, r1 / N = ,因为第一位数字必然会引领一个单调区间。显然, rn / N = ,因为最后一位数字不可能引领一个新的单调区间。那么,对于其他的 ri / N 呢?注意到,第 i – 位、第 i 位和第 i + 位的大小关系一共可能有以下四种情况: 其中,只有第三种情况和第四种情况下,第 i 位才会成为一个新的单调区间的开始。为了计算这两种情况发生的概率,我们只需要算出情况 和情况 发生的概率,再用 来减即可。情况 发生的概率有多大呢?三位数字串一共有 · 个(第一位有 种选择,后面的每一位都只有 种选择,因为要跟前一位不一样)。为了得到递增的数字串,我们只需要选出三个不同的数字,然后把它们从小到大排列即可,这一共有 C(, ) 种方法。因此,情况 的发生概率就是 C(, ) / ( · ) = / 。同理,情况 的发生概率也是 / ,两者加起来就是 / ;反过来,情况 和情况 出现的概率就是 – / = / 了。
因此,我们最终要求的答案就是 + / + / + … + / + = + (n – ) · / 。 这个结论还会引出很多有意思的问题。在一个 位数当中,平均会产生 个单调区间。我们似乎发现了一个很不合理的地方:这岂不意味着,平均每个单调区间的长度只有 / = 1.45 个数字吗?考虑到单调区间的长度不可能恰好是 1.45 个数字,为了得到 1.45 这个平均长度,一定有些区间的长度比 1.45 小,有些区间的长度比 1.45 大。有些区间的长度比 1.45 小,这不就意味着这些区间的长度为 吗?而一个区间的长度显然是不可能为 的。怎么回事?
其实, / = 1.45 这个算式是错的。在这 个单调区间中,除了最后一个区间以外,每一个区间的最后一个数与下一个区间的第一个数都是公共的。因此,这个 位数当中,有 个数被重复使用了。所以,在一个 位数当中,单调区间的平均长度应该是 ( + ) / = 2.4 。
类似的, n 位数的单调区间的平均长度为 (n + (/)(n – )) / ( + (/)(n – )) = (46n – ) / (19n – ) = ( – /n) / ( – /n) 。当 n 无穷大时,其极限为 / 。 参考资料:Ross Honsberger, From Erdos to Kiev: Problems of Olympiad Caliber, pp. -

代码:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <set>
#include <map>
#include <queue>
#include <string>
#define LL long long int main()
{
//freopen("test.in", "r", stdin);
int T, n;
double ans;
scanf("%d", &T);
for (int times = ; times <= T; ++times)
{
printf("Case #%d:\n", times);
scanf("%d", &n);
ans = + (n-2.0)*/;
printf("%lf ", ans);
ans = (-38.0/n) / (-11.0/n);
printf("%lf\n", ans);
}
return ;
}

ACM学习历程—BestCoder 2015百度之星资格赛1006 单调区间(组合数学)的更多相关文章

  1. ACM学习历程—BestCoder 2015百度之星资格赛1004 放盘子(策略 && 计算几何)

    Problem Description 小度熊喜欢恶作剧.今天他向来访者们提出一个恶俗的游戏.他和来访者们轮流往一个正多边形内放盘子.最后放盘子的是获胜者,会赢得失败者的一个吻.玩了两次以后,小度熊发 ...

  2. ACM学习历程—BestCoder 2015百度之星资格赛1001 大搬家(递推 && 组合数学)

    Problem Description 近期B厂组织了一次大搬家,所有人都要按照指示换到指定的座位上.指示的内容是坐在位置i 上的人要搬到位置j 上.现在B厂有N 个人,一对一到N 个位置上.搬家之后 ...

  3. ACM学习历程—BestCoder 2015百度之星资格赛1002 列变位法解密(vector容器)

    Problem Description 列变位法是古典密码算法中变位加密的一种方法,具体过程如下 将明文字符分割成个数固定的分组(如5个一组,5即为密钥),按一组一行的次序整齐排列,最后不足一组不放置 ...

  4. ACM学习历程—BestCoder 2015百度之星资格赛1003 IP聚合(set容器)

    Problem Description 当今世界,网络已经无处不在了,小度熊由于犯了错误,当上了度度公司的网络管理员,他手上有大量的 IP列表,小度熊想知道在某个固定的子网掩码下,有多少个网络地址.网 ...

  5. 模拟 2015百度之星资格赛 1003 IP聚合

    题目传送门 /* 模拟水题,排序后找出重复的ip就可以了 */ #include <cstdio> #include <iostream> #include <algor ...

  6. ACM学习历程—BestCoder Round #75

    1001:King's Cake(数论) http://acm.hdu.edu.cn/showproblem.php?pid=5640 这题有点辗转相除的意思.基本没有什么坑点. 代码: #inclu ...

  7. ACM学习历程—SNNUOJ 1116 A Simple Problem(递推 && 逆元 && 组合数学 && 快速幂)(2015陕西省大学生程序设计竞赛K题)

    Description Assuming a finite – radius “ball” which is on an N dimension is cut with a “knife” of N- ...

  8. ACM学习历程—NPU1045 2015年陕西省程序设计竞赛网络预赛(热身赛)C题 Graph Theory(递推 && 组合数学 && 大数)

    Description In graph theory, a matching or independent edge set in a graph G = (V , E) is a set of e ...

  9. ACM学习历程—NPU 2015年陕西省程序设计竞赛网络预赛(正式赛)F题 和谐的比赛(递推)

    Description 今天西工大举办了一场比赛总共有m+n人,但是有m人比较懒没带电脑,另外的n个人带了电脑.不幸的是,今天机房的电脑全坏了只能用带的电脑,一台电脑最多两人公用,确保n>=m. ...

随机推荐

  1. \\s+ split替换

    出自: http://www.tuicool.com/articles/vy2ymm 详解 "\\s+" 正则表达式中\s匹配任何空白字符,包括空格.制表符.换页符等等, 等价于[ ...

  2. 【翻译自mos文章】执行utlpwdmg.sql之后报ORA-28003, ORA-20001, ORA-20002, ORA-20003, ORA-20004 错误

    执行utlpwdmg.sql之后报ORA-28003, ORA-20001, ORA-20002, ORA-20003, ORA-20004 错误. 适用于: Oracle Server - Ente ...

  3. caffe-ubuntu1604-gtx850m-i7-4710hq----VGG_ILSVRC_16_layers.caffemodel

    c++调用vgg16: ./build/install/bin/classification \ /media/whale/wsWin10/wsCaffe/model-zoo/VGG16//deplo ...

  4. Linux内核源码分析方法_转

    Linux内核源码分析方法 转自:http://www.cnblogs.com/fanzhidongyzby/archive/2013/03/20/2970624.html 一.内核源码之我见 Lin ...

  5. Intellj IDEA光标替insert状态,back键无法删除内容

    Intellj IDEA光标为insert状态,无法删除内容导入项目后,发现打开java文件的光标是win系统下按了insert键后的那种宽的光标,并且还无法删除内容,且按删除(delete)键也只见 ...

  6. 【转】安卓apk反编译、修改、重新打包、签名全过程

    首先明确,反编译别人apk是一件不厚道的事情.代码是程序员辛苦工作的成果,想通过这种手段不劳而获,是不对的.这也说明,代码混淆是非常重要的.本文抱着学习的态度,研究在一些特殊的情况下如果有需要,该怎么 ...

  7. Unity3D研究院之Inspector视图中的get/set使用

    get set 使用起来很方便,但是编辑时在Inspector视图中问题就来了,因为get/set的属性即使是public了,但是在Inspector视图中依然不显示..谷歌一下估计就是下面这样的答案 ...

  8. Java 学习 day03

    01-语句(while) 02-语句(do while) 03-语句(for) 04-语句(for和while的区别) 05-语句(循环语句的其他特点) 06-语句(for语句练习-累加&计数 ...

  9. 九度OJ 1030:毕业bg (01背包、DP)

    时间限制:1 秒 内存限制:32 兆 特殊判题:否 提交:1814 解决:798 题目描述:     每年毕业的季节都会有大量毕业生发起狂欢,好朋友们相约吃散伙饭,网络上称为"bg" ...

  10. 洛谷 4568 [JLOI2011] 飞行路线

    题目戳这里 一句话题意: 有n个点,m条边的有向图,最多可以把k条边变为0,求从起点到终点最短距离. Solution 首先看到这题目,感觉贼难,看起来像DP,貌似也有大佬这么做,但鉴于本蒟蒻思维能力 ...