首页 > 其他分享 >闲话 6.30 -JL 引理

闲话 6.30 -JL 引理

时间:2024-06-30 17:43:20浏览次数:21  
标签:infty par le frac la eps JL 6.30 引理

参考了 https://spaces.ac.cn/archives/8679/comment-page-1,有一些增删。

JL 引理

首先下面需要应用马尔可夫不等式的另一个形式:

\[\newcommand \E{\mathbb E} P(x\ge a)=P(e^{\lambda x}\ge e^{\lambda a})(\lambda >0)\le \min_{\lambda >0}e^{\lambda a}\E[e^{\lambda x}] \]

单位模引理

单位模引理:对于 \(u\in \R^n\),每一维从 \(N(0,1/n)\) 随机采样,则

\[\newcommand\eps{\epsilon}\newcommand\la{\lambda}\newcommand\par{\parallel} \forall \eps \in (0,1)\\ P(\big|\par u\par^2-1\big|>\eps)\le 2\exp(-\eps^2n/8) \]

绝对值两边大概是差不多的(?),先看看 \(\par u\par ^2-1>\eps\)。

套用上述结论有:

\[P(\par u\par ^2-1>\eps)\le \min_{\la >0}e^{-\la (\eps+1)}\E[e^{\la\par u\par^2}]\\ \]

化简 \(\E[e^{\la\par u\par^2}]\),将 \(u\) 正交分解为 \(u_{1:n}\)。

\[\E[e^{\la\par u\par^2}]=\prod_{i=1}^n \E[e^{\la u_i^2}]=\E[e^{\la u_i^2}]^n\\ =\left(\int_{-\infty}^{\infty}P(u_i=t)e^{\la t^2}dt\right)^n\\ =\left(\int _{-\infty}^{\infty}\frac{1}{\sqrt{2n\pi}}e^{-\frac{nt^2}{2}}e^{\la t^2}dt\right)^n=\left(\int _{-\infty}^{\infty}\frac{1}{\sqrt{2n\pi}}e^{-(\frac n2-\la)\frac{t^2}2}dt\right)^n \]

考虑正态分布 \(N(0,\sigma^2)\) 满足 \(\frac{1}{2\sigma^2}=\frac{n-2\la}{2}\),即 \(\sigma=\sqrt{\dfrac{1}{n-2\la}}\)。比较

\[\int_{-\infty}^{\infty}\frac{1}{\sqrt{2\pi}}e^{-\frac{x^2}{2\sigma^2}}dt=1 \]

上式的系数,得到:

\[\E[e^{\la u_i^2}]=\int _{-\infty}^{\infty}\frac{1}{\sqrt{2n\pi}}e^{-(\frac n2-\la)\frac{t^2}2}dt=\frac{\sqrt{n}}{\sqrt{n-2\la}}=\sqrt{\frac{n}{n-2\la}} \]

那么

\[P(\par u\par ^2-1>\eps)\le \min_{\la >0}e^{-\la (\eps+1)}\left(\frac n{n-2\la}\right)^{n/2}\\ \]

设 \(f(\la)=e^{-\la (\eps+1)}\left(\dfrac n{n-2\la}\right)^{n/2}\)。则

