首页 > 其他分享 >「魔怔」鲜花 2023/10/9

「魔怔」鲜花 2023/10/9

时间:2023-10-10 16:33:07浏览次数:37  
标签:10 right frac 鲜花 varphi overline cdots 随机 2023

碎碎念。

嘿嘿,我的嘿嘿。今天怎么这么摆。

什么?你以为这篇文章是鲜花?这是压缩诈骗。

建议大家都学点新东西,不然就会像我一样在 NOIP 模拟赛中被 T2 放 Powerful Number 筛的出题人创死。

我尝试把开头写长一点让预览看不到下面的内容,这样你就要点进来看了,嘿嘿,我是不是很机智。

让我猜猜你随机读鲜花的停时的期望是多少呢?(记住这个问题,下面要考。)

其实这是一篇正经的科普文章。

一些约定

  • 随机过程:对一个参数集随机的一组变量集合,参数集通常是时间 \(T\),若 \(\forall t\in T\),\(X_t\) 为随机变量,那么 \(\{X_t|t\in T\}\) 就是一个随机过程。
  • 条件概率:\(P(A|B)\),即 \(A\) 事件在 \(B\) 事件条件下发生的概率。
  • 条件期望:\(E(A|B)\),类比条件概率,即 \(A\) 事件在 \(B\) 事件条件下的期望。
  • 几乎一定:事件 \(A\) 几乎一定发生,当且仅当 \(P(A)=1\)。\(A\) 不发生的情况集合可能非空,但它们的概率为 \(0\)。样本空间有限时,“几乎一定”和“一定”通常没有区别;但是当样本空间无限时,就有区别了。比如 \([0,1]\) 中均匀随机一个实数 \(x\),则 \(P(x>1)=1\),可以说选出的 \(x\) 几乎一定大于 \(0\),但理论上确实可以存在 \(x=0\) 的样本。

离散时间鞅

一个离散时间鞅为一个以时间为参数集合的随机过程,即随机过程 \(\{X_0,X_1,\cdots \}\),满足对于任意时刻 \(n\in \mathbb{N}\),均有:

  • \(E(|X_n|)<∞\)。
  • \(E(X_{n+1}-X_n|X_{0},X_{1},\cdots ,X_{n})=0\)。

第二个条件看起来不知所云,其实意思就是当 \(X_{0},X_{1},\cdots ,X_{n}\) 已经确定时,\(X_{n+1}\) 的期望等于 \(X_n\)。类似一个公平赌博游戏,若我们已经得到了前 \(n\) 场赌博后你所拥有的财产(观测值),那么要求第 \(n+1\) 场的观测值的期望等于第 \(n\) 场的观测值(不亏不赚),才能满足公平赌博游戏的定义(所以公平赌博游戏其实差不多就是鞅www)。

同时不难发现鞅的线性加减、增减常数都不影响其鞅的性质。

停时

关于随机过程 \(\{X_{0},X_{1},\cdots\}\) 的停时为一个非负的随机变量 \(T\)(可能为 \(∞\),即操作无限进行),满足对任意时刻 \(n\),你可以通过观测 \(X_0,X_1,\cdots,X_{n}\) 的取值得出 \(n\) 和 \(T\) 的大小关系/等量关系。即 \([n=T],[n<T],[n>T]\) 三个随机变量的取值仅与 \(X_{0},X_{1},\cdots ,X_{n}\) 有关。

对于随机过程 \(\{X_0,X_1,\cdots\}\),对应带停时为 \(T\) 的随机过程 \(\{\overline {X_0},\overline {X_1},\cdots\}\) 的定义为:

\[\overline {X_n}=\begin{cases}X_n\quad n\le T\\X_T\quad n>T\end{cases} \]

即 \(\{\overline {X_{1}},\overline {X_{2}},\cdots \}\) 为将 \(\{X_0,X_1,\cdots \}\) 在 \(T\) 之后的项全部覆盖为 \(T\) 后的随机过程,可以看作在时刻 \(T\) 终止随机操作。(所以停时也可以理解为操作的停止时间。)

