Alias采样算法
2021/8/31 1:06:19
本文主要是介绍Alias采样算法,对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!
Alias采样过程分析
学习来源:
1、时间复杂度O(1)的离散采样算法——Alias method
2、Alias采样算法
程序实现
import numpy as np def alias_setup(probs): K = len(probs) q = np.zeros(K) J = np.zeros(K, dtype=np.int) smaller = [] larger = [] for kk, prob in enumerate(probs): q[kk] = K * prob # 保存面积块 if q[kk] < 1.0: smaller.append(kk) # 面积小于1的元素下标 else: larger.append(kk) # 面积大于1的元素下标 while len(smaller) > 0 and len(larger) > 0: small = smaller.pop() large = larger.pop() # J 记录的是填充这一列的下标 J[small] = large # q 记录的是原来事件的占比 q[large] = q[large] - (1 - q[small]) if q[large] < 1.0: smaller.append(large) else: larger.append(large) return J, q def alias_draw(J, q): K = len(J) kk = int(np.floor(np.random.rand() * K)) # 随机选一列 if np.random.rand() < q[kk]: return kk else: return J[kk] if __name__ == "__main__": J, q = alias_setup([1/2, 1/3, 1/12, 1/12]) print(J, q) # [0 0 0 1] [1. 0.66666667 0.33333333 0.33333333] ret = alias_draw(J, q)
构造表的时间复杂度:O(n)
采样的时间复杂度:O(1)
这篇关于Alias采样算法的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!
- 2024-05-15鸿蒙生态设备数量超8亿台
- 2024-05-13TiDB + ES:转转业财系统亿级数据存储优化实践
- 2024-05-09“2024鸿蒙零基础快速实战-仿抖音App开发(ArkTS版)”实战课程已上线
- 2024-05-09聊聊如何通过arthas-tunnel-server来远程管理所有需要arthas监控的应用
- 2024-05-09log4j2这么配就对了
- 2024-05-09nginx修改Content-Type
- 2024-05-09Redis多数据源,看这篇就够了
- 2024-05-09Google Chrome驱动程序 124.0.6367.62(正式版本)去哪下载?
- 2024-05-09有没有大佬知道这种数据应该怎么抓取呀?
- 2024-05-09这种运行结果里的10.100000001,怎么能最快改成10.1?