2017年6月9日金曜日

170609

Ruby


The Kac-Wakimoto Conjecture(1)

The 1-2-3 of Modular Forms には、
定理(!)の片割れだけ載っている。
この定理に従い、値を求めてみた。

def A16(n)
  m = 2 * n + 4
  s = 0
  1.step(m, 2){|b|
    (b + 2).step(m, 2){|a|
      1.step(m, 2){|x|
        1.step(m, 2){|y|
          s += a * b * (a * a - b * b) ** 2 if a * x + b * y == m
        }
      }
    }
  }
  s / 192
end

def A36(n)
  m = 2 * n + 9
  s = 0
  1.step(m, 2){|c|
    (c + 2).step(m, 2){|b|
      (b + 2).step(m, 2){|a|
        1.step(m, 2){|x|
          1.step(m, 2){|y|
            1.step(m, 2){|z|
              s += a * b * c * ((a * a - b * b) * (a * a - c * c) * (b * b - c * c)) ** 2 if a * x + b * y + c * z == m
            }
          }
        }
      }
    }
  }
  s / 141557760
end

n = 10
p (0..n).map{|i| A16(i)}
p (0..n).map{|i| A36(i)}

出力結果
[1, 16, 120, 576, 2060, 6048, 15424, 35200, 73518, 143280, 263584]
[1, 36, 630, 7176, 60165, 398412, 2184078, 10255320, 42321942, 156590980, 527649912]

0 件のコメント:

コメントを投稿

注: コメントを投稿できるのは、このブログのメンバーだけです。