A387522 Index of first term in A386482 that is divisible by the n-th prime.
2, 4, 8, 10, 23, 21, 28, 30, 37, 47, 56, 67, 63, 61, 75, 94, 88, 86, 80, 119, 117, 135, 131, 174, 166, 162, 160, 156, 154, 150, 200, 235, 229, 227, 217, 215, 209, 270, 266, 260, 254, 252, 242, 240, 297, 295, 314, 354, 350, 348, 344, 338, 336, 326, 428, 422, 416, 414, 408, 404, 402, 392, 378, 374, 372, 478, 464, 458, 608, 606, 602
Offset: 1
Links
- Michael De Vlieger, Table of n, a(n) for n = 1..65536 (first 841 terms from N. J. A. Sloane)
Programs
-
Mathematica
Block[{p, s}, p[A386482/b386482.txt%22,%20%22Data%22%5D%5B%5BAll,%20-1%5D%5D;%20Do%5BMap%5BIf%5Bp%5B%23%5D%20==%200,%20Set%5Bp%5B%23%5D,%20n%5D%5D%20&,%20FactorInteger%5Bs%5B%5Bn%5D%5D%5D%5B%5B;;%20,%201%5D%5D%5D,%20%7Bn,%20Length%5Bs%5D%7D%5D;%20TakeWhile%5BArray%5Bp%5BPrime%5B%23%5D%5D%20&,%20120%5D,%20%23%20%3E%200%20&%5D%20%5D%20(*%20_Michael%20De%20Vlieger">] := 0; s = Import["https://oeis.org/A386482/b386482.txt", "Data"][[All, -1]]; Do[Map[If[p[#] == 0, Set[p[#], n]] &, FactorInteger[s[[n]]][[;; , 1]]], {n, Length[s]}]; TakeWhile[Array[p[Prime[#]] &, 120], # > 0 &] ] (* _Michael De Vlieger, Sep 03 2025 *)