基底関数をあつめているBasis Set Exchangeが微妙に便利でにわかクオンタムケミストの僕などは軽く重宝している。

あと、SOAPのサービスがついていて、これは！と思ったんだけど、メソッドがややこしげで、結局どう使えばいいのかわからんまま今に至る。

僕とU隊長の10年来の共通の友人で、脈がないのに脈ありと勝手に盛り上がって妄想が最高潮に達しては自爆するという事を繰り返す子がおるんやけど、最近また豪快に自爆したらしい。

なんで、学習せーへんのかなー？いい年やのになー

なんて会話を毎年のようにしているんだけど、突然思いつきましたよ。 レコーディングダイエットのように脈ありと感じたことを日記につけていけばいいんじゃないかと。

で、それをU隊長に語ったら、

そんなんblogでヒトに見せな効果ないやろ

と返された。なるほど。でも、それだったら小町でええやん。

processing周りをさがしていて偶然見つけた。

そういえばこれ数年前にメタモで見たわ。

javascript周りをjQueryに変えようとごそごそいじってたのだけど、途中でレイアウトいじりだしてしまった。段落の始めにインデントが入ってないのが前から気になっていたのだけどそれもついでになおした。

勝間和代さんの本に論理パズルがとかあったので、ちょっといくつか買っておいたのだ。

論理パズル101―推理の楽しさ、ひらめきの快感 (ブルーバックス)
小野田 博一
講談社 / ￥ 987 ()
通常24時間以内に発送

史上最強の論理パズル―ポイントを見抜く力を養う60問 (ブルーバックス)
小野田 博一
講談社 / ￥ 788 ()
通常24時間以内に発送

101問とかいうやつは紙と鉛筆使ってグリッドぐりぐり書くタイプっぽいので60問のほうがお手軽。まぁ矛盾をつきながらあーだこーだと考えを巡らせるのは楽しい。そういえば小学生の頃はこういうの夢中になってやったよな（頭の体操だっけな）とフラッシュバック。

先週の大阪出張の行き帰りの新幹線の中で30問くらいは解いたので、そのままの勢いで60問全て終了。

論理的な思考力ってのはタスクが複雑に絡み合った時などに特に必要とされますな。普段は実験とかさくっとこなして、言うてる事もまともなんだけど、ちょっと複雑になるとパンクしてしまうヒトっていますよな。でも、社歴が長いからといって論理力が高まるかっていうとそういうもんでもないんだよな（英語と一緒で）。

感覚的には理解している事だけど、本読んで頭の中が整理されると何をやるべきやらないべきが自分の中で固まるのでこういうビジネス本を読むのもいいことだなぁと思った。

そういえば、「この酒はおとこざけですから〜」とか言われて?となったので、帰ってきて調べたのだけど、硬水を使って短期発酵を強めたとかまぁそんな程度の説明しか見つからなかった。

もうちっとプロセスを詳しく！

とか思ったので、後で図書館にでもいって調べてくるかと思った次第。

何人か紹介してもらい、色々面白い話を聞けたが、時間が足りないなぁ。

• useR!2008の話

  • 日本人の参加数は10人いかない位
  • bioinformaticsな感じ?
  • 内容くらいは押さえておかないといけないかも

• RのGUIないの？っていわれたので、R Commandarの話をしといた。僕は使った事ないけども

• SOAPめんどい?

  • pylons-soaplibで超ラクチンとか言ってみたけど伝わらなかったかも。後でなんか書くかも

• Knime

  • Tavernaもありますよねーと言いたかったが、名前が出てこなかった。あとUI的にはYahoo! Pipesのほうが洗練されてていい感じ。

• Plutoそろそろちゃんとやる。

R Commanderハンドブック―A Basic-Statistics GUI for R
舟尾 暢男
九天社 / ￥ 3,360 ()

連続する素数リスト和が素数になるもののうち、100万までのなかからその並びが最も大きいものを探す。

2,3,5...と和をとって100万までの素数の最大値を超えたところが並びの最大。ここからwindowをずらして和をもとめ素数が見つかったらそれが最大。見つからなかったらwindowサイズを一つ小さくしながら同じ事を繰り返す。

def make_primes(n):
    primes = [2]
    nums = range(3,n+1,2)
    while 1:
        prime = nums[0]
        if prime ** 2 > n:
            primes += nums
            return primes
        else:
            primes.append(prime)
            nums = [x for x in nums if x%prime != 0]
    return primes

def c_sum(i,n):
    return reduce(lambda x,y:x+y,primes[i:i+n+1])

def check_max():
    for i in range(1,len(primes)):
        if c_sum(0,i) > primes[-1]:
            return len(primes[:i])

if __name__ == "__main__":
    primes = make_primes(1000000)
    print "primes are %d" % len(primes)

    cm = check_max()
    print "max windowsize %d" % cm

    max = 0
    pl = cm
    for x in range(1,pl):
        window = pl - x
        print "checking length %d..." % window
        for j in reversed(range(x)):
            num = c_sum(j,window)
            if num in primes:
                max = num
                print num
                print primes[j:j+window+1]
                if max > 0: break
        if max > 0: break

風邪ひいてしんどいので、早退。

メタミドホスの基準値って結構ばらついてますな。

だから基準値の何倍とかいう議論ってあれだなぁとか思うのと、最近よくある、モラールごときでヒトが死ぬ状況が作られやすい社会ってなんだかなぁと。

生体系を考えても負のフィードバックが働くからホメオスタシスが保たれるのであって、正のフィードバックなんて一部の危機的状況で働くから有効なんじゃないの？

常に正のフィードバックが働きやすいいまの社会系って異常だよなぁ、強い揺り戻しがどこかでおこるのかなぁとふとよぎった。

problem 63

n乗したらn桁がなるような数字はいくつあるか?

10以上は考えなくてよいのだけど、一応100まで数えてみた。

count = 0
for i in range(1,100):
    if len(str(i)) > 1: break

    j = 1
    while(1):
        n = i ** j
        if len(str(n)) < j: break
        # print n,i,j
        j += 1
        count += 1

print count

problem 99

対数変換して比較するだけなのでlog10で。

import math
f = open("base_exp.txt")

n    = 0
max  = 0
maxn = 0
for l in f:
    n += 1
    a,b = l[:-2].split(',')
    x = math.log(int(a)) * int(b)
    if max < x:
        max = x
        maxn = n
print maxn,max

これで、level2の六面体に昇進