seg木 このへんから高速化とかそういう技術が必要になってきます logNは定数。(要出典) 一般的なの RMQ。ある区間内の最大(最小)値、一か所の変更ができる。 int N; int seg[400005]; void init(){ N=1; while(N<=n)N*=2; rep(i,N)seg[i+N]=-inf; rep(i,n)se…

AOJ 1601 Short PhraseのShort Coding - cookies.txt .scr cookiesくんがショートコーディングをしていたのでそれに対抗してみる。C++のやつ。206byte #include<bits/stdc++.h> main(){ for(char s[50],i,j,p,q;p=atoi(gets(s));printf("%d\n",i)){ for(i=0;i<p;s[i++]=strlen(gets(s+i))); for(i=j=0;j<7;) for(p=i++,j=q=1;q>0;) for(q=j++></p;s[i++]=strlen(gets(s+i)));></bits/stdc++.h>…

要約 木DPをすればよいです。以下では、各変数(ノード)はK種類の離散状態をとるものとし、ノード(変数)はN個あるものとします。 あと、以下のコードは、pythonっぽい疑似コードもどきです。 ある変数xが状態aを取る確率p(x=a)を求める モデルが一本鎖のとき …

おわび ガウス過程のところの話が分からなかったのでその手前までです。ベイズファンのみなさんすみません。以下、6.40のとこのコードです。なんか適当にいじって遊んだってください。 #coding:utf-8 import random import math import numpy as np print "i…