-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathprog.sf
80 lines (61 loc) · 1.42 KB
/
prog.sf
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
#!/usr/bin/ruby
# a(n) is the smallest number which can be represented as the product of n distinct integers > 1 in exactly n ways.
# https://oeis.org/A360590
# Known terms:
# 2, 12, 60, 420, 3456, 60060
# New terms found:
# a(7) = 155520
func almost_prime_numbers(n, k, primes, callback, squarefree = false) {
var sqf = (squarefree ? 1 : 0)
var factors = primes.sort.uniq
var factors_end = factors.end
if (k == 0) {
callback(1)
return nil
}
func f(m, k, i=0) {
if (k == 1) {
var L = idiv(n,m)
for j in (i..factors_end) {
with (factors[j]) {|q|
q > L && break
callback(m*q)
}
}
return nil
}
var L = idiv(n,m).iroot(k)
for j in (i..factors_end) {
with (factors[j]) { |q|
q > L && break
__FUNC__(m*q, k-1, j + sqf)
}
}
}(1, k)
f = nil # to prevent memory leak due to circular references
return nil
}
func isok(k, n) {
var count = 0
almost_prime_numbers(k, n, k.divisors.grep { _ > 1 }, {|t|
++count if (t == k)
}, true)
count == n
}
func a(n) {
for k in (2..Inf) {
if (isok(k, n)) {
return k
}
}
}
for n in (2..100) {
say [n, a(n)]
}
__END__
[2, 12]
[3, 60]
[4, 420]
[5, 3456]
[6, 60060]
[7, 155520]