06/03/11 00:31:55
>>464で、x番目に調べるまでの距離の総和をf(x)とすると、
f(2^n-1)
= 2*Σ{i=1..(n-1)}(2^i-1)
= 2*((2^n-2)-(n-1))
<= 2*2^n
2^(n-1) <= x <= 2^n-1のとき、
f(x) <= f(2^n-1)
<= 2*2^n
= 4*2^(n-1)
<= 4*x
よってO(x)
469:デフォルトの名無しさん
06/03/11 01:10:36
＞＞４６８
なるほど。
助かりました。
ありがとうございました。
470:デフォルトの名無しさん
06/03/14 22:27:00
(0,1)における実数の集合が、可算無限集合ではないことを背理法と対角線論法を使って証明するやつだけど、
いまいち何やってるか分からないんだよね。
分かった気にはなるけど、どうもしっくりこないっつうかなんつうか。
他に証明方法とかあるのかね？
471:デフォルトの名無しさん
06/03/14 22:38:27
>>470
有理数列と実数の部分集合に1対1の対応が作れる
有理数列は自然数→有理数の写像とみなせる
自然数→有理数の写像の濃度は アレフ0^アレフ0 ＝ 2^アレフ0
アレフ0 ＜ 2^アレフ0
（ある集合の濃度が N のとき、その冪集合の濃度は 2^N、
ある集合とその冪集合の間にはどうやっても全単写が作れない）
なので、可算濃度＜実数の濃度
472:デフォルトの名無しさん
06/03/16 13:32:51
冪集合がわからんヒトがいると見た。
473:デフォルトの名無しさん
06/03/18 08:50:19
数学板から誘導されて来ました。
QRを解析するプログラムを作ろうと思っているのですが、
誤り訂正複合のリード・ソロモン符号の複合の仕方が分かりません。
R=(r0,r1,r2,r3...,r25)
R(x)=(r0+r1x+r2x^2+...+r25x^25)
ここにri(i=0～25)は、GF(2^8)の元とする。
とあり、r0-r25に間違っていないデータ（0-255）を代入しているのですが、次のシンドロームを求める式でシンドロームが0になってくれません。
シンドロームSiを求める式
S0=R(1)=r0+r1+r2....+r25
S1=R(a)=r0+r1a+r2a^2....+r25a^25
データをr0-r25に代入して、それらをGF(2^8)の法100011101（a^8+a^4+a^3+a^2+1）で割っているのですが0になりません。
もう数日煮詰まっています。お願いします。何方かご教授して下さい。
474:デフォルトの名無しさん
06/03/18 09:42:13
>>473
誰かのコードを参考にするといいよ
例：
URLﾘﾝｸ(isw3.kankyo-u.ac.jp)
475:デフォルトの名無しさん
06/03/18 11:04:44
>>474さん、ありがとうございます。
ありがたく拝見させてもらいましたが、残念ながら書き込みの処理しか載っていませんでした。
（書き忘れて申し訳ございませんが、473で書いた部分はQRを読み取る部分の処理です）
他にもあるかもと思って「誤り訂正 QR」などで検索を掛けてみましたが、書き込みに関する部分ばかりで、
読み取り解析に関するサイトは見つかりませんでした。
JISの企画書もシンドロームを求めるとしか書いてありませんし、もう、完全に行き詰っています。
どんな些細な事でも結構ですので、アドバイスを頂けないでしょうか。
お願いします。
476:デフォルトの名無しさん
06/03/18 11:27:53
>>475
その部分のコード晒してみて、
+ は XORの事とか 掛算は足算の事とか、そこらへんは判ってるんだよね？
477:デフォルトの名無しさん
06/03/18 11:37:40
>>475
じゃあココのは見たの？ オレは中身はみてないけど
URLﾘﾝｸ(sourceforge.jp)
478:デフォルトの名無しさん
06/03/18 12:37:08
>>476さん、ありがとうございます。
＋は足し算で、掛け算は掛け算で計算していたので、直してみましたが、やはり０になってくれません。
コードはこのようになっています。
//求められたデータ語
int[]Data={
32, 65, 205, 69, 41, 220, 46, 128, 236, 42, 159, 74, 221,
244 ,169, 239, 150, 138, 70, 237, 85, 224, 96, 74, 219, 61 };
//a指数→整数
int[] a_int={1,2,4…（以下続いています）};
public void Syndrome(){
//求めるシンドロームの数を決定
int[]S=new int[8];
for(int i=0;i<S.length;i++){
for(int d=0;d<Data.length;d++){
S[i]=S[i]^(Data[d]+a_int[(d*i)%255]);
}
}
}
>>477さん、
このプログラムは知りませんでした。早速今から落として見て見たいと思います。
ありがとうございます。
479:デフォルトの名無しさん
06/03/18 13:24:07
>>478 やっぱり計算式が違うよ >>474のcalculate.javaをよく見てみて
480:http://www.vector.co.jp/soft/win95/util/se072729.html
06/03/18 19:44:49
TextSS のWindowsXP(Professional)64bit化おながいします
もしくは64bitにネイティブ対応したテキスト置換ソフトありますか？
481:デフォルトの名無しさん
06/03/19 23:37:58
ラジアンをぶつかった壁に対して面対称で
XだけYだけ90°ひっくり返すにはどうすればいいの？
482:デフォルトの名無しさん
06/03/19 23:56:28
いや言ってる意味が良くわかりませんから
483:デフォルトの名無しさん
06/03/20 00:24:24
ビリヤードか
484:デフォルトの名無しさん
06/03/20 00:45:38
vy = -vy;
vx = vx + 90/180*PI
485:デフォルトの名無しさん
06/03/20 11:03:21
>>483 ビリヤードなら90度というのがどこから来るのか・・・・
486:デフォルトの名無しさん
06/03/20 16:05:35
必ず90度で跳ね返るPONGでも作ってみるか
487:デフォルトの名無しさん
06/03/20 16:42:10
その場合、速度はどうなるのん？
488:デフォルトの名無しさん
06/03/20 17:22:36
必ず180で反射するコーナーリフレクターは実在するけど
必ず90度ってどうやるんだ？
489:デフォルトの名無しさん
06/03/20 17:24:03
三角関数？簡単なブロック崩しだったらそんなの使わんぞ
つーか、玉の軌道をＸ軸かY軸で反転すれば問題なし。
490:デフォルトの名無しさん
06/03/20 18:44:58
さて、ここで唐突に「三角関数」とか出てきたわけだが
491:デフォルトの名無しさん
06/03/20 18:52:43
「プログラマが今まで「これはすげえ」と思ったｹﾞｰﾑ」
の149からの流れがここに？
ｽﾚﾘﾝｸ(prog板:149-番)
492:デフォルトの名無しさん
06/03/20 19:24:06
URLﾘﾝｸ(www.dango-itimi.com)
斜面への衝突判定と反射1
493:デフォルトの名無しさん
06/03/23 07:46:12
学校の宿題なのですが・・・。
HはHASH関数で、p は(p-1)/2も素数であるような素数。
aは、1<a<p-1も満たす整数。
g=a^2(mod p)
H(x,y,z,t)=g^(xy+zt) (mod p)
このとき、Hはone-way関数であるが、collision-free関数ではないことを示せ。
なのですが、どなたか分かる方、助けてください。
494:デフォルトの名無しさん
06/03/23 12:28:39
定義は知らんがアバウトに考えて、
剰余が絡んでる段階で多対一関数だから1方向だろうし
collision-freeではなさそうだわな。
証明は2つの値が実際に1個の値になる例を計算すればいいんじゃね？
495:494
06/03/24 16:59:01
エエカゲンに書いたのに誰も突っ込まない…。
ほ、放置プレイ？
496:デフォルトの名無しさん
06/03/24 18:30:11
定義がないのでなんともはや
497:デフォルトの名無しさん
06/03/28 08:18:12
sin(x)/x
Ｑ１ って関数に名前を付けたいけど適切な名前は？
Ｑ２ 浮動小数点なので計算方法は単純に
abs（x)<0.0001 の時は1-x*x/6 でなければ そのまま計算でいいよね？
498:デフォルトの名無しさん
06/03/28 08:53:02
キャラの座標は
左上と真ん中
どちらのポイントを保持して使うべきですか
499:デフォルトの名無しさん
06/03/28 08:55:22
>>497
sine_x_per_x()
500:デフォルトの名無しさん
06/03/28 09:01:48
>>498
いまどきどこでもいい。
処理によって必要な座標は違う。例えばサイドビューのジャンプアクションなら足元座標も使う。
必要な座標を計算するメソッドがあればそれでよい。
数学関係ない。
501:デフォルトの名無しさん
06/03/28 12:19:04
>>497
それは一般的にsinc関数と呼ばれている。
502:497
06/03/28 12:41:21
>>501 ありがとう
503:デフォルトの名無しさん
06/03/28 17:19:04
>>498
RPGの座標には、マップ上の座標と表示する座標があるが、
マップ上の座標で保存する。
表示する時には、表示する座標に変換して表示する。
504:デフォルトの名無しさん
06/03/28 17:58:28
>>503
キャラといえばRPGしか連想できないゲーマーですか?
505:デフォルトの名無しさん
06/03/28 23:55:29
すまん。
でも数ピクセル単位で動かす事になるとおもうから、
その動かせる単位の位置を保存して、
表示する位置に変換するって感じで使えるかも。
506:デフォルトの名無しさん
06/03/29 10:10:24
座標は数学でしょ
507:デフォルトの名無しさん
06/03/29 17:34:46
２Ｄゲームを作っているのですが、
640×480の画面に32×32の画像をしきつめようとすると、
右端と下端がきれいにそろわず、ちょっとはみだしたりしまいます。
これって例えば何が原因で起こるのでしょうか。
508:デフォルトの名無しさん
06/03/29 17:45:27
■■■■■■■■■■■■■■■■■■■■ どういう状況なの？
■■■■■■■■■■■■■■■■■■■■ 20x15 で普通に敷き詰められるよね？
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
509:デフォルトの名無しさん
06/03/29 18:10:03
>507
まさか浮動小数点演算が介在してたりせんだろうな？
510:デフォルトの名無しさん
06/03/29 18:33:53
>>507
>例えば何が原因で
・計算間違い
・勘違い
・お門違い
511:デフォルトの名無しさん
06/03/29 18:39:10
基地（ｒｙ
512:デフォルトの名無しさん
06/03/29 20:09:24
>>507
右端は揃うんじゃﾏｲｶ？
マジレスすると、敷き詰める画像のサイズが640と480の公約数である必要がある。
513:デフォルトの名無しさん
06/03/29 20:17:27
どうしても正方形が好きならね。
514:デフォルトの名無しさん
06/03/29 22:22:06
480って32で……あれぇ？うぅー？
515:デフォルトの名無しさん
06/03/29 23:42:36
勘でいうとAdjustWindowRectしてないから
516:デフォルトの名無しさん
06/03/30 09:24:33
テクスチャで貼っててピクセルがずれている悪寒。
517:デフォルトの名無しさん
06/03/31 01:31:23
>>498
>>507
お前さ
いい加減あちこちにマルチポストするの止めろよ
もしくは解決したなら解決したと全てのスレに書いて回れ
518:デフォルトの名無しさん
06/03/31 23:45:54
地球の下側は反対に回っているってどういう意味ですか？
上も下も繋がっているんだから回る方向は同じじゃないの？
519:デフォルトの名無しさん
06/03/32 02:22:08
>>518
オーストラリアの世界地図 でググってみるといいかも
方向という単語の意味の問題だから、ここはスレ違いかも
520:デフォルトの名無しさん
06/03/32 03:18:09
北半球の人間はオーストラリアが下だろと認識している
南半球の人間はロシアが下だろと認識している
521:デフォルトの名無しさん
06/03/32 04:01:57
南半球では西から日が昇るんだよね?
522:デフォルトの名無しさん
06/03/32 11:13:29
右からだよ
523:デフォルトの名無しさん
06/04/02 08:07:38
ラグナロクみたいなオンラインゲームでステータスを表示したり、
所持アイテムを表示したりする、ウィンドウについて。
これって、どうやって実装してるんでしょうか？
ただ単に固定されてるウィンドウをつくるなら、ウィンドウのテクスチャ描いて、
それを張って、座標が決まっているのでそこに文字やらなんやら表示すれば良いと
思うんですけど、
オンラインゲームに存在するウィンドウって、windowsのウィンドウみたいに
マウスでクリックすると自由に動かすことができますよね？
これって、APIで一つ一つ子ウィンドウを作って、そこに表示させてるんでしょうか？
もしくは、他に方法とかあるんでしょうか、助言をよろしくお願いします。
524:デフォルトの名無しさん
06/04/02 08:27:33
ｽﾚﾘﾝｸ(gamedev板:543番)
525:デフォルトの名無しさん
06/04/02 08:28:25
>>523
つ[DirectX]
後スレ違い。
526:デフォルトの名無しさん
06/04/02 14:00:59
DirectXは直接的にはかんけいないだろ
527:デフォルトの名無しさん
06/04/03 04:30:12
>>523
>>517
528:デフォルトの名無しさん
06/04/07 13:54:31
>526
「数学と算数」よりは関係あるけどな。
529:デフォルトの名無しさん
06/04/24 11:15:56
>528
そうか？
530:デフォルトの名無しさん
06/04/24 13:10:28
>>528
クリッピングや論理演算は数学の範疇だと思うんだけどな
531:デフォルトの名無しさん
06/04/25 12:28:25
ﾊｧ?
DirectXはスレ違いだし、クリッピングの話なんてどこにも出てないが？
532:デフォルトの名無しさん
06/04/25 14:21:17
DirectX 描画者たち
窓の中の昴～♪
533:ルカ
06/04/26 16:58:08
sgn←これなんて読むの？だれかがシグネーチャー？って読んでたんだけどそれでオッケイなんですか？
てかこれどういう意味？符号関数ってなんですか・・・
すいませんここはプログラムとかそういうサイトっぽいけどたぶんここしか頼れる所がないので・・・
大学の数学で出てきたんですけどもしよろしければ誰か教えてください（＞＿＜）
534:デフォルトの名無しさん
06/04/26 17:02:02
もうシグネチャーでいいと思うよ。
まあ、マジレスすると、符号関数でぐぐればいくらでも説明出てくる。
535:デフォルトの名無しさん
06/04/26 18:13:04
現在、固有値および固有ベクトルを求めるプログラムを作成しているのですが
固有ベクトルの求め方で質問させていただきたいことがあります。
固有値はハウスホルダー変換と２分法で求めました。
その後、逆反復法を用いようと思っているのですが、
初期ベクトルの決め方と、繰り返し回数がわからず困っています。
どなたか教えていただけないでしょうか？
536:デフォルトの名無しさん
06/04/26 20:23:48
sgnと関数でぐぐるだけでも出てきそうなもんだが試した訳ではない。
537:デフォルトの名無しさん
06/04/27 01:34:48
sgnでググってみたら一番上がFortranの(だと思う)sgn関数の説明だった。
sgnの読み方はともかく、意味は分かるだろう。符号関数とも書いてるから
探しているものという確認も取れるし。
そしてシグネーチャーでオッケイかを知りたければシグネーチャから思いつく
綴りを全部辞書で調べろ。
538:デフォルトの名無しさん
06/04/27 02:01:54
>>535
初期ベクトルは、求めたい固有ベクトルの成分がゼロ、
もしくは、よほどゼロに近くない限りなんでもいいんじゃない？
ランダムなベクトルにするとか。
例え、偶然にも求めたい固有ベクトルの成分がゼロの初期ベクトルを
選んでしまったとしても、なんだかんだで誤差が出てくるので、なんとかなるかも。
収束は遅くなるだろうけど。
繰り返し回数は、求めたい精度に達するまでやればいいと思う。
Aを行列、aを固有値、xを求めた固有ベクトル、として、
(x,Ax)/(x,x)を計算すれば、固有値aの近似値が求まるから、
その精度を、既に求まっている固有値aと比較するとか。
539:535
06/04/27 11:26:58
なるほど、ありがとうございました。
参考にさせていただきます。
540:デフォルトの名無しさん
06/05/03 12:47:46
行列を使って非斉次の連立方程式を解くんだが、
行列の成分の
90%が0
8%が1
2%がその他
という行列を解くには
どんなアルゴリズムを使えば高速に解ける？
こういう行列に限り高速にできるとか
そんなアルゴリズムないかなと思ったんだが。
541:デフォルトの名無しさん
06/05/03 14:35:25
どこが非ゼロなのか、という情報がないとなんともいえない。
非ゼロの位置に特徴が無い場合は、行列をリスト表現して計算するくらいしか。
542:デフォルトの名無しさん
06/05/03 17:01:03
>>541
サンクス
非ゼロの位置に規則性はないわ。
あきらめとく。
543:デフォルトの名無しさん
06/05/03 18:30:42
"sparse matrix"でｸﾞｸﾞれば英語ばかりだが参考になりそうなサイトが出てくる
544:デフォルトの名無しさん
06/05/05 01:27:02
1:30～
たけしのコマネチ大学数学科
ｽﾚﾘﾝｸ(livecx板)
545:デフォルトの名無しさん
06/05/07 17:26:16
>542
非零が疎な行列だと反復法系のアルゴリズムが案外早いぞ。
546:デフォルトの名無しさん
06/05/23 08:07:13
朝倉から「コンピュータ代数ハンドブック」、定評のあるModern Computer Algebra 2nd ed.の
待望の翻訳!!、なる案内がきたが、定評と待望についてよろしく。
しかし\31,500とはずいぶんだな。
547:デフォルトの名無しさん
06/06/15 10:42:40
保守
548:デフォルトの名無しさん
06/07/16 09:13:35
保守
549:マイク ◆yrBrqfF1Ew
06/07/19 14:31:54
おい以下の問題がわからないので中卒の俺でもわかるようにPerlのコードをまじえて教えろ。
n個の区間Ii = [ai, bi](i=1,2,...,n)が与えられる時重なりを持つ区間の対を
全て列挙したい。なお[ai, bi]は実数の集合{x∈R | ai <= x <= bi}を表し、
二つの区間IiとIjが重なりを持つとは[ai, bi]∩[aj, bj]が空集合でないことを
意味する。全てのi = 1,2,..,nに対してaiとbiは整数でai<=biを満たし、
また任意のiとj(i != j)に対してai != ajを仮定する。区間のデータは
端の値aiとbiが配列で与えられており2つの数の大小比較や四則演算などの基本操作は
全てO(1)時間で可能とする。
(i) 区間対全てに対してそれぞれ重なりの有無を調べて
該当するものを列挙する方法が要する時間量を述べよ。
(ii) 重なりを持つ区間対の総数をkとする時、そのような区間対を列挙する
O(nlogn+k)時間のアルゴリズムを与えよ。
(iii) 重なりを持つ区間対を列挙するのではなく、その総数kのみを出力する
O(nlogn)時間のアルゴリズムを与えよ。
550:デフォルトの名無しさん
06/07/19 15:07:18
Perl 以外でもいいですか？
551:デフォルトの名無しさん
06/07/19 15:23:21
>550
>549に代わってお願いします。
ぜひC++で（ry
552:デフォルトの名無しさん
06/07/19 15:23:34
>>糞ｺﾃ
そういう態度だから何処へ行っても嫌われる。
553:デフォルトの名無しさん
06/07/19 21:59:50
>549
死に晒せｳﾞｫｹ
554:デフォルトの名無しさん
06/07/20 07:11:30
550 じゃないが
>>551
URLﾘﾝｸ(kansai2channeler.hp.infoseek.co.jp)
方針は以下：
(1) は自明．
(2) はソートして左から数える．ソートしたおかげで単調性が得られ，
一度交わらなくなったらそれより先を調べる必要がなくなる．
(3) は (2) でどこまで調べないといけないかを二分探索を行う．
555:デフォルトの名無しさん
06/07/20 08:09:59
拙いPerlですが。
URLﾘﾝｸ(sourcepost.sytes.net)
(i) 区間対の個数はO(nn)、重なりの有無の判定はO(1)だから、全体でO(nn)。
(ii) >>554に同じ。
(iii) 与えられた区間の始点と終点を列挙し、ソートし、各点に何個の始点・終点が重なっているか調べる。
区間のネストの数を把握しながらこの列を走査して重なりの総数を得る。
実際のコードでは始点・終点の個数の計算と最後の操作を一つのループで行っている。
556:マイク ◆yrBrqfF1Ew
06/07/20 19:23:08
>>554-555
･･･
２ｃｈネラにこんなことを言うのもあれだが
お前らすげーな。
どこの金子勇的風貌の京大院生だ？
557:デフォルトの名無しさん
06/07/20 22:30:46
大学で情報やれば大体やるぞ。
今からでも勉強してみたらどうでい。
558:デフォルトの名無しさん
06/07/30 08:15:10
二次元上の多角形の面積を求めるには外積で一発ですが、
三次元上の多面体の体積を求めるには何か方法はありませんか？
559:デフォルトの名無しさん
06/07/30 10:47:51
O = (0,0,0), A = (ax,ay,az), B = (bx,by,bz), C = (cx,cy,cz)
４面体OABCの体積は
1/6( ax by cz + ay bz cx + az bx cy - az by cx - ax bz cy - ay bx cz )
560:デフォルトの名無しさん
06/07/30 10:49:24
外積っていうより行列式に見える……。外積と行列式って関係あったっけ？
561:デフォルトの名無しさん
06/07/30 11:00:04
そりゃあるさ
562:デフォルトの名無しさん
06/07/30 12:06:07
>>558
外積
563:デフォルトの名無しさん
06/07/30 16:12:50
n次元体積を求めるにはベクトル出して行列式で一発じゃないのか？
564:デフォルトの名無しさん
06/07/30 16:14:36
3次元なら体積は (a×b)・c = det(a b c)だし
565:デフォルトの名無しさん
06/07/30 18:27:31
精度や対応できる形状をどこまでの範囲にするか決めないとな
体積を求めたい多面体が凸でない場合は４面体の和に分割するところが面倒かもよ
凸多面体でも分け方によって誤差が変わってくるからそんなに簡単ではないかもよ
566:デフォルトの名無しさん
06/08/06 16:49:43
おのおのの面に対して、原点とその面でつくられる錐体の符号付体積を足せばいいんでないの？
567:デフォルトの名無しさん
06/08/06 18:06:17
向き付けはどうすんのよ？
568:デフォルトの名無しさん
06/08/06 20:54:07
>566
対象が凸多面体で原点を多面体の内部か表面にとればそれは正しい
569:デフォルトの名無しさん
06/08/06 21:45:34
>>568
符号付き、というのはたぶん
原点から見て裏面が見えてるときはマイナスの体積、ということだろうから
そうすると原点に関して星状であるという条件はいらなくなるよ
Ａが星状というのは原点とＡの任意の点とを結ぶ線分がＡに含まれることね
570:デフォルトの名無しさん
06/08/06 23:55:55
あ、ほんとだ
571:デフォルトの名無しさん
06/08/07 12:48:29
NURBS曲面同士の衝突判定がしたいんだけど
便利なライブラリ or 公式とかある？
もしもなかったら、
衝突判定がしやすい曲面の表現方法教えてもらえませんか。
572:デフォルトの名無しさん
06/08/07 21:29:19
どんなものを表現したいのかね
573:デフォルトの名無しさん
06/08/07 21:41:17
NURBSやベジェのように、制御点を移動することで
形状を簡単に変えることができるものであればいいんですが、
球も表現したかったのでNURBSを考えていました。
574:デフォルトの名無しさん
06/08/08 11:21:05
>>571
NURBS曲面かどうか調べてないが、Newtonとか言うGAME向けの物理演算ライブラリがある
575:デフォルトの名無しさん
06/10/01 18:50:02
当選小計率 （数．数％） ＝ 当選小計 ÷ 販売額 × １００倍％
当選口指数 ｛数．数｝ ＝ 当選口数 ÷ 販売口数 × １００万倍
URLﾘﾝｸ(sea.advenbbs.net)
576:デフォルトの名無しさん
06/10/04 20:30:01
中学あたりからの学習要綱を網羅した数学ドリルサイトみたいなのないですか？
577:デフォルトの名無しさん
06/10/12 11:02:11
サンプルコードが理解出来ないレベルなら、素直に諦めて勉強し直した方がいい。
578:デフォルトの名無しさん
06/10/12 11:02:53
ｺﾞﾊﾞｰｸ
579:デフォルトの名無しさん
06/11/03 10:21:15
正射影について教えてください。
マウス（2次元）座標をゲーム画面（3次元）の座標に
変換したいんだけれど、
やり方がわからないです。
[やりたいこと]
大航海時代Onlineのように、
「キャタクターの行先（目的地）」を「地面をマウスクリックする」ことにより
指定し、その後、キャラクターが動き出す。
ということをやりたいと考えています。
580:デフォルトの名無しさん
06/11/03 10:29:29
>>579
DirectXを使っているならDirectX隔離スレへ。
その他3D表示ライブラリを使っているならそれらのスレへ。
そうでないなら自分で勝手に計算汁。
581:579
06/11/03 10:37:51
>>580
スレ違いスマソ。
DirectX隔離スレ逝ってきます。
582:デフォルトの名無しさん
06/11/23 14:10:24
グレゴリウスの定義を用いて入力した日が何曜日であるか、を表示するプログラムのアルゴリズムってどんな感じになるかわかりますか？
ちなみに
・グレゴリウスの定義
・西暦１年１月１日は月曜日
・１年３６５日だが、閏年は３６６日
・４または４００で割り切れる年は閏年
・１００で割り切れるのは閏年ではない
もしわかった方いましたらお願いします
583:デフォルトの名無しさん
06/11/23 14:15:04
>>582
C++の擬似コードだと、こんな感じかなぁ・・・
main(){Youbi 曜日; カレンダーを見る(新聞の日付, &曜日);
if (時計の時間<新聞配達時間) 曜日.進める;
return 曜日;}
584:デフォルトの名無しさん
06/11/23 15:01:22
>>582
普通はAPIを使う。
こんな所で質問してる奴が組んだ物は怖くて信用おけないから。
関係無いけど、ANSI C準拠のtime()関数系の1970/01/01 00:00:00からの通算秒で
有効桁31bitの場合は、全部閏年で換算しておｋ。
585:デフォルトの名無しさん
06/11/23 15:06:49
int zeller(int year, int month, int day)
{
if (month < 3) {
year--; month += 12;
}
return (year + year / 4 - year / 100 + year / 400 + (13 * month + 8) / 5 + day) % 7;
}
586:デフォルトの名無しさん
06/12/10 12:20:07
URLﾘﾝｸ(bal4u.dip.jp)
にあるprime-test()で多倍長整数の場合、最初に判定を間違える数って何？
587:デフォルトの名無しさん
06/12/10 13:01:52
>>582
マジレスすると西暦1年1月1日にはまだグレゴリオ暦は発明されていなかった罠。
588:デフォルトの名無しさん
06/12/10 15:35:33
>>587
遡及的に定義できる。
589:デフォルトの名無しさん
06/12/10 17:39:40
>>586
ミラーラビンテストでぐぐったら強擬素数ってのが出てきた。
あとは頑張れ。
590:デフォルトの名無しさん
06/12/10 20:48:11
>>586 底を何にするかで変わるから、一般的にパッと見つける方法はない。なので一つ一つ調べなさい。
>>589 強擬素数はフェルマーテスト
591:デフォルトの名無しさん
06/12/11 02:24:03
フェルマーテストで誤判定される数は、擬素数とかカーマイケル数。
適当いったらだめなんじゃね。
592:デフォルトの名無しさん
06/12/17 00:24:50
すいません、スレ違いかも知れませんが
他に適当なスレが見つからなくて。
私は会社の事務員なのですが、
手間のかかる手作業がありまして、パソコンにやらせたいのです。
職場の人員のスキル的に見て、自分でやるしかない状況です。
エクセルのVBAで何度かプログラムを組んだことのあるレベルです。
以下が問題の概要です。
1から100までの整数の並びがある。
そのうちランダムに20個程度の数が欠番となっている。
この状態から、残った数を並べ替え、
隣り合う数同士の差の絶対値が全体でなるべく等しくなるように、
つまり等間隔に並ぶようにしたい。
ここ一ヶ月くらい悩んでみましたが、アイディアが浮かびません。
何か既存のアルゴリズムで、応用できそうなものを紹介してもらえませんか。
593:デフォルトの名無しさん
06/12/17 00:31:57
>>592
判りにくいな。例えば、1, 2, 90, 91があったら1, 90, 2, 91に並べ替えるってことか?
#前者は隣同士の差が1, 88, 1で、後者は89, -88, 89だから後者の方が差の絶対値はほぼ均等になる。
594:デフォルトの名無しさん
06/12/17 00:39:09
>>592
総当りで調べるとかかなぁ。
「隣同士の差」の分散辺りを評価基準にして、全部の組み合わせ試して一番分散の小さいのを選ぶ。
595:592
06/12/17 00:52:03
>>593
わかりにくくてすみません。
実例を出すべきでした。
まさに、その後者のように並べたいのです。
今は人間がエクセル上で手作業であれこれと
数値を入れ替えし、納得いく結果になるまで
３日はかかってしまいます。
>>594
分散、とは統計学の用語ですか。
調べてみます。
596:デフォルトの名無しさん
06/12/17 06:08:03
80の総当りだと約10^119通り、３日じゃ終わらん。
単純に、小さい値と大きい値で２つに分けて、ソート後交互に配置じゃだめなん？
597:デフォルトの名無しさん
06/12/17 08:18:01
インターリーブ決定みたいなもんかな。
>596でよさそうだけど。
598:デフォルトの名無しさん
06/12/17 12:29:10
たとえば
1, 22, 23, 24, 25, 26, ..., 99, 100
というデータに対して>>596をすると
1, 61, 22, 62, 23, 63, 24, 64, ..., 59, 99, 60, 100
となって、「隣同士の差」のレンジは20にもなってしまう
しかしたとえば
1, 22, 42, 23, 43, 24, 44, 25, 45, ..., 39, 59, 40, 60, 41, 61,
81, 100, 80, 99, 79, 98, 78, 97, ..., 64, 84, 63, 83, 62, 82
というほうほうならレンジは2
599:デフォルトの名無しさん
06/12/17 16:18:54
>596
それならソートしただけのデータでもｲﾝｼﾞｬﾈ?
相対的な距離に関する最適解は総当りじゃないと求まらないんじゃない？
しかも「分散の小さい物」という条件だから足切りも出来ない。
600:デフォルトの名無しさん
06/12/17 22:37:12
判定関数を>>598のいうようなレンジとかにすれば枝切りもできるようになるね
601:デフォルトの名無しさん
06/12/21 21:47:49
RSA暗号を実装しようとしてて、2048ビットの鍵長の元に暗号化・復号を行っています。
整数乗算についてなんですが、2048ビット程度だと、fftかまして乗算しても恩恵はない。
と聞きました。実際のところはどうなんでしょうか？どなたか検証サイトなどご存知ではありませんか？
602:デフォルトの名無しさん
06/12/22 01:05:13
FFTの計算量評価の式から実際に求めるしかないね。
乗算が加算よりどれだけ遅いかも評価に入れなきゃいけないだろうし。
603:デフォルトの名無しさん
06/12/22 01:45:31
FFTの時間は実装レベル次第だから、実際に組んで比較するのが一番。
604:デフォルトの名無しさん
06/12/22 12:31:04
>>601
RSA の乗算＋剰余演算の場合、FFT による方法じゃなくて、
Montgomery multiplication っていう高速化手法があるよ。
検証したことはないけど、2048 ビット程度では FFT の効果はない
あるいは逆効果ってのはそうだと思う。
せめて1万ワードくらいの長さはないと。
605:デフォルトの名無しさん
06/12/22 12:37:35
ああ、あと、ワード長が2のべきになってる場合、
上位・下位半分ずつ、再帰的に計算することで
乗算回数減らす方法もある。
下位桁の乗算結果使って、上位桁の乗算をサボる。
昔、RSA のプログラムの最適化の仕事したことあるけど、
「単純に畳み込みで乗算 ＆ 引き放し法的な剰余」を
「再帰乗算 ＆ Montgomery 乗算法」に変えたら
鍵長256bitでも動作速度を250倍以上速くできた。
606:デフォルトの名無しさん
06/12/22 12:41:24
補足：
速度250倍達成は、あと、
バイト単位で多倍長整数演算してたのを、
1ワード32ビット単位に変えたのも含めてだわ。
607:デフォルトの名無しさん
06/12/27 13:26:09
おながいします(´･ω･｀)
CADのプログラムですが、CADの掲示板で質問してるけど、３次元の幾何計算の話なので、あまり回答が返ってきません。
誰か助けてぇえええ！！
URLﾘﾝｸ(www.mura.sh)
608:デフォルトの名無しさん
06/12/27 13:52:38
回答した人に失礼な話だ
609:デフォルトの名無しさん
06/12/27 15:42:41
>608
そうはいってもあそこはCADユーザーの掲示板で、
プログラム関係はほとんどレスがつかないんですよ。(´･ω･｀)
餅は餅屋ということで、こっちにも質問してみる。
610:デフォルトの名無しさん
06/12/27 16:54:00
>607,609
マルチポストうぜぇよ。
見たら解答ちゃんと出てんじゃねぇか。
それを理解出来ないお前の頭までは救えねぇよｳﾞｫｹ。
611:デフォルトの名無しさん
06/12/27 16:58:17
>610
横のソリッドがX軸方向に伸びてるならいい。
しかしそれ以外の方向ならお手上げなんだけど。
612:デフォルトの名無しさん
06/12/27 22:31:43
知らないがそれなら全体に回転行列かけて
回転してＸ軸に来るようにしたらいいんじゃね
613:デフォルトの名無しさん
07/01/06 04:41:04
保守アゲ
614:デフォルトの名無しさん
07/01/06 18:10:16
二次元のピーク検出のやり方を教えて下さい。
例えば天気図の気圧の図があったとして、
高気圧が何個あって、どの座標に中心があるかを調べたい
というような問題です。
ピークがない場合やたくさんある場合があります。
こういう計算はどういうジャンルの学問になるんでしょうか。
615:デフォルトの名無しさん
07/01/06 21:05:18
山登り法
616:デフォルトの名無しさん
07/01/06 22:14:53
>>615
いいヒントをありがとうございました。
「山登り法」というキーワードでいろいろ探せそうです。
617:デフォルトの名無しさん
07/01/07 22:14:23
多次元の最大値、最小値問題か
618:デフォルトの名無しさん
07/01/07 23:02:49
最大最小というより極値やね
619:デフォルトの名無しさん
07/01/09 11:29:24
z = f(x, y)
∂z
----
∂x
∂z
----
∂y
620:デフォルトの名無しさん
07/01/09 11:36:15
URLﾘﾝｸ(www.heisei-u.ac.jp)
621:デフォルトの名無しさん
07/01/19 14:01:03
リードソロモン符号でこんなことってできますか？
例えば1MB固定のデータファイルが n個あります。
これに任意の m個のエラー訂正用のファイルを作ります（エラー訂正というより、補完をしたい）。
たとえば、n=8、m=2 で、全部で 10個のファイルを作ります。
この 10個のファイルのうち、どれか n個だけ揃えば元の n個のファイルを復元できる、
みたいなことをやりたいのです。
ここで、n や m の取り方に、なにか制限はあったりするでしょうか？
例えば n は偶数じゃないといけないとか、m は 2のn乗じゃないといけないとか、
そういうなんらかの制限です。
できれば n も m も、その都度自由に設定したいのです。
あと、元データを復元するために必要なファイルの数は、n個だけで充分なのでしょうか。
それとも、n個よりも少し多く必要だったりするでしょうか。
今こんなことをやりたいがためにリードソロモン符号についての勉強を始めたのですが、
今結論だけを端的に知りたいです。
あるいは、他の符号化についての候補などあったら知りたいです。
よろしくです。
622:デフォルトの名無しさん
07/01/19 15:31:54
おれは学生時代、数学をきちんと勉強してなかったので後悔してる…
やっぱりプログラム、コンピュータの基盤となるのは数学だから
数学的な知識があるとないじゃ、雲泥の差だという事を、つくづく思い知らされたよ。
とりあえず、ちょっとでもスキルアップ出来る様に、勉強し直そうと思うよ…
623:デフォルトの名無しさん
07/01/19 16:56:36
最近、分数の割り算をググった事は秘密だ orz
624:デフォルトの名無しさん
07/01/19 17:04:45
>>623
最近は分数が出来ない文系の大学生も多いらしいからねぇ。
書く言うおれも、中学レベルの数学が既にあやしい…orz
625:デフォルトの名無しさん
07/01/19 17:49:25
>624
にほんごもあやしいようだね
626:デフォルトの名無しさん
07/01/19 19:49:50
>>621
Luigi Rizzo氏のHPにReed-Solomon符号を使ったデータ復元プログラムのソースが置いてあったような...
確か、制限はn+m<2^8 またはn+m<2^16。
復元もn+m個中ｎ個を受信すれば大丈夫だったはず。検索して調べてみてください。
627:624
07/01/19 21:30:15
>>625
全てを基本からやり直す必要がありそうだ…orz
628:621
07/01/20 12:46:41
>>626
情報ありがとうございます。
検索の結果それっぽいページにたどり着いたものの、いかんせん英語が不得手でして･･･
実装例は自分が行き詰ったときに見ることにして、とりあえずやりたいことはできるのだろうということはわかりました。
n+m<2^8 とか n+m<2^16 とかの制限は、多項式の係数を何ビットで計算するかによるのでしょうね。
8bit で計算すればエンディアンの問題が無くて都合が良さそうですが、できるだけ多ビットで計算した方が計算回数が
少なくて高速になりそうですね。
# 実は 8bit しか念頭に置いていませんでしたが、実験として 32bit くらいでやってみようかと思いました
629:デフォルトの名無しさん
07/01/21 23:08:41
当方受験生なのですが、数学出来ません。物理できません。化学の計算できません。
国語や英語はできます。論理性はある方だとは思います。理系の教科は勉強し
ても全然伸びませんでした。理系の科目が出来ないのは、恐らく、数学的思考が出来ていない
為だと思われます。
やってるうちに色々自分の欠点に気づき、その欠点を克服したいと思っています。
プログラミングが数学のスキルアップによいと聞き、それが欠点の克服に繋がると、
思い、やってみたいのですが、なにから始めたら良いでしょうか?
何か他にもアドバイスがあったら教えて下さいませ。
630:デフォルトの名無しさん
07/01/21 23:20:09
LOGO でタートルグラフィックでもやってみてはどうだろう
631:デフォルトの名無しさん
07/01/22 07:54:17
>>629
今受験生ってことは、もう残り1,2ヶ月？
プログラミングに費やす時間あったら、他のことにまわしたほうがいいぞ
長い目で見ればプログラミングに手を出すのも悪いことではないけど、一朝一夕ではどうかなぁ
632:デフォルトの名無しさん
07/01/22 08:05:35
プログラミング自体は文系のひとでもかなり出来るひとはいる
文系か理系かで分けるんじゃなくて
論理的に考えられるかどうかが重要
それが出来ないと何やってもだめ
633:デフォルトの名無しさん
07/01/22 08:17:07
現代のプログラミングは数学のスキルアップにはならない
何を作ろうとしてもライブラリ。
六法全書以上に分厚い資料との格闘。 記憶力や事務処理能力の方が要求される分野かもね
上で出たLOGOやpascalレベルなら別だろうけど、実用プログラミングはね
634:デフォルトの名無しさん
07/01/22 09:46:34
Mathematicaとかならともかく
現代とかどうとか関係なくそもそもプログラミングと数学は全然別物だと思うが
635:629
07/01/22 15:01:04
今年は無理なんで、来年受けようと思っています。
このスレの上の方でBASiC等が数学のスキルアップに良い。と書かれてますが、
実際どうなんでしょう?
英語は外人と喋ったり、チャット・メールをして上達しています。
数学も、こんな感じで実用面からアプローチすれば上達するかな?
と、思っています。
>>630 一度調べてみます。有難う御座います。
636:デフォルトの名無しさん
07/01/22 15:17:21
プログラマと数学者との関係をほかのもので例えるとしたら、作家と言語学者は別物というのが近い感じだと思う。
それにいくら外人と話してしても、日常会話は上達するだろうが、英語学者にはなれないだろう。
これも近い例だと思う。
637:デフォルトの名無しさん
07/01/22 15:17:47
>このスレの上の方でBASiC等が数学のスキルアップに良い。と書かれてますが、
どこのこと？
>>135ならセンター試験のプログラムの問題の話だが。
>数学も、こんな感じで実用面からアプローチすれば上達するかな?
プログラミング(の大部分)は高校数学の応用ではない。
確かに数学とプログラミングは関連しているし、考え方も共通性がない訳じゃないけど、
受験勉強としてやるには効率が悪すぎると思う。
現代国語の対策として小説の書き方を勉強するようなものじゃなかろうか。
638:デフォルトの名無しさん
07/01/22 23:20:20
今年のセンター試験のプログラミングは難化したな
片手間でやろうとした奴はそうとう焦っただろ
二分法がどういうものか分かってないと話にならない
まあ数値計算の章では普通に教科書にも出てるから
本当にやってる人ならなんでもなかったのだろうけど
統計の問題もそうだけど今年は「冷やかしお断り」っていう
センターの意図がすごい感じられたよ
639:デフォルトの名無しさん
07/01/22 23:28:19
へー
640:629
07/01/23 00:02:36
LOGOやpascalを肯定してる人も居ますが、これはどうなんでしょう?
やはり地道にがんばるのがいいんでしょうか?
何とかして数学的思考を出来る様にならなきゃいけないんですが・・・
641:デフォルトの名無しさん
07/01/23 00:13:28
今の世の中、数学のためにコンピュータがあるのではなく、
コンピュータのために数学を学ばせているような気がする。
642:デフォルトの名無しさん
07/01/23 00:13:44
「数学的思考」についてお前の思うところを述べよ
643:629
07/01/23 00:34:40
数学的思考とは、物事を順序立てて考える事であり、具体的な文字を使わず、抽象的な文字を用いる
思考です。又、余計な言葉は入れず必要な事だけを述べて行く・書き下して
行く様にする思考です。
644:629
07/01/23 00:36:24
この思考が出来ない事が、数学・物理・化学の計算を出来なくしている大きな要因
だと思っています。
645:デフォルトの名無しさん
07/01/23 01:50:51
具体的なものを使わず、じゃない。
数字や記号や式を具体的なものとして捉える力のほうが大事。
俺の持論だけど。
646:デフォルトの名無しさん
07/01/23 02:41:21
F = ma
なんのことか分からん罠
647:デフォルトの名無しさん
07/01/23 06:25:56
>>643
おまいさんは 脚の総本数 = 4 * 亀の個体数 + 2 * 鶴の個体数 から遣り直した方がいい。
648:デフォルトの名無しさん
07/01/23 07:19:32
ツルとカメと、あわせて８ひきいます。
足の数はぜんぶで、２６本。
ツルは何匹、カメは何匹？
649:デフォルトの名無しさん
07/01/23 07:21:05
URLﾘﾝｸ(www.morinogakko.com)
650:デフォルトの名無しさん
07/01/23 07:30:39
ｽﾚﾘﾝｸ(tech板:10番)
【 課題 】プログラミング演習
1)N桁の数字を各桁毎に分解し、その数字の並びの組み合わせで得られる
最大の数から最小の数を引くという演算を、結果が収束または循環するまで繰り返す。
例：
a)1111 のとき 1111 - 1111 = 0 収束
b)2006 のとき 6200 - 0026 = 6174 さらに 7641 - 1467 = 6174 収束
2)上の演算をN桁で構成される全ての数字について行い、
その桁数における収束パターンの一覧を結果として表示する。
例：
N = 4 のとき 1000 ～ 9999 までの数字で確認
3)Nを1～10程度まで可変にしてそれぞれの結果を出力。
【 形態 】Javaアプリケーション
【 GUI 】Swing(CUIでも可)
【 Ver 】1.5以上
【 期限 】月曜１７：００
【 補足 】
できれば数字の組み合わせが重複するケースや
計算の途中で既知の値(系列)に到達した場合は
それ以降の計算を省略する等して高速化を図ること
651:デフォルトの名無しさん
07/01/23 23:26:57
>629
受験数学が駄目なのか、微積分とか代数幾何とかの
教科書レベルの基本的概念を掴むところから駄目なのか
想定している大学のレベルはどれくらいか
たとえば、物理の物体の運動を数式の暗記ではなくて、
微積分の概念を用いて理解するなら数値積分の結果を
グラフにしたりするのがよくあるパターンかと思うが、
それならプログラミングなんてする必要はなく、excelで
も十分な罠。
652:デフォルトの名無しさん
07/01/24 02:21:07
>>650
でかい配列作って計算結果の一覧表を作ればおｋ
653:629
07/01/24 22:15:45
>>650 これは私宛てでしょうか?これはいったい何なんでしょうか?
>>651 受験数学は駄目ですね。これは基本概念が掴めていない所為だと
思います。志望は阪大・大阪市立の医学部です。
excelですか。一度調べてみます。有難う御座います。
654:デフォルトの名無しさん
07/01/25 10:38:57
>629
>論理性はある方だとは思います。
>数学的思考が出来ていない
矛盾してないか？
655:デフォルトの名無しさん
07/01/26 03:07:02
国語における論理（≒常識的論理）と、数学やプログラムのそれとは、
若干のズレがあるからな。
たとえば、「明日の花見は、雨ならば中止します。」と幹事が言ったとして、
「では晴れたらどうするんですか？」なんて質問をしたら、馬鹿呼ばわりされること必至。
（日常会話の「ならば」は、必要十分条件を意図して使われていることが多い。）
数学（プログラム）的には、この質問は至極自然かつ妥当なものなのだが。
656:デフォルトの名無しさん
07/01/26 10:25:47
>653
基本概念が掴めていなければ，応用もできんわな
そんなヤツに医者になってほしくない・・・
657:デフォルトの名無しさん
07/01/26 11:03:26
>>655
でも、その国語と数学のずれを認識できない奴のことを、
はたして論理性があるといえるのか？
658:デフォルトの名無しさん
07/01/26 22:47:39
それを論理の違いとして捉えているところがすでにダメダメだろ
「数学では論理を正確にするが，日常ではルーズにする」とかな
そうじゃなくて言葉の意味が違うだけ
力学で言う「仕事」が普通に言う「仕事」と違うっていうようなもの
659:629
07/01/26 23:54:56
>>654 いえ、私は、
数学的思考とは、物事を順序立てて考える事であり、具体的な文字を使わず、
抽象的な文字を用いる 思考です。又、余計な言葉は入れず必要な事だけを述べて行く・書き下して
行く様にする思考です。
こう言いました。
この私が呼んでる数学的思考という能力が劣っていても、論理性が高いという事は
ありえます。
先ず、具体的な文字を使わず、抽象的な文字を使う。という事と論理性とは
なんの関係もありません。(まぁ、厳密に調べて言ったら、多少なりとも関係
しているでしょうが)
次に、余計なことを述べず必要な物だけを述べる。と言うのも、関係有りません。
別に不必要な物を付け足しても、論理的な展開はでき、論理力が高いという事はあります。
最後に、物事を順序立てて考える。と言うことですが、これは正直よく反論出来ません。
別に物事を順序だてて考えなくても、論理力が高いと言う事はある。脳みそは複雑でバラバラな
事でも統合的に処理できている。というぐらいしか言えません。
また、頭の良い人から論理性が高い。と言われますし、実際国語の点数も良いです。
なので、私の論理性が高くない、と言う事はあまり考えるべきでないと思います。
660:629
07/01/27 00:04:12
>>655 国語における論理（≒常識的論理）と、数学やプログラムの論理は別に
考えているつもりです。
>>657 >>658 何処でどう私の事を、その国語と数学のずれを認識できない奴と、
思ったのですか?
もし、論理性の欠如から来る、ミスなら気をつけて下さい。
話がいらん、何も有益が無い方向に曲がってしまいます。
661:デフォルトの名無しさん
07/01/27 00:11:08
論理がどうこうよりも、急がば回れで、小中学生レベルの算数・数学の問題を解きまくったほうが、計算はできるようになると思う。
662:デフォルトの名無しさん
07/01/27 01:18:26
数板かどっかに行って欲しいなあ。
本人は論理的なつもりなのかも知らんが、誤った前提を元に論理も糞も無いだろうに。
663:629
07/01/27 02:19:15
>>662 うだうだ文句垂れるだけならサルでも出来る。
>>661 演算は出来ます。というか誰も演算できないなんて言ってません。
まぁ、大体このスレでの有益な意見は出たので、これで失礼します。
一部の人が怒りっぽく、荒れてきましたし。
私が言ってないことを言ったかの様に述べたり、確定してない事を確
定してると思い込むのは、恐らく、頭脳から来る頭の悪さでは無く、
人格の欠陥から来る頭の悪さでしょう。
人格の欠陥から来る頭の悪さは大概自己愛の欠落からきています。
自己が自分や他人(特に両親)から肯定されていないんですよ。
一度両親の元へ行って、愛を与えて貰いに行くのがいいと思います。
もしくは、両親との不仲を直し、互いに愛し合う関係になったほうが
良いです。それで、もしかしたら直ると思いますよ。ではノシ
664:デフォルトの名無しさん
07/01/27 02:30:09
怒りっぽくなったり荒れたりしてるのは、お前さんだけのように見えるが。
論理性があって計算もできる奴が、どうして
＞数学出来ません。物理できません。化学の計算できません。
なんてなるのかが不思議でしょうがない。
665:デフォルトの名無しさん
07/01/27 03:07:16
663見りゃ分かるとおり、コイツは勝手に論理飛躍した上で他人の人格欠陥を指摘するタイプだ。
どこも論理的じゃないから、不思議でもなんでもない。
666:デフォルトの名無しさん
07/01/27 10:23:34
>659
「数学的思考=抽象化」ということか
オレは「数学的思考=論理的思考」だと思ってたんだが
そもそも，抽象化は考える前の「単純作業」だと思ってたし
その辺，みんなはどうよ？
>別に物事を順序だてて考えなくても、論理力が高いと言う事はある
これはさすがに違うんじゃねえか？
もう見てないのかな
667:デフォルトの名無しさん
07/01/27 10:26:52
自分の考えが否定されたり，お望みの回答が得られなかったら人格攻撃をするなんてのは
アレな人の典型のようにも思える
668:デフォルトの名無しさん
07/01/27 13:53:01
>>666
呼びたいのか？ｗ
俺はもう帰ってきて欲しくないが
669:デフォルトの名無しさん
07/01/27 15:25:34
>>659
日本語の勉強を、中学校辺りから遣り直してください。
仮令国語の点数がよくても、このような文章としてなっていないものを書いているようではお郷が知れます。
670:デフォルトの名無しさん
07/01/27 15:30:12
>>668
大丈夫、この手の香具師は自意識過剰だから自分の話が出ないか気になって必ずヲチしている。
プライドが充分に低いか高いかしなければ二度と書き込まないでくれるはずなのだが……
671:デフォルトの名無しさん
07/01/27 17:05:37
最初は普通の受験生かと思ってたのに、こんな奴だったのかよｗ
672:デフォルトの名無しさん
07/01/27 21:17:51
それが論理的かどうかは別だけど、
物事をシンプルにまとめる力がないと、プログラムでは苦労するだろうな、
と、経験的に思う。
673:629
07/01/27 21:19:32
呼ばれたんで出て来ました!w
>>670さんは中々鋭いですね。私はナルシストです。
では本題に、
>>666 実際、多くの受験生が苦手としている科目が数学と物理です。
これは人が抽象的な思考をするのが苦手だからだと思います。
より具体的に説明された方が理解しやすいのは経験からも判ると思います。
ので、抽象化を'単なる'単純作業と捕らえるのは違うと思います。
しかし、出来る人からすると、単なる単純作業なんでしょうね。
>別に物事を順序だてて考えなくても、論理力が高いと言う事はある
これはさすがに違うんじゃねえか？
これについては、一例を出しますと、
ある繋がった1～10迄の事があるとします。順序だてて考えて行くとは、
1.2.3.4.5.6.7.8.9.10と一から十迄規則的に考えて行く思考法の事です。
一方の私がそうでなくてもいける、と言ったのは、
2.5.1.4.6.3.7.9.10.8.と不規則ながらもこれを頭の中で正確に関連づける思考法が
人には出来るという事です。不規則な順序で進めても、それらを正確に関連づけて
最終的に上手く統合出来る。そしてそれは確固とした論理性であり論理力である。
と言うことが言いたかった訳です。
674:デフォルトの名無しさん
07/01/27 22:01:54
>>673
>670によると、プライドが高いか低いかどっちかだそうだけどどっち?
675:デフォルトの名無しさん
07/01/27 22:05:17
＞多くの受験生が苦手としている科目が数学と物理です。
これは非常に簡単な理由がある。
こららは所謂理数系と言われる教科なのに、教えている教師は
所謂文系と言う数学物理を苦手とする傾向のある連中が多いからだ。
676:デフォルトの名無しさん
07/01/27 22:17:17
>>673
×ナルシスト
○ナルシシスト
677:デフォルトの名無しさん
07/01/27 22:52:52
高校数学物理で抽象化して捉える能力が重要だと思ったことはないなぁ
大学以降は抽象化されたものを扱う機会が増えたが、どちらかというとそういう考え方は苦手だった。
そんな自分でも高校のときは数学や物理には苦労しなかったので、
受験レベルの内容に抽象的な思考はそれほど必要ないと考えている。
>>629には数学や物理ができない原因が本当に抽象的な思考の部分にあるのかを考えてみてほしいな。
あと元々はプログラミングが役に立つかという話だったけど、
プログラミングに手を出したコストに対して得られる対価が小さいかなと思う。
もっと楽にスキルを上げる方法はありそう。
少しでも興味を持ったのなら、受験は抜きにして趣味でやってみるのもいいかもしれない。
嫌だったら辞めればいいだけの話だし。
678:デフォルトの名無しさん
07/01/27 22:57:28
青チャートやれよ
それが終わったら赤に行け
679:661
07/01/27 23:40:15
>>663
そうか、すまん。
お前が論理性があるというので、学校の数学から論理性を取り除いて残るものと言ったら計算しか思いつかなかった。
それだけのことだ。
680:デフォルトの名無しさん
07/01/28 11:10:15
ナルシス
→ナルシシズム、ナルシシスト
ナルス
→ナルシズム、ナルシスト
URLﾘﾝｸ(ja.wikipedia.org)
お前同郷かｗ 頑張れよｗ
抽象化じゃなくて一般化という用語を使えばもっと楽だろうに。
～ist の用法を一般化したらナルシシストになりそうなことは想像できるだろう。
ナルシスを知らなければしょうがないが。
681:デフォルトの名無しさん
07/01/28 14:02:55
誤用とはいえ広く浸透してる言葉の揚げ足取っても仕方ないやね
682:デフォルトの名無しさん
07/01/28 15:52:17
-istの一般的用法だと、人名の後ろにはistでなくてianをつけないか？
○machiavellian
×machiavellist
683:デフォルトの名無しさん
07/01/28 15:56:56
暮らし安心ナルシアン
684:デフォルトの名無しさん
07/01/30 05:08:18
大学の時に数値計算法の本をたくさん勉強した。
かなり役に立った。特に微積に関してはかなり深く勉強できたと思う。
でも今はCOBOLのプログラマだ・・・orz
全く数値計算法なんて役に立たねえっす。
685:629
07/01/31 05:51:03
数値計算法の本って具体的にどんなですか?
良かったら教えてくださいまし。
686:デフォルトの名無しさん
07/01/31 10:30:54
最小二乗法ライブラリはどこにありまつか？
687:デフォルトの名無しさん
07/01/31 11:03:43
>>686
URLﾘﾝｸ(www.netlib.org)
688:686
07/01/31 11:13:39
ｻﾝｸｽ>>687
LAPACKでググったら、直線のサンプルもすぐ見つかりました。
URLﾘﾝｸ(www.sip.eee.yamaguchi-u.ac.jp)
実際にやりたいのは真円のフィッティングだから、もうちょっと調査しないと。
689:デフォルトの名無しさん
07/01/31 11:33:59
円だと最小２乗では普通にやったら直線のように解けないよ。
数値解を繰り返しで求めるか
URLﾘﾝｸ(www.tensyo.com)
中心からの距離の２乗の誤差の２乗で代用すれば解けるそうだ
690:デフォルトの名無しさん
07/01/31 11:40:36
>>685
ニューメリカルレシピ・イン・シー
691:686
07/01/31 11:42:04
あ、そうですか。では、そのやり方でやります。
やりたいのは４点（真円になるか”？”）から円の中心と半径を決めたいです。
＞中心からの距離の２乗の誤差の２乗で代用
ちょっと難？
関係ないけど今回の場合、
URLﾘﾝｸ(okwave.jp)
は参考にならないような気がしました。
692:デフォルトの名無しさん
07/01/31 12:02:12
円は３点で求まる。 >>689にも ”３点を通る円” で公式が書かれているね
４点なら最小２乗を使うほどにはデータ量が多くない。
4点から３点を取り出しては公式で中心を求めて
４つ求まった中心の平均を出したらどうだろ？
その中心から４点への平均距離を求めて、それを半径の代用としたら？
693:デフォルトの名無しさん
07/01/31 12:09:42
＞＞中心からの距離の２乗の誤差の２乗で代用
＞ちょっと難？
これはデータの特性にもよるだろな。 使った感じでは十分だよ。
もっとも画像からのデータでデータ量は十分あったが
何より式一発で求まるのは非常にありがたい
694:691
07/01/31 12:19:57
>>693
ごめんなさい、書き方が悪かったです。
そのやり方が悪いという意味じゃなくて、文章難しくて、
その文章どおりにするには何をどうすれば良いのか考えなきゃという意味です。
>>692
その通りやってみたんですが、トンデモない円が出来て（３点だと上手く行く）。。。
どういう算数でやろう？
最小二乗法は？
「中心からの距離の２乗の誤差の２乗で代用」ってどういう意味？ ← 今ここ
といった感じです。
695:デフォルトの名無しさん
07/01/31 12:22:31
>>692の方法でとんでもない円になるようなデータなら
最小２乗を使ってもとんでもない結果しか出ないよ
これがデータ１００点あって１個だけ飛び離れてるような状態ならマシだけどさ
696:691
07/01/31 12:29:27
>>695
その通りなんですが、
未知データの解析なので、
「最小二乗法を使っても円にならない」ならそれもひとつの結論、
みたいな。
697:デフォルトの名無しさん
07/01/31 13:38:16
何も考えずにたった4点で最小二乗法をやるのは馬鹿げてる。
もし論文なんかで出したら一笑されるぞ。
少数観測点からデータを復元する話は GPS などの位置計測で
よくある話で、そういう場合にまともな推定をしようとすれば
観測点の状態が分かってることが必要になる。
たとえば、真の値からのずれがある確率分布に従うと仮定を
置けば、最尤推定などによって最もそれっぽい場所が出る。
正規分布に従うと仮定してよい場合は最尤推定は最小自乗法に
一致するため、多くの GPS ではこの仮定を（知ってか知らずか）置いている。
多量のデータで最小自乗法が有効なのは、中心極限定理によって
分布が正規分布に漸近するため。
698:691
07/01/31 13:40:21
あ、このページ
＞★最小２乗法による円弧推定
ってそのまんまのタイトルがあった。
３点から円を求めるのもこのページから見つけたのに気付かなかったOTL
699:691
07/01/31 13:43:33
>>697
＞何も考えずにたった4点で最小二乗法をやるのは馬鹿げてる。
＞もし論文なんかで出したら一笑されるぞ。
了解、納得でつ。
必ずしも円になるのか分からないデータの解析なので、
逆にデータ解析結果から観測点が円にならないという結論を出しますが。
ある部分は円、ある部分は円にならず従って応力受けてる、みたいな。
700:デフォルトの名無しさん
07/01/31 13:48:36
>逆にデータ解析結果から観測点が円にならないという結論を出しますが。
これはいえない。最小自乗法では円にならなくても、ほかの方法では
円になる場合があるかもしれない。
701:699
07/01/31 13:57:36
>>700
じゃあ、ほかの方法教えて。
ところで、
＞★最小２乗法による円弧推定
を読んだら、
＞この３つの方程式を解けば良い．
で終わってる。
解かないといけないのか。
702:デフォルトの名無しさん
07/01/31 14:12:46
>>701
上で述べられていた任意の三点の中心の平均や重心は有力。
また、697 で述べたように、最尤推定はひとつの別の方法。
ほかにも最小自乗や最尤推定を部分に含むが、適当な
ペナルティ関数を設定して最小化問題として解く手もある。
最尤推定やペナルティ法は、円になるかどうか分からないものの
判定なら、おそらく最小自乗を用いるよりも良い結果が得られる。
（最小自乗では、一直線上のデータは半径が（ほぼ）無限大の円と
判断されるが、通常の応用ではそれを弾きたいことが多い）
703:701
07/01/31 14:33:09
＞最尤推定やペナルティ法は
これらはライブラリにもなってなさそうだね。
実装難しそう...
704:701
07/01/31 18:42:23
＞2a( X*X/N -XX)+ 2b(X*Y/N -XY) =X(XX+YY)/N-XXX-XYY ----1)
＞2a( X*Y/N -XY)+ 2b(Y*Y/N -YY) =Y(YY+XX)/N-YYY-XXY ----2)
＞r*r=(XX+YY-2.0(a*X+b*Y) )/N+a*a+b*b; ----3)
＞この３つの方程式を解けば良い．
やっぱこれ何度読んでも分からないんだけど、
４点をどうやって代入するわけでしょうか？
解答が目の前にありながら使いこなせないorz
705:デフォルトの名無しさん
07/02/01 12:41:01
>>704
X → ∑xi XX → ∑xi*xi XXX → ∑xi*xi*xi
Y → ∑yi YY → ∑yi*yi YYY → ∑yi*yi*yi
以下略
の意じゃね？ xi, yi が配列要素で直線上にない３点が最低限必要と。
706:デフォルトの名無しさん
07/02/05 17:23:56
>>704
後は、1) 2)式から a,b２変数の連立一次方程式を解けば中心が求まるでしょ？
それを3)に代入すれば半径も求まる
707:デフォルトの名無しさん
07/02/14 10:09:26
たぶん、質問者の目的には、
4点中 ３点で求まる円から、 残った1点までの円からの距離2乗（|中心からの距離-半径|)を出して
それを４つ加算した結果、
またはそれを半径で割って正規化したもの
でいいんじゃないのか？
708:デフォルトの名無しさん
07/02/28 21:55:01
保守アゲと質問募集
709:デフォルトの名無しさん
07/03/25 15:02:16
幅 W、高さ H の楕円上に N個の点を等間隔に打とうと思いました。
for(i = 0; i < N; i++) {
int x = cx + sin(i * 2 * PI / N) * (W / 2);
int y = cy + cos(i * 2 * PI / N) * (H / 2);
plot(x, y);
}
ところがというか当然というか、これだと W = H の円の時にしか点が等間隔になりません。
どうしたら良いでしょうか？
知っていなくてはいけない公式などありましたら教えてくださいませ。
710:デフォルトの名無しさん
07/03/25 15:24:58
少し上のほうにあった，放物線を速さ一定で進むというのと同じ話だよ
711:デフォルトの名無しさん
07/03/25 15:26:04
あれ，と思ったけどそれは他スレか
712:デフォルトの名無しさん
07/03/25 15:44:51
楕円の円周の長さは 楕円積分という厄介な分野で
それを等間隔に分割というのは、その厄介な問題に挑戦するという事になる
数式では基本的に解けないから、数値積分で頑張るしかないよ
713:デフォルトの名無しさん
07/03/25 15:46:43
W > H > 0になるようにして
int x = cx + W * cos(i * 2 * PI / N);
int y = cy + H * sin(i * 2 * PI / N);
でどうでっしゃろ？
714:713
07/03/25 15:47:20
あ、無理だったごめん
715:デフォルトの名無しさん
07/03/25 16:06:21
楕円関数の近似式をｋの値を適当に決めたルックアップテーブルみたいな形で用意しとけばいいんじゃね？
それつかって円周を求めた(ここではLとおく)として
始点(仮に中心からx軸方向に幅/2移動した点としよう)から
円周にそってl ( = n×L, nは適当な実数)分移動したときの座標(x,y)ってのは楽に求まるの？
精度気にしないのならこれでいけると思うんだけど、どうでしょうか？
716:デフォルトの名無しさん
07/03/25 16:20:25
近似でいいのなら、
Nの個数を30倍くらいに増やし一度テーブルを作成して
補間で、近似値を求めるのが一番簡単だろう
717:デフォルトの名無しさん
07/03/25 16:40:25
もう少し詳しく。
そのテーブルの中味と、
何のデータを元にして、
何の式で補間するか書いてください
お願いします(-∧-)
718:709
07/03/25 16:51:16
ありゃ、思いのほか厄介な問題だったんですね・・・
>>715 円周に沿って l の座標を求める方法がよくわからないです。すみません。
>>716 で言う方法かはわかりませんが、自分ならこうやりますがどうでしょう？
1. 30N の全ての点の座標を求たテーブルを作る
2. 隣り合う点との距離を三平方の定理で求めることを全ての点について行い、円周の近時値を計算する
3. 円周を N等分した長さを求め、これを当初の N個の頂点間の距離 d とする
4. 30N のテーブル中の頂点の中から、それぞれの点の円周上の位置（始点からの距離）に近いものを選んで点を打つ
719:デフォルトの名無しさん
07/03/25 16:53:40
補間は、とりあえず線形補間とすれば、
テーブルの中身：
double ax[N*30];
double ay[N*30];
double aL[N*30];
double L=0;
ax[i] = cx + W * cos(i * 2 * PI / N);
ay[i] = cy + H * sin(i * 2 * PI / N);
if( i != 0 ) L+= hypot(cx-ocx , cy-ocy);
aL[i]=L;
ocx=cx;
ocy=cy;
のテーブルを作るでしょ？
Lを周長の代用として、
c= L*i/N で aL[k]<=c && aL[k+1]<c になる k を見つけて比例で分割して
w = ((c-aL[k]) /(aL[k+1]-aL[k]) + k) :2*Pi/N
を角度にするという感じ
720:デフォルトの名無しさん
07/03/25 17:04:31
>>718
＞円周に沿って l の座標を求める方法
楕円積分の逆関数だから、楕円関数使って求まる気がする。
もちろん、楕円関数の値は数値的に計算するものだけど。
721:719
07/03/25 17:13:07
ゴメンミス テーブルは aL だけでいい。
double aL[N*30];
double L=0;
for( i=0;i<30*N;i++}{
ax = cx + W * cos(i * 2 * PI / N);
ay = cy + H * sin(i * 2 * PI / N);
if( i != 0 ) L+= hypot(ax-oax , ay-oay);
aL[i]=L;
oax=ax;
oay=ay;
}
次のループは
for( i=0;i<N;i++}{
で while(aL[k]>c) k++;
722:719
07/03/25 17:16:37
ありゃ、元の式は 違うのか、 上の
ax = cx + sin(i * 2 * PI / N) * (W / 2);
ay = cy + cos(i * 2 * PI / N) * (H / 2);
に訂正、
続き、
w = ((c-aL[k]) /(aL[k+1]-aL[k]) + k)*2*PI/N;
int x = cx + sin(w) * (W / 2);
int y = cy + cos(w) * (H / 2);
plot(x, y);
}
723:709
07/03/25 17:30:02
>>722
丁寧にどうもありがとうございます。
よくわかりました。
>>720
まずは楕円積分というのを勉強しなくてはいけなそうです。
おいおいスキルアップしたいと思います。
>>710
よろしければそのスレを教えてもらえないでしょうか。
放物線を速さ一定で進むってのも興味あります。
724:デフォルトの名無しさん
07/03/25 17:44:59
楕円の周長なら
URLﾘﾝｸ(www.tensyo.com)
の後ろの方に計算方法が書いてあったけど
コレは今回の問題には応用できないな
725:しょう
07/03/25 20:12:39
ある学校の生徒数は １年生が全体の三分の一であり、２年生と ３年生の生徒数の比は５：６である。１年生の生徒数をa人、２年生の生徒数をb人とするとき、bをaの式で表せ。 この問題誰か解いてください
726:デフォルトの名無しさん
07/03/25 20:19:03
プログラミング関係ないな。
宿題は他をあたりなさい。
727:デフォルトの名無しさん
07/03/25 20:22:52
1) a*3=b+c
2) b*6=c*5
1' a*3*5=b*5+c*5
a*3*5=b*5+b*6=b*11
b=a*3*5/11
728:デフォルトの名無しさん
07/03/25 20:23:34
うわーなんてむずかしいもんだいなんだー
729:しょう
07/03/25 20:27:12
解らないですか？
730:デフォルトの名無しさん
07/03/25 20:28:17
ヒントください
せめて何学年まであるか、一学年何クラスかぐらいはわからないと・・・
731:しょう
07/03/25 20:31:33
高校受験の問題です。問題用紙のままカキコしました。
732:デフォルトの名無しさん
07/03/25 20:34:09
……('д`)帰れリア厨
733:デフォルトの名無しさん
07/03/25 20:41:40
SPI並の難しさだな
中学生じゃ解けないよ
大学入って掃き出し法を習うまでとっておきなさい
734:デフォルトの名無しさん
07/03/25 21:04:47
厨房でも溶けるだろ
>727間違ってる
1)a*2=b+c
だな。
で、2)はb:c=5:6から出した式ですね
あとは連立方程式
735:デフォルトの名無しさん
07/03/25 21:31:43
このスレの上の方で文系は数学から逃げ～とか言ってたけど
俺数学選択で早稲田受かったよ
736:デフォルトの名無しさん
07/03/25 21:38:14
おめでと(^O^)/
737:デフォルトの名無しさん
07/03/25 21:42:19
ありがと(^0^)／
まぁ政経だけどねｗショボいけど
738:デフォルトの名無しさん
07/03/25 22:12:11
理系で数学から逃げた俺は退学して工場労働者やってるお(^o^)
739:デフォルトの名無しさん
07/03/26 00:05:29
>>723
これだった
ｽﾚﾘﾝｸ(gamedev板)
740:デフォルトの名無しさん
07/03/26 00:10:22
>>734
それを言うなら、 a = (a + b + c) * 1/3 だろ。実際のところ、3学年しかないという前提はどこにもないわけだが。
741:デフォルトの名無しさん
07/03/26 00:26:09
>>723
楕円の周長計算、
「使えるだけでいい」という姿勢なら
単に数値計算ライブラリからソースコピればいける気が。
742:デフォルトの名無しさん
07/03/27 01:25:31
テスト
URLﾘﾝｸ(www.youtube.suppa.jp)
743:デフォルトの名無しさん
07/03/27 02:01:51
>>742はコチラへ移動しました
★☆YouTubeのCMを連続動画に☆★
//ame.x0.com/main/070327015730.html
URLﾘﾝｸ(ame.x0.com)
744:デフォルトの名無しさん
07/03/28 11:26:38
円と曲線がどれだけズレてるかの面積？みたいなものはどうやって計算すればよいですか？
745:デフォルトの名無しさん
07/03/28 11:31:26
曲線が定義済なら、数値積分すればいいじゃない
746:デフォルトの名無しさん
07/03/28 11:37:36
あ、曲線は自由曲線（ドットというか。。。）です。
747:デフォルトの名無しさん
07/03/28 11:43:15
５０レスほど戻れば、似たような話をしているよ
>>689のリンク先のページは見た？
748:デフォルトの名無しさん
07/03/28 13:40:09
>>746
その離散点使って、数値積分的なことすれば？
Σ（点から直線への距離 × 傾き）
みたいなのを。
749:デフォルトの名無しさん
07/03/28 19:00:37
曲率の分散とか？
750:デフォルトの名無しさん
07/03/28 23:08:36
>>746
>あ、曲線は自由曲線（ドットというか。。。）です。
ドットを数えろ。
751:デフォルトの名無しさん
07/03/29 08:38:46
ﾗｼﾞｬ>>750
752:デフォルトの名無しさん
07/03/29 08:57:37
円と離散点群がどれだけズレているかなら、
案１、 (点と中心との距離-半径)^2 の平均 ÷ 半径^2
案２、 ( 点と中心との距離^2 - 半径^2 )^2 の平均 ÷ 半径^4
あたりだろう。
753:デフォルトの名無しさん
07/03/29 14:03:24
「ズレている」がきちんと定義されんとなんとも。
基準円よりも広いところで円をなす点たちと
基準円上のある点にのみ集中する点たちで
どっちが「円からズレている」かは一概には言えない。
754:デフォルトの名無しさん
07/03/29 14:22:10
そこで残差の二乗の総和ですよ
755:デフォルトの名無しさん
07/03/29 15:48:26
何が「そこで」なの？
残差自乗みたいな簡単な尺度では「（概形が）円からズレてない」
みたいな位相的な構造はとても捕まえられないはずだけども。
756:デフォルトの名無しさん
07/03/29 16:07:10
>>753 それは、 ズレを最小にする円があるかどうかの問題になるんじゃないのか？
指定された円とのズレという量があればそれを最小にする半径、中心も求められるわけで・・・
それとも何か素晴らしいアイデアをお持ちで？
757:デフォルトの名無しさん
07/03/29 16:33:07
>>756
だから「ズレている」を定義してくれと言ってるんだけどな。
一点に集中してても「ズレてない」とするなら残差自乗で十分だし、
そうでないならより輪郭線抽出などの手法が要るかもしれない。
758:デフォルトの名無しさん
07/03/29 16:37:32
欠点はあるが 簡単な定義は>>752くらいしか無いだろう？
もう少しややこしくするなら、点同士がどれだけ中心からの角度で分散しているかの数値を入れるかい？
759:デフォルトの名無しさん
07/03/29 16:44:25
なんとなく想像だけど、手書きの丸と円のずれ具合を定量化したいんじゃないのかな?
だとすれば>750で充分だと思うのだけど。
#目的も判らずに数学的な意味を見出そうとしても虚しいばかりだ。
760:デフォルトの名無しさん
07/03/29 16:52:34
ドットを数えるってどうやるの？
761:ラジアンの比較で躓いてます
07/03/29 16:57:43
別スレで ｽﾚﾘﾝｸ(tech板:9番) の質問をして、最終的にこちらに誘導されてきました。
質問の回答として、以下のコードを教えてもらいました。
v1 … p1→p2 のベクトル
v2 … p2→p3 のベクトル
struct point { double x, y; };
bool isJustLeft(point v1, piont v2)
{
double corssProd = v1.x * v2.y - v1.y * v2.x; //外積
double norm1 = v1.x * v1.x + v1.y * v1.y; // |v1|^2
double norm2 = v2.x * v2.x + v2.y * v2.y; // |v2|^2
if( corssProd < 0.0 ) return false; // 時計回りはダメ
if( crossProd * crossProd == norm1 * norm2 ) return true; // 直角判定
return false;
}
でも、ベクトルの外積って通常3次元のベクトルを返しますよね。上記だと
> double corssProd = v1.x * v2.y - v1.y * v2.x; //外積
とスカラー値を返しているのですが、今一つやっている意味が判りません。
URLﾘﾝｸ(yuki.to)
ぐぐったらこんな掲示板見つけたけど、回答者の答えがイマイチ判りません。
コードを通して、ベクトルを理解したいのですが、誰か教えてもらえませんか？
762:デフォルトの名無しさん
07/03/29 17:08:28
v1,v2の外積の結果は、その２つに直角な方向ですが、
v1,v2が平面上なのでZ成分のみとなります。 だから省略したのでしょう
763:デフォルトの名無しさん
07/03/29 17:24:43
外積の定義って曖昧というか、人によって違うというか。
1つは、>>762 の言うように、3次元ベクトルの外積の考え方を使って、
(x1, y1, 0) × (x2, y2, 0) = (0, 0, x1 y2 - x2 y1) の z 成分のみを取ったもの。
もう1つは、n 次元のベクトルは、n - 1 本あれば、それらに垂直なベクトルが決まるので、
n - 1 本のベクトル → 1 本の n 次元ベクトルを求める演算を外積と呼ぶ。
こっちの流儀だち、2次元ベクトルの外積は1本→1本の演算になって、
「積」というとちょっと微妙。
まあ、x1 y2 - x2 y1 は、外積の値というか、
符号付面積、あるいは行列式よね。
3次元ベクトルの外積は、その絶対値が符号付面積になってるから、
その類推で、2次元ベクトル2本の貼る平行四辺形の符号付面積を外積と呼ぶのかも。
764:ラジアンの比較で躓いてます
07/03/29 17:36:19
>762
なるほど・・・時計回りだとZ値が下方向
反時計回りだと上方向になる性質を利用して、
Z値だけで判断すればいいという事ですね！！
765:ラジアンの比較で躓いてます
07/03/29 17:39:31
もうひとつ質問ですが、ここは何をやっているのでしょうか？
double norm1 = v1.x * v1.x + v1.y * v1.y; // |v1|^2
double norm2 = v2.x * v2.x + v2.y * v2.y; // |v2|^2
766:デフォルトの名無しさん
07/03/29 17:47:02
>>765
え、えっと、変数名もnormだし、コメントもnormだし、normを計算しているんじゃないかなあ。
767:デフォルトの名無しさん
07/03/29 17:53:54
直角であるかどうかは 内積x1*x2+y1*y2 が0になる事
(v1.x * v2.y - v1.y * v2.x)^2 -(v1.x * v1.x + v1.y * v1.y)*( v2.x * v2.x + v2.y * v2.y)
を変形してくと・・・・
って内積計算した方が計算量少ないかもしれないが
まあ、折角 外積計算したからって所じゃないのかな
768:ラジアンの比較で躓いてます
07/03/29 18:03:57
>766
一瞬正規化？とか思っちゃいましたが、「ノルム」でしたか・・・orz
769:デフォルトの名無しさん
07/03/29 18:19:33
＞「ノルム」
って日本語の数学では何だっけ？
770:デフォルトの名無しさん
07/03/29 18:27:37
世の中には日本語の数学と英語の数学があるらしい。
771:デフォルトの名無しさん
07/03/29 18:37:01
定まった和訳は無く、日本語をあててる本も特に知らないなあ。
参考までに、中国語では「范数」と書くそうな。
772:デフォルトの名無しさん
07/03/29 18:41:31
>>769
ノルムはノルムじゃない？
数学用語としてじゃなくて、一般には基準とか模範って訳すけど。
ノルムに似たので（というか、絶対値の一般化）付値ってのがあるけど、
それは英語でも valuation。
773:デフォルトの名無しさん
07/03/30 08:56:16
則（のり）じゃなかったか？
片方を９０度回転したベクトルで内積をとっても、時計回り判定はできる。
９０度回転操作を (x, y) → (-y, x) とすると、外積とコメントされた式と同じになる。
好きな方で解釈するといい。
774:デフォルトの名無しさん
07/03/30 09:20:27
内積外積を使わなくてもこれは解ける
片方のベクトルが水平(y成分が0)になるように回転変換し、
もう一方のベクトルのx成分が0なら垂直で y成分の符号を見ればいい
でも、それが内積と外積になっちゃうんだけどね
775:デフォルトの名無しさん
07/03/30 11:56:04
プログラミングの学習を先にはじめて、その必要に迫られて
その都度、数学の教養を身に着ける順序でも遅くなくねえ？
問題集をひたすら解くだけの抽象的な数学の本ばかり読んで
いると、生きることの意味がわからなくなってくるよ。
776:デフォルトの名無しさん
07/03/30 11:57:48
既出でしたら、ごめんなさい
半径１０センチの球表面の座標(XYZ)をファイルに出力したいと考えております
点の間隔は0.1センチぐらいでいいかな、と
ファイルに落とす部分は、わかっているんですが
座標を算出するアルゴリズムが、さっぱりわからなくて
お分かりになる方、御教授いただけると助かります
777:デフォルトの名無しさん
07/03/30 12:20:47
>>755
それが許されればな。
>>756
3次元空間なら
x^2 + y^2 + z^2 ＝ 半径^2
を満す実数だったと。
ざっぱになら x と y で for ループ回しながら z の値を算出すすとか。
もしくは 三次元空間の極座標
x = 半径 * sin(th1) * sin(th2)
y = 半径 * sin(th1) * cos(th2)
z = 半径 * cos(th1)
解説： URLﾘﾝｸ(hp.vector.co.jp)
で角度(th1, th2)でループまわすことも考えられる。
角度で回すのなら、球面三角法
URLﾘﾝｸ(ja.wikipedia.org)
も見おくべし。
778:776
07/03/30 12:56:17
>>>777
ありがとうございます。できたっぽいです
たすかりました
779:デフォルトの名無しさん
07/03/30 20:41:06
そのような極座標だと、目の細かい場所と粗い場所ができないか
できてもさしつかえないならいいけど、もしさしつかえあるなら
ユニバーサルメルカトル図法みたいな雰囲気で局所座標系を
とったりするとよさそうな気が
780:デフォルトの名無しさん
07/03/31 01:53:29
等間隔にするならジオデシックスフィア（日本語でなんて言うのか知らん）の頂点として出すとか
781:デフォルトの名無しさん
07/03/31 02:39:25
まあ、やっぱり極付近ほど密になるけど、↓みたいなのはある。
URLﾘﾝｸ(www.cubido.at)
ジオデシックドームの頂点求めるんだたら↓これ？
URLﾘﾝｸ(www2.tokai.or.jp)
782:デフォルトの名無しさん
07/03/31 04:54:50
正２０面体に重心細分を繰り返して得られる多面体とかではどうかね
783:デフォルトの名無しさん
07/04/11 11:07:50
自由な曲線（ベジェ曲線か、折れ線の点列）を円弧のあつまりで近似したいんですが、
ヒントはないでしょうか？
784:デフォルトの名無しさん
07/04/11 11:16:08
円弧は半径が決まっているの？ つまりフライスのようなので削るというような場合？
単純に円弧で近似したいのがどういう状況か判らないのだけど
曲線の場合は、微分が一致するように接続してゆけばいいのだけど
円弧の場合は２点と半径で求まってしまうので、
どうやっても接続点が尖がってしまう
785:デフォルトの名無しさん
07/04/12 00:12:15
円弧と線分ならどうにかそれっぽくなるかも
786:デフォルトの名無しさん
07/04/12 08:40:07
円弧の半径に制限無かったら、無限小の円弧になるだけだろ。
787:デフォルトの名無しさん
07/04/13 04:39:53
無限小の円弧の集合では曲線は近似できないのでは。
788:784
07/04/13 06:06:45
まてよ。
接続点で中心の方向が一致すればいいと解けば
接続点が尖らないように出来るか
789:デフォルトの名無しさん
07/04/13 06:16:25
>>788
）
（
）
（
こういうこと？
790:デフォルトの名無しさん
07/04/13 13:54:27
円弧も極小の長さで繋いでゆけばどんな曲線でも表現出来るし、
直線も半径を限りなく大きな値にすれば可能っぽいね。
適当に曲線から３点抽出して、３点を通る円を求めればいいんじゃね？
必要精度に達していなければ間隔を短くし、足りてれば長くして情報量を落とせばいけそう。
791:デフォルトの名無しさん
07/04/13 20:32:08
半径無限大に飛ばせば曲率ゼロだしな。
792:デフォルトの名無しさん
07/04/14 00:14:38
よくあるフォームの座標系を
0｜
―＋――→x
｜
｜
↓
y
を、
y
↑
｜
｜
―＋――→x
0 |
に変換する行列教えてください。
793:デフォルトの名無しさん
07/04/14 00:26:24
1, 0; 0, -1
794:デフォルトの名無しさん
07/04/14 15:36:47
>「network.standard-url.escape-utf8」を「false」にしてください。
>about:configで「network.standard-url.encode-utf8」を「true」にします。
上記の設定で、無事、日本語になりました。
気になるのは、IE7では、『"』⇒『"』でしたが、
URLﾘﾝｸ(forum.mozilla.gr.jp)'国際化'
URLﾘﾝｸ(forum.mozilla.gr.jp)国際化"
FireFox2では、『"』⇒『%22』になっていました、少々オシイです。
URLﾘﾝｸ(forum.mozilla.gr.jp)'国際化'
URLﾘﾝｸ(forum.mozilla.gr.jp)国際化%22
『%22』を『"』に戻す作業が残ってしまいます。
共通化としてOpera9の国際化URL設定も分かると良いと思います。
795:デフォルトの名無しさん
07/04/18 10:03:30
数学的要素が少ないプログラムの分野は何ですか？
ゲームプログラムは数学的要素満載だと思うのですが。
796:デフォルトの名無しさん
07/04/18 10:46:25
事務web系
797:デフォルトの名無しさん
07/04/18 13:57:03
>>795
どんな分野であっても、ｽﾚﾀｲが読める程度の日本語力は必要。
798:デフォルトの名無しさん
07/04/18 21:51:26
>>796
事務web系って何ですか？
799:デフォルトの名無しさん
07/04/18 22:46:08
全国の支社から出退勤のデータを収集して給与を計算するとか
ネット通販の注文を受け付けて倉庫に発送を指示するとか
800:デフォルトの名無しさん
07/04/19 02:44:55
そういうのって地味だから長く働けそうですね。
801:デフォルトの名無しさん
07/04/19 02:48:30
追記
出来るだけ長く働ける分野が良いんです
802:デフォルトの名無しさん
07/04/19 11:29:46
息子が大学に行く直前まで働ける職種でおk？
803:デフォルトの名無しさん
07/04/19 11:38:45
数学、算数はプログラムに必要か？
ｽﾚﾘﾝｸ(tech板)
を汚した奴だろ？
あちこち汚してゆくような奴が奴隷の待遇の上下を聞いたって嫌われるだけだろうに
804:デフォルトの名無しさん
07/04/19 23:48:04
就職先探すのとかは完全にすれ違い。帰れ。
805:デフォルトの名無しさん
07/04/21 06:19:39
>>802
あなたの息子を舐め舐めしますよ
806:デフォルトの名無しさん
07/04/21 09:01:30
これから数学の勉強始めてみようと思ってるんだが、プログラミングに役立つ
勉強方法が知りたい。大学受験レベルの数学問題をひたすら解くのがいいのか、
公式を理解しつつ大学レベルの数学を勉強するのがいいのか。
どんな風にすればいいか教えてください。。。。。
807:デフォルトの名無しさん
07/04/21 09:31:49
大学レベルの数学もいろいろあって、いまのあなたの
プログラミング・数学の能力と、役立てようとしている
方面によって勉強すべきものは変わる。
目的に近いところを読んでから、不足してそうな部分を
適宜補うってのが順当なごく普通の勉強の仕方。
どんなのがやりたいか言ってくれれば本などは紹介するよ。
因みに受験レベルの数学は単純計算すら覚束ないなら
仕方が無いけれど、そうでなければやる必要なしと思う。
808:デフォルトの名無しさん
07/04/21 16:00:13
> プログラミングに役立つ勉強方法が知りたい。
これは例えるなら「スポーツに役立つ練習方法を知りたい」と言ってるようなもの。
スポーツの種目によって練習方法が違うように、プログラミングの分野によって役立つ
数学の分野も違う。
# プログラミングで飯を食っていきたいなら、数学より簿記でもやった方が潰しがきくような気がしないでもない。
809:806
07/04/22 13:35:06
>>807-808
レスサンクス。
プログラミングはCの基本的な文法を理解している程度で、数学に関しては
高校レベルもｻﾊﾟｰﾘな状況です。
簿記の知識って業務系だとそんなに役に立つのか。。。。
それなら一応日商の簿記一級を持ってるんで、この知識を生かして業務系に
進んだほうがいいのかなと思ってるんだが、
「生涯現役でプログラマなんだぜ？」
な漏れとしては組込み系や制御系のほうがいいのかなと考えたり、、、
今現在はどの分野にいこうか迷ってるところです。
どうやらプログラミングに必要な数学は各分野ごとにまちまちで、漏れの
やりたいこと自体もまだ定まってないので、とりあえず高校レベルの数学を
勉強する。そしてその時の勉強法としては大学受験を目指す感じのやり方
じゃなく、基本的な公式を理解する程度でよいってことでFA？
810:デフォルトの名無しさん
07/04/22 21:45:37
人に分からないことは全く問題ないけど、
全部頼るのは良くない。
自分で決める、たぶんそれが一番大事だ。
811:デフォルトの名無しさん
07/04/22 22:39:13
＞人に分からないことは全く問題ないけど、
日本語でOK。
812:デフォルトの名無しさん
07/04/24 00:45:08
進路相談は完全に板違いなんでとっとと消えろな
毎年毎年現れる上に、人によって状況が違うし、いちいち答えてたらキリが無い
813:全知全能者
07/04/25 05:10:58
いつの世も物を言うのは「力」だ。
原始時代は「筋力」
江戸時代は「家柄」
そして現代は「金」
現代社会では金を持っている人間が強い。
革新的なパラダイムの転換が無い限りこの価値観は変わらない。
814:デフォルトの名無しさん
07/04/25 15:49:14
そして、この先は『人柄』が力となる。
815:デフォルトの名無しさん
07/04/29 00:22:43
ピクセルのフェード計算についての質問をさせてください。
実行したいのは以下の式です。
// dst[0 ～ 255]: 転送先ピクセルの色要素)
// src[0 ～ 255]: 転送元ピクセルの色要素)
// rate[0 ～ 255]; srcの比率 )
// (すべてbyte型です)
//
フェード式
= (dst * (255 - rate) + src * rate) >> 8;
~= dst + (((src - dst) * rate) >> 8);
(~= はﾆｱﾘｰｲｺｰﾙです)
これを実現するために以下のような計算方法がよく使われています。
1) short tmp = (short)src - (short)dst; //< 符号付き２バイト数に拡張します
2) tmp = (short)(tmp * rate); //< 演算結果の下位２バイトを結果として受け取ります
3) tmp = (word)tmp >> 8; //< 無符号型としてシフト()
4) byte result = (byte)(dst + (byte)tmp); //< tmpの下位１バイトのみを足し込みます
これだけ見ると変に複雑に見えますが、
実は計算にはmmxを使っていて４要素まとめて演算します。
そこで、１要素につき２バイトの範囲内で
(src - dst){-255 ～ 255} * rate{0 ～ 255}
の符号付き乗算をしないといけないため、このようなことになっています。
816:815
07/04/29 00:24:22
>>815の続きです^^
例えば、src=0, dst=255, rate=255の場合、
結果としては０が期待されますが、実際に計算してみると
1) tmp = -255 = 0xff01
2) tmp = (short)(0xff01 * 0xff) = (short)0xfff01ff = 0x01ff
3) tmp = 0x01
4) result = (byte)(0xff + 0x01) = 0x00
となり正しい結果が得られます。
また、src=200, dst=225, rate=200の場合は、
204程度が望まれますが、実際に計算してみると
1) tmp = -25 = 0xffe7
2) tmp = (short)(0xffe7 * 0xc8) = (short)0xc7ec78 = 0xec78
3) tmp = 0xec
4) result = (byte)(0xe1 + 0xec) = 0xcd = 205
とほぼ正しい結果が得られます。
上の式は有名なライブラリで使われている式でもあり、
正しいことはほとんど保証されているのですが、
これがなぜ正しいのか証明できる方はいませんか？
少なくとも自分には理解できないです。
文献でもいいです。よろしくお願いします。
817:デフォルトの名無しさん
07/04/29 00:48:54
>>815
どうしてもなにも、定義そのままの計算式じゃん。
> 1) short tmp = (short)src - (short)dst; //< 符号付き２バイト数に拡張します
(src-dst) は、-255～+255 なので、short に納まる。
> 2) tmp = (short)(tmp * rate); //< 演算結果の下位２バイトを結果として受け取ります
(src-dst)*rate は、-65280～+65280 なので、2の補数で17ビット必要→shortだと1bit足りない→後述
> 3) tmp = (word)tmp >> 8; //< 無符号型としてシフト()
((src-dst)*rate)>>8 は -255～+255なので、2の補数で9ビット必要だが、2)の時点で9ビット目の情報は落ちている
> 4) byte result = (byte)(dst + (byte)tmp);
ここで、8ビット整数で計算するのだから、上記3)の所で、結果は8ビットあれば十分。9ビット目の値は要らない。
8ビットで演算するのに「-1(0xff)を足す」のも「255(0xff)を足す」のも、まったく同じ結果になることに気付けばOK。
3)の結果で、-255 が出てきたのを、+1 として処理してもまったく問題ない。
→ -255～-1 を +1～+255 で処理しても結果は同じ。
818:デフォルトの名無しさん
07/04/29 14:31:47
>>817
おお！
８ビットまで精度を保証すればうまくいくってことね^^
どうもありがとう
819:デフォルトの名無しさん
07/04/29 18:19:52
数学の知識はそりゃあった方が良いと思うけど、
一番大事なのは物事を合理的に考える事が出来るか。
820:デフォルトの名無しさん
07/04/29 19:08:45
合理的に考えれば2chで質問するほど無駄なことは無い訳だが
821:デフォルトの名無しさん
07/04/29 19:46:50
自由な曲線（ベジェ曲線か、折れ線の点列）を円弧のあつまりで近似したいんですが、
ヒントはないでしょうか？ と、質問した者です。
点列が１個の円や直線にフィットするかどうかは、最小二乗法などで、
解決できると思います。
したがって、どういう風なグループで円弧や直線にすると、さらに最小になるか、
という問題になるような気がするのですが、そういう問題には、どのような考え方で臨めばよいでしょうか。
822:デフォルトの名無しさん
07/04/29 20:18:42
やりたい事が判らないから困ってしまうのだけど、
円は３点で決まるから、３点毎に円弧を描いても、近似と言い張れば近似になる。
823:デフォルトの名無しさん
07/04/29 20:40:03
４点あるときにどの３点を選んで円弧にするかっていう話？
824:デフォルトの名無しさん
07/04/30 01:30:11
円弧だけじゃなくて、適宜直線も使いたいってことかい？
825:デフォルトの名無しさん
07/05/02 18:06:14
たとえば、N点（３０００点とか）からなる折れ線の図形があったとして、人間は図をかけば、適度にこの部分は円弧だろうとか、直線だろうとか、あてはめることができます。
それを、精度をあたえることで、コンピュータに計算で円弧＋直線に解かせられないでしょうかね。
826:デフォルトの名無しさん
07/05/02 18:27:41
使う円弧の数に制限が無いのなら、>>822の方法で指定点での誤差ゼロで描ける
使う円弧の数を減らしたいという要求があるなら
単純に端から３点でフィッティングして、次の点が誤差の範囲内ならと処理してもいいし
最小２乗円を求めては、誤差の範囲が収まるならとやってもいいとおもう
ただし、誤差だけでは、つなぎ目がガクガクに見えるという事になる
だから、それが嫌なら、点との誤差が幾ら以内で、接続が滑らかであるというような条件を追加しないといけない
827:デフォルトの名無しさん
07/05/02 19:42:29
>>825
用途を書いてくれればもう少し具体的なアドバイスが出てくると思われ
828:デフォルトの名無しさん
07/05/02 20:19:07
>>825
滑らかじゃなくていいなら既に色々レスついてるからそれ参照。誤差ゼロ。
滑らかである必要があるなら、確かに精度(許容誤差)を与える必要はあるだろう。
滑らかって何とか言うならもっと勉強しなさい。
829:デフォルトの名無しさん
07/05/04 00:50:19
ペゾルドにかいてなかったっけ？
830:デフォルトの名無しさん
07/05/15 23:20:02
URLﾘﾝｸ(www.emit.jp)
の高速除算なのですが、
どうしてこれで正しく計算できるのか分かる方いませんか？
あと、割られる数が負数のときでも上手くいくような
つまり
X / D = Q (Qは整数)
のときに
-X / D = -Q
となるような高速除算の方法を知っている方はいますか？
831:デフォルトの名無しさん
07/05/16 01:13:49
830のリンク先はなーんか誤差が出そうなやり方だなー。
定数除算なら、誤差項を１ビット誤差以下に抑えればいいから、こういうことは出来るけど。
32bit同士の乗算がオーバーフローせずに使える場合。
X/3 = [X/3 + 2n/(3*2^32)] (∵ 2n/(3*2^32) < 1bit )
= [ (2^32 + 2)/3 * (n/2^32) ]
= 1431655766 * n / 2^32
= 1431655766 * n >> 32
832:デフォルトの名無しさん
07/05/16 01:16:49
あ、途中でXがnになっちまった……。
nが負なら最後に１足してくんなまし。
833:デフォルトの名無しさん
07/05/16 07:36:56
>>830
Web上にある情報を、どれも正しいと思うな
Q＝trunc(X/D +0.5) と計算したい筈なのに
m = 2*n-------------------1)
R = (2^m - 1) / D ---------2)
Q = (R * X + 2^n) >> m ------3) から>>m を２のベキに変更して
Q = (R * X + 2^n)/2^m Rに代入して
Q = ((2^m - 1) / D * X + 2^n)/2^m
Q = (2^m - 1) /2^m * X/ D + 2^(n-m)
~~~~~~~~~~~~~~~~ ~~~~~~~
(2^m-1)/2^m は１ではない
2^(n-m) は 0.5 ではない
これは単に使いたかった範囲で巧くいっただけだろ
834:デフォルトの名無しさん
07/05/16 11:48:36
>>831,832
ありがとうございます。
一応、プログラムを組んで実験してみたのですが、
違う結果しか得られませんでした。
自分の理解が足りていないのかもしれません。
チェック用プログラムをアップしますので、確認していただければ幸いです。
URLﾘﾝｸ(kissho.xii.jp)
>>833
整数における割り算なので、やりたいのは四捨五入ではなく切り捨て、
つまりQ＝trunc(X/D) です。
あと、使える値の範囲も>>830のページ上に書いてありますよ。
835:デフォルトの名無しさん
07/05/16 16:28:10
>>834
831に書いた方法は、
R=(2^32 + F)/Dのとき、F=D-1では、多くのDで誤差が出る。
大事なのはRが割り切れることと、F/2^32が1bitより小さい正の値になること。
これが満たされない場合は、
R=(2^33+F)/Dとか、R=(2^34+T)/D-2^32なんかを使う必要がある。
(おまけに最後の例は最後に2^32分の補正が必要)
以上のようにR の決定がそんなに単純じゃないんだー。
だから定数除算って書いた。
っていうかハッカーのたのしみって本を買えばこういうことが書いてある。
興味あるならオススメ。
836:デフォルトの名無しさん
07/05/16 17:04:05
>>830
URLﾘﾝｸ(swox.com)
> Division by Invariant Integers using Multiplication
>
> 4 Unsigned division
> 5 Signed division, quotient rounded towards 0
> 6 Signed division, quotient rounded towards -∞
837:デフォルトの名無しさん
07/05/16 18:06:56
ターゲットのチップにも寄るだろうけど、Intel系でって話なら普通に割り算命令使った方が速いんじゃね？
838:デフォルトの名無しさん
07/05/16 18:13:07
>>837
実測してみないとなんともだけど割り算の結果を遅延評価するといいかもね。
839:デフォルトの名無しさん
07/05/16 23:59:10
>>835
そうなんですか^^
ハッカーの教科書は持っていたのですがこんな話題があるとは知りませんでした。
ちょっと研究してみようと思います。
>>836
そのページを見てみたのですが、正確にやるのは結構大変なんですね^^
よく使う定数の除算のみを最適化して
あとはテーブルなりなんなりで処理することにしました。
>>837
ｽﾚﾘﾝｸ(tech板:99-番)
↑を見て欲しいのですが、x86のdiv命令は
おそらくもっとも効率化がめざましい命令の一種です。
最新の石ならばそのままが一番速いと思うのですが、
昔のCPUを考慮するとなかなかそういうわけにもいきません。
840:デフォルトの名無しさん
07/05/17 00:08:18
>>839
この手の高速化は高級言語で記述してもコンパイラが糞なら無意味だし、そもそも可読性が落ちるので歓迎されない。
そうなるとアセンブラでの記述に限定されるわけだが、肝心のターゲットが書かれていないので議論するだけ無駄。
というかそもそもこのスレ向きの話題じゃない。