由于 \(\overline{X_n}\) 中 \(n\le T\) 时 \(\overline{X_n}=X_n\) 期望不变,而 \(n>T\) 时我们钦定它期望不变,所以满足鞅期望观测值不变的性质,所以 \(\{\overline{X_1},\overline{X_2},\cdots\}\) 也是一个鞅。

鞅的停时定理

设 \(T\) 为离散时间鞅 \(\{X_0,X_1,\cdots\}\) 的停时,且 \(T\) 几乎一定有限,当下面三个条件之一成立时,有 \(E(X_T)=E(X_0)\)

  • \(T\) 几乎一定有界,即存在常数 \(K\),\(P(T\le K)=1\)。
  • \(E(|X_{i+1}-X_i|)\) 几乎一定有界,且 \(E(T)\) 有限,即存在常数 \(K\) 使得 \(P(E(|X_{i+1}-X_i|)\le K)=1\)
  • \(\overline{X_i}\) 几乎一定有界,即存在常数 \(K\) 使得 \(P(|\overline{X_i}|\le K)=1\)

为了帮助理解,给出下面若干例子:

  • 数轴上从 \(0\) 开始随机向右游走,每次走 \(1\) 或 \(2\) 的距离,走过一个常数 \(K\) 后停止。那么 \(T\) 一定有界,满足第一个条件。
  • 从 \([0,1]\) 中随机选择实数,直到选出非 \(0\) 数为止。\(T\) 几乎一定有界,满足第一个条件。
  • 数轴上从 \(0\) 开始随机左右游走,每次走 \(1\) 的距离,走到区间 \([-K,K]\)(\(K\) 为常数) 边界后停止。由于 \(\overline{X_n}\) 有界 \([-K,K]\),满足第三个条件。
  • 我们回到开头那个问题(还记得吗?)。你在读这篇文章,假设你每次随机向目前阅读位置的后面跳一个位置然后开始读,读完为止。那么很显然 \(\overline{X_n}\) 有界(文章长度),也满足第三个条件。
  • 第二个条件是大多数题目所包含的(比如让你求停时的期望,每次操作带有的变化量有界等等)。

势能分析

其实前面这么一大堆定义、性质,统统不用记

考虑经典模型:

假设现在有随机过程 \(\{X_0,X_1,\cdots \}\),\(T\) 为其停时,给定初始状态 \(X_0\),终止状态 \(\overline{X_T}\),求 \(E(T)\)。

考虑构造势能函数 \(\varphi(X)\),同时描述了时间状态两个变量:

  • \(E(\varphi(X_{n+1})-\varphi(X_{n})|X_{0},X_{1},\cdots,X_{n})=-1\)
  • \(\varphi(X_{T})\) 为常数,且 \(\nexists i<T,\varphi(X_{i})=\varphi(X_T)\)

第一个限制保证了每个时刻,势能的变化量期望为 \(-1\),那么 \(E(\varphi(X_0))-E(\varphi(X_n))=n\);第二个限制保证了末状态势能唯一,也就是保证末状态为第一个满足终止条件的状态的需要。

构造 \(Y_i=\varphi(X_i)+i\),根据定义,\(\{Y_0,Y_1,\cdots\}\) 为一个离散时间鞅。那么根据停时定理,可以得到 \(E(Y_T)=E(Y_0)\),即 \(E(\varphi(X_T)+T)=E(\varphi(X_0))\),即 \(E(T)=E(\varphi(X_0))-E(\varphi(X_T))\),根据定义右边两项都是常数值,于是就能方便地求出 \(E(T)\)。

如何准确地构造势能函数,一般构造关于题面中状态有关量为自变量的函数,然后根据第一条限制解函数方程/递推。

你可能在高中物理课听说过势能是相对的,这允许我们对于一些常数 \(C\),钦定 \(f(C)\) 的值。

