計算アルゴリズム【�U】

1 名前：デフォルトの名無しさん mailto:sage [2005/10/15(土) 20:42:23 ] 前スレ創設者 FeaturesOfTheGod ◆UdoWOLrsDM の言葉 >プログラム板の皆さん、こんにちは。 >無謀にもこんなスレを立ててみました。 >四則演算、初等関数、その他の関数の関数値を求めるアルゴリズムについての話をしましょう。 >人間にとって計算しやすい方法についても別途語ることにしましょう。 前スレ↓ pc8.2ch.net/test/read.cgi/tech/1090227743/ 152 名前：デフォルトの名無しさん mailto:sage [2005/11/13(日) 17:16:00 ] >>151 入力データの仕様をください 153 名前：デフォルトの名無しさん mailto:sage [2005/11/13(日) 18:28:11 ] >>147 入力の仕様がわからなかったので，以下の仕様に従うことにした．標準入力から読み込む． 「1行目: 都市の数，2行目から: 都市1 都市2 その間の距離，終端: EOF」 #include <iostream> #include <vector> #include <algorithm> using namespace std; const int inf = 1 << 30; // brute force using recursion ( O(n!) ) int solve(int pos, int start, vector< vector<int> >& adj, vector<bool>& used) {   used[pos] = true;   int d = inf;    if (find(used.begin(), used.end(), false) == used.end()) d = adj[pos][start];   for (int i = 0; i < used.size(); ++i)      if (!used[i]) d = min(d, adj[pos][i] + solve(i, start, adj, used));   used[pos] = false;   return d; } int main() {   int n; cin >> n;   vector< vector<int> > adj(n, vector<int>(n, inf));   int a, b, d;   while (cin >> a >> b >> d) adj[a][b] = adj[b][a] = d;   vector<bool> used(n, false);   int start = 0; used[start] = true;   cout << solve(start, start, adj, used) << endl; } 154 名前：147 mailto:sage [2005/11/14(月) 06:25:03 ] >>153 ありがとうございます。仕様はそのとおりでいきましょう。 走らせるために #define min(a,b) (((a)<(b))?(a):(b)) と return 0; を追加しましたが入力の仕方がいまいち分かりません。 都市a bがint型ですよね…。 -------------------------------------------------- 3 a b 3 -1073741824 Press any key to continue -------------------------------------------------- これ↓ではエラーが出て駄目ですし… -------------------------------------------------- 3 1 2 3 4 5 6 Press any key to continue -------------------------------------------------- どう入力するんでしょうか？m(__)m 155 名前：147 mailto:sage [2005/11/14(月) 11:24:19 ] >>153 ダメダメな自分でもようやく分かりました。 Ctrl+Qで終了というのが渋いですね。 少し改良しました。 -------------------------------------------------- Number of Cities: 4 Source Destination Distance: 0 1 2 0 2 4 0 3 3 1 2 3 1 3 6 2 3 1 ^Q 9 Press any key to continue -------------------------------------------------- 後はもう少しコメントなどを入れて改良してみます。(当然n数は10を超えない程度で) ありがとうございました！ 156 名前：147 mailto:sage [2005/11/14(月) 11:58:32 ] 仕様を少し変更しました。パスの数をユーザーから求めて、その数だけ入力がされたら終了、としました。 #include <iostream> #include <vector> #include <algorithm> #define min(a,b) (((a)<(b))?(a):(b)) using namespace std; const int inf = 1 << 30; // brute force using recursion ( O(n!) ) int solve(int pos, int start, vector< vector<int> >& adj, vector<bool>& visited) { visited[pos] = true; int d = inf; if (find(visited.begin(), visited.end(), false) == visited.end()) d = adj[pos][start]; for (int i = 0; i < visited.size(); ++i) if (!visited[i]) d = min(d, adj[pos][i] + solve(i, start, adj, visited)); //←ここで何か記録しないといけないんでしょうが… visited[pos] = false; return d; } void main() { int cities, paths; cout << "Number of Cities: "; cin >> cities; cout << "Number of Paths: "; cin >> paths; vector< vector<int> > adj(cities, vector<int>(cities, inf)); int a, b, d; cout << " S D Distance" << endl; for (int i = 0; i < paths; ++i) { cout << "Path" << i+1 << ": "; cin >> a >> b >> d; adj[a][b] = adj[b][a] = d; } vector<bool> visited(cities, false); int start = 0; visited[start] = true; cout << "The minimum distance is " <<solve(start, start, adj, visited) << endl; } あの、で、最短のルートを表示するにはどうしたらいいんでしょうか？本当にしつこくてすみません。m(__)m 157 名前：デフォルトの名無しさん mailto:sage [2005/11/14(月) 17:58:34 ] いくつか手はあるけど，手っ取り早いのは「各 visited の状態から次何を選んだか」を グローバルの map<vector<bool>, int> にでも覚えさせておくこと． 表示するときは適当な visited からスタートして，visited = 全部 true になるまで覚えさせた map を辿っていく． 具体的には下みたいな感じにやる． kansai2channeler.hp.infoseek.co.jp/cgi-bin/joyful/img/1082.cpp 158 名前：157 mailto:sage [2005/11/15(火) 06:58:17 ] solve(i, start, adj, visited) を二回も呼んでるのが非効率なのでなんとかしといてください 159 名前：147 mailto:sage [2005/11/15(火) 15:43:57 ] >>157 ありがとうございます！m(__)mﾊﾊｧｰ　m(__)mﾊﾊｧｰ　m(__)mﾊﾊｧｰ　m(__)mﾊﾊｧｰ　m(__)mﾊﾊｧｰ 再帰の途中で覚えさせる方法が分からなくて困っていました。 なるほど、ああやって割り込むんですね。本当に勉強になります。 solve(i, start, adj, visited)を二回呼んでるのは後でなんとかします。 ああ、もっと賢くなりたい。←自分 160 名前：フローチャート mailto:fh [2005/11/23(水) 20:25:34 ] 価格Ａ，Ｂの商品をそれぞれＸ，Ｙ個、購入する場合の 支払い金額を計算するフローチャートを示せ。ただし、値引き率８％と 消費税５％を考慮すること 161 名前：デフォルトの名無しさん mailto:sage [2005/11/23(水) 20:28:11 ] >>160 そりは難しい課題だねえ。 まずはフローチャートを描くルーチンを作成しないと。 環境は何？　まさか、線を引くルーチンから作らなくちゃいけない？ 162 名前：デフォルトの名無しさん mailto:sage [2005/11/26(土) 17:42:03 ] その前に有理数体での加減乗除を定義しないと。 また、有理数の円未満の端数処理をどうするか検討しないとな。 切捨て？切りage？四捨五入？設定で選択できるようにしないと 実用に耐えないかもよ。 163 名前：デフォルトの名無しさん mailto:sage [2005/11/26(土) 20:15:28 ] 意地悪なやつらめ。そんなおまえらが好きではあるけど。 164 名前：デフォルトの名無しさん mailto:sage [2005/11/26(土) 21:21:57 ] まあ 137 あたりの経緯があるからな 165 名前：デフォルトの名無しさん mailto:sage [2005/11/27(日) 18:06:06 ] >>162 有理数体上の加減乗除は自明じゃねーの？ 166 名前：デフォルトの名無しさん mailto:sage [2005/11/27(日) 21:26:38 ] ていうかコンピュータ上では有理数全てを表現できないし 167 名前：デフォルトの名無しさん mailto:sage [2005/11/27(日) 21:35:56 ] ディジタルコンピュータ上ではできないね。 168 名前：デフォルトの名無しさん mailto:sage [2005/11/27(日) 21:45:26 ] 全ては表現できなくても適当な位相さえ入れば十分でしょ 169 名前：デフォルトの名無しさん mailto:sage [2005/12/01(木) 12:06:34 ] てか、実数から見たら表現できる数値なんて穴だらけだよな。 170 名前：デフォルトの名無しさん mailto:sage [2005/12/01(木) 18:13:16 ] つ 区間演算 171 名前：デフォルトの名無しさん mailto:sage [2005/12/04(日) 01:41:03 ] 配列nに対してn-1個の値があり、残りひとつはNULLでも0でもいいのでデータとしては何も入ってない状態。 その配列nの中のn-1個の数列を配列nの中だけでソートしたいのですが何かいいアイデアないでしょうか？ つまりデータ構造やポインタは勿論、 二分木や多分木を使ってヒープも駄目ということです。 172 名前：デフォルトの名無しさん mailto:sage [2005/12/04(日) 02:00:23 ] バブルソートなりクイックソートなり 173 名前：デフォルトの名無しさん mailto:sage [2005/12/04(日) 03:54:54 ] >171 状況がよくわからん。 具体的に。 174 名前：デフォルトの名無しさん mailto:age [2005/12/05(月) 16:33:28 ] カタカナの「ノ」の形の離散データの極率を求める数値計算法ってありますか？ 書物で調べたりググったりしていますが、見つかりません。 175 名前：174 mailto:age [2005/12/05(月) 16:35:59 ] 極率　ではなく　曲率　でした。 176 名前：デフォルトの名無しさん mailto:sage [2005/12/05(月) 16:40:06 ] 二次回帰分析とかって話? 177 名前：174 mailto:age [2005/12/05(月) 16:58:41 ] >>176 レスありがとうございます。 二次回帰分析ですと、y=a(x^2)+bx+c のa,b,cを決定するとこになると思われます。 今やりたいことは、例えば y = exp(-3*p)+q^2 に対応する離散データから (p, q)を求めるというものです。 exp(x)は２次関数ではありませんが、例えばテイラー展開等で近似をすれば、 適応出来る物なのでしょうか？ 178 名前：デフォルトの名無しさん mailto:sage [2005/12/05(月) 19:46:30 ] 最小二乗フィッティングだね． 理論的には，s_1, ..., s_n をパラメータとする関数 f(x; s_1, ..., s_n) を データ列 (x_1, y_1), ..., (x_N, y_N) にフィッティングするときは誤差の二乗和 φ(s_1, ..., s_n) = Σ(f(x_i; s_1, ..., s_n) - y_i)^2 を最小化する s_1, ..., s_n を求めてやればいい． で，それは結構難しくて，凸関数とか条件付かないと最良解が出る保証は無い． ある程度の解でいいなら，「勾配法 + アニーリング」くらいで，それなりに求まってくれるはず． 179 名前：デフォルトの名無しさん [2005/12/05(月) 19:50:09 ] 画像調整のガンマをスクロールバーで実装したいのですが、 イマイチ正確な定義が分かりません。 それか、組み込めるライブラリどこかに落ちてないでしょうか？ 180 名前：174 mailto:age [2005/12/05(月) 20:15:07 ] >>178 詳しいレスありがとうございます。 今現在では、範囲と精度を与えて、全部の(p, q)について調べ、その誤差の二乗和が 最小になる(p, q)の組を解としています。 勾配法、アニーリングを調べてみて、検討してみます。 181 名前：デフォルトの名無しさん mailto:sage [2005/12/05(月) 20:17:32 ] >>179 質問の意味がわからないけど多分スレ違い。 ・言語と開発ツールは何を使ってるのか ・具体的にやりたいことは何か ・今自分は何がわからないのか この三点をはっきりさせて、適切なスレに行って聞いてください。 182 名前：デフォルトの名無しさん [2005/12/05(月) 21:10:48 ] あるデータ列に Y = (A(X^3)+B(X^2))/(CX+D) の最小自乗法を適応するばあい、 残差の２乗和 S に対して A, B, C, Dそれぞれで偏微分したもの = 0 の連立方程式を解いて、 A = Σを含む式 B = ... C = ... D = ... とすれば良いのでしょうか？ 簡単な式の最小自乗法しか使っていなかったので、 そのまま適応して良いものかどうかが分かりません。 お願いします。 183 名前：デフォルトの名無しさん mailto:sage [2005/12/05(月) 21:18:39 ] >>182 OK． ただ偏微分したもの = 0 が陽に解ける保証はどこにもないし，解けたとしても解は一般に複数存在するので 出した解が本当に良いフィッティングになってるかは十分議論する必要がある． 184 名前：182 [2005/12/05(月) 21:30:11 ] >>183 ありがとうございます。勉強になります。 185 名前：アルゴリズム mailto:fh [2005/12/06(火) 10:20:52 ] N人分のデータ（氏名、体重、身長、年齢）がＤＡＴＡ文で入力されているプログラムが ある。これを用いて次のプログラムをＢＡＳＩＣで作成しなさい 年齢が３０歳以下の人の、体重と身長の平均値を計算し表示する 186 名前：デフォルトの名無しさん mailto:sage [2005/12/06(火) 10:59:37 ] select average(weight), average(height) from DATA where (age <= 30); 187 名前：デフォルトの名無しさん mailto:sage [2005/12/06(火) 12:50:40 ] >185 スレタイくらい嫁。 くだ質スレか宿題スレ池や。 188 名前：デフォルトの名無しさん mailto:sage [2005/12/06(火) 20:36:58 ] >177 エクセルにデータを叩き込んでexpの形で近似曲線出すのが一番手軽な希ガス。 ところでexp（）の形になると判ってるなら、logを取ったらどうなる？ 189 名前：デフォルトの名無しさん mailto:sage [2005/12/06(火) 21:51:46 ] プチ統計学講座だなｗ 190 名前：デフォルトの名無しさん mailto:sage [2005/12/06(火) 22:09:54 ] >>188 ヒント：Windowsとは限らない 191 名前：デフォルトの名無しさん mailto:sage [2005/12/06(火) 22:28:03 ] ひんと：Excel以外にも統計処理できるソフトはあるし、そもそもMacでもExcelは動く。 192 名前：デフォルトの名無しさん mailto:sage [2005/12/06(火) 22:34:45 ] >>190-191 ヒント：そもそも汎用機とも限らない　例：宇宙計算用スパコン だからこそのアルゴリズムスレだろ？ 193 名前：デフォルトの名無しさん mailto:sage [2005/12/06(火) 23:14:51 ] だからlogとったらどうなるかと聞いてみたんだけど 194 名前：デフォルトの名無しさん mailto:sage [2005/12/06(火) 23:30:44 ] だから前半にだけ突っこんだんだけど 195 名前：デフォルトの名無しさん mailto:sage [2005/12/07(水) 19:22:24 ] まぁまぁ、次のネタを待ちましょう 196 名前：デフォルトの名無しさん [2005/12/07(水) 19:55:37 ] 上にも最小自乗法の話題が出ていますが･･･ y=(A(x^3)+B(x))/(C(x^2)+D) の関数に最小自乗法を適応することは出来ますか？ 残差の２乗和が最小になる(A,B,C,D)を見つけるわけですが、 例えばCで偏微分して=0となるCは、単なる極地のCであり、 それが極小かつ最小となるわけではありませんよね？ 197 名前：196 [2005/12/07(水) 19:57:51 ] あれ！？ よく見たら>>182とまったく同じ式でした･･･。 当方、学生ではないので課題とかでもないのですが。 同じ書物を読んでいるのかしら。 >>183を参考にします。 無視してください。 198 名前：デフォルトの名無しさん mailto:sage [2005/12/07(水) 20:29:05 ] フローチャートしか知らない 199 名前：デフォルトの名無しさん [2005/12/08(木) 10:14:07 ] 二次元配列と、 その中の３つの点で定義される三角形があって、 ある点がその中か外か判定するアルゴリズム教えて下さい。 200 名前：デフォルトの名無しさん mailto:sage [2005/12/08(木) 11:58:18 ] >>199 (0,0), (0,2), (2,0) という要素数3個の2次元配列が与えられた場合、 (1,1) という点がその三角形の中にあるかどうかということだよな？ 1)与えられた3点が三角形を構成するかどうかを吟味。 2)その3点から各々2点を取り出し(3通り)、2点を通る直線の式を求める。 3)その3本について、直線で座標平面を分割した場合、残りの頂点を含む領域を不等式(3通り)であらわす。 4)中か外かを判定したい点について、その3つの不等式を同時に満たせば中、そうでなければ外。 201 名前：デフォルトの名無しさん [2005/12/08(木) 12:26:41 ] そのまんまじゃん 202 名前：デフォルトの名無しさん mailto:sage [2005/12/08(木) 13:04:30 ] ナップザックを背負ったバックパッカー*がいるとする。ナップザックの中は、 できるだけ使い勝手の良いもので満たさなければならない。ナップザックの大きさは限られているので、 どのような荷物で満たすかが非常に大事な問題となる。ナップザック問題とは、ある目的のコストが 最大となるように限られたスペースにアイテム(item)を配置する問題のことである。次の問いに答えよ *リュックサックに全ての荷物を入れて旅行する人のこと。 問1 0/1 ナップザック問題とは何か？ 問2 全てのアイテムが(1)同じコストあるいは(2)同じサイズであればナップザック問題はどん な問題となるのだろうか？ 問3 アイテムのコストが大きさに比例する場合は、ナップザック問題はどんな問題となるの か？ 問4 ナップザックの大きさがアイテムに比べて非常に大きい場合は、一体何が起こるか？参考 書にあるプログラムを使って調べ、現象から生じた問題を指摘し、それが何故生じるのか 記述せよ． 203 名前：デフォルトの名無しさん mailto:sage [2005/12/08(木) 15:35:02 ] >>202 ここは君が来るべきスレでは無いよ。 204 名前：デフォルトの名無しさん mailto:sage [2005/12/08(木) 15:59:27 ] >>203 いいじゃんべつに。 お前の個人的な感情は出さなくてもいいよ。 そんなもの読んでいる人間には全く価値がないから。 文句があるならどのスレにいけ、ぐらいは言えよ。 あ、頭悪スレとかに誘導するのは全く無意味だからやめてね。 205 名前：デフォルトの名無しさん mailto:sage [2005/12/08(木) 18:15:33 ] >>204 つまんないﾈﾀ書くのも自由だし それを批判するのも自由。 206 名前：デフォルトの名無しさん [2005/12/08(木) 18:37:18 ] 座標３つ（どれがどの位置か不明）から、３角形の面積を計算する方法はありますか？ 207 名前：デフォルトの名無しさん mailto:sage [2005/12/08(木) 18:40:49 ] >>205 自重しろよ。 自由ではないんだよ？あまりひどいとアク禁だ。 208 名前：デフォルトの名無しさん mailto:sage [2005/12/08(木) 19:52:21 ] >>206 3点から辺の長さ求めてHeronの公式、とか。 209 名前：デフォルトの名無しさん mailto:sage [2005/12/08(木) 20:01:43 ] 三点のベクトルが v1, v2, v3 のとき [v1-v2, v3-v2]/2 の絶対値が面積になるね。 [x,y] は x と y の外積ね。 210 名前：デフォルトの名無しさん mailto:sage [2005/12/08(木) 20:02:07 ] >>206 ある点を原点に移動させるように全体を並行移動させて、 (0,0),(a,b),(c,d)の三角形の面積が(1/2)*|a*d - b*c| であることを 使った方が簡単かも。 ３点が三角形が作られない位置に無いことを確かめてからね。 211 名前：デフォルトの名無しさん mailto:sage [2005/12/08(木) 20:03:50 ] 実は （>>182の式）≠（>>196の式） である件について 212 名前：デフォルトの名無しさん mailto:sage [2005/12/08(木) 20:55:54 ] >>206 www.tensyo.com/urame/prog/linealgo.htm の多角形の面積でやれば？ 213 名前：デフォルトの名無しさん mailto:sage [2005/12/08(木) 23:07:36 ] >>199 >>99,>>107-108 >>204 >>187 てか、質問する前にスレ内くらいは読もうぜ？ 214 名前：デフォルトの名無しさん mailto:sage [2005/12/13(火) 18:24:45 ] >>206 でも、3点とも座標がわかってて位置不明ってどゆこと？ 215 名前：デフォルトの名無しさん mailto:sage [2005/12/13(火) 19:24:25 ] >214 「固定座標じゃない」って意味だろ。 確かに意味は解らないが、気持ちは伝わってきた。 216 名前：デフォルトの名無しさん mailto:sage [2005/12/13(火) 20:30:04 ] ヘロンの公式のことをたずねているの？ 217 名前：デフォルトの名無しさん mailto:sage [2005/12/13(火) 20:39:02 ] 辺の長さが与えられてるわけじゃなく、頂点座標だから ヘロンの公式をわざわざ使うより 素直に台形公式による閉曲線積分 218 名前：デフォルトの名無しさん mailto:sage [2005/12/13(火) 20:56:04 ] >>217 頂点座標？閉曲線積分？言葉はちゃんと使おうぜ 219 名前：デフォルトの名無しさん mailto:sage [2005/12/14(水) 01:18:24 ] >>214 任意の三点ってことだろ。 220 名前：デフォルトの名無しさん mailto:sage [2005/12/14(水) 11:53:04 ] 三角形なんだから、3点を P, Q, R として、 2つのベクトル a ＝ PQ、b ＝ PR を定義。 (1/2) |a × b| ただし、×は3次元ベクトルの外積、|| はベクトルの絶対値。 っつーか、ぐぐれ。 homepage2.nifty.com/PAF00305/math/triangle/node3.html 221 名前：BASIC mailto:fh [2005/12/14(水) 15:57:56 ] 次のプログラムをＢＡＳＩＣで作成しなさい データを１件ずつ、ユーザ定義関数を利用して 「公共料金」を計算させることにします ＜水道料金の計算＞ 基本料金：1250円 定量制料金単価は ・使用量が10立方米以下であれば80円／立方米 ・ 〃 10立方米超20立方米以下であれば185円／立方米 ・ 〃 20立方米超50立方米以下であれば205円／立方米 ・ 〃 50立方米超100立方米以下であれば240円／立方米 ・ 〃 100立方米超200立方米以下であれば275円／立方米 ・ 〃 200立方米超であれば310円／立方米 （計算例） 使用量が8立方米であれば、1250＋8×80 使用量が15立方米であれば、1250＋10×80＋(15−10)×185 ＝2050＋(15−10)×185 使用量が40立方米であれば、1250＋10×80＋10×185＋(40−20)×205 ＝3900＋(40−20)×205 使用量が75立方米であれば、1250＋10×80＋10×185＋30×205＋ (75−50)×240＝10050＋(75−50)×240 使用量が120立方米であれば、1250＋10×80＋10×185＋30×205 ＋50×240＋(120−100)×275＝22050＋(120−100)×275 使用量が230立方米であれば、1250＋10×80＋10×185＋30×205 ＋50×240 ＋100×275＋(230−200)×310 ＝49550＋(230−200)×310 222 名前：デフォルトの名無しさん mailto:sage [2005/12/14(水) 15:58:42 ] 作成しなさい、だってさ↓ 223 名前：デフォルトの名無しさん mailto:sage [2005/12/14(水) 16:01:00 ] 　　　　　　　　　　　｜　　　　　∧＿∧ 　　　　　　　　　　　｜.　　　 　(･ω･` )ｽﾙｰ 　　　　　　　　　　　｜ｽｽ…　/J　　 J 　　　　　　　　　　　↓　　　,,, し―-J 224 名前：デフォルトの名無しさん mailto:sage [2005/12/14(水) 16:51:10 ] アルゴリズムじゃねえな 225 名前：デフォルトの名無しさん mailto:sage [2005/12/14(水) 19:15:46 ] アゴリズムだな 226 名前：デフォルトの名無しさん mailto:sage [2005/12/15(木) 00:36:36 ] アルゴリズムじゃねえって！ 227 名前：デフォルトの名無しさん mailto:sage [2005/12/15(木) 04:40:03 ] アルゴリズム体操！アルゴリズム体操！ 228 名前：デフォルトの名無しさん mailto:sage [2005/12/16(金) 15:42:28 ] 立方米って一般的なのか。 m^3と書かないの。 229 名前：デフォルトの名無しさん mailto:sage [2005/12/16(金) 15:51:32 ] >>228 メートル=米、だからね。 他にも糎とか粍とか。 でも、普通は立方米と書かずに立米(りゅうべい)と書く希ガス。 #平米は割と一般的でしょ。 230 名前：デフォルトの名無しさん mailto:sage [2005/12/17(土) 12:32:02 ] 口頭ではよく「へーべー」と言うな。 俺の父親が。 231 名前：デフォルトの名無しさん mailto:sage [2005/12/17(土) 12:34:10 ] なんでメートルが米なの？？その由来は？ 232 名前：デフォルトの名無しさん mailto:sage [2005/12/17(土) 12:58:32 ] 亜米利加　と同じじゃないの？　メの元の漢字だと思ってたけど ttp://www.fudemojiya.com/syo/katakana.htm によると　女　らしい 233 名前：デフォルトの名無しさん mailto:sage [2005/12/17(土) 13:05:23 ] 明治時代か江戸末期に「メートル(或いはメーター)」という音に似ているから当てたんじゃないの? 他にもグラムやインチなどにも当てられた漢字があることだし。 234 名前：デフォルトの名無しさん mailto:sage [2005/12/17(土) 19:27:50 ] ttp://freett.com/ToshiN/kanji/kanji02.html 発案した人も途中からヤケクソになっていったのだと推測 235 名前：デフォルトの名無しさん mailto:sage [2005/12/17(土) 22:42:29 ] >>234 そこまでして漢字にせんでもええやん…とおもた。 236 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 03:23:45 ] 当時ハ送リ仮名ガ片仮名デアッタタメ 外来語ヲ片仮名ニ割リ当テルトイウ 発想ハ無カッタンダロウナ 237 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 04:08:19 ] 質問なのですが、例えば、 1、２、３、４、５、４０、５０、６０、７０　といった代表値である数字があるとします。 で、例えばXがその代表値のどれに一番近いかを求めるとき すべてを探索すれば答えはわかると思いますが、できるだけ高速に 求めたいと考えたとき、代表地を木にすればいいのではと私は考えました、 例えばXが6だとするとまず、5が一番近いことがわかります これをどうやって木にしたらいいのかがわかりません。 また１，２、３・・・・といっていましたが （X,Y,Z)と三次元に増えた場合、一番近い点を見つけるにはどうしたらいいのでしょう？ 一致するのを探すのではなく一番近いものを探すので苦労しています。 近さの判定としてはユークリッド距離を用います。 どうぞよろしくお願いいたします 238 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 04:13:08 ] 代表値の数列が整列されている前提ならバイナリサーチでいいのでは？ 239 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 04:32:56 ] >>238 ありがとうございます。 バイナリサーチですと完全に一致する場合は問題ありませんが 完全に一致せず一番近いものを選択するときに問題がおきます 237の例で6に一番近いものをバイナリサーチで探そうとすると 40になると思われます 240 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 04:46:58 ] バイナリサーチで最終的に 40 がヒットするわけだけど、 この 40 と X を比較することによって 5 < X < 40 という関係はすぐにわかる。 あとはどちらに近いかを計算するだけじゃないのかな。 241 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 04:54:28 ] >239 最終段でヒットしなかった時、left/rightどちらの方が近いかチェックしろ。 n次の場合は、a=f(X,Y,...)で、aの値でソート。 242 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 05:01:13 ] >>240 >>241 ありがとうございます。 ｎ次の場合はうまくいかないように感じるのですが できれば具体例を出していただけないでしょうか？ 241さんのを自分なりに解釈すると、まずはXについてどれが一番近いか しらべて次はY、Z・・・というふうなことだと感じているのですが 例えばXがかなり近くてもZ,Yがかなり遠いと総合的な距離は遠いわけで Xが多少遠くてもY,Zが近ければ総合的な距離は近くなるので あまりうまくいかなようなきがします。 間違っていたらすみません 243 名前：241 mailto:sage [2005/12/18(日) 05:37:33 ] >242 スマソ、完全に寝ぼけてた。 1次以上の距離で探したいのなら、総当りで検索するしかないんじゃね？ 244 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 05:40:16 ] >>243 ありがとうございます。 総当りですか・・・それですと時間がかかるので 何かいい方法がないかなと思いまして 245 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 07:16:28 ] >>241 ２次元（平面）ならボロノイ図・近傍で似たような問題がある ３次元以上だとできるかどうか知らない 総当りにしたら 246 名前：245 mailto:sage [2005/12/18(日) 07:17:31 ] >>245 の >>241 は >>237 の間違いです ｽﾏｿ 247 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 07:55:24 ] >>245 ありがとうございます。 ボロノイ図については知っていまして、三次元でも可能のようです。 ですが、これをどう木にすればいいのかよくわからなかったもので 何かいい方法があればと思って相談してみました、。 248 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 09:03:12 ] >>237 空間データベースって分野のお話だね． 望んだ形でデータを持ってていいなら，k-d tree を使うのが手軽かな． 空間に一点が追加されるたびに，その一点を通る超平面を引き，その超平面のどちら側に 点があるかで二分木を作る．探索は木を手繰りながら範囲を絞り込んで， その中で全探索をする．最悪 O(n) だけど，適当にばらけてれば O(log n)． ほかにも SR-tree とか SS-tree とか色々あるので，参考になりそうなのを挙げとく： www.kecl.ntt.co.jp/csl/sirg/people/kani/db.html vision.unige.ch/~moennen/publi/moennen_tech0505.pdf H. Samet, The Design and Analysis of Spatial Data Structures, Addison-Wesley 249 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 10:00:28 ] >>248 ありがとうございます 望んだ形でデータを持っていても大丈夫です。 K−ｄ　tree　についてなのですが、これでは一番近い代表値に近似することはできないと 考えられます。 例えば www.kecl.ntt.co.jp/csl/sirg/people/kani/db.html のｋ−ｄ　ｔｒｅｅ　のところの図で説明しますと 代表地がＰと考えると、Ｐ４とＰ５の線の中に点があるのですがＰ９に 一番近い点（一番近い代表地に近似したい点）があったとすると これはＰ４かＰ５に近似されることになります。 実際はＰ９に一番近いのにこれには近似されないとおいうことです 木を作る際にＰ７からでなくＰ６からはじめたとしても同様の問題はおきると考えられます しかしＲ−ｔｒｅｅというのは使えそうな気がします ちょっと読んでみますがいかんせん説明が少ないので理解できるかわかりませんが なんとかがんばってみます ありがとうございました 250 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 10:04:36 ] www.kecl.ntt.co.jp/csl/sirg/people/kani/db.html のＲ−ｔｒｅｅを読んでみましたがこれでも正確に一番近い代表地に 近似されないように感じます。 うーんどうしよう 251 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 10:05:18 ] 一番近い点は一回では求まらないだろう。 　求めた点までの距離を半径とした点で総当りをその後ですればいいだろう。 　たとえば2次元なら　x 座標がその範囲にある条件、次にy座標がその範囲にあるという条件で絞りこんで 　　矩形で絞り込めばいいでしょ 252 名前：デフォルトの名無しさん mailto:sage [2005/12/18(日) 10:07:59 ] あと、この問題は、最悪の場合は、全数検索が必要になる。 たとえば円周上に点が並んでいて、中心に近似しようとした場合など

---

---

[ 続きを読む ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef