nth prime number contain n as their substring within 100000000

82 篇文章 0 订阅
49 篇文章 0 订阅

A small VBA code to get these numbers:

 

  1. Sub Main()
  2.     Dim i&, j&, k&, n&, num&, a(100000000) As Byte, p(1 To 10000000) As Long
  3.         num = 1
  4.      p(num) = 2    'The 1st prime
  5.     n = 100000000    '10^8
  6.     k = Int(Sqr(n)) 'sqrare root of n
  7.     For i = 3 To k Step 2
  8.         If a(i) = 0 Then
  9.             num = num + 1
  10.             p(num) = i
  11.             If InStr(i, num) Then Debug.Print "P(" & num & ")=" & i
  12.             For j = i * i To n Step 2 * i    'Eractosthenes
  13.                 a(j) = 100    'Not prime number
  14.             Next
  15.         End If
  16.     Next
  17.     For i = k + 1 To n Step 2    'List all prime numbers to array p()
  18.         If a(i) = 0 Then
  19.             num = num + 1
  20.             'p(num) = i
  21.             If InStr(i, num) Then Debug.Print "P(" & num & ")=" & i
  22.         End If
  23.     Next
  24. End Sub

 

It returns:

 

P(7)=17
P(6455)=64553
P(6456)=64567
P(6457)=64577
P(6459)=64591
P(6460)=64601
P(6466)=64661
P(9551)=99551
P(303027)=4303027
P(440999)=6440999
P(968819)=14968819
P(5517973)=95517973

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值