CF1025G Company Acquisitions

定义整个局面 \(X\) 的势能函数为 \(\varphi(X)\),由于我们只关心每个点跟在屁股后面的点的个数,定义单个点 \(u\) 的势能为 \(f(c_u)\),\(c_u\) 为局面 \(X\) 下跟在 \(u\) 后的点的个数,\(f\) 为关于 \(c\) 的函数,令 \(\varphi(X)=\sum\limits_{u}f(c_u)\)。

则 \(X\to X'\) 的势能变化量只和每次随机选出来的两个点 \(u,v\) 的 \(c\) 有关,令 \(c_u=x,c_v=y\):

\[\begin{aligned}E(\varphi(X'))-E(\varphi(X))&=-1\\E(\varphi(X'))&=E(\varphi(X))-1\\\frac{1}{2}(f(x+1)+yf(0))+\frac{1}{2}(f(y+1)+xf(0))&=(f(x)+f(y))-1\end{aligned} \]

我们令 \(f(0)=0\):

\[\frac{1}{2}f(x+1)+\frac{1}{2}f(y+1)=f(x)+f(y)-1 \]

要求对任意 \(x,y\) 均成立,那么一个简单的想法是直接拆两半:

\[\frac{1}{2}f(x+1)=f(x)-\frac{1}{2} \]

结合 \(f(0)=0\) 得到:

\[f(x)=1-2^x \]

那么答案就是 \(\left(\sum\limits_{u}f(c_u)\right)- f(n)\),复杂度 \(O(n)\)。

CF1479E School Clubs

我们只关心每个组内的人数,所以依旧令局面 \(X\) 的势能值 \(\varphi(X)=\sum\limits_{i=1}^{m}f(a_i)\) 为各个组势能之和,一组的势能定义为其人数 \(a_i\) 的 \(f\) 值。

每次操作,将随机发电的那个人所在的组拎出来,设为第 \(i\) 组:

  • 这个人新开了一个组:\(\varphi(X')=\varphi(X)-f(a_i)+f(a_i-1)+f(1)\)。
  • 这个人回到了原来的组:\(\varphi(X')=\varphi(X)\)。
  • 这个人跑到了其他组,设为 \(j\) 组:\(\varphi(X')=\varphi(X)-f(a_i)-f(a_j)+f(a_i-1)+f(a_j+1)\)

然后就是一些 Dirty Work 了:

\[\begin{aligned}&E(\varphi(X')-\varphi(X)|X)\\=&\sum\limits_{i=1}^m\frac{a_i}{2n}\left(-f(a_i)+f(a_i-1)+f(1)+\frac{a_i}{n}\varphi(X)+\sum\limits_{j\neq i}\frac{a_j}{n}\left(-f(a_i)-f(b_i)+f(a_i-1)+f(b_i+1)\right)\right)\\=&\sum\limits_{i=1}^m\frac{a_i}{2n}\left(-\left(2-\frac{a_i}{n}\right)f(a_i)+\left(2-\frac{a_i}{n}\right)f(a_i-1)+f(1)+\sum\limits_{j\neq i}(f(a_j+1)-f(a_j))\right)\\=&\frac{1}{2}f(1)+\sum\limits_{i=1}^n\frac{a_i}{2n}\left(-\left(3-\frac{2a_i}{n}\right)f(a_i)+\left(2-\frac{a_i}{n}\right)f(a_i-1)+\left(1-\frac{a_i}{n}\right)f(a_i+1)\right)\\=&-1\end{aligned} \]

不妨设 \(f(1)=-2\):

\[\sum\limits_{i=1}^n\frac{a_i}{2n^2}\left(-(3n-2a_i)f(a_i)+(2n-a_i)f(a_i-1)+(n-a_i)f(a_i+1)\right)=0 \]

\[-(3n-2a)f(a)+(2n-a)f(a-1)+(n-a)f(a+1)=0 \]

然后 \(O(n)\) 递推 \(f\) 的值即可。这是能过的。

标签:10,right,frac,鲜花,varphi,overline,cdots,随机,2023
From: https://www.cnblogs.com/Ender32k/p/17755046.html

相关文章

  • 【2023年10月10日】STF60_docker_Day01(下午)
     STF60_docker_Day01(下午)容器运行先导入镜像 dockerload</home/centos-lamp.tar 给导入的镜像命名 dockertag0b8dnickistre/centos-lamp.tar 交互式运行容器一般就是临时用用,看看配置文件等等dockerrun-it镜像id或镜像名:tag版本/bin/bash或bash......
  • 单机10万TCP连接测试记录
    转自:https://www.cnblogs.com/fuhua/p/16904864.html单机10万TCP连接测试记录 目录前言准备工作安装DotNet6环境服务端代码客户端代码编译测试记录失败尝试1(Linux可用端口范围限制)解决Linux端口范围限制查看端口范围修改端口范围失败尝试2(可用端口耗尽......
  • 看雪2023CTF
    文章目录Game-第一关:-.--/---/..-/.--/../-.Game-第二关:二维码Game-第三关:错误的MD5Game-第四关:盲文Game-第五关:看雪的历史Game-第六关:凯撒留下了什么?Game-第七关110米要跨几个栏Web-签到题Web-2023签到Web-[强网杯2019]随便注Web-动不了Web......
  • GBU810-ASEMI高性能整流桥GBU810
    编辑:llGBU810-ASEMI高性能整流桥GBU810型号:GBU810品牌:ASEMI封装:GBU-4恢复时间:>50ns正向电流:8A反向耐压:1000V芯片个数:4引脚数量:4类型:整流桥、功率整流器件特性:功率整流器件、高性能整流桥浪涌电流:200A正向压降:1.10V封装尺寸:如图工作温度:-55°C~150°CGBU810特性超......
  • 2023年高二10月月考21题空间直角坐标系
    ......
  • 2023年最全得软件测试工程师 学习知识架构体系
    一、Python编程入门到精通 二、接口自动化项目实战 三、Web自动化项目实战 四、App自动化项目实战 五、一线大厂简历 六、测试开发DevOps体系 七、常用自动化测试工具 八、JMeter性能测试 只有不断超越自己的勇气,才能让梦想破茧而出......
  • 10月9日总结
    一.今天做了什么上午上工程实训课,老师先讲了许多铁道的知识,接着讲解了如何为高铁办理发车,和如何开高铁。然后我们就实际上手操作。这节课收获很多。下午上java课,写数据库,感觉很难好在最后写出来了,以后有时间写个教程。二.遇到的问题,如何解决无......
  • Redis淘汰策略-231005
    Redis的内存淘汰策略有哪些:noeviction:当内存不足以容纳新写入数据时,新写入操作会报错;allkeys-lru:当内存不足以容纳新写入数据时,在键空间中,移除最近最少使用的key。(这个是最常用的);allkeys-random:当内存不足以容纳新写入数据时,在键空间中,随机移除某个key。设置过期时间的键空间......
  • 【2023年10月10日】STF61_docker_Day01(上午)
     STF61_docker_Day01(上午)1. 什么是docker?docker类似于VMware软件,也能虚拟出来很多的系统,虚拟出来的系统不叫虚拟机,叫容器。docker:linux系统上的虚拟机2. docker和传统虚拟机的区别VM:使用VMware提供虚拟机的运行平台,管理每个VM中操作系统的运行。每个VM都有自己......
  • 10.10算法
    爬楼梯假设你正在爬楼梯。需要n 阶你才能到达楼顶。每次你可以爬1或2个台阶。你有多少种不同的方法可以爬到楼顶呢? 示例1:输入:n=2输出:2解释:有两种方法可以爬到楼顶。1.1阶+1阶2.2阶示例2:输入:n=3输出:3解释:有三种方法可以爬到楼顶。1.1阶+1阶+1阶2.......