This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
def s(k, m, n)
s = 0
(1..n).each{|i| s += i if n % i == 0 && i % k == m}
s
end
def A007325(n)
ary = [1]
a = [0] + (1..n).map{|i| s(5, 1, i) + s(5, 4, i) - s(5, 2, i) - s(5, 3, i)}
(1..n).each{|i| ary << (1..i).inject(0){|s, j| s - a[j] * ary[-j]} / i}
ary
end
def mul(f_ary, b_ary, m)
s1, s2 = f_ary.size, b_ary.size
ary = Array.new(s1 + s2 - 1, 0)
(0..s1 - 1).each{|i|
(0..s2 - 1).each{|j|
ary[i + j] += f_ary[i] * b_ary[j]
}
}
ary[0..m]
end
def A285441(n)
ary1 = A007325(n)
ary2 = Array.new(n + 1, 0)
(0..n / 2).each{|i| ary2[i * 2] = ary1[i]}
ary = [-1] + mul(ary1, mul(ary2, ary2, n), n)[0..-2]
mul(ary2, (0..n).map{|i| -ary[i]}, n)
end
p A285441(100)
Comments