matlab有序生成整数,matlab – 有效地生成唯一的整数对

悠悠 2023-01-15 08:08 59阅读 0赞

容易,吝啬,以正确的方式观看.

您希望生成n对整数\[p,q\],使得p和q位于间隔\[1,m\]中,p

有多少对可能？总数只是m \*(m-1)/ 2. (即,从1到m-1的数字之和)

所以我们可以在\[1,m \*(m-1)/ 2\]的范围内生成n个随机整数.兰德尔做得很好. (旧的matlab版本不允许第二个参数为randperm.)

k = randperm(m/2\*(m-1),n);

(注意,我用一个有趣的方式写了这个表达式,除以2可能是一个奇怪的地方,这避免了在上限附近的某些值的精度问题.)

现在,如果我们将每个可能的对\[p,q\]与k中的一个整数相关联,我们可以从k中生成的整数向后工作到一对\[p,q\].因此,该列表中的前几个对是：

\{\[1,2\], \[1,3\], \[2,3\], \[1,4\], \[2,4\], \[3,4\], ..., \[m-1,m\]\}

我们可以把它们看作是m的严格上三角阵列中的元素,这些元素在主对角线之上.

q = floor(sqrt(8\*(k-1) + 1)/2 + 1/2);

p = k - q.\*(q-1)/2;

看到这些公式从k中展开的元素中恢复p和q.我们可以说服自己确实有效,但也许这里只是一个简单的方法：

k = 1:21;

q = floor(sqrt(8\*(k-1) + 1)/2 + 3/2);

p = k - (q-1).\*(q-2)/2;

\[k;p;q\]'

ans =

1 1 2

2 1 3

3 2 3

4 1 4

5 2 4

6 3 4

7 1 5

8 2 5

9 3 5

10 4 5

11 1 6

12 2 6

13 3 6

14 4 6

15 5 6

16 1 7

17 2 7

18 3 7

19 4 7

20 5 7

21 6 7

测试它的另一种方式是显示所有的对都是为小案例生成的.

m = 5;

n = 10;

k = randperm(m/2\*(m-1),n);

q = floor(sqrt(8\*(k-1) + 1)/2 + 3/2);

p = k - (q-1).\*(q-2)/2;

sortrows(\[p;q\]',\[2 1\])

ans =

1 2

1 3

2 3

1 4

2 4

3 4

1 5

2 5

3 5

4 5

嗯,看起来一切都很完美.现在尝试一些大量的m和n来测试所用的时间.

tic

m = 1e6;

n = 100000;

k = randperm(m/2\*(m-1),n);

q = floor(sqrt(8\*(k-1) + 1)/2 + 3/2);

p = k - (q-1).\*(q-2)/2;

toc

Elapsed time is 0.014689 seconds.

该方案对于m大致为1e8,在由于双精度的精度误差而失败之前将起作用.在m / 2 \*(m-1)超过2 ^ 53之前,精确极限应该不大于134217728.一个很好的功能是,不需要对重复对进行拒绝.

发表评论取消回复

表情：

评论列表（有 0 条评论，59人围观）

还没有评论，来说两句吧...

相关阅读

相关 matlab代码生成150个1到250的随机整数

您可以使用以下代码生成150个在1到250之间的随机整数: rng(0); % 设置随机数生成器种子，保证每次生成的随机数相同 rand_ints = ran

男娘i/ 2024年03月25日 17:25/ 0 赞/ 41 阅读

相关 Matlab生成遮蔽Rician随机数

Matlab生成遮蔽Rician随机数在通信和图像处理等领域中，Rician分布是一种常见的概率分布模型。在一些实际应用场景中，可能需要使用Rician随机数来模拟某些物理

迈不过友情╰/ 2024年03月17日 14:00/ 0 赞/ 84 阅读

相关 GetTimeNumber()当前时间生成的整数

' ============================================ ' GetTimeNumber()当前时间生成的整数 ' ========

迈不过友情╰/ 2024年02月17日 22:35/ 0 赞/ 25 阅读

相关材料生成算法 Matlab

材料生成算法 Matlab 近年来，材料生成算法在材料科学领域引起了广泛的关注。这些算法利用计算机模拟和机器学习技术，能够预测和生成具有特定性质的新材料。在本文中，我们将介绍

ゞ浴缸里的玫瑰/ 2023年10月15日 23:34/ 0 赞/ 45 阅读

相关 Matlab 整数线性规划问题模型代码

![1489774-20190917102243668-148038964.png][] 整数线性规划问题的基本内容 > 整数线性规划解决的是自变量在一定的线性约束条件

偏执的太偏执、/ 2023年08月17日 16:12/ 0 赞/ 152 阅读

相关 matlab有序生成整数,matlab – 有效地生成唯一的整数对

容易,吝啬,以正确的方式观看. 您希望生成n对整数\[p,q\],使得p和q位于间隔\[1,m\]中,p 有多少对可能？总数只是m \(m-1)/ 2. (即,从1到m-1

悠悠/ 2023年01月15日 08:08/ 0 赞/ 60 阅读

相关 c 调用matlab生成的dll,vs2008中调用matlab生成的dll

matlab版本为:matlab7.11.0(R2010b) vs2008 1.在matlab中生成dll 在m文件中实现如下功能：获得N个魔方矩阵的秩。m文件的内容

妖狐艹你老母/ 2023年01月15日 02:28/ 0 赞/ 217 阅读

相关 Matlab 如何生成一个[a,b]范围内随机整数的2种方法【已经解决】

如何使用MATLAB生成一个\[a,b\]范围内的随机整数。比如：随机生成\[9,13\]范围内的一个(或多个)整数~ 首先感谢@slandarer的指正，现已经更改。

灰太狼/ 2022年09月03日 10:25/ 0 赞/ 27 阅读

相关生成伪随机整数

生成伪随机整数在头文件<stdlib.h>中，声明了函数rand(); <span style="font-family:SimSun;">const int r

水深无声/ 2022年08月04日 10:46/ 0 赞/ 281 阅读

相关 matlab 保留固定长度的整数位

z为某一位整数,需要保留2位整数需要输出固定长度的整数位是: fprintf('%02d',z) 输出是0z 如果在字符串中需要用到是 sprintf('%0

深碍√TFBOYSˉ_/ 2022年03月31日 10:54/ 0 赞/ 637 阅读