http://w.atwiki.jp/imolib/ Imoのアルゴリズムライブラリ ja 2010-06-26T17:18:46+09:00 1277540326 https://w.atwiki.jp/imolib/pages/27.html 2010-06-26T17:18:46+09:00 1277540326 https://w.atwiki.jp/imolib/pages/26.html *レポートの体裁について レポートは指定のクラスファイル（icsreport.cls(202)）を用いて pLaTeX2e で作成し，A4用紙に印刷して提出すること．尚，クラスファイルの使用例は/home/exp/exppub/exp-material/PC2000 に sample-icsreport.tex というファイル名で置いてあるので参考にせよ． 表紙には，クラスファイルの使用例に従い，科目名，課題名，担当教官名，提出者氏名，学籍番号，コース，学年，メールアドレスを明記すること． レポートは内容ごとに節（section）に分割し，さらに必要に応じて小節（subsection）に分割すること． レポートは左端の2箇所をホッチキスでとめて提出すること．クリップによる提出は認めない． プログラムの作成を要求した課題に対するレポートには，課題内容，プログラムの仕様，主要処理の概要，アルゴリズム，プログラムの動作テストの方法，プログラムの動作結果，工夫点，感想，プログラムリストを含めること．プログラム作成以外の課題の場合は，課題内容と課題に対する答え，および感想を含めること．ただし，10行程度のシェルスクリプトの作成はここでいうプログラムの作成とはみなさないが，作成したスクリプトの説明は必ず記述すること． プログラム（短いシェルスクリプトを除く）を作成した場合は，プログラムリストをレポートの最後に付録として添付すること．印刷用紙の節約のため，A4用紙に2ページ分を印刷することが望ましい． レポートは課題ごとに特に指定のない限り，表紙とプログラムリストを除き全体で10ページ以内で作成すること． *レポートの内容について レポート本文は，簡潔にかつ分かりやすく書くこと．第三者が読むことを常に念頭に入れて書くこと． レポート本文は，あくまで報告書であることから客観的な記述の仕方に徹するべきである．従って，感想以外の部分で，個人の感想や主観的な意見，コメントを書くべきではない．また，過去にチャット調の文章を書いていた者もいるが，そのようなレポートは課題の達成度がいかに高かろうと，報告書としての価値は皆無であるとみなされるので注意せよ． 課題内容は，レポート課題の問題文をまる写ししてはならない． レポートに書くべ 2008-05-26T23:13:13+09:00 1211811193 https://w.atwiki.jp/imolib/pages/25.html *Max Power **問題(Code: MXP) a1, a2, …, anとb1, b2, …, bnの二つの数列があった時、 akのbk乗が最大になるkを求めよ。 &italic(){引用: [[Open 2008 (SPOJ)>https://www.spoj.pl/ZFUN08/problems/MXP/]]} **入力 The first line of input contains a positive integer n, not greater than 10000. In the second line you are given a set of positive integers ai separated by spaces, and in the third line – integers bi. All numbers in both sequences are not greater than 10000. It is guaranteed that all power values are different. **計算量 O(N) **ソースコード #include <stdio.h> #include <math.h> int n, o, i; double m, q, a[10001], b[10001]; int main(void) { scanf("%d", &n); { for (i = 1; i <= n; ++i) scanf("%lf", &a[i]); for (i = 1; i <= n; ++i) scanf("%lf", &b[i]); m = log(a[1]) * b[1]; o = 1; for (i = 2; i <= n; ++i) { q = log(a[i]) * b[i]; if (m < q) m = q, o = i; } } printf("%d\n", o); return 0; } **確認 済み //**参考 //***全点対間最短距離 //#list_by_tag(ライブラリ,全点対間最短距離,100,sort=pagena 2008-05-04T23:27:10+09:00 1209911230 https://w.atwiki.jp/imolib/pages/24.html 2008-04-24T00:47:02+09:00 1208965622 https://w.atwiki.jp/imolib/pages/23.html *Page Hopping **問題 Given a graph in which all nodes can be reached from any starting point, your job is to find the average shortest path length between arbitrary pairs of nodes. &italic(){引用: [[UVa 821>http://icpcres.ecs.baylor.edu/onlinejudge/external/8/821.html]]} **概要 グラフが与えられた時、辿り着きうる任意の二点を選んだ時の平均の最短距離を求める **計算量 O(N&u(){3}) **ソースコード // Page Hopping #include <stdio.h> #include <string.h> int main(void) { int a, b, p = 0, m[100][100]; while (scanf("%d%d", &a, &b) && a && b) { memset(m, 1, sizeof(m)); p++; do { m[a - 1][b - 1] = 1; } while (scanf("%d%d", &a, &b) && a && b); // cf. Floyd Warshall for (int k = 0; k < 100; k++) for (int i = 0; i < 100; i++) for (int j = 0; j < 100; j++) if (m[i][k] + m[k][j] < m[i][j]) m[i][j] = m[i][k] + m[k][j]; for (int i = 0; i < 100; i++) for (int j = 0; j < 100; j++) if (i != j && m[i][j] < 100) { a += m[i][j]; b++; } printf("Case 2008-04-17T00:18:58+09:00 1208359138 https://w.atwiki.jp/imolib/pages/22.html *問題 #list_by_tag(問題,100,sort=pagename,sort2=asc,errmsg=該当する問題はありません) 2008-04-16T23:57:56+09:00 1208357876 https://w.atwiki.jp/imolib/pages/21.html *素因数分解 **説明 素因数の発見は順番に探索していく単純な方法。 n=p&u(){k}q&u(){m}p&u(){n}…の約数の個数は(k+1)(m+1)(n+1)…であることを利用して解いている。 **計算量 O(N&u(){0.5}) **使い方 divisor_count(n)でnの約数の個数を返す。 **必要なライブラリ //#list_by_tag(map,1,sort=pagename,errmsg=map) -なし **ソースコード 約数の個数 int divisor_count(int n) { int j, c = 1; for (int i = 2; i * i <= n; i++, c *= j) for (j = 1; !(n % i); j++) n /= i; return (1 < n) ? c * 2 : c; } **確認 なし 2008-04-12T22:40:47+09:00 1208007647 https://w.atwiki.jp/imolib/pages/20.html *素因数分解 **説明 順番に探索していく単純な素因数分解。 **計算量 O(N&u(){0.5}) **使い方 factorize(n)でnを素因数分解した結果をmap<素因数,回数>で返す。 **必要なライブラリ #list_by_tag(map,1,sort=pagename,errmsg=map) **ソースコード 素因数分解 map<int, int> fractorize(int n) { map<int, int> m; for (int i = 2; i * i <= n; i++) for (; !(n % i); m[i]++) n /= i; if (1 < n) m[n]++; return m; } **確認 なし 2008-04-12T22:40:04+09:00 1208007604 https://w.atwiki.jp/imolib/pages/19.html *素数列挙 **説明 エラトステネスの篩および単純な探索。 **計算量 O(N&u(){1.5}) **使い方 primes(n)でnまでの素数を列挙して返す。 **必要なライブラリ #list_by_tag(vector,1,sort=pagename,errmsg=vector) **ソースコード 単純探索 vector<int> primes(int n) { vector<int> p; p.reserve(1.26 * n / log(n)); p.pb(2); int s = 0, t = 4, f = 0; for (int i = 3; i < n; i++) { f = 1; rep(j, s) if (f = i % p[j], !f) break; if (f) if (t - i) p.pb(i); else t = p[++s] * p[s]; } return p; } エラトステネスの篩 vector<int> primes2(int n) { vector<int> p; p.reserve(1.26 * n / log(n)); p.pb(2); int *isd = new int[n]; memset(isd, 0, sizeof(int) * n); for (int i = 3; i < n; i += 2) if (!isd[i]) { p.pb(i); for (int j = i * i; j < n; j += i) isd[j] = 1; } delete [] isd; return p; } **確認 なし **参考 ***素数関連 #list_by_tag(素数,100,sort=pagename,sort2=asc) 2010-06-26T17:10:33+09:00 1277539833 https://w.atwiki.jp/imolib/pages/18.html -実験ー - Imo 2008-04-09 23:17:20 &dummy(1207750640.39) -間違いやら希望があったらここへー - Imo 2008-04-09 23:17:49 &dummy(1207750669.24) 2008-04-09T23:17:49+09:00 1207750669