読者です 読者をやめる 読者になる 読者になる

Tamflexの貯蔵庫

やる気のない備忘録

aoj DPL_1_A

Coin Changing Problem | Aizu Online Judge コイン問題 貪欲法では解けないのでDPを使う. dp[i][j] := (i枚目までのコインを使った時の金額jの最小実現枚数) とすれば計算量O(nm)で解ける #include <bits/stdc++.h> #define REP(i,n) for(int i=0;i<n;i++) #define FOR(i,a,b) for(int i=a;i<=b;i++) #define INF 1 << 30 using namespace std; int main() { int n, m; cin >> n …</n;i++)></bits/stdc++.h>

aoj ALDS1_14_D

Multiple String Matching | Aizu Online Judge 文字列Sの接尾辞配列を構築すると文字列Tの検索が二分探索で行えるためO(|T|log|S|)のオーダーで行える. 接尾辞配列の構築にかかるコストはO(|S|(log|S|)^2)なので何度もTの検索を行う場合有用である. #includ…

C++でのBM法の実装

文字列検索のアルゴリズムの中で比較的高速なものとして有名なBM法 文字列検索を行った時に役に立ったので書いたソースを載せるなお以下の文献を参考にさせていただきました. http://www.comp.cs.gunma-u.ac.jp/~koichi/ALG2/9beamerBM.pdf #include <bits/stdc++.h> using </bits/stdc++.h>…

aoj ALDS1_14_C

Pattern Search | Aizu Online Judge二次元配列の効率的なハッシュの精製方法について配列 から以下のような行方向のハッシュの配列 を生成する. この時 のように計算すればO(h)のオーダーでハッシュが計算できる.同様に以下のような列方向のハッシュの配列 …

aoj ALDS1_14_B

String Search | Aizu Online Judge 逐一比較だとTLEしてしまったがsubstrなるものを使うと通る. アルゴリズム的には同じオーダなのだけど内部処理で何かが違うのだろう. std::string::substr - C++入門 #include <bits/stdc++.h> #define REP(i,n) for(int i=0;i<n;i++) #define FOR(i,a,b) for(int i=a;i<=b;i++) using namespace std; int main() { string t,p; cin>>t>>p; int</n;i++)></bits/stdc++.h>…

aoj ALDS1_14_A

