最終更新: 2011-03-02T05:24+0900
何も考えずにコーディングしただけ。一瞬 CPUが考え込みます。
generators = [ lambda{ n = 0 lambda{ n+=1; n*(n+1)/2 } }.call, lambda{ n = 0 lambda{ n+=1; n*n } }.call, lambda{ n = 0 lambda{ n+=1; n*(3*n-1)/2 } }.call, lambda{ n = 0 lambda{ n+=1; n*(2*n-1) } }.call, lambda{ n = 0 lambda{ n+=1; n*(5*n-3)/2 } }.call, lambda{ n = 0 lambda{ n+=1; n*(3*n-2) } }.call, ] # 数を準備 d4polynumbers = generators.map{|g| () while (p = g.call) < 1000 a = [p] a.push(p) while (p = g.call) < 10000 a } # 端緒(の集まり) bunch_of_chain = d4polynumbers[d4polynumbers.size-1].map{|p| [[p, d4polynumbers.size-1]] } # 端緒を伸ばすもの extender = lambda{|chain, pool| xx = chain.last.first.to_s[-2,2] ( (0...(pool.size)).to_a - chain.map{|_| _.last } ).map{|i| [i, pool[i]] }.map{|i, nums| nums.find_all{|num| num.to_s[0,2] == xx }.map{|num| chain + [[num, i]] } }.inject(&:+) } # 伸ばしていく (d4polynumbers.size-1).times{ bunch_of_chain = bunch_of_chain.map{|chain| extender[chain, d4polynumbers] }.inject(&:+) } # 輪っか? bunch_of_cyclic_chain = bunch_of_chain.reject{|chain| chain.first.first.to_s[0,2] != chain.last.first.to_s[-2,2] } # 出力 bunch_of_cyclic_chain.each{|chain| puts chain.map{|a,_| a }.join("\t") puts chain.map{|_,b| "P#{b+3}" }.join("\t") puts "sum: #{chain.map{|a,_| a }.inject(&:+)}" }
先は長いのにもう失速してる。「良いもの。悪いもの。: Project Eulerを100問解いてみた」テトレーションとか聞いたこともない単語なんだけど……。
中学生の時に 3^{50}
の一の位は何かという問題が出た。でも Problem 188は何乗したらいいかもわからない。下手の考え休むに似たりっていうけどどうしたもんかなあ。ない知恵を絞るのも悪くないと思うんだけど。