\[f'(\la)=e^{-\la(\eps+1)}\left(\dfrac n{n-2\la}\right)^{n/2-1}\left(\frac{n^2}{(n-2\la)^2}-\frac{n(\eps+1)}{n-2\la}{}\right) \]

所以不妨取

\[\la=\frac{n\eps}{2(\eps+1)} \]

所以:

\[P(\par u\par ^2-1>\eps)\le e^{-n\eps/2}e^{n\ln(\eps+1)/2}=e^{n(\ln(1+\eps)-\eps)/2}\\ \]

而设 \(g(\eps)=\ln(\eps+1)-\eps+\eps^2/4\)。

\[g'(\eps)=\frac1{1+\eps}-1+\frac\eps2\le 0(0<\eps<1) \]

所以 \(g(\eps)\le g(0)=0\),所以:

\[P(\par u\par ^2-1>\eps)\le \exp(-n\eps^2/8) \]

JL 引理

对于 \(N\) 个向量 \(v_{1:N}\in \R^m,n>\frac{24\ln N}{\eps^2}\),随机矩阵 \(A\in R^{n\times m}\) 采样于 \(N(0,1/n)\),\(\eps\in (0,1)\),则至少有 \(\frac{N-1}{N}\) 的概率满足:

\[\forall i\neq j,\\ (1-\eps)\par v_i-v_j\par ^2\le \par Av_i-Av_j\par ^2 \le(1+\eps)\par v_i-v_j\par \]

证明:

容易发现,\(\forall u\in \R^m\),\((Au)_i\) 服从 \(N(0,1/n)\)。

使用 Union Bound,带入 \(u=\frac{v_i-v_j}{\par v_i-v_j\par}\)。那么:

\[P\left(\exists i\neq j,\left| \frac{A(v_i-v_j)}{\par v_i-v_j\par}-1\right|\ge \eps\right)\le 2\binom N2\exp(-\eps^2n/8) \]

若 \(n>\dfrac{24\ln N}{\eps^2}\),有:

\[1-N(N-1)\exp(-\eps^2n/8)>1-\frac{N(N-1)}{N^3}>\frac{N-1}{N} \]

结束。

但是在计算中,\(24\ln N/\eps^2\) 并不是非常小的值,如果不使用一些更加高级的降维算法的话,差不多只有在 \(\eps>0.5\) 的时候才能发挥作用了(悲)。

标签:infty,par,le,frac,la,eps,JL,6.30,引理
From: https://www.cnblogs.com/british-union/p/18276712/JL

相关文章

  • 【抽代复习笔记】21-群(十五):循环群引理及定义
    例4:证明,如果σ=(i1i2…ik)是Sn中的一个k-循环,而r∈Sn,则rσr^(-1)也是一个k-循环,且rσr^(-1)=(r(i1),r(i2),…,r(ik))。证:①设σ=(i1i2…ik)=(i1ik)(i1ik-1)…(i1i2),则rσr^(-1)=r(i1i2…ik)r^(-1)=r(i1ik)(i1ik-1)…(i1i2)r^(-1)=r(i1ik)[r^(-1)r](i1ik-1)[......
  • 老玩家BJL百家三珠路打法及技巧----实战技巧总结篇
    更多技巧可移步围脖—老晨谈赌​​三珠路打法由来已久,因其简单实用,很多人都会用到。什么叫三珠路,就是把大路庄闲,按三个一组进行划分,找一个路单图片,按水平和垂直排列。把牌路按照三珠路排列后,样式如下图:每一列为一个图形,会出现8种图形:如果按照前面的统计数据,那么以下......
  • 老晨谈赌详解AG百家和BJL下三路的实战技巧打法个人经验
    更多技巧可移步围脖【老晨谈赌】技术可以通过学习来获得,经验可以通过实战来得到,心态可以通过调节来增强。每一个人都不是生来都无比强大的,我也是如此,也是通过无数个黑夜的煎熬最后才研究出来的,所以如果说幸运,我们都幸运,如果说不幸运,我们其实都一样。可以不设置止盈点,但是你一......
  • P3267 JLOI2016/SHOI2016 侦察守卫
    P3267JLOI2016/SHOI2016侦察守卫互相赋值的双dp思路设\(f[u][i]\)表示包括\(u\)子树内所有关键点都被覆盖(包括\(u\)),且至少还可以向\(u\)的父亲方向覆盖\(i\)层的最小代价。设\(g[u][i]\)表示向下距离大于等于\(i\)的点全部被覆盖,剩下距离小于\(i\)的点随意......
  • [JLU]校园网上网攻略汇总与补充
    前言如题,陆爻齐为了汇集一下觉得比较有用的JLU校园网相关的资源,同时对于一些比较重要但比较少被提及的地方做点补充而写本文。希望能对吉林大学(也许)的各位有所帮助如果陆爻齐再次重装系统,或许也用的到罢:)正文现有攻略推介文章官方老文如果你会搜索下“吉林大学校园网”,......
  • P4568 [JLOI2011] 飞行路线
    题目P4568[JLOI2011]飞行路线要求找到在最多可以免费乘坐k条航线的情况下,从城市s到城市t的最少花费。这是一个典型的分层图问题。分层图的建模1.建立层次将原图分成k+1层,表示在0到k次免费乘坐的情况下的状态。第i层表示已经使用了i次免费乘坐机会的状态。2.建立节点和边......
  • OpenCloudOS Kernel SIG 月度动态:发布 OCK 6.6.30-4 版本,新增特性支持
    作为OpenCloudOS核心SIG之一,KernelSIG负责OpenCloudOS内核的路线规划、方案讨论、技术能力输出,为OpenCloudOS提供更加稳定、可靠的内核基座。一、整体进展1、发布OCK6.6.30-4版本,合入Intel、龙芯、Phytium、兆芯等特性支持;2、OpenCloudOSKernel文档内容更新,实......
  • LGV引理
    在一张有向无环图DAG中,有边权,给定起点点集A,终点点集B,且A,B中的点数一致。定义P表示DAG中的一条路径。定义w(P)表示路径P上的边权乘积。定义e(a,b)表示a到b的所有路径的边权乘积之和,即\(e(a,b)=\sum_{P_i\in(a\tob)}w(P_i)\)定义一组A到B的不相交路......
  • LGV 引理学习笔记
    \(\text{LGV}\)引理学习笔记\(\text{LGV}\)引理一般用于求解有向无环图中多条不相交路径的方案数,引理内容如下。引理定义\(w(P)\)指的是路径\(P\)上所有边权的乘积(在路径计数问题中认为所有边权均为\(1\)即可),\(e(A,B)\)指的是\(A\toB\)的所有路径的\(w\)和。对......
  • P3270 [JLOI2016] 成绩比较
    记\(a_i=N-R_i,n=N-1\)。先不考虑有多少人被碾压,计算出符合排名限制的方案数。枚举每门课和B神在每门课中的得分,选出\(a_i\)个人得分小于等于他,即:\[\prod\limits_{i=1}^m\dbinom{n}{a_i}\sum\limits_{j=1}^{U_i}j^{a_i}(U_i-j)^{n-a_i}\]设\(s(x)=\sum\limits_{j=1}^{......