Tamflexの貯蔵庫

やる気のない備忘録

c++

複数プロジェクトが併存できるC++用の汎用的なMakefile

c++

汎用的に使えるMakefileを作ってみた github.com詳細はこちらから qiita.com

codeIQ 1621

なんとなくlower_boundでやりたかった. vector<ll> p; void init() { p.resize(10000); fill(p.begin(),p.end(),1LL<<30); int m = 0; p[m++] = 2; for(int i=3;i<=100000;i+=2) { bool flag = true; for(int j=0;p[j]*p[j]<=i;j++) { if(i%p[j]==0) { flag = fa</ll>…

codeIQ 1630

abc029にも似た, というかほぼ同じな問題があった. 念の為long long intを用いて解いた. void solve(ll n) { n++; ll d = 1, t, ans = 0; while(n>d) { t = d; d *= 10; ans += (n/d)*t+(n%d>t*7?(n%d>=t*8?t:(n%d)%t):0); } cout << ans << "\n"; } int mai…

codeIQ 2717

dfsで全探索する. int n; vector<int> a; int dfs(int k, int sum) { if(k==a.size()) return sum; int x = dfs(k+1,sum+a[k]); int y = dfs(k+1,sum); return abs(n-x)<abs(n-y)? x : y; } int main() { int t; cin >> n; while(cin >> t) a.push_back(t); sort(a.begin(),a.end()); cout << abs(n-dfs(0,0)) <</abs(n-y)?></int>…

abc 039

abc039.contest.atcoder.jp白の周囲を白にして、その後黒の周囲を黒にして元に戻るか判定する. int h,w; int a[100][100]; int b[100][100]; int c[100][100]; bool range(int x, int y) { return x >= 0 && x < h && y >= 0 && y < w; } int main() { cin >…

abc 037

abc037.contest.atcoder.jp DPの問題. 点を始点とする場合の数はその点の周囲のマス+1になる.あとは再帰的に計算してメモ化すればできる. ll T[1010][1010]; ll dp[1010][1010]; int h, w; ll f(int x, int y, ll p) { if(x<0||x>=h||y<0||y>=w) return 0; i…

c++11でグラフ理論

競プロで頻繁にグラフが絡む問題が出題されます. 隣接リストvector<vector<Edge>>を継承してグラフの様々な値を求めることができるクラスを作ってみました. 使い方 ・宣言 graph<T> g((int)n): 重みの型がTの,頂点の数がn個の隣接グラフの宣言. ・辺の挿入 g.direct(src,dst,</t></vector<edge>…

aoj DPL_3_B

最大長方形 | 動的計画法 | Aizu Online Judge最大長方形を求めさせる問題. こちらを参考に応用する. ヒストグラムの中の最大長方形を求める問題はlog(n)のオーダで求まる. あとはそれを行について行えばいいので計算量はlog(w*h)のオーダーで行ける. struct…

codeIQ 2749

オセロのAIと同じ要領でできる. 初手判定に要注意 int a[25]; int dx[] = {1,1, 1, 0,-1,-1,-1,0}; int dy[] = {1,0,-1,-1,-1, 0, 1,1}; vector<string> split(const string &str, char delim) { istringstream iss(str); string tmp; vector<string> res; while(getline(iss</string></string>…

aoj DPL_1_E

Edit Distance (Levenshtein Distance) | Aizu Online Judge レーベンシュタイン距離を求めるアルゴリズム. 下が詳しい レーベンシュタイン距離 - Wikipedia dpでいける. int solve(string s, string t) { int sn = s.size(); int tn = t.size(); vector<vector<int>> dp</vector<int>…

aoj DPL_1_D

Longest Increasing Subsequence | Aizu Online Judge まずはTLEする回答 dp[i] := (i番目の数まで使った時の最長部分増加列の大きさ) と定義しnについて2重ループを回せばO(n^2)で解ける. すなわちが以前の数より大きければ のような漸化式が書ける int mai…

aoj DPL_1_C

Knapsack Problem | Aizu Online Judge 0/1ナップサック問題とループの向きが単に逆になったパターン. 今回は中間結果の情報は破棄しているが問題によっては必要になる. #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 DOWN(i,b,a) for(int i=b;i>=a;i--) using namespace std; int main() {</n;i++)></bits/stdc++.h>…

aoj DPL_1_B

0-1 Knapsack Problem | Aizu Online Judge 前にも書いたが0/1ナップサック問題には大きく分けて3つの解法がありそのうちdpである2つの解法で解いてみた dp[i][j] := (i番目の品物まで使った時の合計重さjを実現できる最大合計価値) dp[i][j] := (i番目の品…

C++でのBM法の実装

c++

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

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

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 = …