Naive String Search | Aizu Online Judge やるだけ a,b=$<.map{|s|s.chomp};(c=a.size).times{|v|a[v..c-1]=~/^#{b}/&&p(v)}

aoj 0099

ワカサギ釣り | Aizu Online Judge priority_queueは最大値を一番上に持ってくるSTLでこれを使えばとても便利. データごとにpqに挿入し現在の魚の数を集計. 現在の状態と一致するものの中で最大のものを出力し残りはpopする. こうすることによって古いデータ…

aoj 0097

整数の和 | Aizu Online Judge 動的計画法を用いてやる. dp[i][j] := (i個の異なる0〜100の数を用いて和jを実現する場合の数) あとは0〜100まで順に更新する. ここでポイントはdp[i][j]の更新の順序で必ずiについて逆順で更新しなければならない. なぜならば…

aoj 0096

4つの整数の和 | Aizu Online Judge 動的計画法を使ってやる dp[i][j] = (i+1個の数を使って和jを実現する場合の数) あとは適当に1000までという条件を組み込めばOK #include <bits/stdc++.h> #define REP(i,n) for(int i=0;i</bits/stdc++.h>

aoj 0095

ワカサギ釣り | Aizu Online Judge ホントは$> gets;x,y=$<.max_by{|s|eval s.split*'*-1+101*'}.split;puts x+' '+y

aoj 0094

土地面積 | Aizu Online Judge 一文字のリテラル/文字列リテラルは?をつけることで実現できる(ex 'd'=?d) リテラル (Ruby 2.0.0) evalをつければ簡単に計算できる. p eval(gets.split*?*)/3.30578

aoj 0092

正方形探索 | Aizu Online Judge 動的計画法の問題. dp(i,j)を(i,j)を右下の頂点とするの最大サイズと定義し以下のような漸化式を計算すれば良い. dp(i,j) = min(dp(i-1,j),dp(i-1,j-1),dp(i,j-1))+1 あとは最大値を求めればよい. int main() { string s; in…

aoj 0091

にじみ | Aizu Online Judge dfsを使うところまでは良かったが探索の仕方がまずかった. まず最初に試したのは全ての探索木の節/葉において全部(8*8)を探索する方法でこれだとTLEしてしまった. 自力では解けなかったので先人達の答えを拝見させていただいた. …

aoj 0090

シールの重なり | Aizu Online Judge 二円の交点がどのくらい重なっているか数えれば良い. ついでに座標を扱うクラスを拡張した. class P { public: double x,y; P(){};P(double x,double y):x(x),y(y){}; P operator+(const P&q){P t;t.x=x+q.x;t.y=y+q.y;r…

aoj 0089

菱形上の最短経路 | Aizu Online Judge 添字の扱いが厄介 vector<string> split(const string &str, char delim) { istringstream iss(str); string tmp; vector<string> res; while(getline(iss, tmp, delim)) res.push_back(tmp); return res; } int main() { vector<int> v; st</int></string></string>…

aoj 0087

逆ポーランド記法 | Aizu Online Judge やるだけ vector<string> split(const string &str, char delim) { istringstream iss(str); string tmp; vector<string> res; while(getline(iss, tmp, delim)) res.push_back(tmp); return res; } int main() { string s; while(getl</string></string>…

aoj 0086

パトロール | Aizu Online Judge 一筆書き可能か判定する問題, すなわちグラフGが与えられた時そのグラフがオイラー路であるか判定すればよい. オイラー路となる条件は, すべての頂点の次数が偶数⇔オイラー閉路(オイラーグラフ) 奇数次数の頂点が2つで残りが…

codingame - Skynet: the Virus

Ambush Finish the 4th test of the "Skynet: the virus" puzzle with 50 or more remaining links and get a 100% score. (訳) 待ち伏せ "Skynet: the virus"の4番目のテストケースを50辺以上残して完了し,100%のスコアを獲得せよ. というAcheivementがあ…

codingame - Horse-racing Duals

bashをプログラミング言語と読んでいいのかよくわからないけどbash縛りのAchivementがあったので解いてみた. read N;s=($(for((i=0;i

codingame - Mayan Calculation

連想配列を使えば楽に解ける long long を使うことに注意 #include <bits/stdc++.h> #define REP(i,n) for(int i=0;i<n;i++) using namespace std; typedef long long ll; int main() { string s; ll l,h,ans,x[2]; cin>>l>>h; vector<string> a(20,""); map<string,ll> b; REP(i,l) { getline(cin,s); if(s==""){i--;continue;} REP(j,s.size())…</string,ll></string></n;i++)></bits/stdc++.h>

codingame - The Gift

貪欲的に小さい順にとけば良い. #include <bits/stdc++.h> #define REP(i,n) for(int i=0;i<n;i++) using namespace std; typedef long long ll; int main() { ll n,c; cin>>n>>c; vector<ll> a(n,0); REP(i,n) cin >> a[i]; string ans = ""; sort(a.begin(),a.end()); REP(i,n) { int t=(a[i]</ll></n;i++)></bits/stdc++.h>

codingame - Scrabble

割と面倒だった(225byte) y='';$><<(a=$<.to_a.map{|s|s.chomp.split y})[0,n=a.size-1].select{|v|a[n].inject(v*y){|s,b|s=s.sub(b,y)}==y}.max_by{|s|eval({qz:10,jx:8,k:5,fhvwy:4,bcmp:3,dg:2,'a-z'=>1}.inject(s*'+'){|a,v|a=a.gsub(/[#{v[0]}]/,v[1].…

codingame - Bender, a depressed robot

ループ判定が面倒 自分は同じ所を同じ向きで通った場合にループ判定を出すようにした. しかしこれだけでは最後のテストケースが落ちるので20回以上同じ所を同じ向きで通った場合判定を出すようにした. struct P { int x,y; P(){}; P(int x, int y): x(x), y(…

codingame - Dwarfs standing on the shoulders of giants

メモ再帰によって簡単に解ける template<typename T> struct Edge { int dst; T w; Edge() {}; Edge(int dst): dst(dst) {w=1;}; Edge(int dst, T w): dst(dst), w(w) {}; bool operator < (const Edge &e){return w != e.w ? w > e.w : dst < e.dst;} }; typedef vector<Edge<ll></edge<ll></typename>…

codingame - Telephone Numbers

例によってrubyで挑戦してみた(73byte) gets;p$<.inject([]){|a,s|[*1..s.size-1].each{|v|a<

codingame - Conway Sequence

気合の131byte l,r=$<.read.split;$><<[*2..r.to_i].inject([l.to_i]){|s|y=0;n=s[0];t=[];(s+['']).each{|v|(v==n)?(y+=1):(t+=[y,n];n=v;y=1)};s=t;}*' '

codingame - Network Cabling

x軸またはy軸に平行なメインケーブルから各家庭に枝を伸ばせば良い. 中央値は絶対偏差の和を最小にすることが知られているので, メインケーブルのy座標/x座標は中央値を用いると良さそう. 証明はこちらなど↓ 連続関数に関してはこちら↓ http://web.uvic.ca/~…

codingame - Teads Sponsored Contest

グラフの直径を求めれば良い 以下が参考になる Spaghetti Source - 木の直径前に解いたabc019のdouble sweepのアルゴリズムを思い出した. 計算量はO(E) この手の問題は多いのでそのうちグラフをクラス化して保存しておきたい. 直径,最小全域木等をg.diameter…

codingame - Indiana - Level 1

やるだけ int main() { int W; // number of columns. int H; // number of rows. cin >> W >> H; cin.ignore(); vector<vector<int>> T(W,vector<int>(H,0)); REP(i,H)REP(j,W) cin >> T[i][j]; int EX; // the coordinate along the X axis of the exit (not useful for thi</int></vector<int>…

codingame - The Paranoid Android

そこまで難しくはないが例外処理に手間取った. struct Info { int floor; int pos; Info(){}; Info(int f,int p):floor(f),pos(p){}; bool operator==(const Info &a) { return floor==a.floor && pos==a.pos; } }; const int SPACE = 0; const int BLOCK = …

codingame - Stock Exchange Losses

int main() { ll n; cin >> n; ll ans = 0, m = LLINF; vector<ll> a(n); REP(i,n) cin >> a[i]; DOWN(i,n-1,0) { m = min(a[i],m); ans = min(m-a[i],ans); } cout << ans << endl; return 0; }</ll>

aoj 0085

ヨセフのおイモ | Aizu Online Judge 愚直にやるとこうなる main(n,m,i,c,p){while(scanf("%d %d",&n,&m),n){int a[1010]={};c=p=0;while(c++!=n){for(i=0;i++

aoj 0084

検索エンジン | Aizu Online Judge splitで指定文字を除くのがミソ $>Errorになるのでこれで妥協 puts gets.split(/[ ,.]/).select{|v|v[2]&&!v[6]}*' '

aoj 0083

西暦と和暦の変換 | Aizu Online Judge 無理やり感ある #!ruby -nrdate a=Date;y,m,d=$_.split.map &:to_i s=((b=a.new(y,m,d))

aoj 0082

メリーゴーランド | Aizu Online Judge 愚直そのもの int main() { int a[] = {1,2,1,2,1,4,1,4}; int p[8],i=0; while(cin >> p[(i++)%8]) { if(i%8==0) { int ans=0, md = 1e5, dist; REP(j,8) { dist = 0; REP(k,8) dist+=(p[k]>a[(k+j)%8])?p[k]-a[(k+j)…

aoj 0081

線対称の位置にある点 | Aizu Online Judge回転行列を利用してといてみた. class P { public: double x,y; P(){};P(double x,double y):x(x),y(y){}; P operator+(const P&q){P t;t.x=x+q.x;t.y=y+q.y;return t;} P operator+=(const P&q){x+=q.x;y+=q.y;ret…

aoj 0080

3乗根 | Aizu Online Judge やるだけ void solve(const double q) { double a = q*0.5; while(fabs(a*a*a-q)>=1e-5*q) a -= (a*a*a-q)/(3*a*a); printf("%lf\n",a); } int main() { double q; while(cin>>q,q>0) solve(q); return 0; }

aoj 0078

魔方陣 | Aizu Online Judge やるだけ int main() { int n; while(cin>>n,n) { vector<vector<int> > T(n,vector<int>(n,0)); int k=1, x=n/2+1, y=n/2; while(k<=n*n) { while(T[x][y]!=0) { x = (x+1)%n; y = (y-1)%n; y = y<0? y+n : y; } T[x][y] = k; x = (x+1)%n; y = (</int></vector<int>…

aoj 0077

ランレングス | Aizu Online Judge やるだけ int main() { string s, t = ""; while(cin >> s) { int i = 0; while(i < s.size()) { if(s[i]=='@'){REP(j,s[i+1]-'0')t+=s[i+2];i+=3;} else{t+=s[i];i++;} } cout << t << endl; t = ""; } return 0; }

aoj 0076

宝探し | Aizu Online Judge 久々の投稿 複素数を使うと以下のような漸化式になる あとはやるだけ int main() { int n; complex<double> T[1001]; T[0] = 0.0; T[1] = 1.0; FOR(i,2,1000) T[i] = complex<double>(1.,1./sqrt(i-1))*T[i-1]; while(cin >> n, n!=-1) { printf(</double></double>…

aoj 0075

BMI | Aizu Online Judge (64byte) $<.map{|s|a,b,c=s.split(',').map &:to_f;p a.to_i if b/c**2>=25} すごい答え(42byte) #!ruby -n i,w,h=eval"[#$_]" w/h/h<25||p(i) $_について module Kernel (Ruby 2.0.0) オートスプリットを利用した答え(52byte) #!r…

aoj 0074

ビデオテープ | Aizu Online Judge 更に短いコードのためもう一度cから勉強したほうがいいかもしれない. #include <bits/stdc++.h> #define A printf("%02d:%02d:%02d\n",a/e,a/d%d,a%d); main(){int h,m,s,a,d=60,e=d*d;for(;scanf("%d%d%d",&h,&m,&s),~h;){a=e*2-h*e-m*d-</bits/stdc++.h>…

aoj 0073

四角すいの表面積 | Aizu Online Judge #include <bits/stdc++.h> main(){int x,h;while(scanf("%d%d",&x,&h),x)printf("%f\n",x*(x+sqrt(4*h*h+x*x)));} $<.map(&:to_i).each_slice(2){|a|puts (x=a[0]**2)+Math.sqrt(x*(x+4*a[1]**2))if a[0]!=0}</bits/stdc++.h>

aoj 0072

灯篭 | Aizu Online Judge 最小全域木のコストを求めさせる問題. プリムのアルゴリズムを用いる. Spaghetti Source - 最小全域木 (Prim) typedef int Weight; struct Edge { int src, dst; Weight weight; Edge(int src, int dst, Weight weight) : src(src)…

aoj 0071

爆弾の連鎖 | Aizu Online Judge x,y座標の変換に注意. int a[8][8]; void f(int x, int y) { a[x][y]=0; for(int i=1;i<=3;i++) { if(x+i>7) break; if(a[x+i][y]) f(x+i,y); } for(int i=1;i<=3;i++) { if(x-i<0) break; if(a[x-i][y]) f(x-i,y); } for(in…

aoj 0070

組み合わせの個数 | Aizu Online JudgedfsでやるとTLEしてしまった... int a[10]; int n, s, cnt; void dfs(int m,int r) { if(r > s) return; if(m==n+1) { if(r==s) cnt++; return; } FOR(i,0,9) { if(a[i]) { a[i] = 0; dfs(m+1,r+m*i); a[i] = 1; } } } …

vector<int>::size_typeについて

c++のvectorのsize()を使うときに注意すること. using namespace std; int main() { vector<int> a(5); cout << 4-a.size() << endl; return 0; } このプログラムの出力は-1ではなく以下のようになる. 18446744073709551615原因はsize()の返り値がunsigned intで</int>…

aoj 0068

輪ゴムで囲んだ釘 | Aizu Online Judge 凸包を求めるアルゴリズムが蟻本に乗っていたのでそのまま書いてみた. class P { private: double add(double a,double b){if(abs(a+b)<1e-10*(abs(a)+abs(b)))return 0;else return a+b;} public: double x,y; P(){};…

aoj 0067

島の数 | Aizu Online Judgeやるだけ. int a[144]; void flip(int x, int y) { if(x>11||y>11||x<0||y<0) return; if(a[x+y*12]) { a[x+y*12] = 0; flip(x+1,y); flip(x-1,y); flip(x,y+1); flip(x,y-1); } } void solve() { int cnt = 0; REP(i,144) if(a[i…

aoj 0066

三目並べ | Aizu Online Judge やるだけ. char solve(string s) { REP(i,3) { if(s[i*3]==s[i*3+1]&&s[i*3]==s[i*3+2]&&s[i*3]!='s') return s[i*3]; if(s[i]==s[i+3]&&s[i]==s[i+6]&&s[i]!='s') return s[i]; } if(s[0]==s[4]&&s[4]==s[8]&&s[4]!='s') retu…