プログラミングの為の数学と算数 vol.2
レス数が1000を超えています。これ以上書き込みはできません。
0001デフォルトの名無しさん
NGNG語りましょう。TIPS/Q&Aスレです。
0002デフォルトの名無しさん
NGNG∧ ∧
(,,゚Д゚) <ハァ?
/ |
(,,_/
/
0003デフォルトの名無しさん
NGNGhttp://pc5.2ch.net/test/read.cgi/tech/997150743/
プログラミングに最も必要なのは数学力
http://pc5.2ch.net/test/read.cgi/tech/1045224980/
ゲームのための物理(力学)
http://pc5.2ch.net/test/read.cgi/gamedev/1018255106/
【Collision Detection】
http://pc5.2ch.net/test/read.cgi/gamedev/1015484028/
0004デフォルトの名無しさん
NGNGフーリエ変換について教えてください
http://piza2.2ch.net/test/read.cgi?bbs=tech&key=996929748
交差判定アルゴリズム
http://piza2.2ch.net/test/read.cgi?bbs=tech&key=996157997
サウンドプログラミング
http://piza2.2ch.net/test/read.cgi?bbs=tech&key=996171508
最速の素数判定アルゴリズム
http://piza2.2ch.net/test/read.cgi?bbs=tech&key=993457354
公開鍵の高速な生成
http://piza2.2ch.net/test/read.cgi?bbs=tech&key=991992977
0005デフォルトの名無しさん
NGNGとりあえず。
グラフ理論とかも、知っておいたほうがよさげ
論理的・抽象的に考える力を養うために数学を勉強するのが最も大切なわけだが。
0007デフォルトの名無しさん
NGNG0008デフォルトの名無しさん
NGNGソフトハウスに就職した。研修始まった頃、
0〜n−1がn個であることがわかんねー、とか、16進のことで、
アルファベットがなんで数字?、とかいっていた。
んでも、配属されて開発やってんだってさ。どんな仕事やってんだろ。
0009デフォルトの名無しさん
NGNG0010デフォルトの名無しさん
NGNG0011デフォルトの名無しさん
NGNG0012デフォルトの名無しさん
NGNGネットの普及で、情報へアクセスしたり、論壇に上る特権がなくなってしまった。
よほど専門的な分野でない限り、文殊の知恵のほうが強い。
0013デフォルトの名無しさん
NGNG0014デフォルトの名無しさん
NGNG文殊は智慧をつかさどる菩薩だから、強いのは当たり前だな。
0015デフォルトの名無しさん
NGNG文系のヤツの方が英語苦手だし。
0016デフォルトの名無しさん
NGNG0017デフォルトの名無しさん
NGNG専門分野によってわけたのが理系と文系。
もともとあまり明確な区分は無いけれど、
おもにの物質的なもの・ことや論理などから研究するのが理系で、
人文をもとに研究するのが文系じゃないかなぁ…
人というのは間違いを侵すものだから、人文から結論を出す文系は
いつも間違えだらけ。
0018デフォルトの名無しさん
NGNG資料から結論を出すのが文系
だけど、資料が多くなると計算して答えを求める必要があるから、結局両方理系。
0019デフォルトの名無しさん
NGNG0020デフォルトの名無しさん
NGNG0021デフォルトの名無しさん
NGNG0022デフォルトの名無しさん
NGNGところが学生に教える教員が純文系で算数ダメなので終了してる
0023デフォルトの名無しさん
NGNG0024デフォルトの名無しさん
NGNG0025デフォルトの名無しさん
NGNGで、便器に座ってふと横を見たら
「トイレットペーパー以外の物は詰まりの原因になりますので流さないでください」
と書いてあった
だから俺はウンコは便器の横でして、尻を拭いたペーパーだけを流してしまった
そしたら店員が来て、大きなピンセットみたいな物でウンコを一粒一粒拾いながら
「おまえわざとか?」と言ってきた
俺は悪くないのにイヤミ言われてアタッマきた
0026デフォルトの名無しさん
NGNG確かに論理的には合ってるが…
0027デフォルトの名無しさん
NGNGしかし流すなと、そういうことだったんだろう。
流すなら紙だけ流せと。
はっきり書いておかない店側が悪いな。
0028マイク ◆yrBrqfF1Ew
NGNG「ウンコを便器の中にして流す」というのが通常の行為だと思ったから
わざわざ便所に行ったんだろ?
今まで便器外にそれらを配置してくるのがお前にとって当たり前のことだったか?
便所は汚物を排泄する場所であり、
例えば包茎の皮を剥いて他の男と亀頭をこすり合わせて遊ぶような場所ではないと
通常わかるなら便器の横にモノを外してOKだと思えないだろうな。
ウンコ2ちゃんネラと言われても致し方ない
0030デフォルトの名無しさん
NGNG0031デフォルトの名無しさん
NGNG0032デフォルトの名無しさん
NGNG最高。使える!
0033マイク ◆yrBrqfF1Ew
NGNGと書いてあったのであって「うんこを流すな」とそのまま明確には書いていない。
そして常識で考えればこういう表現は「汚物とトイレットペーパー以外」を意味する。
一般に便所でのこういった注意は汚物のことを当然と考えて省略している
論理としてだけただ見れば間違いだが
この社会で生活していくにはその程度の知識は必要不可欠であろう。
0034デフォルトの名無しさん
NGNGX1,X2,X3,…,X7
の…をX4,X5,X6
という風に認識するのも暗黙の了解。
0035マイク ◆yrBrqfF1Ew
NGNG@numbers = 1..100;
の .. だろうな。
0036デフォルトの名無しさん
NGNGそれはちゃんとした定義があるだろう
細かいことになるけど
0037デフォルトの名無しさん
NGNG俺も使えると思った。
0038デフォルトの名無しさん
NGNG0039デフォルトの名無しさん
NGNGしかし、うんこしてたら店員に見られるトイレもいやだ。
0040デフォルトの名無しさん
NGNG日本だと数学から逃げられるからダメすぎですよね。
他の国はどこも数学必須でしょ。
まあ、俺も逃げましたが。(だめじゃん)
0041デフォルトの名無しさん
NGNG0042デフォルトの名無しさん
NGNGこれはウンコだったんだね
0043デフォルトの名無しさん
NGNG何も出来ない奴はとりあえず文系に進むって言うのが問題なんだと思う。
ちゃんとした文系の人がかわいそう。
まあ、そんなことこのスレ的にはどうでもいいしさ、
誰か俺と数学を語り合おうぜ。
0044デフォルトの名無しさん
NGNGあなたが自分のことを俺とよぶメガネの似合う髪がさらさらのきれいなお姉さんであれば、考えんでもない
0045デフォルトの名無しさん
NGNG0046デフォルトの名無しさん
NGNG何か他にいい求めかたないですか?
0047デフォルトの名無しさん
NGNG数学板とかで聞けば?
0048デフォルトの名無しさん
NGNG収束悪いって?
時間がかかるって言うことを言っているの?
それとも精度?
0049デフォルトの名無しさん
NGNGlog(x)/2 = (x-1)/(x+1) + ((x-1)/(x+1))^3/3 + ... +((x-1)/(x+1))^(2k-1)/(2k-1)+...
とか有名ですね.他にもあると思うので適当に数値計算の本を当たってみてください.
0050デフォルトの名無しさん
NGNGわかりました明日あたり図書館行ってきます
>>48
時間(?)と収束半径ともにです
0051デフォルトの名無しさん
NGNG3 名前:1 投稿日:2001/08/07(火) 11:28
早速だけど、
非力なCPUで対数を計算させたいのです。
log(a)を求めるのに、a=b*2^n として
bを0.5〜の範囲に調整し
x = 1-b として(xは0〜0.5の範囲)
ln(b) = -(x+x^2/2+x^3/3+x^3/4+x^5/5+x^6/6+x^7/7)
さらに高速なのは
x=(1-b)/(1+b) として(xは0〜0.333の範囲)
ln(b) = -2*(x+x^3/3+x^5/5+x^7/7+x^7/9)
log(b) =-0.86859*(x+x^3/3+x^5/5+x^7/7+x^7/9)
6 名前:デフォルトの名無しさん 投稿日:2001/08/07(火) 14:05
>>3
除算がそれほど苦じゃないなら 0.5〜1にした後で
√0.5以下なら√0.5で割ってから級数展開すれば収束が早くなるよ
後でlog(√0.5)を足せばいい。
あるいは平方根(lこれはニュートン法3回程度で十分)出してから
級数展開してもいいし
10 名前:1 投稿日:2001/08/07(火) 18:55
log(a)を求める方法ですが、a=b*2^n として
bを0.5〜1の範囲に調整し
log(a)= -0.933*b*b + 2.775*b -1.8437 + log(2)*n
で求める精度が出ました
0052デフォルトの名無しさん
NGNGlog(x)/2 = (x-1)/(x+1) + ((x-1)/(x+1))^3/3 + うんこ + ((x-1)/(x+1))^(2k-1)/(2k-1) + うんこ
0053デフォルトの名無しさん
NGNGありがとうございます参考にします
0054デフォルトの名無しさん
NGNG逆双曲線関数を使うなら、CORDICも使えそうだね
0055デフォルトの名無しさん
NGNGfloatangle_s = (2.0f * PI / (float)segment);
// リンググループを生成する
for(DWORD i=0; i<ring; i++){
floatr0 = sinf((float)(i + 0) * angle_r);
floatr1 = sinf((float)(i + 1) * angle_r);
floaty0 = cosf((float)(i + 0) * angle_r);
floaty1 = cosf((float)(i + 1) * angle_r);
// 現在のリングについてセグメントグループを生成する
for(DWORD j=0; j<segment+1; j++){
floatx0 = sinf(j * angle_s);
floatz0 = cosf(j * angle_s);
floatx1 = sinf(j * angle_s);
floatz1 = cosf(j * angle_s);
// 頂点データを入力する
vert->vPos = D3DXVECTOR3(r0 * x0, y0, r0 * z0) * radius;
vert->vNormal = D3DXVECTOR3(x0, y0, z0);
vert->fU = -((float)j) / (float)segment;
vert->fV = (float)(i + 0) / (float)ring;
vert++;
vert->vPos = D3DXVECTOR3(r1 * x1, y1, r1 * z1) * radius;
vert->vNormal = D3DXVECTOR3(x1, y1, z1);
vert->fU = -((float)j) / (float)segment;
vert->fV = (float)(i + 1) / (float)ring;
vert++;
}
}
頂点シェーダで球を作るものですが、
rがsinでyがcosの理由がよく分からないのですが。
長文コピペすみません。
0056デフォルトの名無しさん
NGNG円柱座標系なんじゃないの?
0057デフォルトの名無しさん
NGNG0058デフォルトの名無しさん
NGNG0059デフォルトの名無しさん
NGNG円周率も3だし
わかりやすいDQN判別法として有効じゃん
0060デフォルトの名無しさん
NGNG0061デフォルトの名無しさん
NGNG0062デフォルトの名無しさん
NGNG0063デフォルトの名無しさん
NGNG>>61はxが小さいときに成り立つ公式
いきなり出てくる意図はわからん。
0064デフォルトの名無しさん
NGNG(1 + x)^n = 1 + n・x
log_(1 + x) (1 + x)^n = log_(1 + x) (1 + n・x)
n = log_(1 + x) n・(1/n + x)
n = n・log_(1 + x) (1/n + x)
log_(1 + x) (1/n + x) = 1
1/n + x = 1 + x
n = 1
0065デフォルトの名無しさん
NGNGまちがってはいないね。
>>61は公式というか近似式だな。
0066デフォルトの名無しさん
NGNG実際に間違えてないか?
0067デフォルトの名無しさん
NGNGガウス消去法って、どなたかご存じないですか?
そのプログラムを組みたいと思ったのです。
3x+9y-6z = -6 …α
4x+10y-2z = 16 …β
2x+3y-z = -1 …γ
ってのがあったとします。そこで、はじめ第1番目方程式の
第1番要素の係数を1にするみたいです。
つまり、3x+9y-6z = -6 を3でわって『 x+3y-2z = -2 』にするみたいですね。
そのあとβ,γの第1要素にある係数を掛けて0(ゼロ)にするらしいですが、
ゼロにするために掛ける、係数の求め方を忘れてしまいました。
適当だったら、γに-2を掛けてβと減じればゼロになるのでしょうが、法則性がわからなければ
プログラムが組めません。
法則性があったような、なかったような…
ご存知の方がいらっしゃれば、情報提供お願いいたします。。。
0068デフォルトの名無しさん
NGNGガウス消去法ってのは知らんけど、
各式に x の係数の逆数をかければ係数全部1になる。
後はただの引き算。
0069デフォルトの名無しさん
NGNG検索くらいしる。
ttp://www.fuka.info.waseda.ac.jp/~kozo/suuchi/simple_equation/simple_equation_1.html
0070デフォルトの名無しさん
NGNGガウス 消去 #include
0071デフォルトの名無しさん
NGNGどうも、アリガトウございました^^
お陰でやっとわかりました。一応検索したんですが、、、むぅ。
>>68
すごいっす。そんなやり方もあったのですね。
貴重な情報どうもです!
0072デフォルトの名無しさん
NGNGサウンドプログラミングのスレにいったんですが、
微妙に違う感じなんで、こっちで質問させてください。
最終的には2つの異なる音声信号にかけたいと思っていますが、
とりあえず最初のうちは文字列比較などの一般的な話から
理解につとめたいと考えております。
スレ違いだったら御指摘くださいませ。
0073デフォルトの名無しさん
NGNG0074デフォルトの名無しさん
NGNGようは
uzeeezou と
uzeeeezo を比較するのに 順にマッチさせて、一つ先とのマッチが大きければ先にすすめるだけ
0075デフォルトの名無しさん
NGNGhttp://hp.vector.co.jp/authors/VA013937/google.html
0076デフォルトの名無しさん
NGNG>>28
トリビアの泉だが。
0077デフォルトの名無しさん
NGNG「音声認識の基礎(上)」(NTTアドバンステクノロジ)がくどいほど詳しく載ってるよ。
0078デフォルトの名無しさん
NGNGレスありがとう参考にしてみます。
いちおうウェブ検索しましたが、
私にわかるほど噛み砕いた解説がなかったのです。
0079デフォルトの名無しさん
NGNGこの数列の規則を見つけてくれ。
(その道の専門家ならば、一目見ればすぐ分かるかも?)。
1 1 2 1 2 1 3 2 2 1 3 1 2 2 4 1 3 1 3 2 2 1 4 2 2 3 3 1 3 1 5 2 2 2 4 1 2 2 4 1
3 1 3 3 2 1 5 2 3 2 3 1 4 2 4 2 2 1 4 1 2 3 6 2 3 1 3 2 3 1 5 1 2 3 3 2 3 1 5 4
2 1 4 2 2 2 4 1 4 2 3 2 2 2 6 1 3 3 4 1 3 1 4 3 2 1 5 1 3 2 5 1 3 2 3 3 2 2 5 2
2 2 3 3 4 1 7 2 3 1 4 2 2 4 4 1 3 1 4 2 2 2 6 2 2 3 3 1 4 1 4 3 3 2 4 1 2 2 6 2
5 1 3 3 2 1 5 2 3 3 3 1 3 3 5 2 2 1 5 1 3 2 4 2 3 2 3 4 3 1 7 1 2 3 4 1 4 1 5 2
2 2 4 2 2 3 5 2 4 1 3 2 2 2 6 2 2 2 4 2 3 1 6 4 2 1 4 1 3 3 4 1 4 2 3 2 3 1 6 1
3 5 3 3 3 2 4 2 4 1 5 2 2 3 8 1 3 2 4 3 2 1 5 2 3 2 3 1 5 1 5 3 2 3 4 1 2 3 5 1
3 1 3 3 3 2 7 2 3 2 3 1 4 2 4 4 2 2 5 2 2 2 5 2 4 1 4 2 3 1 5 1 2 4 3 1 3 2 7 2
3 2 6 3 2 2 4 2 4 1 3 3 2 2 6 1 3 2 4 2 4 3 4 3 2 1 4 1 4 4 6 1 3 2 3 3 2 1 6 2
2 3 4 2 3 1 5 3 3 2 4 1 3 4 4 2 5 1 4 2 2 1 8 3 2 3 3 1 4 2 5 2 2 2 5 1 2 3 6 1
3 2 3 5 3 2 5 1 3 2 3 2 4 2 6 2 3 1 5 1 2 3 4 3 3 2 3 3 3 1 7 1 3 3 3 2 3 1 5 4
3 1 4 2 2 2 7 1 5 2 3 2 2 3 5 1 2 4 4 1 4 1 5 3 2 1 5 2 3 2 4 2 3 3 4 3 2 1 7 2
2 3 4 2 6 1 4 2 4 1 4 2 3 4 5 2 3 1 5 2 2 1 6 2 3 3 3 1 4 2 9 4 2 2 4 2 3 2 5 1
4 1 3 4 2 2 6 2 3 3 4 2 3 2 4 2 2 3 6 1 2 2 6 2 4 1 3 3 4 2 5 2 2 3 3 1 4 2 6 3
2 1 4 2 2 5 4 1 4 1 4 2 3 3 8 1 3 2 4 2 3 2 4 4 2 1 5 2 3 2 5 1 5 3 3 2 3 1 6 1
3 3 3 3 3 1 6 3 3 2 5 1 2 3 5 1 3 1 4 4 2 2 6 4 2 3 3 2 5 1 4 2 2 2 4 3 3 3 8 1
3 1 4 3 3 1 7 2 4 3 3 1 3 2 5 3 3 1 5 1 2 3 4 3 4 2 3 2 3 2 7 1 2 5 4 1 3 2 5 2
3 1 5 2 4 2 5 2 4 1 3 4 2 2 5 2 2 2 5 1 5 2 7 3 2 2 4 1 3 3 4 2 4 3 3 2 2 1 7 2
3 2 3 3 4 1 5 6 3 2 4 1 2 4 6 2 4 1 4 3 3 1 5 2 2 3 4 2 5 1 5 2 3 2 6 1 2 3 5 1
3 2 3 4 2 2 9 1 4 2 3 1 4 3 4 3 2 2 5 2 3 4 6 2 3 1 3 2 3 2 6 2 2 3 3 1 4 2 7 3
0080デフォルトの名無しさん
NGNG0081デフォルトの名無しさん
NGNG0082デフォルトの名無しさん
NGNG21101111121201211222002
01202210211212220221222
008379
NGNGもう少し待った方がいい?
間違いなく擬似乱数なんかではないよ。
008479
NGNG今は全て一桁の数だけど、
そのうち2桁や3桁の数も出てくる筈。
0085デフォルトの名無しさん
NGNG0087デフォルトの名無しさん
NGNG0が無い
0088デフォルトの名無しさん
NGNG1
1
2
1 2 1
3
2 2 1 3 1 2 2
4
1 3 1 3 2 2 1 4 2 2 3 3 1 3 1
と、初出の数nが出たら、そこから次にそのnが出る位置と反対側に次の数n+1が出てるような希ガス。
で、上で抜き出した数列の項数は順番に1,3,7,15,31,…
階差を取ると2,4,8,16,…
というのがちゃんと続いているのかどうかは2行目までしか確かめてません
008979
NGNG2=2[=>1], 3=3[=>1], 4=2*2[=>2], 5=5[=>1], 6=2*3[=>2], 7=7[=>1], 8=2*2*2[=>3],
9=3*3[=>2], 10=2*5[=>2], 11=11[=>1], 12=2*2*3[=>3], 13=13[=>1], 14=2*7[=>2],
15=3*5[=>2], 16=2*2*2*2[=>4], 17=17[=>1], 18=2*3*3[=>3], 19=19[=>1],
20=2*2*5[=>3], 21=3*7[=>2], 22=2*11[=>2], 23=23[=>1], 24=2*2*2*3[=>4],
25=5*5[=>2], 26=2*13[=>2], 27=3*3*3[=>3], 28=2*2*7[=>3], 29=29[=>1],
30=2*3*5[=>3], ...
やっぱり数の世界は奥が深いね・・・
0090デフォルトの名無しさん
NGNG0091デフォルトの名無しさん
NGNGというかこの場合、1=1[=>0]もあってしかるべきではないのか。
0092デフォルトの名無しさん
NGNG0093デフォルトの名無しさん
NGNG0094デフォルトの名無しさん
NGNG1が素数でないことと、この数列から省略されることは関係ないだろ。
数列の値が「素因数分解して出てくる素因数の数」だとするなら
1には0という正当な解が存在している以上、
省略されているのは不当な難易度の上昇だと思うんだが。
0095デフォルトの名無しさん
NGNG0096デフォルトの名無しさん
NGNG素数に馴染んでない普段の行いが悪いんだよ。多分。
0097デフォルトの名無しさん
NGNG問題としていまいち面白みに欠ける。
0098デフォルトの名無しさん
NGNG0099デフォルトの名無しさん
NGNG0100デフォルトの名無しさん
NGNG0101デフォルトの名無しさん
NGNGフェルマーの大/小定理と素数と約数の数は全て関連がある
0102デフォルトの名無しさん
NGNGお勧めの本はありますか?
0103デフォルトの名無しさん
NGNGしかしまあ、ニュートン法とかは大抵の入門書のレベルだと思うが
0104デフォルトの名無しさん
NGNGCによる理工学問題の解法
薄いし、大きいし、見やすい。
0105デフォルトの名無しさん
NGNGhttp://www.amazon.co.jp/exec/obidos/ASIN/4526036323/
これですね。
他にも何かありましたら、推薦お願いします。(Fortran, Pascal などでかかれたものはパス)
0106デフォルトの名無しさん
NGNG本じゃなくてウェブサイトがあれば教えてほしいです。
やっぱりただ読みが一番。
0107デフォルトの名無しさん
NGNGFortran も Pascal のも簡単だから覚えておいたほうがいいよ。
数値計算の本は Fortran で書かれているものが結構あるよ。
0108デフォルトの名無しさん
NGNGwebサイトではあまり良いサイトはないな。
タダだと内容までチープだよ。
0109デフォルトの名無しさん
NGNG著者/訳者名 : 本田大地
出版社名 : 秀和システム
発売予定日 : 2004年06月下旬
サイズ : B5判・460頁
販売価格 : 4,410円(税込)
最新のゲームでは当たり前に使われるようになった、物理シミュレーションプログラミングの書籍。
ゲーム中のキャラクタや背景の動きを、物理的な法則に則って表現するためのプログラミング技法を解説する。
物理法則の基本から、プログラミング方法、更に実際にゲームで利用するための方法などを収録。
付属CDには物理シミュレーションプログラミングを利用したキャラクタを操作して動かせるデモなどを収録予定
これ、どーなったんだ???
0110デフォルトの名無しさん
NGNGそれは糞本
0111デフォルトの名無しさん
NGNGなんで喪前らが論理を求めるかってーと
喪前らの顔って黄金比率が狂っているからだろうな
つまり気持ち悪がられる
よって 友達いなーーーーーーーーーーーーーーい
そうだろ?wwwwwwwwwwwww
HA-HA-HA-HA-HA-HAWWWWWWWW
0112デフォルトの名無しさん
NGNG0113デフォルトの名無しさん
NGNG俺もいない
0114デフォルトの名無しさん
NGNG0115デフォルトの名無しさん
NGNGどうすれば人間になれるでしょうか?
数学的な対象を頭の中に構成する能力は
生まれたときに既に決まっているような気がします。
0116マイク ◆yrBrqfF1Ew
NGNG何故だめなんだ?
0117デフォルトの名無しさん
NGNG黙れ基地外
0118デフォルトの名無しさん
NGNG0119デフォルトの名無しさん
NGNG0120デフォルトの名無しさん
NGNG0121デフォルトの名無しさん
NGNG発散するか振動するかだ・・・・orz
0122デフォルトの名無しさん
NGNG0123デフォルトの名無しさん
NGNG禿同。同じように語学や音楽の才能も生まれたときに決まってる気がする。
努力である程度まではカバーできてしまうところから多くの誤解が生まれるけど。
0124デフォルトの名無しさん
NGNG0125デフォルトの名無しさん
NGNGということが言いたいんだとESP。
才能があっても努力しなければ伸びないのは当然だから
乳幼児期の教育、環境が充実していることが前提での話であることは明らか。
0126デフォルトの名無しさん
NGNG0127マイク ◆yrBrqfF1Ew
NGNG0128デフォルトの名無しさん
NGNG0129デフォルトの名無しさん
NGNGシャノンも厳しいと思いますが
0130デフォルトの名無しさん
NGNG0131デフォルトの名無しさん
NGNG0132デフォルトの名無しさん
NGNG思うけど、
こういうプログラム書いてるととか、こういう言語
使ってるとこれこれの分野の数学がわかるように
なりますよ、ってのはないのですか。
一番下の娘の大学受験が来月から始まるので。
もう手遅れか・・・。
0133デフォルトの名無しさん
NGNGマセマティカ買ってやれ。
0134デフォルトの名無しさん
NGNGAHAHAHAHAHA!!!!!!
0135デフォルトの名無しさん
NGNGBASIC(VBにあらず)を勉強するとセンター試験の数学で有利だ。
卑怯なくらい有利だ。
他に使途がない上にもう手遅れっぽいけどな。
0136デフォルトの名無しさん
NGNGBASICの問題は検算ができないからお勧めしないぞ。
確かに圧倒的に簡単なのは簡単なんだけど。
0137デフォルトの名無しさん
NGNG0138デフォルトの名無しさん
NGNG情報論とか
記号論理とか
数学基礎論とか
形式言語論とか
プログラムの意味論とか
プログラミング言語やってると
すんなり入りやすい(上におもしろい)
数学はいろいろあるが
どれもこれも大学入ってからの科目で
受験には役立ちそうもないな。
0139デフォルトの名無しさん
NGNGAO入試では役に立つかもよ。
0140デフォルトの名無しさん
NGNGデジタル信号処理関係の本を連続して読んでる(谷萩・辻井先生の本あたり)んだけど
線形予測分析付近から
R(0) R(1)とか、E[f(t)^2] とかがどの本見ても具体的説明が無く出て来て困惑しています。
R(i-j)=E[s(t-i)*s(t-j)] とあり、自己相関関数を表しているとありますので
この事からR(0)は R(0) = Σs(t)^2 R(1) = Σ[s(t)*s(t-1)]付近なんだろなと思うのあけど
Q E[式] Eは期待値を意味しており・・・・という説明はあるのですが、具体的には
平均か、Σとほぼ同じ意味だろうと想像してるのですが、どちらが正しいでしょうか?
どちらでもあっても、矛盾は無い感じは感じなんだけど・・・
0141デフォルトの名無しさん
NGNG普通は平均。
どっちでも単に定数倍になっただけだからあんまりその後読み進めていくうえで支障はなさそう。
0143デフォルトの名無しさん
05/01/31 22:46:49どうすればいいでしょうか?
一次関数の場合は公式を本で見つけましたが、それ以外についてはやり方が分かりません。
ごり押しの計算でもいいので求め方を教えてください。
0144デフォルトの名無しさん
05/02/01 01:45:24漏れ、プログラマーでは無い(実験職)けど
方法は関数電卓のマニュアルに載ってるよ。
対数関数や指数関数くらいは。
もっと詳しく知りたければ大学理系学部レベルの数理統計学の教科書や線形代数学、物理学大辞典でも調べなさい。
非線形関数の取扱いとかも乗ってますよ。
0145デフォルトの名無しさん
05/02/01 07:44:28有限範囲に最小点があるなら、微分するとゼロにる点のどれかが最小だから
R'(各パラメータ群)=0
これを満たす各パラメータから、Rを計算して実際に最小の点になる点を探す
0146デフォルトの名無しさん
05/02/01 08:29:260147デフォルトの名無しさん
05/02/01 08:54:430148デフォルトの名無しさん
05/02/01 09:25:020149デフォルトの名無しさん
05/02/01 09:33:47ttp://www.eli.hokkai-s-u.ac.jp/~kikuchi/ma2/chap08.html
ttp://aoki2.si.gunma-u.ac.jp/LaTeX/nonlinear.pdf
納得いかないときはここ嫁
ttp://aoki2.si.gunma-u.ac.jp/lecture/mb-arc/arc015/399.html
0150デフォルトの名無しさん
05/02/03 18:44:36サンプル点は、x0, x1, x2, ... xm のm+1点あり、m≧nとする
近似すべき式は、以下の行列式になる
( 1 x0 x0^2 x0^3 ... x0^n ) ( a0 ) ( f(x0) )
( 1 x1 x1^2 x1^3 ... x1^n ) ( a1 ) ( f(x1) )
( 1 x2 x2^2 x2^3 ... x2^n ) ( a2 ) = ( f(x2) )
( 1 x3 x3^2 x3^3 ... x3^n ) ( a3 ) ( f(x3) )
( | ) ( | ) ( | )
( 1 xm xm^2 xm^3 ... xm^n ) ( an ) ( f(xm) )
上の式の行列を前から順に、A X Yとおくと
X A = Y
Xの転置行列をtXとすると
tX X A = tX Y
tX X は正方行列なので、通常は逆行列が存在するそれをinv(tX X)とおくと
inv(tX X) tX X A = inv(tX X) tX Y
↓
A = inv(tX X) tX Y
Aがわかれば、多項式の係数がわかるので、近似式がわかる
とまあ、こんな感じだったと思うぞ
0151デフォルトの名無しさん
05/02/03 18:46:45( 1 x0 x0^2 x0^3 ... x0^n ) ( a0 ) ( f(x0) )
( 1 x1 x1^2 x1^3 ... x1^n ) ( a1 ) ( f(x1) )
( 1 x2 x2^2 x2^3 ... x2^n ) ( a2 ) = ( f(x2) )
( 1 x3 x3^2 x3^3 ... x3^n ) ( a3 ) ( f(x3) )
( | ) ( | ) ( | )
( 1 xm xm^2 xm^3 ... xm^n ) ( an ) ( f(xm) )
0152デフォルトの名無しさん
05/02/04 00:32:11特殊な場合を除いて特異値分解して解く
つーかNumerical Recipes嫁 載ってるコードは糞らしいが
0153デフォルトの名無しさん
05/02/14 18:56:310154デフォルトの名無しさん
05/02/14 18:57:33マルチポストには教えられません。
0155デフォルトの名無しさん
05/03/18 00:31:280156デフォルトの名無しさん
05/03/18 00:34:150157デフォルトの名無しさん
05/03/18 00:47:02n + 1 のべき乗とか文字みるとクラクラしてくるんです!!
0158デフォルトの名無しさん
05/03/18 00:51:53数学ができなくても生きていけるよ。
馬鹿は無理をせずに馬鹿を貫け。
0159デフォルトの名無しさん
05/03/18 00:57:560160デフォルトの名無しさん
05/03/18 04:16:180161デフォルトの名無しさん
05/03/19 12:18:57ぐぐれ。
0162デフォルトの名無しさん
05/03/20 18:57:430163デフォルトの名無しさん
2005/03/25(金) 09:25:380164デフォルトの名無しさん
2005/03/28(月) 06:03:35俺の脳のスタックが足りねえこと痛感よ
0165デフォルトの名無しさん
2005/03/31(木) 08:56:580166デフォルトの名無しさん
2005/04/11(月) 19:27:12どうなんでしょう
0167デフォルトの名無しさん
2005/04/11(月) 19:34:52sage
0168デフォルトの名無しさん
2005/04/11(月) 19:35:26http://133.86.144.64/
http://133.86.144.64/~ss.jpg
http://dawn3/
http://dawn3/~ss.jpg
0169デフォルトの名無しさん
2005/04/11(月) 19:36:200170デフォルトの名無しさん
2005/04/11(月) 20:45:27atanの定義はイクナイ
0171デフォルトの名無しさん
2005/04/12(火) 05:18:260172デフォルトの名無しさん
2005/04/12(火) 11:00:18C言語のライブラリに atan は要らないってことよね?
それなら禿同。
0173デフォルトの名無しさん
2005/04/12(火) 19:00:030174デフォルトの名無しさん
2005/04/12(火) 22:03:00まあ、それは俺も。
っていうかね、紙の上で微分方程式解いたりするよりも暗算の方が苦手。
0175デフォルトの名無しさん
2005/04/12(火) 22:54:07もまえらこれやってみろ。安産力付くぞ
ttp://r_akky.at.infoseek.co.jp/mental.html
漏れは10000越えると死ぬ
0176デフォルトの名無しさん
2005/04/14(木) 13:33:16全頂点の入次数が1以下であることを確認するよりも効率的な方法はありますか?
0177デフォルトの名無しさん
2005/04/14(木) 15:32:31の意味が判らないけど、これで十分条件になるの?
A→B, B→C, C→A みたいな循環構造が木と判定されない?
0178デフォルトの名無しさん
2005/04/27(水) 19:42:54なんか言えよ雑魚
が
0179デフォルトの名無しさん
2005/04/28(木) 01:01:100180デフォルトの名無しさん
2005/04/28(木) 03:17:27なぜ曜日が求まるのか、詳しく教えていただけませんか・・・?
0181デフォルトの名無しさん
2005/04/28(木) 03:25:56とりあえずぐぐってみて一番正確かつ詳しく説明してそうなページを探し、
その上でどこがわからないかを文章にしてみたら?
0182デフォルトの名無しさん
2005/04/28(木) 05:14:03http://www.amazon.co.jp/exec/obidos/ASIN/4797329734/ref=ase_2217800-22/249-3298830-5914740
0183デフォルトの名無しさん
2005/04/28(木) 06:53:53ただし、7で割った余りのみが問題になるので、何日たったかを正確に求める必要は無く
7の倍数だけずれた日数を求められれば、それを7で割った余りで曜日を求めるには十分である
一年は365日なので、うるう年を無視するとn年間で日数は 365 * n 日
しかし途中にうるう年が入るとややこしくなる、うるう年で増える一日を
一年の最後に持っていくと都合がいい
そうしないと、うるう年となる年の曜日を求める際に、年の途中ではいる2月29日以降は
一日だけ日数を増やす必要がでてくる
したがって、1月は前年の13月、2月は前年の14月として扱う
この変形グレゴリウス暦の西暦0年3月1日からY年間で何日になるかは、
うるう年を考慮に入れて、以下のようになる
365 * Y + [Y / 4] - [Y / 100] + [Y / 400]
= 364 * Y + Y + [Y / 4] - [Y / 100] + [Y / 400]
曜日を考える分には7の倍数は無視できる
364 * Y = 7 * 52 * Yで、7の倍数だから無視すると
Y + [Y / 4] - [Y / 100] + [Y / 400]
0184デフォルトの名無しさん
2005/04/28(木) 06:55:03100 * c + y + [(100 * c + y) / 4] - [(100 * c + y) / 100] + [(100 * c + y) / 400]
= 100 * c + y + 100/4 * c + [y / 4] - 100/100 * c - [y / 100] + [100/400 * c + y / 400]
= (100 + 25 - 1) * c + y + [y / 4] - [y / 100] + [c / 4 + y / 400]
ここで、0 <= y < 100より、[y / 100] = 0, [c / 4 + y / 400] = [c / 4]なので
= 124 * c + [c / 4] + y + [y / 4]
= 119 * c + 5 * c + [c / 4] + y + [y / 4]
119 * c = 7 * 17 * cで、7の倍数なので無視して
5 * c + [c / 4] + y + [y / 4]
= [21/4 * c] + [5/4 * y]
[21/4 * c] + [5/4 * y]を西暦2005年3月1日(火曜日)について計算すると、111で
111を7で割った余りは、6だから、6が火曜日となる
ツェラーの公式に合わせて6が土曜日になるには、3を足せばよい
[21/4 * c] + [5/4 * y] + 3
これを、7で割った余りから、各年の3月1日の曜日は求められる
0185デフォルトの名無しさん
2005/04/28(木) 06:55:56m月1日(3 <= m <= 14)が3月1日から何日たっているかを求めればよい
その日数を f(m) であらわすとすると
[21/4 * c] + [5/4 * y] + 3 + f(m) を求め、それを7で割った余りから各月の一日の曜日は求められる
なお、うるう年に増える一日は年の最後の日であるので、f(m) の値ははうるう年の影響を受けない
f(3) = 0 は自明、3月は31日あるので、f(4) = f(3) + 31 = 31 だが
例によって7の倍数分は無視しても曜日には影響は無いので28引いて、f(4) = 3 としても差し支えない
つまり、f(m) を求めるに際しては、3月は3日しかないと考えても結果は変わらない
この考えに従って、各月の日数を3月から順に列挙すると
3, 2, 3, 2, 3, 3, 2, 3, 2, 3, 3, 0 or 1となる
1日(ついたち)の曜日を考える上では最後の月である14月が何日であるかは影響が無い
したがってこの数列は 3, 2, 3, 2, 3, 3, 2, 3, 2, 3, 3 と考えてよい
これをよく見ると、3,2,3,2,3というパターンの単純な繰り返しとなっている
このことから、5ヶ月は13日(3+2+3+2+3=13)であると考えることができるのである
したがってこのことから、f(m) は、[13/5 * m]に定数項を足したもので表現できそうである
[13/5 * m]の値を m = 3 から m = 14 について列挙してみよう
7, 10, 13, 15, 18, 20, 23, 26, 28, 31, 33, 36
この数列の隣同士の値の差が各月の日数と合致していればよい
実際に求めると 3, 3, 2, 3, 2, 3, 3, 2, 3, 2, 3
これは、上述の数列とよく似ているが一つ横にずれている(残念)
0186デフォルトの名無しさん
2005/04/28(木) 06:57:32[13/5 * m]ではなく[13/5 * (m + 1)]の値を m = 3 から m = 14 について列挙してみよう
まあこれは、再計算する必要も無く、さっきの数列が左に一つずれただけだ
10, 13, 15, 18, 20, 23, 26, 28, 31, 33, 36, 39
この数列の隣同士の値の差は、やはり左に一つずれて
3, 2, 3, 2, 3, 3, 2, 3, 2, 3, 3となり、これで目的の数列と一致する
つまり、 f(m) は、[13/5 * (m + 1)]に定数項を足したもので表現できる
ここで[13/5 * (m + 1)]の値を m = 3 で求めると、[13/5 * 4] = 10
7の倍数は無視できるので、この値は 3 と同じとみなせる
これは目的の値 f(3) = 0 より 3 多い
したがって f(m) = [13/5 * (m + 1] - 3 となる
各月の1日(ついたち)の曜日は
[21/4 * c] + [5/4 * y] + 3 + f(m)
= [21/4 * c] + [5/4 * y] + [13/5 * (m + 1]
を7で割った余りから求められる
2日(ふつか)以降の曜日については、一日ごとに曜日が一個ずれるから
m月d日について求めるときは、上式に d - 1 を加えればよい
[21/4 * c] + [5/4 * y] + [13/5 * (m + 1] + d - 1
この値を7で割った余りが、100 * c + y年m月d日の曜日となる
0187デフォルトの名無しさん
2005/04/28(木) 06:58:30しかしこれであってるのだろうか?
0188マイク ◆yrBrqfF1Ew
2005/04/28(木) 15:56:08nが自然数でn^2+1が素数になる数がいくつあるのか中卒の僕に優しく教えてくださいペコリー
0189デフォルトの名無しさん
2005/04/28(木) 18:57:37http://ja.wikipedia.org/wiki/%E7%B4%A0%E6%95%B0
ここの「素数に関連する未解決の問題」の欄を見るといい。
0190デフォルトの名無しさん
2005/04/28(木) 20:55:41やるな・・・
0191デフォルトの名無しさん
2005/04/29(金) 13:33:400193デフォルトの名無しさん
2005/05/04(水) 02:14:140194デフォルトの名無しさん
2005/05/04(水) 02:43:37「トイレットペーパー以外の物は詰まりの原因になりますので流さないでください」
と書いてあったが、俺はウンコは流しても問題ないよなと思った
だから俺はウンコをし、尻を拭いたペーパーも流してしまった
そしたらトイレが詰まってしまい、怖いお兄さんが飛んできて後始末をしながら
「トイレットペーパー以外流さないで下さいって書いてあるだろ」と言ってきた
結局殴られる変わりに、一万円払って何とか逃げ出しました
以上、ぼったくりトイレは怖いと思い知った一日でした
0195デフォルトの名無しさん
2005/05/11(水) 15:41:210196デフォルトの名無しさん
2005/05/30(月) 13:26:46> >>188
> http://ja.wikipedia.org/wiki/%E7%B4%A0%E6%95%B0
> ここの「素数に関連する未解決の問題」の欄を見るといい。
のぞいてみたら、紛らわしい記述があった。
>素数は無限に存在する。
>
> * エウクレイデスによる証明
>
> 背理法による。
エウクレイデス(ユークリッド)のは背理法ではないよ。
0197デフォルトの名無しさん
2005/05/30(月) 17:15:51何故、背理法ではないと思うの?
0198デフォルトの名無しさん
2005/05/30(月) 17:35:050199デフォルトの名無しさん
2005/05/30(月) 17:45:110200デフォルトの名無しさん
2005/05/30(月) 20:55:230201デフォルトの名無しさん
2005/05/30(月) 21:27:480202デフォルトの名無しさん
2005/05/31(火) 04:39:00矛盾が生じる。
0203デフォルトの名無しさん
2005/05/31(火) 05:11:54実は単なる対偶法、ということはよくあるけどね
0204デフォルトの名無しさん
2005/05/31(火) 08:16:120205デフォルトの名無しさん
2005/05/31(火) 13:42:44あるわけで、2からPまでの全部の積に1を加えた数を考えて
みよう。これは明らかにPより大きい数だが、Pを整除するには?
0206デフォルトの名無しさん
2005/05/31(火) 17:21:161以外ではできない あれ、これは素数ではないか
→すると仮定に矛盾する。
0207デフォルトの名無しさん
2005/05/31(火) 17:48:36君は原論第9巻の20を読んだことはないということだ。そこには
「素数の個数はいかなる定められた素数の個数よりも多い」とある。
0208デフォルトの名無しさん
2005/05/31(火) 18:13:33メール:[email protected]
>先週はフィリピンで旧日本兵の方が発見されたとかいう未確認情報も
>流れて結構話題になっていたりしていますが、それにしてもフィリピンと
>言えば前の戦争で50万人も日本人が死んでいるんですよね。
>しかもこれ、ほとんどが飢え死にだそうで。話題になっている第三十師団も
>二週間分の食料しか持たずにジャングルを50日間行軍して6.7割の人が
>餓死したのだそうですが(インパール作戦みたいなものですね)、
>いやはやなんと馬鹿みたいな戦争をしたものです。
>しかも結局何も得るものも無かったわけですから、死んだ兵士はまさに犬死に
>としか言いようがありませんね。まったく二度と戦争なんて馬鹿げたことはしないで
>欲しいものです。とりあえず日本には九条があるので大丈夫かとは思いますけれどね。
0209デフォルトの名無しさん
2005/05/31(火) 18:53:02ここに貼ってどうする
政治経済方面に貼りまくるんだ
0210デフォルトの名無しさん
2005/05/31(火) 20:12:29だね。
0211デフォルトの名無しさん
2005/05/31(火) 20:13:010212デフォルトの名無しさん
2005/06/01(水) 04:14:300213デフォルトの名無しさん
2005/06/02(木) 05:35:26原論は昔読んだ事があるんだが・・・
ごめん、あなたの方が正しい気がしてきた。もう一度勉強しなおしてくる。
0214デフォルトの名無しさん
2005/06/06(月) 09:52:300215デフォルトの名無しさん
2005/06/06(月) 09:57:55話がでてるけど余談
エルデシュ(Erdoes) はこの証明を
"The BOOK" すなわち彼が信仰し存在することを
前提としている「究極のエレガントな証明集」の
代表例といっているそうだ
0216デフォルトの名無しさん
2005/06/06(月) 10:03:28この証明,何かだまされたような気がする人がいるのは
背理法,もっというと排中律を暗黙の仮定としているからだろう.
実際にかつて,排中律を使わずに数学を構成しよう,
という立場の人があらわれた.
より詳しい話は数学の歴史の本を参照してください.
0217デフォルトの名無しさん
2005/06/06(月) 16:02:550218デフォルトの名無しさん
2005/06/06(月) 20:05:18つまり廃中立を認めない立場では
素数は有限しかないことないのか
0219デフォルトの名無しさん
2005/06/06(月) 20:48:03背理法使わずに証明できるよ
0220デフォルトの名無しさん
2005/06/06(月) 21:19:05ひたすら直観主義論理で算術を構成してくっていうやつで
その最後の問題が素数の非有限性だったな
俺は証明できなかったけど
0221デフォルトの名無しさん
2005/06/06(月) 23:22:03てゆーかあれで証明なのかどうか理解できにくい
0222デフォルトの名無しさん
2005/06/06(月) 23:56:45あれは自然数の公理のうちに入ってる
0223デフォルトの名無しさん
2005/06/07(火) 22:46:31それを公理に加えちまうってのが、
いかにも後付け的でインチキくせーってことでしょ。
0224デフォルトの名無しさん
2005/06/07(火) 23:32:50まあ、無限集合を定義するのにどの道必要だし。
0225デフォルトの名無しさん
2005/06/08(水) 11:36:090226デフォルトの名無しさん
2005/06/08(水) 21:57:09それはつまり、「実数全体からなる集合は、
自然数全体からなる集合と1対1対応が取れる」
ということですか?
0227デフォルトの名無しさん
2005/06/08(水) 22:15:030228デフォルトの名無しさん
2005/06/08(水) 22:24:40流石にスレ違いのような
0229デフォルトの名無しさん
2005/06/09(木) 09:13:05それはウソ
実数は可算ではないことが証明できる
有名な Cantor の対角線論法
0230デフォルトの名無しさん
2005/06/09(木) 09:17:40可算の濃度と実数の濃度に間の
濃度を持つ集合は存在するか?
0231デフォルトの名無しさん
2005/06/09(木) 10:48:21それは「どっちでもいい」で解決したはず
0232デフォルトの名無しさん
2005/06/09(木) 14:27:31上に同じ
0233デフォルトの名無しさん
2005/06/10(金) 00:02:39正確には、ZFC公理系の範囲では否定も肯定もできない。
もうちょっと公理を足せば否定も肯定もできる。
0234デフォルトの名無しさん
2005/06/14(火) 10:17:37ZF公理系の無矛盾性を証明せよ
0235デフォルトの名無しさん
2005/06/14(火) 20:52:41ネタフリ乙!
ttp://ja.wikipedia.org/wiki/%E3%82%AF%E3%83%AB%E3%83%88%E3%83%BB%E3%82%B2%E3%83%BC%E3%83%87%E3%83%AB
0236デフォルトの名無しさん
2005/06/15(水) 00:34:45(2)仮定により矛盾が導かれる。
(3)故に仮定(1)は否定される。
0237デフォルトの名無しさん
2005/07/16(土) 12:47:120238デフォルトの名無しさん
2005/08/15(月) 13:33:360239デフォルトの名無しさん
2005/08/20(土) 13:56:07まさに再帰アルゴリズムだよね。
0240デフォルトの名無しさん
2005/08/20(土) 13:59:47背理法
0241デフォルトの名無しさん
2005/08/20(土) 14:26:530242デフォルトの名無しさん
2005/08/20(土) 14:44:34Googleの「日本海」の表記が韓国からの抗議により「東海」に変更
↓↓↓詳しくはこちら↓↓↓
【日本海】(;^ω^)VIP vs VANKlt;`∀´*gt;part8【海戦】
http://ex11.2ch.net/test/read.cgi/news4vip/1124499103/l50
今日の17時、20時、23時に 何 か を 「 撃 ち 」 ま す 。
抗議メール等、ご協力お願いいたします。
~~
↑
重要
0243デフォルトの名無しさん
2005/08/20(土) 16:06:16大学に入れば勝ち
それだけで
仕事で使えもしないのに
おおいばりできるからな
0244デフォルトの名無しさん
2005/08/23(火) 09:45:03違うだろ
それでやっと「スタートライン」だろ
0245デフォルトの名無しさん
2005/08/24(水) 00:58:46http://www.amazon.co.jp/exec/obidos/ASIN/4320026683/
一応、既出ではないよね?
コンピュータの名著100にも選ばれた超定番。
スタンフォード大学で使われてる教科書だけど、
ハノイの塔とかイメージしやすいところから始まって高度なところまで
自然に導いてくれるから、思ったほど敷居は高くない。
日本の教科書みたいにトップダウン式の記述ではなくかなり親切に書かれてるから、
腰をすえて読めばきっと満足できるはず。
値段的にアレだけど…。
0246デフォルトの名無しさん
2005/08/24(水) 09:43:17URL長すぎ
もっと短くしる
0247デフォルトの名無しさん
2005/08/24(水) 12:35:21http://www.amazon.co.jp/o/ASIN/4320026683
と出来る。
まだ長い?
0248デフォルトの名無しさん
2005/08/24(水) 12:48:08こう?
0249デフォルトの名無しさん
2005/08/24(水) 16:24:06www.amazon.co.jp/o/ASIN/4320026683
0250デフォルトの名無しさん
2005/08/24(水) 16:29:12">>245"
0251デフォルトの名無しさん
2005/08/24(水) 22:57:500252デフォルトの名無しさん
2005/08/25(木) 10:14:14> コンピュータの数学 グレアム著
この本はマジでいいよ
値段がアレだけど
Knuth の例の The Art... から
おいしいとこだけもってきて
きれいに系統だてたかんじ
TeX による組版のできもいい
数式は普通はイタリックを使うのが慣例だけど
別の書体を使ってるのがおもしろい
そもそもイタリックを使うのは普通の文章の
ローマンと区別するためのもので
区別させつけば数式がイタリックである必然性はない
0253デフォルトの名無しさん
2005/08/25(木) 16:13:28俺は洋書で買って、すぐ単元ごとに切り分けてしまったよ。
定評あるコンピュータ本って豪華な装丁だったり、
ページ数が膨大だったりするから、ホントはいつでも気軽に参照したいのに
(気持ち的に)アクセスが重いってことがある。
自分の血肉とするためには、本の体裁など気にせず
とことん自分仕様に改良するのがよい。
同じ本の違うページを同時に開いて見れるしね。
0254デフォルトの名無しさん
2005/08/26(金) 21:52:02わかった!
さっそく切り刻んでくる!
0255デフォルトの名無しさん
2005/08/26(金) 23:03:010256デフォルトの名無しさん
2005/08/30(火) 01:26:450257デフォルトの名無しさん
2005/08/30(火) 09:44:390258デフォルトの名無しさん
2005/08/30(火) 13:08:36>>6
イタリア語で85000(8500だったかな・・)は
「 オッタッタチンコミーレ !」
0259デフォルトの名無しさん
2005/08/30(火) 22:30:49http://www.geocities.jp/supermisosan/jacobi.html
b配列に代入されている初期値の意味は何ですか?
0260デフォルトの名無しさん
2005/08/30(火) 22:50:212.0*x + 5.0*y + -8.0*z = 8.0
1.0*x + -4.0*y + 7.0*z = 2.0
を解くプログラムで、a[n][n]が左辺の係数、b[n]が右辺ではないかい。
0261デフォルトの名無しさん
2005/08/31(水) 01:42:08すばらしいぃ!!
疑問が晴れました。ありがとうございます。
0262デフォルトの名無しさん
2005/09/13(火) 01:54:25>>245 の「コンピュータの数学」グレアム著、
の初歩を紹介した感じの内容だね。
なんつーか、この人わかりやすい比喩を考えるのがうまい。
www.amazon.co.jp/o/ASIN/4797329734/
0263デフォルトの名無しさん
2005/10/02(日) 01:33:08非定常の差分方程式を解くプログラムをCで作りたいのですが
おすすめの参考書とかありますか?
Tm+1 + Tm-1 + 2Tm / �凅^2 = 1/α δTp/δt
の方程式です
0264デフォルトの名無しさん
2005/10/05(水) 11:38:180265デフォルトの名無しさん
2005/10/05(水) 11:56:16本当に基礎的な証明でも膨大な時間をかけて講義してたらしい
0266デフォルトの名無しさん
2005/10/05(水) 12:11:37そのレベルなら、俺が師匠の手伝いしておごってもらった時に書いた本でも作ったな
てか、初歩の数値解析関係で、差分方程式のサンプルにまず100%含まれてる
後は、境界条件と制約条件をどうするかだけしかないな
0267デフォルトの名無しさん
2005/10/09(日) 02:17:46例で書いた数式1次元ですね。すいません
んー、3次元の非定常のプログラムがなかなかなくて・・・
参考書はほとんど二次元で終わってまして
0268デフォルトの名無しさん
2005/10/12(水) 17:21:19x^2/a^2 + y^2/b^2 = 1とネットで調べると出てくるのですが、
円の中心(本来はそういわないかもしれない)が(0,0)で
ない場合の式は、この式じゃないような気がするんですが。
0269デフォルトの名無しさん
2005/10/12(水) 17:26:540270デフォルトの名無しさん
2005/10/12(水) 17:34:17さら平行移動させたり、回転させたりすると、当然複雑になるけど、複雑にしてもメリットないでしょ?
実際に楕円描く時はその式では不便で、sin/cosの媒介変数型使うだろうし
・・・・いや、実際に楕円描く場合はベジェ曲線で近似するだろうな
0271デフォルトの名無しさん
2005/10/12(水) 18:28:16応用が効かない言い方をするならば、
中心(c,d)なら(x-c)^2/a^2+(y-d)/b^2=1と考えていい。
0272デフォルトの名無しさん
2005/10/12(水) 19:36:15曲線f(x,y)=0を x軸方向に+a, y軸方向に+b したら f(x-a,y-b)=0
x=f(t), y=g(t) だったら x=f(t)+a, y=g(t)+b だし
0273デフォルトの名無しさん
2005/10/12(水) 21:49:370274デフォルトの名無しさん
2005/10/15(土) 20:36:140275デフォルトの名無しさん
2005/10/16(日) 22:27:040276デフォルトの名無しさん
2005/10/19(水) 23:05:17一応言っておくと「楕円方程式」と「楕円の方程式」は別物ですよ
ここで言ってるのは前者ですね
0277デフォルトの名無しさん
2005/10/19(水) 23:06:430278デフォルトの名無しさん
2005/10/21(金) 20:14:160279デフォルトの名無しさん
2005/10/21(金) 23:53:520280デフォルトの名無しさん
2005/10/22(土) 00:56:25前者に決まっとる。
0281デフォルトの名無しさん
2005/10/22(土) 03:32:24どっちもどっちw
0282デフォルトの名無しさん
2005/10/24(月) 13:41:540283デフォルトの名無しさん
2005/10/24(月) 21:16:45参考:>>268
0284デフォルトの名無しさん
2005/10/26(水) 14:03:530285デフォルトの名無しさん
2005/10/26(水) 23:53:36嘘つきになる要素が・・・
0286デフォルトの名無しさん
2005/10/27(木) 13:49:130287デフォルトの名無しさん
2005/10/27(木) 14:31:050288デフォルトの名無しさん
2005/10/27(木) 15:00:160289デフォルトの名無しさん
2005/10/27(木) 22:38:310290デフォルトの名無しさん
2005/10/28(金) 00:58:510291デフォルトの名無しさん
2005/10/28(金) 01:04:290292291
2005/10/28(金) 01:05:140293デフォルトの名無しさん
2005/10/28(金) 01:10:03素数を生成するコードを書いて落ち着くんだ…
0294デフォルトの名無しさん
2005/10/28(金) 01:21:5420 N = N + 1
30 K=2
40 IF N = K THEN LPRINT N
50 IF N MOD K = 0 THEN GOTO 20
60 K = K + 1
70 GOTO 40
0295デフォルトの名無しさん
2005/10/28(金) 10:58:220296デフォルトの名無しさん
2005/11/04(金) 02:06:550297デフォルトの名無しさん
2005/11/04(金) 22:23:21「だぜ」に突っ込むか・・・。
「やれやれだぜ」自体が漫画の台詞。
ジョジョ。
0298デフォルトの名無しさん
2005/11/04(金) 23:27:220300& ◆BSGaog1Z0Q
2005/11/08(火) 14:18:52//>>294
std::ostream& dumpPrimes(std::ostream& os, int max) {
for(int n=2; n <= max; ++n){
for(int k=2; ; ++k){
if(n==k) {
os << n << "\n";
}
if(n % k == 0) {
break;
}
}
}
return(os);
}
0301デフォルトの名無しさん
2005/11/08(火) 14:36:01if(n==k){
os << n << "\n";
}
return((n % k == 0)? os: testDiv(os, n, k+1));
}
std::ostream& testNum(std::ostream& os, int max, int n){
return((n > max)? os: testNum(testDiv(os, n, 2), max, n+1));
}
std::ostream& dumpPrimes(std::ostream& os, int max) {
return(testNum(os, max, 2));
}
0302デフォルトの名無しさん
2005/11/16(水) 22:10:39> 素数は無限に存在する。
>
> エウクレイデスによる証明
> 背理法による。
> 素数が有限個しかないと仮定し、それらを次のようにおく。
>
> pi, i <= n
> ただし n は定数。
>
> q = p1p2p3...pn + 1
>
> を考えよう。q は合成数であるか素数であるかのいずれかである。
> q が合成数だとすると q は pi のいずれかを用いて積の形に表されるはずである。その一方で q は pi のいずれで割っても 1 があまり、矛盾する。
> 素数だとすると、これは pi のいずれとも異なるから素数が有限個しかないことに反する。
> Q.E.D.
この証明変だよね。
0303デフォルトの名無しさん
2005/11/16(水) 22:18:18どこで引っかかった?何が変?
0304デフォルトの名無しさん
2005/11/16(水) 22:22:32「qが合成数」も「素数が有限」も否定していいんだよ。
0305302
2005/11/16(水) 22:29:52そうとは限らない
2 + 1 =
3 is a prime.
2 x 3 + 1 =
7 is a prime.
2 x 3 x 5 + 1 =
31 is a prime.
2 x 3 x 5 x 7 + 1 =
211 is a prime.
2 x 3 x 5 x 7 x 11 + 1 =
2311 is a prime.
2 x 3 x 5 x 7 x 11 x 13 + 1 =
30031 is 509 x 59
2 x 3 x 5 x 7 x 11 x 13 x 17 + 1 =
510511 is 97 x 5263
2 x 3 x 5 x 7 x 11 x 13 x 17 x 19 + 1 =
9699691 is 27953 x 347
2 x 3 x 5 x 7 x 11 x 13 x 17 x 19 x 23 + 1 =
223092871 is 703763 x 317
0306デフォルトの名無しさん
2005/11/16(水) 22:37:59> 2 + 1 =
>3 is a prime.
2<3
> 2 x 3 + 1 =
>7 is a prime.
3<7
> 2 x 3 x 5 + 1 =
>31 is a prime.
5<31
> 2 x 3 x 5 x 7 + 1 =
>211 is a prime.
7<211
> 2 x 3 x 5 x 7 x 11 + 1 =
>2311 is a prime.
11<2311
> 2 x 3 x 5 x 7 x 11 x 13 + 1 =
>30031 is 509 x 59
13<59
> 2 x 3 x 5 x 7 x 11 x 13 x 17 + 1 =
>510511 is 97 x 5263
17<97
> 2 x 3 x 5 x 7 x 11 x 13 x 17 x 19 + 1 =
>9699691 is 27953 x 347
19<347
> 2 x 3 x 5 x 7 x 11 x 13 x 17 x 19 x 23 + 1 =
>223092871 is 703763 x 317
23<317
論外なところから値を引っ張ってきて論じようとしている。
23が「最大」の素数であれば、317は素数でも合成数でもないのでどちらにしても矛盾がしょうじる。
0307デフォルトの名無しさん
2005/11/16(水) 22:41:27この推論は「素数が有限個しかない」という偽の前提があるから、反例を持ち出してもしょうがない。
> q が合成数だとすると q は pi のいずれかを用いて積の形に表されるはずである。
これは、「合成数は有限個の素数の積で表される」という定理と「素数はpiしかない」からの帰結。
0308302
2005/11/16(水) 22:49:50> q が合成数だとすると q は pi のいずれかを用いて積の形に表されるはずである
そうでない場合もあるとは思いませんか?
そして、その場合についても論ずるべきでは?
0309デフォルトの名無しさん
2005/11/16(水) 22:57:11だから、>>307に書いたように、
> q が合成数だとすると q は pi のいずれかを用いて積の形に表されるはずである
は「素数がpiしかない」という仮定の下で証明できる。
証明できるからには、これが偽である可能性を検討する必要はない。
0310302
2005/11/16(水) 23:37:43>「合成数は有限個の素数の積で表される」という定理
を前提としているわけですか?
0311デフォルトの名無しさん
2005/11/16(水) 23:45:05群論か環論の本を読んでくださいな、素元、既約元のあたりです。あとできれば背理法についてもね。
0312デフォルトの名無しさん
2005/11/16(水) 23:57:48細かい違いはあるかも知れない(例えば、その定理の代わりに「合成数は素因数を持つ」でも十分)けど、
基本的にはその通りだと思う。
0313302
2005/11/17(木) 00:16:30> を考えよう。q は合成数であるか素数であるかのいずれかである。
> q が合成数だとすると q は pi のいずれかを用いて積の形に表されるはずである。その一方で q は pi のいずれで割っても 1 があまり、矛盾する。
> 素数だとすると、これは pi のいずれとも異なるから素数が有限個しかないことに反する。
こう変えてもいいですか?
を考えよう。q は合成数であるか素数であるかのいずれかである。
q が合成数だとすると q は pi のいずれかを用いて積の形に表されるとき、q は pi のいずれで割っても 1 があまり、矛盾する。
また q は pi のいずれかを用いて積の形に表されないとき、q は pi を超える素数を持ち、矛盾する。
素数だとすると、これは pi のいずれとも異なるから素数が有限個しかないことに反する。
0314デフォルトの名無しさん
2005/11/17(木) 00:19:56128bit暗号を解読できるくらいのを
0315デフォルトの名無しさん
2005/11/17(木) 00:31:12ちょっとおかしいところ(「piを超える素数」じゃなくて「pi以外の素因数」というべき)があるけど、大筋としてはそれでもいい。
でも、そう書き換えたところで推論の構造はほとんど変わっていないし、
「合成数は素因数をもつ」ことに依存しているのも変わらない。
0316302
2005/11/17(木) 00:46:27自分なりに納得できて助かりました。
ありがとうございました。
0317デフォルトの名無しさん
2005/11/17(木) 01:16:310318デフォルトの名無しさん
2005/11/17(木) 11:36:57いや、真っ先に背理法を理解すべきでは?
0319デフォルトの名無しさん
2005/11/17(木) 12:27:310320デフォルトの名無しさん
2005/11/17(木) 13:05:040321302
2005/11/17(木) 16:06:10素数を順番にかけるとint型では、23までしかいけなかったので
double型でやったら、途中からおかしな答えになりました。
下から2番目なんですが、答えが1大きくなりました。
なんでなんでしょう?
2 x 3 x 5 x 7 x 11 x 13 x 17 x 19 x 23 x 29 + 1 =
6469693231 is 331 x 19545901
2 x 3 x 5 x 7 x 11 x 13 x 17 x 19 x 23 x 29 x 31 + 1 =
200560490131 is a prime.
2 x 3 x 5 x 7 x 11 x 13 x 17 x 19 x 23 x 29 x 31 x 37 + 1 =
7420738134811 is 181 x 40998553231
2 x 3 x 5 x 7 x 11 x 13 x 17 x 19 x 23 x 29 x 31 x 37 x 41 + 1 =
304250263527211 is 61 x 4987709238151
2 x 3 x 5 x 7 x 11 x 13 x 17 x 19 x 23 x 29 x 31 x 37 x 41 x 43 + 1 =
13082761331670032 is 2 x 6541380665835016
2 x 3 x 5 x 7 x 11 x 13 x 17 x 19 x 23 x 29 x 31 x 37 x 41 x 43 x 47 + 1 =
614889782588491390 is 2 x 307444891294245700
0322デフォルトの名無しさん
2005/11/17(木) 16:08:500323デフォルトの名無しさん
2005/11/17(木) 17:01:31浮動小数点表現の仮数部が 52 ビットのとき、精度は 15.6 桁。
…って、計算機の初歩の初歩だろうに。
0324302
2005/11/17(木) 17:11:45レスありがとうございます。
2進数の表し方がわからなかったので、
16進数で見てますが、所々おかしなところがあり、
手直ししています。
>>323
そうでしたか、出直してきます。
0325デフォルトの名無しさん
2005/11/17(木) 21:07:29巡回セールスマン問題解くプログラム作ってるんですけど
このスレでいいっすか?
それとも数学板に行ったほうがいいですかね?
http://www.watanabe.nuie.nagoya-u.ac.jp/member/jien/lectures/tokuron/toku5.pdf
↑でとりあえず理屈はわかったんだけど
摂動行列の決め方がわからん。
ていうか
あるルート候補から次のステップのルート候補作るときに
全部シャッフルしたらさすがに効率悪そうだし
隣同士入れ替えただけだともろ局所解に引っかかりそうだし。
なんか巡回セールスマン向けのうまい方法
(あるルートXnからX(n+1)を作るアルゴリズム)
だれかしらないですか?
0326デフォルトの名無しさん
2005/11/17(木) 22:09:18アニーリング関数が適当なら、その隣同士とかでやったとして
局所的なくぼ地にはまってしまうような場合であっても
少ない確率ではあるけどそのくぼ地を乗りこえることかできるから
結局それなりにうまくいくんじゃないのかな
0327デフォルトの名無しさん
2005/11/18(金) 00:37:31つ[long long]
0328デフォルトの名無しさん
2005/11/19(土) 09:22:25つ[多倍長演算]
0329デフォルトの名無しさん
2005/11/19(土) 10:35:170330デフォルトの名無しさん
2005/11/19(土) 11:18:12OK!
0331デフォルトの名無しさん
2005/11/19(土) 11:26:320332デフォルトの名無しさん
2005/11/19(土) 11:29:020333デフォルトの名無しさん
2005/11/19(土) 13:23:14サンクス。
変換しても出てこないので違うのかなと思っていた。
0334デフォルトの名無しさん
2005/11/19(土) 13:51:31かな漢字変換辞書が全ての技術語を網羅しているわけがないだろう。
それともあれか、藻前は「へんかんできないからそんな字ない」とか言っちゃう小学生か?
0335デフォルトの名無しさん
2005/11/19(土) 19:32:310336デフォルトの名無しさん
2005/11/19(土) 20:10:17そんなにキツく言わんでも・・。
0337デフォルトの名無しさん
2005/11/20(日) 22:49:35そのタプルが大きさ N の全順列中において辞書順で何番目かを返すアルゴリズムを
実装したいんですが,どういうアルゴリズムが考えられますか?
例えば, {1, 2, 0} を与えたら 4 を返すような感じのアルゴリズムです.
0338デフォルトの名無しさん
2005/11/20(日) 22:53:590339デフォルトの名無しさん
2005/11/20(日) 23:01:25tuple
【名】 組、《コ》タプル
0340デフォルトの名無しさん
2005/11/20(日) 23:08:44それの例のまま、 {1, 2, 0} で考えると、
最初が 1 → このタプルより前には {0, *, *} がある → 前に 2! 個 = 2個 ある
次が 2 → でも1は出てる → これより前で、{0, *, *} 以外のものは {1, 0, *} → 1個
だから、{1, 2, 0} より前には3個のタプルがあるから、
{1, 2, 0} は4番目。って感じ?
1つ数字増やして {2, 3, 1, 0} とかでいくと、
最初、{0, *, *, *} と {1, *, *, *} 数えて、 3! = 6個
次、{2, 0, *, *} と {2, 1, * *} で 2×2! = 4個
次、{2, 3, 0, *} があるから 1個
計11個前にあるから {2, 3, 1, 0} は12番目。
0341デフォルトの名無しさん
2005/11/20(日) 23:10:03で、マッチングをとる。
0342デフォルトの名無しさん
2005/11/20(日) 23:17:15最後の結果も18だわ。
>>341
それだと順列生成自体が糞重たいし、コーディング結構難しくない?
0343デフォルトの名無しさん
2005/11/21(月) 00:32:410344デフォルトの名無しさん
2005/11/21(月) 01:58:190345デフォルトの名無しさん
2005/11/21(月) 11:28:06std::vector<int> src, arr;
src.push_back(1);
src.push_back(2);
src.push_back(0);
arr = src;
size_t ct = 0;
do {
#ifdef _DEBUG
std::cout << "{";
for (int i=0; i<arr.size(); i++)
std::cout << (i? ", ": "") << arr[i];
std::cout << "}" << std::endl;
#endif // _DEBUG
} while (ct++, std::prev_permutation(arr.begin(), arr.end()));
std::cout << ct << std::endl;
結果:
{1, 2, 0}
{1, 0, 2}
{0, 2, 1}
{0, 1, 2}
4
0346デフォルトの名無しさん
2005/11/21(月) 15:45:08素直に説明を読むと、
「1からNまでのN個の順列組み合わせを辞書式順序で並べたとき
各組合わせは何番目?」
という問だと思うのだが、いきなり組の要素に0があったりするのはなぜ?
0347337
2005/11/21(月) 18:58:04素直にやるなら,前から順に読んでいって現れた数字を記録しないとならないですよね.
前から順に読んで現れた数字より大きい数字を全てデクリメントしていくやり方も
あるのかな?参考になります.
>>345
これって 340 さんの考え方じゃなくてむしろ 341 さんの考え方ですよね?
>>346
すいません.最初,説明を「0 から N-1」と書いてて後で直したんですが
例を書き直すのを忘れてました.337 の正しい例は {2, 3, 1} です.
0348デフォルトの名無しさん
2005/11/22(火) 04:02:57それでいいんでね?
(1)result=0で初期化
*ここで、1〜Nの整数が、適当に並べ替えられて並んでいたとして、
(2)N=1ならresult+1を返し終了。
(3)左端の数aに対し (a-1)*(N-1)! を計算し、resultに加算。
(4)列の中で、aより大きい要素は全て1だけ減らす(小さい要素はいじらない)。
(5)先頭のaを切り落とし(N=N-1)、残った列に対し*に戻る。
0349340
2005/11/22(火) 12:02:43それ見てて思ったけど、
(3)(階乗計算)とか(4)の計算に掛かる時間を考えると、
>>345 のコードよりも計算量大きくなるかも。
>>345 だと標準ライブラリの機能そのまま使っててコード量的にも少ないし。
順列生成のアルゴリズムって、>>345 の結果見てる感じ、
バブルソートを1ステップずつやるようなものなのね。
0350デフォルトの名無しさん
2005/11/22(火) 16:31:371 から N までの N 個の異なる自然数による順序付タブル をランダムに出力せよ
0351デフォルトの名無しさん
2005/11/22(火) 18:36:000352337
2005/11/23(水) 22:03:52345 さんのコードだと平均・最悪で O(N!) な一方で
348 さんのだと平均・最悪で(階乗計算も含めて) O(N^2) ですから,
特別な事情がなければ 348 さんのを選択するのが良いと思います.
というか 348 さんのアルゴリズムで実装できました.ありがとうございました.
0353デフォルトの名無しさん
2005/11/28(月) 23:10:17フーリエ変換が趣味です。
http://hobby8.2ch.net/test/read.cgi/hobby/1133005088/
0354デフォルトの名無しさん
2005/11/29(火) 19:35:47べての方程式を解くことは困難なので、そのプログラムで解ける範囲を指定し記述せよ。』
この問題のアルゴリズムと、ソースを書いてくれませんか?
ここで悩んでいるので是非とも力を貸してください!!お願いします。
0355デフォルトの名無しさん
2005/11/29(火) 19:36:390356デフォルトの名無しさん
2005/11/29(火) 19:36:55○答え教えろ
0357デフォルトの名無しさん
2005/11/29(火) 20:08:28とりあえずオイラー法を用意して、それで解ける範囲を明示しておけばレポート通るんじゃね?
0358デフォルトの名無しさん
2005/11/30(水) 01:23:230359デフォルトの名無しさん
2005/12/07(水) 23:46:24上記のサイトで、Sinをテイラー展開で求めているようですが、
sin(0.00) = NaN (非数値)
sin(0.01) = NaN (非数値)
sin(0.02) = NaN (非数値)
・
・
となってしまいます。どこが間違っているのでしょうか?
0360デフォルトの名無しさん
2005/12/08(木) 00:11:24double Sin(double x)
{
double xx = -x * x;
double fact = 1;
double sin = x;
for(int i=1; i<100; )
{
fact *= ++i; fact *= ++i;
x *= xx;
sin += x / fact;
}
return sin;
}
0361デフォルトの名無しさん
2005/12/08(木) 00:16:080362359
2005/12/08(木) 00:39:53なるほど! 数学関連の話で、何が合っていてどこから手をつけていいものか
困り果てているところでした。非常にお早いレス感服いたします。
ありがとうございました。
>>361
すみません。いくつか間違いがみられ、一応メールアドレスを探してみた
のですがわかりやすい位置になかったので、すぐに断念してしまいました。
0363デフォルトの名無しさん
2005/12/08(木) 00:42:210364デフォルトの名無しさん
2005/12/08(木) 01:32:540365340
2005/12/08(木) 11:29:34トップページたどれって。
全ページに連絡先入れといた方がいいのか?
>>363
それいいだすとさ、
π/2 の範囲に収まるように mod 2π とか三角関数の公式使って変換して、
sin と cos のどちらか選んで計算するとか、凝りだすときりない。
餅は餅屋と言う奴で、そういうのは数値計算屋さんにお任せしないと。
0366デフォルトの名無しさん
2005/12/08(木) 22:43:03ウェブページの作り方的な話なら、
全ページにaddress要素は入れるのが基本
address要素内には連絡先として通常メールアドレスが入る
0367デフォルトの名無しさん
2005/12/08(木) 23:13:51それ、確かに推奨はされてるけど、
ボットにアドレス収集されるのが嫌で、実際やってる人いなくない?
0368デフォルトの名無しさん
2005/12/08(木) 23:25:27俺はやってるよ
あるいは自働収集が嫌って言って
画像を使ったり余計な文字をはさんだりしてる人もいる
要は形式上の問題じゃなくて、
>>362みたいな状況にならないようにっていうのが
目的なんだし。
0369 【凶】 【308円】
2006/01/01(日) 07:54:110371デフォルトの名無しさん
2006/01/09(月) 00:38:04マージソートについて勉強してます
マージソートの配列の分割についてなんですけど
例えば8だけの場合は
4 4
2 2 2 2
1 1 1 1 1 1 1 1
とだけ分割できると思います。
ですが例えば要素が34の場合
17 17
と一回だけ分割すると素数になってしまいます。
更に極端にしますと
7や23などの要素の個数の場合はどう分割するのでしょうか・・
まだ勉強初めて1日だけなんでアレですが参考のページでもあると嬉しいです
0372デフォルトの名無しさん
2006/01/09(月) 00:49:45そりゃ、8と9に割ればいいじゃん。
0373デフォルトの名無しさん
2006/01/09(月) 00:54:24しかしその場合併合(マージ)しようとするとどんな風になるのでしょうか
グーグル先生のイメージ検索はエレメントの個数が偶数ばかりのもので・・
勉強してきます(´・ω・`)
0374デフォルトの名無しさん
2006/01/09(月) 00:56:370375デフォルトの名無しさん
2006/01/09(月) 03:02:59べつに偶数のときとかわらん。
違う数に分割されたからといって
コードの何処も変える必要なんかない。
そんな必要があったら
そもそもマージソートじゃない
0376デフォルトの名無しさん
2006/01/09(月) 10:30:58例えば9の場合、
4, 5
(2, 2), (2, 3)
((1, 1), (1, 1)), ((1, 1), (1, 2))
((1, 1), (1, 1)), ((1, 1), (1, (1, 1)))
と分割すればいい。
てか図解すると並列に見えるけど、実際の処理は再帰的に行われるから、
部分的に階層が深くなっても問題ない。
ttp://oku.edu.mie-u.ac.jp/~okumura/algo/archive/algo.lzh
の mergsort.c 読めば解るよ。
0377デフォルトの名無しさん
2006/01/09(月) 13:56:17つーか、聞く前に一度自分で実際にやってみ。
マージって全然個数の違う2つの列相手にもできるから。
例えば極端な話 1, 16 とかで分けてもマージできる。
0378マイク ◆yrBrqfF1Ew
2006/01/23(月) 00:20:53堀江由衣の処女
0379デフォルトの名無しさん
2006/01/23(月) 04:26:20/ 、゙ヽ、 ‐-'´ ヽ‐- / / ヽ
,/´ .., ヽ,,l_),' '、,ト/ / ヽ
/ ヽ,r' l .,、z:ュ、,_. ,、=,゙.-〈__r,'、 ヽ_
_.l ヽ」 , l. ´ ,r'ャ、`' i'rャ;| ゙‐ヽ、_,, /l
,l l| −'´ll ,ll ヽ. ''`¨¨´ ヽ | .,//゙l //\
l`l| l|ヽ ヽ 入 ,ィ _. ', l |l // }
l \ l| ,l ヽ_. ' `'゙`'‐'i゙ ,' |,l // l バカジャネーノ
/ '\ l|`l l ヽ`'. ,∠.ニフ / l ヽ // |
,l '\ l| .lヽ_l ` 、 、 い.... ,' /___/ | ∨/ ,}
| ヽl | ,| .ヽ', .ヽ`二フ.,' ヽ ,| ,l
| l ,l ' .、,_`,.ィ l \ /
ヽ | \. ヽ/ l ヽ /j
\ / ヽ ヽ | l /
゙l\.. / ヽ ヽj | , /
ヾ ヽ ヽ ヽ / ,l
ヽ、 ヽ l } / ,r
ヽ ヽ | /′ ,,...''
`'':..、 ___ ___,..-.. |, ,l , :..-‐'"´
 ̄ /lr‐‐‐'--、_..... l_,..-'''""'- "
0380デフォルトの名無しさん
2006/01/30(月) 23:54:04中心点から角度θの軸を考え、横2W縦2Hの長方形の中心がこの軸を移動するとして
長方形の頂点が円周上にくる点の座標は?
0382デフォルトの名無しさん
2006/01/31(火) 00:41:23人にモノ聞く態度じゃねぇな。
出直して来い。
0383デフォルトの名無しさん
2006/01/31(火) 09:44:41もしかして日本語が不自由な方なんでわ?
0384デフォルトの名無しさん
2006/01/31(火) 10:41:190385デフォルトの名無しさん
2006/01/31(火) 18:07:460386デフォルトの名無しさん
2006/01/31(火) 18:42:49r=-(Wcosθ+Hsinθ)+√(R^2-(Wsinθ-Hcosθ)^2) として
(rcosθ+W, rsinθ+H) …(1)
θの範囲は0≦θ<π/2だけどね。他の範囲は(1)式にある符号が変化するだけのはず。
0387デフォルトの名無しさん
2006/02/01(水) 08:22:33それだと、 W = H = 0 の時にしか、長方形の頂点が円周上にこないんじゃない?
「点の座標」ってのが長方形の中心だとしても頂点だとしても。
多分やつは問題を理解できてないんだろうなあ。
2Dゲーム作ってて、円形フィールドからはみ出さないようにしたい、と推測したがどうか。
0388デフォルトの名無しさん
2006/02/01(水) 17:20:04それだと浮動小数点でやってる場合、
点の座標を求めて=で判別したらイタイ目見そうだ。
(丁度等しくならなかった時に判定をスリ抜ける。)
不等式で「内側」の度合い評価できる基準を見つけないと。
0389デフォルトの名無しさん
2006/02/01(水) 17:29:09たぶん長方形の外接円で判定してごまかすのが楽なんじゃないかなぁ。
中心点同士の距離と半径だけで判定できる。
2乗ノルムでやれば平方根もいらないし。
もし精密にやりたければ上記の方法で簡易判定して明らかに接触しない事例を除いた上で
拙速するかもしれないケースだけ長方形を構成する線分が交差するかどうかを判定すれば
いいのかな。
0390デフォルトの名無しさん
2006/02/01(水) 17:31:22軌道は具体的に長方形の位置を計算する以外は無関係だね。
0391マイク ◆yrBrqfF1Ew
2006/02/04(土) 00:33:27を満たすxの値の範囲は
a < x <= b
このaとbと導かれる過程がわからん。
包茎じゃなくて東京六大学に在学中か卒業した人だけ答えてくださーい。
0392デフォルトの名無しさん
2006/02/04(土) 11:10:140393デフォルトの名無しさん
2006/02/05(日) 13:36:53| ̄| /> _ _ / _| | | />
| |// \ \ __ / / \ \. | |//
| \ \ \/ \/ / 〉 ヽ. | \
| |\ \ \ /\ / | ̄ ̄ ノ | |\ \
/::::/'"  ̄ヾi /  ̄ ̄ \ /::::/'" ̄ ̄ヾi /  ̄ ヽ
|:::::::| ,,,,,_ ,,,,,,| | ^ ^ | |:::::::| ,,,,,_ ,,,,,,| | ^ ^ |
|r-==( 。);( 。) | >ノ(、_, )ヽ、.| |r-==( 。);( 。) | >ノ(、_, )ヽ、.|
( ヽ :::__)..:: } ! ! -=ニ=- ノ ! ( ヽ :::__)..:: } ! ! -=ニ=- ノ !
ヽ ー== ; \ `ニニ´ / ヽ ー== ; \ `ニニ´ /
\___ !  ̄ ̄ \___ !  ̄ ̄ ̄
0394デフォルトの名無しさん
2006/02/06(月) 04:30:36あんなかでまともな偏差値誇ってんのは
東京 慶應 早稲田 だけだろ
あとの三つは……やめた。
0395デフォルトの名無しさん
2006/02/06(月) 11:34:15名大出のしがない漏れは
頑張って答えようとしたりせずにニヤニヤしてれば好いわけね?
0396デフォルトの名無しさん
2006/02/06(月) 19:59:060397デフォルトの名無しさん
2006/02/07(火) 22:59:21しかし>391ってまるっきり高校レベルじゃないか
0398デフォルトの名無しさん
2006/02/08(水) 12:53:00シッ! 皆わかっててニヤ(・∀・)ニヤしてるんだから!
0399マイク ◆yrBrqfF1Ew
2006/02/08(水) 18:49:120400デフォルトの名無しさん
2006/02/09(木) 10:40:54イヤー、オヂサン参ったな。ナント院卒だが答えられないぞ。
>391が出した問題の最後に示された条件:
> 包茎じゃなくて東京六大学に在学中か卒業した人
を満たすことができないからなー。何せ漏れは重度の冠頭包茎だし、
名古屋生まれの名古屋育ちで大学も院も名古屋大学だからな。
漏れがこの問題に答えるためには
東京6大学(慶応、上智、東京、明治、立教、早稲田)のどっかに入らないと
イカンわけだろ?
いやー回答者の資格に制限がある問題は難しいわ、ホント難しい、
全然答えられそうにないよwww
最近のセンター試験ってこんな答えるのが難しい問題出すの?(・∀・)
0401デフォルトの名無しさん
2006/02/09(木) 10:43:17相変わらず莫迦なコテだな。
0402400
2006/02/09(木) 11:08:11My包皮も切らないとイカンのかー。やっぱ難しいわ。
漏れ包皮ついてるほうが好きだし、
臆病者だから体の一部切るなんてコワいしなー。
0404400
2006/02/09(木) 17:45:26思うんだが学歴を今後も問題の条件に入れ続ける気なら
学歴板でも行ったらどうかと思う。
ここはプログラム板の数学スレだから学歴は関係ない。
そしてプログラムに関係ないセンター試験の問題について
受験生向け解説依頼を受け付けるスレでもない。
受験生ならオトナシク受験板でも行ってたほうがいいんじゃない?
大学受験板
http://etc4.2ch.net/kouri/
学歴板
http://tmp6.2ch.net/joke/
0405マイク ◆yrBrqfF1Ew
2006/02/09(木) 18:10:56答えなくなければね;)
0406デフォルトの名無しさん
2006/02/09(木) 18:16:130407デフォルトの名無しさん
2006/02/09(木) 18:39:45今後スルーが相応と思うがどうか?
0408デフォルトの名無しさん
2006/02/09(木) 19:39:20ー 再開 ー
0409マイク ◆yrBrqfF1Ew
2006/02/09(木) 21:36:08毎度スルースルーと言ってる割には毎度我慢できなくなってレスしちゃうんだね。
子供?;)
0410デフォルトの名無しさん
2006/02/10(金) 00:01:50近似的に解くべきだなw
0411デフォルトの名無しさん
2006/02/10(金) 02:34:30えーと、どこで笑えばいいのかな?
0412デフォルトの名無しさん
2006/02/10(金) 16:15:10じゃぁ、両端を探すために二分探索でもするか?
0413デフォルトの名無しさん
2006/03/01(水) 22:02:54m9(^Д^)
じゃなくて
;)
を無理して使うとこ
0414デフォルトの名無しさん
2006/03/03(金) 14:58:110415デフォルトの名無しさん
2006/03/09(木) 09:58:52とかがかわいくて好きでつ
0416デフォルトの名無しさん
2006/03/09(木) 11:39:330417デフォルトの名無しさん
2006/03/09(木) 11:56:571+1=
0418デフォルトの名無しさん
2006/03/09(木) 12:14:220419デフォルトの名無しさん
2006/03/09(木) 12:24:330420デフォルトの名無しさん
2006/03/09(木) 12:34:320421デフォルトの名無しさん
2006/03/09(木) 13:24:560422デフォルトの名無しさん
2006/03/09(木) 23:36:010423デフォルトの名無しさん
2006/03/10(金) 00:06:04a÷b の結果を四捨五入して ret に取得する処理を以下のようにしました。
[バージョンA]
// 除数と被除数の符号チェック
if ((a ^ b) < 0) {
// a, bが異符号
// ret = (a / b) - (1 / 2)
ret = (2 * a - b) / (2 * b);
}
else {
// a, bが同符号
// ret = (a / b) + (1 / 2)
ret = (2 * a + b) / (2 * b)
}
この方法だと正の場合0.5→1、負の場合-0.5→-1となり
数値0に対して正負の結果が対称になります。
(続く・・・)
0424423
2006/03/10(金) 00:07:03除算結果の整数部: n、小数部: s (s >= 0) としたとき
四捨五入後の結果
・s < 0.5 のとき n
・s >= 0.5 のとき n + 1
[実例]
・-0.6 = -1 + 0.4 = -1 + 0 = -1
・-0.5 = -1 + 0.5 = -1 + 1 = 0
・-0.4 = -1 + 0.6 = -1 + 1 = 0
・0.4 = 0 + 0.4 = 0 + 0 = 0
・0.5 = 0 + 0.5 = 0 + 1 = 1
としたいのですが、上手い方法が見つかりません。
一応、自分なりに考えて以下のように実装したら上手くいきました。
[バージョンB]
if ((a ^ b) < 0) {
// 正数にして計算を行う
a = abs(a); b = abs(b);
// ・整数除算の結果を -1 したもの
// ・小数部を割合化したもの(?)である (b - (a % b)) / b + (1 / 2) を四捨五入したもの
// を加えて求める。
ret = -(a + b) / b + (2 * (b - (a % b)) + b) / (2 * b);
}
else {
ret = (2 * a + b) / (2 * b);
}
除数、被除数の符号チェックをしたりしてスマートではないので
もっとシンプルにできる整数演算での上手い方法はあるのでしょうか?
よろしくお願いします。
0425デフォルトの名無しさん
2006/03/10(金) 00:11:330426423
2006/03/10(金) 01:16:10画像処理で使用するため浮動小数点は、できるだけ使用しないようにしています。
ちなみにバージョンA,Bともに四捨五入をするときは
ret = (a / b) + 0.5
= (a / b) + (1 / 2)
= (2 * a + b) / (2 * b) ←通分(だったっけ?)
のように0.5を足すようにしています。
0427423
2006/03/10(金) 01:32:01単純に0.5を足して切り捨てると除算結果が負数の場合に問題があるのです。
(たとえば、結果が-2のときは -2.0+0.5 → -1.5 → -1になってしまう)
そのためにバージョンAでは、除算の結果が正負で場合わけをして
+0.5か-0.5を切り分けることにしました。
0428デフォルトの名無しさん
2006/03/10(金) 02:46:54そもそもここに間違いがあると思うだけどな
0429デフォルトの名無しさん
2006/03/10(金) 11:44:35方法はそれしかない。可搬性を確保したい場合符号チェックは必然。
挙げられた例題は結果の精度として整数値しか必要でない
(小数点以下0ビットの精度)場合の固定小数点演算と看做すことができる。
固定小数点演算とは例えば小数点以下に2ビットの精度が必要な場合に
3ビット下駄を履かせて1→8、0.5→4とするなどして整数演算によって
一定精度の実数演算を行う方法だ。
その場合結局四捨五入の処理も必要になる。
最下位ビットが0か1かを決めるために剰余を使うのもまさに例題と同じだ。
(固定小数点演算という枠組みで考える理由は精度が異なる場合も
同じ考え方で統一的に考えられるというだけだ。)
そして符号の処理も結局必要になる。
ただ、符号付整数除算のハードウェア仕様としてはAもBもありえて、
ハードウェアの仕様を調べてそれに依存するなら処理を省略できる可能性はある。
通常の整数は2の補数表示をすることで正数に対する処理を転用して
負数を扱っているで0に対して表現が元々対称でない。
だからBバージョンが目当てなら見込みは割とある。
ハードウェアに依存しちゃうけどね。
0430デフォルトの名無しさん
2006/03/10(金) 11:58:452ビットの精度が必要な場合→3ビットの精度が必要な場合
精度は悪化するけど2bitの精度が必要な場合に3bit取って
剰余は見ないで最下位ビットだけ見て四捨五入って手はあるけどね。
0431423
2006/03/10(金) 14:08:49やはり、符号チェックは必要なのですね・・・
浮動小数を使いたくないというのは、参考にしているライブラリの処理速度を計測したところ
その結果から浮動少数は使っていないと思われるためです。
ちなみに、そのライブラリは Win32API のウィンドウとビューポート間の座標変換処理で
比較対照としているものは LPtoDP() という関数です。
こいつが結構くせもので、整数部を n 、小数部を s としたとき
だいたい s >= 0.47 で四捨五入して n + s → n + 1 としているのです。
(負数の場合も -1.53 = -2.0 + 0.47 = -1 [入]、-1.54 = -2.0 + 0.46 = -2 [捨])
上記のように、四捨五入の仕様は>>424のバージョンBと同じです。
整数演算で 0.47 くらいで四捨五入なんて特殊なことをしているので
何か整数演算独自の四捨五入の方式があるのかと思い質問させていただきました。
個人的にはバージョンBの除数、被除数の符号が異なる場合に
除算を2回行うというのに満足できないので、もう少し紙とペンで色々と考えてみます。
0432デフォルトの名無しさん
2006/03/10(金) 15:50:22その方が解析しやすい。
0433デフォルトの名無しさん
2006/03/10(金) 15:55:04小数点以下の数を10進表記するとそれだけで誤差が含まれる。
だからこの場合、実数に換算しないで考えたほうがよい。
つまり計算させてる事例に関わる整数の比がないと
何が行われているか正確なことはわからない。
LPtoDP()ってことは窓の大きさとディスプレイの解像度が絡むんだろ?
0434デフォルトの名無しさん
2006/03/10(金) 16:18:42ret = (2 * a + b*sgn) / (2 * b);
として sgn を 1か-1にすればいい
あるいは
ret = (2 * a + b*(sgn+1)-b) / (2 * b);
とすれば sgn+1 は 0か2なので 0か-1の変数fを使い
ret = (2 * a-b + 2*b&f) / (2 * b);
xor結果の最上位で fを-1か0にすればいい
・右へのビット幅だけシフト
・インラインアセンブラを使って符号拡張命令
して、符号ビットを埋めて ゴチャゴチャやる方法があるけど、そんなの使いたい?
0435デフォルトの名無しさん
2006/03/10(金) 16:27:53ret = (2 * a + (b^t)) / (2 * b);
として tを 0か-1とする方法もある
でもたぶん
ret = ( a + ((b^t)>>1) ) / b;
あたりでやってんじゃないかな
0436423
2006/03/10(金) 18:01:03>>432 >>433
座標変換のための設定は以下の通りです。
// マッピングモード設定
::SetMapMode(hDc, MM_ANISOTROPIC);
// ウィンドウ領域 (0, 0) - (1000, 1000)
::SetWindowExtEx(hDc, 1000, 1000, NULL);
::SetWindowOrgEx(hDc, 0, 0, NULL);
// ビューポート設定 (0, 0) - (10, 10)
::SetViewportExtEx(hDc, 10, 10, NULL);
::SetViewportOrgEx(hDc, 0, 0, NULL);
単純にウィンドウからビューポートへ(1/100)倍する変換です。
ウィンドウ、ビューポートのx座標をそれぞれ wx, vx として
-1000 <= wx <= 1000 の範囲で変換してます。
変換式は憶測ですが
vx = wx * (10 / 1000)
で求めていると思われます。
四捨五入の「入」、「捨」の境界は以下のとおり 0.46〜0.47 です。
(これ以外の146, 246, ...、-154, -254, ... でも同様の結果です)
wx = -54 → vx = -1 (-54/100 → -0.54 → -1 + 0.46 → -1)
wx = -53 → vx = 0 (-53/100 → -0.53 → -1 + 0.47 → 0)
wx = 46 → vx = 0 (46/100 → 0.46 → 0 + 0.46 → 0)
wx = 47 → vx = 1 (47/100 → 0.47 → 0 + 0.47 → 1)
>>434
その方法だと>>423のバージョンAと同じになってしまうんです。
今は>>424のバージョンB方式の四捨五入の実装方法で迷ってるんです。
0437423
2006/03/10(金) 18:03:48MSDNで調べたところウィンドウ領域は32bitを保証してビューポートは27bitしか保証しないと
明記されてます。
残り5bitを小数部とした固定小数点で計算とかをしているんですかね?
0438434
2006/03/10(金) 18:21:37折角書いてやったんだから、ちゃんと読め!
いいか その>>424のバージョンB ってのは a/bの符号によって
符号負 ret = (2 * a - b) / (2 * b);
符号正 ret = (2 * a + b) / (2 * b)
としたいわけだろ?
符号を sgn +1/-1 なら
ret = (2 * a - b*sgn) / (2 * b);
だろが!
0439デフォルトの名無しさん
2006/03/10(金) 18:57:40論外
ポイントになるのは、0.1刻みとして(-2.5〜-1.6), (-1.5〜-0.6), (-0.5〜0.4), (0.5〜1.4)
をどうやって同じグループにするかということ
0440デフォルトの名無しさん
2006/03/10(金) 19:21:540441434
2006/03/10(金) 19:27:10単に >>435の符号を入替えて
int div(int x,int y){
int sgn=x^y;
sgn=sgn>>31;
return (x-((-y^sgn)>>1))/y;
}
x/y = div(x,y)
-20/10= -2
-19/10= -2
-18/10= -2
-17/10= -2
-16/10= -2
-15/10= -1 -5/10= 0 5/10= 1 15/10= 2
-14/10= -1 -4/10= 0 6/10= 1 16/10= 2
-13/10= -1 -3/10= 0 7/10= 1 17/10= 2
-12/10= -1 -2/10= 0 8/10= 1 18/10= 2
-11/10= -1 -1/10= 0 9/10= 1 19/10= 2
-10/10= -1 0/10= 0 10/10= 1 20/10= 2
-9/10= -1 1/10= 0 11/10= 1
-8/10= -1 2/10= 0 12/10= 1
-7/10= -1 3/10= 0 13/10= 1
-6/10= -1 4/10= 0 14/10= 1
これでいいんだろ?
0442434
2006/03/10(金) 19:32:15ホントは四捨五入でret = (2 * a + (b^t)) / (2 * b) を使いたかったけど2つある2倍が嫌なんで
ret = ( a + ((b^t)>>1) ) / b; としたら、プラス側が6で変化したんで
ret = ( a - ((-b^t)>>1) ) / b; として、まあマイナス側に-6で変化したっていいやで 計算量優先にしたんだろ
0443434
2006/03/10(金) 19:59:13abs(x) の代わりに (x>>31)^x のようなのを使ってて出た誤差かな
0444デフォルトの名無しさん
2006/03/10(金) 20:05:540445デフォルトの名無しさん
2006/03/10(金) 20:28:42またおまえか。
0446423
2006/03/10(金) 21:07:51とても参考になりました。
特に>>434さんの方法には脱帽しました。
異符号の場合に -1 と XOR して1の補数を用いるなんて思いもつきませんでした。
と言っても、まだ完全には理解できてはいないのですが
先にお礼を言っておきたかったので。
本当にありがとうございました。
0447434
2006/03/10(金) 21:28:04#include <stdlib.h> に div という関数がある 除算とあまりを出す関数だ
div_t d=div(x+y/2,y);
if(d.rem<0) d.quot--; で d.quot を使えばいい
条件判断を無くしたいなら d.quot+d.rem>>31 でいい
たぶんコレが正解だろう
0448434
2006/03/10(金) 21:30:52普通の四捨五入をやりたいって事にやっと気付いた。 すまんな。
0449434
2006/03/10(金) 21:40:25#include <stdlib.h>
int divd(int x,int y){
div_t d=div(x+y/2,y);
return d.quot+(d.rem>>31);
}
int divd(int x,int y){
div_t d=div(x*2+y,y*2);
return d.quot+(d.rem>>31);
}
結果はどっちも >>441 と y=10では同じになる
0450434
2006/03/10(金) 22:11:31>単純に0.5を足して切り捨てると除算結果が負数の場合に問題があるのです
に騙されてしまった。
単純に0.5を足して切り捨てるのをやりたかったのだろう。
ただ、X86では除算の結果が負数になる場合は余りも負数になる。
a/b= n余りsなら
a = n*b + s = s+b+(n-1)*b となる修正をすればいい
アセンブラで書けば、
cdq
idiv
sqr edx,#31
add eax,edx
と4命令
0451デフォルトの名無しさん
2006/03/10(金) 22:13:38かけ算はいってるけど
int func(int a, int b){
int absa = (a >> 31) ^ a;
return (a + absa*b + (b>>1)) / b - absa;
}
0453434
2006/03/10(金) 22:21:39アセンブラだと4行なのに
使わないと除算とmodを別に計算するか div 関数を使う必要があるのが面倒な所
div関数だと結果も構造体渡しだからメモリアクセスが入って遅くなる
0454デフォルトの名無しさん
2006/03/10(金) 22:27:09あんなの使う気にならないんだけど。
0455434
2006/03/10(金) 22:34:07だったらインラインアセンブラでやるといいよ。
cだけで書くなら
x+=y/2;
int r=x/y;
if( (x % y)<0) r--;
return r;
0456434
2006/03/10(金) 22:54:48原理としては、 余りが負数にならないように巧くオフセットを加えてるわけだよね
巧い方法だけど、 a b が大きい時にオーバフローの問題が起きるね。
abs*a ではなくて
aよりも少しだけ大きい bの倍数 を計算させた方がいいのでは?
この場合 >>436 のように座標計算に使うのだから、 マッピングモード設定 時に予め計算させておけばいい
0457434
2006/03/10(金) 22:56:560458デフォルトの名無しさん
2006/03/10(金) 23:12:33廊下にたっていて、向かい側の壁にはたくさんの開くドア又は開かないドアA,B,C。。。。が無限にあって、
それを開くかどうか確認したい。
スタートはAとーAの間にいる。
。。。|D|C|B|A|−A|−B|ーC|−D|。。。っとドアが続く
最初に地点から一番近い、開くドアを見つけたいが、動く距離をxとして、
距離の総和がO(x)ペースになるように探したい。
例えば、A,−A,B,−B,C,ーCの順番で探していくと、
動く距離が、1、2、3、4、5,...nとなり、距離の総和は1/2*(n)*(n-1)となり、
O(X^2)のペースになるから駄目である。
っていう問題なのですが、何か良い探し法、アルゴリズムありますかね?
0459434
2006/03/10(金) 23:13:40( abs(a/b)+1)*b でどうだろ?
除算が遅いなら | b|をシフトしていって |a| を超えた所でもいいか
0460434
2006/03/10(金) 23:21:43で、開いてるか開いていないかの確率はどの程度なの?
というか確率を仮定して
右方に N1内で調べてみてなければ右側にN2個調べて 見つかる確率を求めてみたら?
右側で M番目に開けば左側でM番目まで調べ調べればいいでしょ
0461デフォルトの名無しさん
2006/03/10(金) 23:29:11それも考えたのですが、例えば、3つずつの固まりで調べていくとして、
C,B,A,-C,-B,-A,
F,E,D,−F,−E,−D
距離を考えていくと、
(3+1+1)+(3+1+1)+
(6+1+1)+(9+1+1)+、、、
となって、総和はどうなるのでしょう。。。
0462デフォルトの名無しさん
2006/03/10(金) 23:40:17初期位置より左側を A , C , E と2個おきに移動 して末端 n で(奇数個偶数個で微調整か?)
.... F D B と戻ってくれば O(x) っぽくならない?
0463434
2006/03/10(金) 23:47:520464デフォルトの名無しさん
2006/03/10(金) 23:51:11左の方を一つ探す「A」
右の方を二つ探す「-A, -B」
左の方を四つ探す「B, C, D, E」
右の方を八つ探す「-C, -D, -E, -F, -G, -H, -I, -J」
この要領でいけないかな。
0465デフォルトの名無しさん
2006/03/10(金) 23:55:52>A,−A,B,−B,C,ーCの順番で探していくと、
>動く距離が、1、2、3、4、5,...nとなり、距離の総和は1/2*(n)*(n-1)となり、
>O(X^2)のペースになるから駄目である。
「調べないけど移動してる」に オーダーのコストかかってる?
0466デフォルトの名無しさん
2006/03/11(土) 00:06:37そうすると、kブロックに区切って、
(1)+(1+1)+(3+1+1+1)+(7+1+1+1+1+1+1+1)+(15+1+1+1+1+1+...)+
= 1+1*2+3*2+7*2+15*2+....+?
0467デフォルトの名無しさん
2006/03/11(土) 00:08:01動く距離=そこで調べるの意味だと思います。
0468デフォルトの名無しさん
2006/03/11(土) 00:31:55f(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)
0469デフォルトの名無しさん
2006/03/11(土) 01:10:36なるほど。
助かりました。
ありがとうございました。
0470デフォルトの名無しさん
2006/03/14(火) 22:27:00いまいち何やってるか分からないんだよね。
分かった気にはなるけど、どうもしっくりこないっつうかなんつうか。
他に証明方法とかあるのかね?
0471デフォルトの名無しさん
2006/03/14(火) 22:38:27有理数列と実数の部分集合に1対1の対応が作れる
有理数列は自然数→有理数の写像とみなせる
自然数→有理数の写像の濃度は アレフ0^アレフ0 = 2^アレフ0
アレフ0 < 2^アレフ0
(ある集合の濃度が N のとき、その冪集合の濃度は 2^N、
ある集合とその冪集合の間にはどうやっても全単写が作れない)
なので、可算濃度<実数の濃度
0472デフォルトの名無しさん
2006/03/16(木) 13:32:510473デフォルトの名無しさん
2006/03/18(土) 08:50:19QRを解析するプログラムを作ろうと思っているのですが、
誤り訂正複合のリード・ソロモン符号の複合の仕方が分かりません。
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になりません。
もう数日煮詰まっています。お願いします。何方かご教授して下さい。
0474デフォルトの名無しさん
2006/03/18(土) 09:42:13誰かのコードを参考にするといいよ
例:
http://isw3.kankyo-u.ac.jp/project/2004/project/1013063.pdf
0475デフォルトの名無しさん
2006/03/18(土) 11:04:44ありがたく拝見させてもらいましたが、残念ながら書き込みの処理しか載っていませんでした。
(書き忘れて申し訳ございませんが、473で書いた部分はQRを読み取る部分の処理です)
他にもあるかもと思って「誤り訂正 QR」などで検索を掛けてみましたが、書き込みに関する部分ばかりで、
読み取り解析に関するサイトは見つかりませんでした。
JISの企画書もシンドロームを求めるとしか書いてありませんし、もう、完全に行き詰っています。
どんな些細な事でも結構ですので、アドバイスを頂けないでしょうか。
お願いします。
0476デフォルトの名無しさん
2006/03/18(土) 11:27:53その部分のコード晒してみて、
+ は XORの事とか 掛算は足算の事とか、そこらへんは判ってるんだよね?
0477デフォルトの名無しさん
2006/03/18(土) 11:37:40じゃあココのは見たの? オレは中身はみてないけど
http://sourceforge.jp/projects/qrcode/
0478デフォルトの名無しさん
2006/03/18(土) 12:37:08+は足し算で、掛け算は掛け算で計算していたので、直してみましたが、やはり0になってくれません。
コードはこのようになっています。
//求められたデータ語
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さん、
このプログラムは知りませんでした。早速今から落として見て見たいと思います。
ありがとうございます。
0479デフォルトの名無しさん
2006/03/18(土) 13:24:070480http://www.vector.co.jp/soft/win95/util/se072729.html
2006/03/18(土) 19:44:49もしくは64bitにネイティブ対応したテキスト置換ソフトありますか?
0481デフォルトの名無しさん
2006/03/19(日) 23:37:58XだけYだけ90°ひっくり返すにはどうすればいいの?
0482デフォルトの名無しさん
2006/03/19(日) 23:56:280483デフォルトの名無しさん
2006/03/20(月) 00:24:240484デフォルトの名無しさん
2006/03/20(月) 00:45:38vx = vx + 90/180*PI
0485デフォルトの名無しさん
2006/03/20(月) 11:03:210486デフォルトの名無しさん
2006/03/20(月) 16:05:350487デフォルトの名無しさん
2006/03/20(月) 16:42:100488デフォルトの名無しさん
2006/03/20(月) 17:22:36必ず90度ってどうやるんだ?
0489デフォルトの名無しさん
2006/03/20(月) 17:24:03つーか、玉の軌道をX軸かY軸で反転すれば問題なし。
0490デフォルトの名無しさん
2006/03/20(月) 18:44:580491デフォルトの名無しさん
2006/03/20(月) 18:52:43の149からの流れがここに?
http://pc8.2ch.net/test/read.cgi/prog/1134791216/149-
0492デフォルトの名無しさん
2006/03/20(月) 19:24:06斜面への衝突判定と反射1
0493デフォルトの名無しさん
2006/03/23(木) 07:46:12Hは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関数ではないことを示せ。
なのですが、どなたか分かる方、助けてください。
0494デフォルトの名無しさん
2006/03/23(木) 12:28:39剰余が絡んでる段階で多対一関数だから1方向だろうし
collision-freeではなさそうだわな。
証明は2つの値が実際に1個の値になる例を計算すればいいんじゃね?
0495494
2006/03/24(金) 16:59:01ほ、放置プレイ?
0496デフォルトの名無しさん
2006/03/24(金) 18:30:110497デフォルトの名無しさん
2006/03/28(火) 08:18:12Q1 って関数に名前を付けたいけど適切な名前は?
Q2 浮動小数点なので計算方法は単純に
abs(x)<0.0001 の時は1-x*x/6 でなければ そのまま計算でいいよね?
0498デフォルトの名無しさん
2006/03/28(火) 08:53:02左上と真ん中
どちらのポイントを保持して使うべきですか
0499デフォルトの名無しさん
2006/03/28(火) 08:55:22sine_x_per_x()
0500デフォルトの名無しさん
2006/03/28(火) 09:01:48いまどきどこでもいい。
処理によって必要な座標は違う。例えばサイドビューのジャンプアクションなら足元座標も使う。
必要な座標を計算するメソッドがあればそれでよい。
数学関係ない。
0501デフォルトの名無しさん
2006/03/28(火) 12:19:04それは一般的にsinc関数と呼ばれている。
0503デフォルトの名無しさん
2006/03/28(火) 17:19:04RPGの座標には、マップ上の座標と表示する座標があるが、
マップ上の座標で保存する。
表示する時には、表示する座標に変換して表示する。
0504デフォルトの名無しさん
2006/03/28(火) 17:58:28キャラといえばRPGしか連想できないゲーマーですか?
0505デフォルトの名無しさん
2006/03/28(火) 23:55:29でも数ピクセル単位で動かす事になるとおもうから、
その動かせる単位の位置を保存して、
表示する位置に変換するって感じで使えるかも。
0506デフォルトの名無しさん
2006/03/29(水) 10:10:240507デフォルトの名無しさん
2006/03/29(水) 17:34:46640×480の画面に32×32の画像をしきつめようとすると、
右端と下端がきれいにそろわず、ちょっとはみだしたりしまいます。
これって例えば何が原因で起こるのでしょうか。
0508デフォルトの名無しさん
2006/03/29(水) 17:45:27■■■■■■■■■■■■■■■■■■■■ 20x15 で普通に敷き詰められるよね?
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
■■■■■■■■■■■■■■■■■■■■
0509デフォルトの名無しさん
2006/03/29(水) 18:10:03まさか浮動小数点演算が介在してたりせんだろうな?
0510デフォルトの名無しさん
2006/03/29(水) 18:33:53>例えば何が原因で
・計算間違い
・勘違い
・お門違い
0511デフォルトの名無しさん
2006/03/29(水) 18:39:100512デフォルトの名無しさん
2006/03/29(水) 20:09:24右端は揃うんじゃマイカ?
マジレスすると、敷き詰める画像のサイズが640と480の公約数である必要がある。
0513デフォルトの名無しさん
2006/03/29(水) 20:17:270514デフォルトの名無しさん
2006/03/29(水) 22:22:060515デフォルトの名無しさん
2006/03/29(水) 23:42:360516デフォルトの名無しさん
2006/03/30(木) 09:24:330517デフォルトの名無しさん
2006/03/31(金) 01:31:23>>507
お前さ
いい加減あちこちにマルチポストするの止めろよ
もしくは解決したなら解決したと全てのスレに書いて回れ
0518デフォルトの名無しさん
2006/03/31(金) 23:45:54上も下も繋がっているんだから回る方向は同じじゃないの?
0519デフォルトの名無しさん
2006/03/32(土) 02:22:08オーストラリアの世界地図 でググってみるといいかも
方向という単語の意味の問題だから、ここはスレ違いかも
0520デフォルトの名無しさん
2006/03/32(土) 03:18:09南半球の人間はロシアが下だろと認識している
0521デフォルトの名無しさん
2006/03/32(土) 04:01:570522デフォルトの名無しさん
2006/03/32(土) 11:13:290523デフォルトの名無しさん
2006/04/02(日) 08:07:38所持アイテムを表示したりする、ウィンドウについて。
これって、どうやって実装してるんでしょうか?
ただ単に固定されてるウィンドウをつくるなら、ウィンドウのテクスチャ描いて、
それを張って、座標が決まっているのでそこに文字やらなんやら表示すれば良いと
思うんですけど、
オンラインゲームに存在するウィンドウって、windowsのウィンドウみたいに
マウスでクリックすると自由に動かすことができますよね?
これって、APIで一つ一つ子ウィンドウを作って、そこに表示させてるんでしょうか?
もしくは、他に方法とかあるんでしょうか、助言をよろしくお願いします。
0524デフォルトの名無しさん
2006/04/02(日) 08:27:330525デフォルトの名無しさん
2006/04/02(日) 08:28:25つ[DirectX]
後スレ違い。
0526デフォルトの名無しさん
2006/04/02(日) 14:00:590527デフォルトの名無しさん
2006/04/03(月) 04:30:12>>517
0528デフォルトの名無しさん
2006/04/07(金) 13:54:31「数学と算数」よりは関係あるけどな。
0529デフォルトの名無しさん
2006/04/24(月) 11:15:56そうか?
0530デフォルトの名無しさん
2006/04/24(月) 13:10:28クリッピングや論理演算は数学の範疇だと思うんだけどな
0531デフォルトの名無しさん
2006/04/25(火) 12:28:25DirectXはスレ違いだし、クリッピングの話なんてどこにも出てないが?
0532デフォルトの名無しさん
2006/04/25(火) 14:21:17窓の中の昴〜♪
0533ルカ
2006/04/26(水) 16:58:08てかこれどういう意味?符号関数ってなんですか・・・
すいませんここはプログラムとかそういうサイトっぽいけどたぶんここしか頼れる所がないので・・・
大学の数学で出てきたんですけどもしよろしければ誰か教えてください(>_<)
0534デフォルトの名無しさん
2006/04/26(水) 17:02:02まあ、マジレスすると、符号関数でぐぐればいくらでも説明出てくる。
0535デフォルトの名無しさん
2006/04/26(水) 18:13:04固有ベクトルの求め方で質問させていただきたいことがあります。
固有値はハウスホルダー変換と2分法で求めました。
その後、逆反復法を用いようと思っているのですが、
初期ベクトルの決め方と、繰り返し回数がわからず困っています。
どなたか教えていただけないでしょうか?
0536デフォルトの名無しさん
2006/04/26(水) 20:23:480537デフォルトの名無しさん
2006/04/27(木) 01:34:48sgnの読み方はともかく、意味は分かるだろう。符号関数とも書いてるから
探しているものという確認も取れるし。
そしてシグネーチャーでオッケイかを知りたければシグネーチャから思いつく
綴りを全部辞書で調べろ。
0538デフォルトの名無しさん
2006/04/27(木) 02:01:54初期ベクトルは、求めたい固有ベクトルの成分がゼロ、
もしくは、よほどゼロに近くない限りなんでもいいんじゃない?
ランダムなベクトルにするとか。
例え、偶然にも求めたい固有ベクトルの成分がゼロの初期ベクトルを
選んでしまったとしても、なんだかんだで誤差が出てくるので、なんとかなるかも。
収束は遅くなるだろうけど。
繰り返し回数は、求めたい精度に達するまでやればいいと思う。
Aを行列、aを固有値、xを求めた固有ベクトル、として、
(x,Ax)/(x,x)を計算すれば、固有値aの近似値が求まるから、
その精度を、既に求まっている固有値aと比較するとか。
0539535
2006/04/27(木) 11:26:58参考にさせていただきます。
0540デフォルトの名無しさん
2006/05/03(水) 12:47:46行列の成分の
90%が0
8%が1
2%がその他
という行列を解くには
どんなアルゴリズムを使えば高速に解ける?
こういう行列に限り高速にできるとか
そんなアルゴリズムないかなと思ったんだが。
0541デフォルトの名無しさん
2006/05/03(水) 14:35:25非ゼロの位置に特徴が無い場合は、行列をリスト表現して計算するくらいしか。
0542デフォルトの名無しさん
2006/05/03(水) 17:01:03サンクス
非ゼロの位置に規則性はないわ。
あきらめとく。
0543デフォルトの名無しさん
2006/05/03(水) 18:30:420544デフォルトの名無しさん
2006/05/05(金) 01:27:02たけしのコマネチ大学数学科
http://live22x.2ch.net/test/read.cgi/livecx/1146759594/
0545デフォルトの名無しさん
2006/05/07(日) 17:26:16非零が疎な行列だと反復法系のアルゴリズムが案外早いぞ。
0546デフォルトの名無しさん
2006/05/23(火) 08:07:13待望の翻訳!!、なる案内がきたが、定評と待望についてよろしく。
しかし\31,500とはずいぶんだな。
0547デフォルトの名無しさん
2006/06/15(木) 10:42:400548デフォルトの名無しさん
2006/07/16(日) 09:13:350549マイク ◆yrBrqfF1Ew
2006/07/19(水) 14:31:54n個の区間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)時間のアルゴリズムを与えよ。
0550デフォルトの名無しさん
2006/07/19(水) 15:07:180551デフォルトの名無しさん
2006/07/19(水) 15:23:21>549に代わってお願いします。
ぜひC++で(ry
0552デフォルトの名無しさん
2006/07/19(水) 15:23:34そういう態度だから何処へ行っても嫌われる。
0553デフォルトの名無しさん
2006/07/19(水) 21:59:50死に晒せヴォケ
0554デフォルトの名無しさん
2006/07/20(木) 07:11:30>>551
http://kansai2channeler.hp.infoseek.co.jp/cgi-bin/joyful/img/2409.cpp
方針は以下:
(1) は自明.
(2) はソートして左から数える.ソートしたおかげで単調性が得られ,
一度交わらなくなったらそれより先を調べる必要がなくなる.
(3) は (2) でどこまで調べないといけないかを二分探索を行う.
0555デフォルトの名無しさん
2006/07/20(木) 08:09:59http://sourcepost.sytes.net/sourcepost/sourceview.aspx?source_id=28099
(i) 区間対の個数はO(nn)、重なりの有無の判定はO(1)だから、全体でO(nn)。
(ii) >>554に同じ。
(iii) 与えられた区間の始点と終点を列挙し、ソートし、各点に何個の始点・終点が重なっているか調べる。
区間のネストの数を把握しながらこの列を走査して重なりの総数を得る。
実際のコードでは始点・終点の個数の計算と最後の操作を一つのループで行っている。
0556マイク ◆yrBrqfF1Ew
2006/07/20(木) 19:23:08・・・
2chネラにこんなことを言うのもあれだが
お前らすげーな。
どこの金子勇的風貌の京大院生だ?
0557デフォルトの名無しさん
2006/07/20(木) 22:30:46今からでも勉強してみたらどうでい。
0558デフォルトの名無しさん
2006/07/30(日) 08:15:10三次元上の多面体の体積を求めるには何か方法はありませんか?
0559デフォルトの名無しさん
2006/07/30(日) 10:47:514面体OABCの体積は
1/6( ax by cz + ay bz cx + az bx cy - az by cx - ax bz cy - ay bx cz )
0560デフォルトの名無しさん
2006/07/30(日) 10:49:240561デフォルトの名無しさん
2006/07/30(日) 11:00:040562デフォルトの名無しさん
2006/07/30(日) 12:06:07外積
0563デフォルトの名無しさん
2006/07/30(日) 16:12:500564デフォルトの名無しさん
2006/07/30(日) 16:14:360565デフォルトの名無しさん
2006/07/30(日) 18:27:31体積を求めたい多面体が凸でない場合は4面体の和に分割するところが面倒かもよ
凸多面体でも分け方によって誤差が変わってくるからそんなに簡単ではないかもよ
0566デフォルトの名無しさん
2006/08/06(日) 16:49:430567デフォルトの名無しさん
2006/08/06(日) 18:06:170568デフォルトの名無しさん
2006/08/06(日) 20:54:07対象が凸多面体で原点を多面体の内部か表面にとればそれは正しい
0569デフォルトの名無しさん
2006/08/06(日) 21:45:34符号付き、というのはたぶん
原点から見て裏面が見えてるときはマイナスの体積、ということだろうから
そうすると原点に関して星状であるという条件はいらなくなるよ
Aが星状というのは原点とAの任意の点とを結ぶ線分がAに含まれることね
0570デフォルトの名無しさん
2006/08/06(日) 23:55:550571デフォルトの名無しさん
2006/08/07(月) 12:48:29便利なライブラリ or 公式とかある?
もしもなかったら、
衝突判定がしやすい曲面の表現方法教えてもらえませんか。
0572デフォルトの名無しさん
2006/08/07(月) 21:29:190573デフォルトの名無しさん
2006/08/07(月) 21:41:17形状を簡単に変えることができるものであればいいんですが、
球も表現したかったのでNURBSを考えていました。
0574デフォルトの名無しさん
2006/08/08(火) 11:21:05NURBS曲面かどうか調べてないが、Newtonとか言うGAME向けの物理演算ライブラリがある
0575デフォルトの名無しさん
2006/10/01(日) 18:50:02当選口指数 {数.数} = 当選口数 ÷ 販売口数 × 100万倍
http://sea.advenbbs.net/bbs/yybbs.cgi?mode=res&no=1593&id=loto66
0576デフォルトの名無しさん
2006/10/04(水) 20:30:010577デフォルトの名無しさん
2006/10/12(木) 11:02:110578デフォルトの名無しさん
2006/10/12(木) 11:02:530579デフォルトの名無しさん
2006/11/03(金) 10:21:15マウス(2次元)座標をゲーム画面(3次元)の座標に
変換したいんだけれど、
やり方がわからないです。
[やりたいこと]
大航海時代Onlineのように、
「キャタクターの行先(目的地)」を「地面をマウスクリックする」ことにより
指定し、その後、キャラクターが動き出す。
ということをやりたいと考えています。
0580デフォルトの名無しさん
2006/11/03(金) 10:29:29DirectXを使っているならDirectX隔離スレへ。
その他3D表示ライブラリを使っているならそれらのスレへ。
そうでないなら自分で勝手に計算汁。
0582デフォルトの名無しさん
2006/11/23(木) 14:10:24ちなみに
・グレゴリウスの定義
・西暦1年1月1日は月曜日
・1年365日だが、閏年は366日
・4または400で割り切れる年は閏年
・100で割り切れるのは閏年ではない
もしわかった方いましたらお願いします
0583デフォルトの名無しさん
2006/11/23(木) 14:15:04C++の擬似コードだと、こんな感じかなぁ・・・
main(){Youbi 曜日; カレンダーを見る(新聞の日付, &曜日);
if (時計の時間<新聞配達時間) 曜日.進める;
return 曜日;}
0584デフォルトの名無しさん
2006/11/23(木) 15:01:22普通はAPIを使う。
こんな所で質問してる奴が組んだ物は怖くて信用おけないから。
関係無いけど、ANSI C準拠のtime()関数系の1970/01/01 00:00:00からの通算秒で
有効桁31bitの場合は、全部閏年で換算しておk。
0585デフォルトの名無しさん
2006/11/23(木) 15:06:49{
if (month < 3) {
year--; month += 12;
}
return (year + year / 4 - year / 100 + year / 400 + (13 * month + 8) / 5 + day) % 7;
}
0586デフォルトの名無しさん
2006/12/10(日) 12:20:07にあるprime-test()で多倍長整数の場合、最初に判定を間違える数って何?
0587デフォルトの名無しさん
2006/12/10(日) 13:01:52マジレスすると西暦1年1月1日にはまだグレゴリオ暦は発明されていなかった罠。
0588デフォルトの名無しさん
2006/12/10(日) 15:35:33遡及的に定義できる。
0589デフォルトの名無しさん
2006/12/10(日) 17:39:40ミラーラビンテストでぐぐったら強擬素数ってのが出てきた。
あとは頑張れ。
0590デフォルトの名無しさん
2006/12/10(日) 20:48:11>>589 強擬素数はフェルマーテスト
0591デフォルトの名無しさん
2006/12/11(月) 02:24:03適当いったらだめなんじゃね。
0592デフォルトの名無しさん
2006/12/17(日) 00:24:50他に適当なスレが見つからなくて。
私は会社の事務員なのですが、
手間のかかる手作業がありまして、パソコンにやらせたいのです。
職場の人員のスキル的に見て、自分でやるしかない状況です。
エクセルのVBAで何度かプログラムを組んだことのあるレベルです。
以下が問題の概要です。
1から100までの整数の並びがある。
そのうちランダムに20個程度の数が欠番となっている。
この状態から、残った数を並べ替え、
隣り合う数同士の差の絶対値が全体でなるべく等しくなるように、
つまり等間隔に並ぶようにしたい。
ここ一ヶ月くらい悩んでみましたが、アイディアが浮かびません。
何か既存のアルゴリズムで、応用できそうなものを紹介してもらえませんか。
0593デフォルトの名無しさん
2006/12/17(日) 00:31:57判りにくいな。例えば、1, 2, 90, 91があったら1, 90, 2, 91に並べ替えるってことか?
#前者は隣同士の差が1, 88, 1で、後者は89, -88, 89だから後者の方が差の絶対値はほぼ均等になる。
0594デフォルトの名無しさん
2006/12/17(日) 00:39:09総当りで調べるとかかなぁ。
「隣同士の差」の分散辺りを評価基準にして、全部の組み合わせ試して一番分散の小さいのを選ぶ。
0595592
2006/12/17(日) 00:52:03わかりにくくてすみません。
実例を出すべきでした。
まさに、その後者のように並べたいのです。
今は人間がエクセル上で手作業であれこれと
数値を入れ替えし、納得いく結果になるまで
3日はかかってしまいます。
>>594
分散、とは統計学の用語ですか。
調べてみます。
0596デフォルトの名無しさん
2006/12/17(日) 06:08:03単純に、小さい値と大きい値で2つに分けて、ソート後交互に配置じゃだめなん?
0597デフォルトの名無しさん
2006/12/17(日) 08:18:01>596でよさそうだけど。
0598デフォルトの名無しさん
2006/12/17(日) 12:29:101, 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
0599デフォルトの名無しさん
2006/12/17(日) 16:18:54それならソートしただけのデータでもインジャネ?
相対的な距離に関する最適解は総当りじゃないと求まらないんじゃない?
しかも「分散の小さい物」という条件だから足切りも出来ない。
0600デフォルトの名無しさん
2006/12/17(日) 22:37:120601デフォルトの名無しさん
2006/12/21(木) 21:47:49整数乗算についてなんですが、2048ビット程度だと、fftかまして乗算しても恩恵はない。
と聞きました。実際のところはどうなんでしょうか?どなたか検証サイトなどご存知ではありませんか?
0602デフォルトの名無しさん
2006/12/22(金) 01:05:13乗算が加算よりどれだけ遅いかも評価に入れなきゃいけないだろうし。
0603デフォルトの名無しさん
2006/12/22(金) 01:45:310604デフォルトの名無しさん
2006/12/22(金) 12:31:04RSA の乗算+剰余演算の場合、FFT による方法じゃなくて、
Montgomery multiplication っていう高速化手法があるよ。
検証したことはないけど、2048 ビット程度では FFT の効果はない
あるいは逆効果ってのはそうだと思う。
せめて1万ワードくらいの長さはないと。
0605デフォルトの名無しさん
2006/12/22(金) 12:37:35上位・下位半分ずつ、再帰的に計算することで
乗算回数減らす方法もある。
下位桁の乗算結果使って、上位桁の乗算をサボる。
昔、RSA のプログラムの最適化の仕事したことあるけど、
「単純に畳み込みで乗算 & 引き放し法的な剰余」を
「再帰乗算 & Montgomery 乗算法」に変えたら
鍵長256bitでも動作速度を250倍以上速くできた。
0606デフォルトの名無しさん
2006/12/22(金) 12:41:24速度250倍達成は、あと、
バイト単位で多倍長整数演算してたのを、
1ワード32ビット単位に変えたのも含めてだわ。
0607デフォルトの名無しさん
2006/12/27(水) 13:26:09CADのプログラムですが、CADの掲示板で質問してるけど、3次元の幾何計算の話なので、あまり回答が返ってきません。
誰か助けてぇえええ!!
http://www.mura.sh/bbs/wwwbbs.cgi?print+200612/06120169.txt
0608デフォルトの名無しさん
2006/12/27(水) 13:52:380609デフォルトの名無しさん
2006/12/27(水) 15:42:41そうはいってもあそこはCADユーザーの掲示板で、
プログラム関係はほとんどレスがつかないんですよ。(´・ω・`)
餅は餅屋ということで、こっちにも質問してみる。
0610デフォルトの名無しさん
2006/12/27(水) 16:54:00マルチポストうぜぇよ。
見たら解答ちゃんと出てんじゃねぇか。
それを理解出来ないお前の頭までは救えねぇよヴォケ。
0611デフォルトの名無しさん
2006/12/27(水) 16:58:17横のソリッドがX軸方向に伸びてるならいい。
しかしそれ以外の方向ならお手上げなんだけど。
0612デフォルトの名無しさん
2006/12/27(水) 22:31:43回転してX軸に来るようにしたらいいんじゃね
0613デフォルトの名無しさん
2007/01/06(土) 04:41:040614デフォルトの名無しさん
2007/01/06(土) 18:10:16例えば天気図の気圧の図があったとして、
高気圧が何個あって、どの座標に中心があるかを調べたい
というような問題です。
ピークがない場合やたくさんある場合があります。
こういう計算はどういうジャンルの学問になるんでしょうか。
0615デフォルトの名無しさん
2007/01/06(土) 21:05:180616デフォルトの名無しさん
2007/01/06(土) 22:14:53いいヒントをありがとうございました。
「山登り法」というキーワードでいろいろ探せそうです。
0617デフォルトの名無しさん
2007/01/07(日) 22:14:230618デフォルトの名無しさん
2007/01/07(日) 23:02:490619デフォルトの名無しさん
2007/01/09(火) 11:29:24∂z
----
∂x
∂z
----
∂y
0620デフォルトの名無しさん
2007/01/09(火) 11:36:150621デフォルトの名無しさん
2007/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個よりも少し多く必要だったりするでしょうか。
今こんなことをやりたいがためにリードソロモン符号についての勉強を始めたのですが、
今結論だけを端的に知りたいです。
あるいは、他の符号化についての候補などあったら知りたいです。
よろしくです。
0622デフォルトの名無しさん
2007/01/19(金) 15:31:54やっぱりプログラム、コンピュータの基盤となるのは数学だから
数学的な知識があるとないじゃ、雲泥の差だという事を、つくづく思い知らされたよ。
とりあえず、ちょっとでもスキルアップ出来る様に、勉強し直そうと思うよ…
0623デフォルトの名無しさん
2007/01/19(金) 16:56:360624デフォルトの名無しさん
2007/01/19(金) 17:04:45最近は分数が出来ない文系の大学生も多いらしいからねぇ。
書く言うおれも、中学レベルの数学が既にあやしい…orz
0625デフォルトの名無しさん
2007/01/19(金) 17:49:25にほんごもあやしいようだね
0626デフォルトの名無しさん
2007/01/19(金) 19:49:50Luigi Rizzo氏のHPにReed-Solomon符号を使ったデータ復元プログラムのソースが置いてあったような...
確か、制限はn+m<2^8 またはn+m<2^16。
復元もn+m個中n個を受信すれば大丈夫だったはず。検索して調べてみてください。
0628621
2007/01/20(土) 12:46:41情報ありがとうございます。
検索の結果それっぽいページにたどり着いたものの、いかんせん英語が不得手でして・・・
実装例は自分が行き詰ったときに見ることにして、とりあえずやりたいことはできるのだろうということはわかりました。
n+m<2^8 とか n+m<2^16 とかの制限は、多項式の係数を何ビットで計算するかによるのでしょうね。
8bit で計算すればエンディアンの問題が無くて都合が良さそうですが、できるだけ多ビットで計算した方が計算回数が
少なくて高速になりそうですね。
# 実は 8bit しか念頭に置いていませんでしたが、実験として 32bit くらいでやってみようかと思いました
0629デフォルトの名無しさん
2007/01/21(日) 23:08:41国語や英語はできます。論理性はある方だとは思います。理系の教科は勉強し
ても全然伸びませんでした。理系の科目が出来ないのは、恐らく、数学的思考が出来ていない
為だと思われます。
やってるうちに色々自分の欠点に気づき、その欠点を克服したいと思っています。
プログラミングが数学のスキルアップによいと聞き、それが欠点の克服に繋がると、
思い、やってみたいのですが、なにから始めたら良いでしょうか?
何か他にもアドバイスがあったら教えて下さいませ。
0630デフォルトの名無しさん
2007/01/21(日) 23:20:090631デフォルトの名無しさん
2007/01/22(月) 07:54:17今受験生ってことは、もう残り1,2ヶ月?
プログラミングに費やす時間あったら、他のことにまわしたほうがいいぞ
長い目で見ればプログラミングに手を出すのも悪いことではないけど、一朝一夕ではどうかなぁ
0632デフォルトの名無しさん
2007/01/22(月) 08:05:35文系か理系かで分けるんじゃなくて
論理的に考えられるかどうかが重要
それが出来ないと何やってもだめ
0633デフォルトの名無しさん
2007/01/22(月) 08:17:07何を作ろうとしてもライブラリ。
六法全書以上に分厚い資料との格闘。 記憶力や事務処理能力の方が要求される分野かもね
上で出たLOGOやpascalレベルなら別だろうけど、実用プログラミングはね
0634デフォルトの名無しさん
2007/01/22(月) 09:46:34現代とかどうとか関係なくそもそもプログラミングと数学は全然別物だと思うが
0635629
2007/01/22(月) 15:01:04このスレの上の方でBASiC等が数学のスキルアップに良い。と書かれてますが、
実際どうなんでしょう?
英語は外人と喋ったり、チャット・メールをして上達しています。
数学も、こんな感じで実用面からアプローチすれば上達するかな?
と、思っています。
>>630 一度調べてみます。有難う御座います。
0636デフォルトの名無しさん
2007/01/22(月) 15:17:21それにいくら外人と話してしても、日常会話は上達するだろうが、英語学者にはなれないだろう。
これも近い例だと思う。
0637デフォルトの名無しさん
2007/01/22(月) 15:17:47どこのこと?
>>135ならセンター試験のプログラムの問題の話だが。
>数学も、こんな感じで実用面からアプローチすれば上達するかな?
プログラミング(の大部分)は高校数学の応用ではない。
確かに数学とプログラミングは関連しているし、考え方も共通性がない訳じゃないけど、
受験勉強としてやるには効率が悪すぎると思う。
現代国語の対策として小説の書き方を勉強するようなものじゃなかろうか。
0638デフォルトの名無しさん
2007/01/22(月) 23:20:20片手間でやろうとした奴はそうとう焦っただろ
二分法がどういうものか分かってないと話にならない
まあ数値計算の章では普通に教科書にも出てるから
本当にやってる人ならなんでもなかったのだろうけど
統計の問題もそうだけど今年は「冷やかしお断り」っていう
センターの意図がすごい感じられたよ
0639デフォルトの名無しさん
2007/01/22(月) 23:28:190640629
2007/01/23(火) 00:02:36やはり地道にがんばるのがいいんでしょうか?
何とかして数学的思考を出来る様にならなきゃいけないんですが・・・
0641デフォルトの名無しさん
2007/01/23(火) 00:13:28コンピュータのために数学を学ばせているような気がする。
0642デフォルトの名無しさん
2007/01/23(火) 00:13:44「数学的思考」についてお前の思うところを述べよ
0643629
2007/01/23(火) 00:34:40思考です。又、余計な言葉は入れず必要な事だけを述べて行く・書き下して
行く様にする思考です。
0644629
2007/01/23(火) 00:36:24だと思っています。
0645デフォルトの名無しさん
2007/01/23(火) 01:50:51数字や記号や式を具体的なものとして捉える力のほうが大事。
俺の持論だけど。
0646デフォルトの名無しさん
2007/01/23(火) 02:41:21なんのことか分からん罠
0647デフォルトの名無しさん
2007/01/23(火) 06:25:56おまいさんは 脚の総本数 = 4 * 亀の個体数 + 2 * 鶴の個体数 から遣り直した方がいい。
0648デフォルトの名無しさん
2007/01/23(火) 07:19:32足の数はぜんぶで、26本。
ツルは何匹、カメは何匹?
0649デフォルトの名無しさん
2007/01/23(火) 07:21:050650デフォルトの名無しさん
2007/01/23(火) 07:30:39【 課題 】プログラミング演習
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以上
【 期限 】月曜17:00
【 補足 】
できれば数字の組み合わせが重複するケースや
計算の途中で既知の値(系列)に到達した場合は
それ以降の計算を省略する等して高速化を図ること
0651デフォルトの名無しさん
2007/01/23(火) 23:26:57受験数学が駄目なのか、微積分とか代数幾何とかの
教科書レベルの基本的概念を掴むところから駄目なのか
想定している大学のレベルはどれくらいか
たとえば、物理の物体の運動を数式の暗記ではなくて、
微積分の概念を用いて理解するなら数値積分の結果を
グラフにしたりするのがよくあるパターンかと思うが、
それならプログラミングなんてする必要はなく、excelで
も十分な罠。
0652デフォルトの名無しさん
2007/01/24(水) 02:21:07でかい配列作って計算結果の一覧表を作ればおk
0653629
2007/01/24(水) 22:15:45>>651 受験数学は駄目ですね。これは基本概念が掴めていない所為だと
思います。志望は阪大・大阪市立の医学部です。
excelですか。一度調べてみます。有難う御座います。
0654デフォルトの名無しさん
2007/01/25(木) 10:38:57>論理性はある方だとは思います。
>数学的思考が出来ていない
矛盾してないか?
0655デフォルトの名無しさん
2007/01/26(金) 03:07:02若干のズレがあるからな。
たとえば、「明日の花見は、雨ならば中止します。」と幹事が言ったとして、
「では晴れたらどうするんですか?」なんて質問をしたら、馬鹿呼ばわりされること必至。
(日常会話の「ならば」は、必要十分条件を意図して使われていることが多い。)
数学(プログラム)的には、この質問は至極自然かつ妥当なものなのだが。
0656デフォルトの名無しさん
2007/01/26(金) 10:25:47基本概念が掴めていなければ,応用もできんわな
そんなヤツに医者になってほしくない・・・
0657デフォルトの名無しさん
2007/01/26(金) 11:03:26でも、その国語と数学のずれを認識できない奴のことを、
はたして論理性があるといえるのか?
0658デフォルトの名無しさん
2007/01/26(金) 22:47:39「数学では論理を正確にするが,日常ではルーズにする」とかな
そうじゃなくて言葉の意味が違うだけ
力学で言う「仕事」が普通に言う「仕事」と違うっていうようなもの
0659629
2007/01/26(金) 23:54:56数学的思考とは、物事を順序立てて考える事であり、具体的な文字を使わず、
抽象的な文字を用いる 思考です。又、余計な言葉は入れず必要な事だけを述べて行く・書き下して
行く様にする思考です。
こう言いました。
この私が呼んでる数学的思考という能力が劣っていても、論理性が高いという事は
ありえます。
先ず、具体的な文字を使わず、抽象的な文字を使う。という事と論理性とは
なんの関係もありません。(まぁ、厳密に調べて言ったら、多少なりとも関係
しているでしょうが)
次に、余計なことを述べず必要な物だけを述べる。と言うのも、関係有りません。
別に不必要な物を付け足しても、論理的な展開はでき、論理力が高いという事はあります。
最後に、物事を順序立てて考える。と言うことですが、これは正直よく反論出来ません。
別に物事を順序だてて考えなくても、論理力が高いと言う事はある。脳みそは複雑でバラバラな
事でも統合的に処理できている。というぐらいしか言えません。
また、頭の良い人から論理性が高い。と言われますし、実際国語の点数も良いです。
なので、私の論理性が高くない、と言う事はあまり考えるべきでないと思います。
0660629
2007/01/27(土) 00:04:12考えているつもりです。
>>657 >>658 何処でどう私の事を、その国語と数学のずれを認識できない奴と、
思ったのですか?
もし、論理性の欠如から来る、ミスなら気をつけて下さい。
話がいらん、何も有益が無い方向に曲がってしまいます。
0661デフォルトの名無しさん
2007/01/27(土) 00:11:080662デフォルトの名無しさん
2007/01/27(土) 01:18:26本人は論理的なつもりなのかも知らんが、誤った前提を元に論理も糞も無いだろうに。
0663629
2007/01/27(土) 02:19:15>>661 演算は出来ます。というか誰も演算できないなんて言ってません。
まぁ、大体このスレでの有益な意見は出たので、これで失礼します。
一部の人が怒りっぽく、荒れてきましたし。
私が言ってないことを言ったかの様に述べたり、確定してない事を確
定してると思い込むのは、恐らく、頭脳から来る頭の悪さでは無く、
人格の欠陥から来る頭の悪さでしょう。
人格の欠陥から来る頭の悪さは大概自己愛の欠落からきています。
自己が自分や他人(特に両親)から肯定されていないんですよ。
一度両親の元へ行って、愛を与えて貰いに行くのがいいと思います。
もしくは、両親との不仲を直し、互いに愛し合う関係になったほうが
良いです。それで、もしかしたら直ると思いますよ。ではノシ
0664デフォルトの名無しさん
2007/01/27(土) 02:30:09論理性があって計算もできる奴が、どうして
>数学出来ません。物理できません。化学の計算できません。
なんてなるのかが不思議でしょうがない。
0665デフォルトの名無しさん
2007/01/27(土) 03:07:16どこも論理的じゃないから、不思議でもなんでもない。
0666デフォルトの名無しさん
2007/01/27(土) 10:23:34「数学的思考=抽象化」ということか
オレは「数学的思考=論理的思考」だと思ってたんだが
そもそも,抽象化は考える前の「単純作業」だと思ってたし
その辺,みんなはどうよ?
>別に物事を順序だてて考えなくても、論理力が高いと言う事はある
これはさすがに違うんじゃねえか?
もう見てないのかな
0667デフォルトの名無しさん
2007/01/27(土) 10:26:52アレな人の典型のようにも思える
0668デフォルトの名無しさん
2007/01/27(土) 13:53:01呼びたいのか?w
俺はもう帰ってきて欲しくないが
0669デフォルトの名無しさん
2007/01/27(土) 15:25:34日本語の勉強を、中学校辺りから遣り直してください。
仮令国語の点数がよくても、このような文章としてなっていないものを書いているようではお郷が知れます。
0670デフォルトの名無しさん
2007/01/27(土) 15:30:12大丈夫、この手の香具師は自意識過剰だから自分の話が出ないか気になって必ずヲチしている。
プライドが充分に低いか高いかしなければ二度と書き込まないでくれるはずなのだが……
0671デフォルトの名無しさん
2007/01/27(土) 17:05:370672デフォルトの名無しさん
2007/01/27(土) 21:17:51物事をシンプルにまとめる力がないと、プログラムでは苦労するだろうな、
と、経験的に思う。
0673629
2007/01/27(土) 21:19:32>>670さんは中々鋭いですね。私はナルシストです。
では本題に、
>>666 実際、多くの受験生が苦手としている科目が数学と物理です。
これは人が抽象的な思考をするのが苦手だからだと思います。
より具体的に説明された方が理解しやすいのは経験からも判ると思います。
ので、抽象化を'単なる'単純作業と捕らえるのは違うと思います。
しかし、出来る人からすると、単なる単純作業なんでしょうね。
>別に物事を順序だてて考えなくても、論理力が高いと言う事はある
これはさすがに違うんじゃねえか?
これについては、一例を出しますと、
ある繋がった1〜10迄の事があるとします。順序だてて考えて行くとは、
1.2.3.4.5.6.7.8.9.10と一から十迄規則的に考えて行く思考法の事です。
一方の私がそうでなくてもいける、と言ったのは、
2.5.1.4.6.3.7.9.10.8.と不規則ながらもこれを頭の中で正確に関連づける思考法が
人には出来るという事です。不規則な順序で進めても、それらを正確に関連づけて
最終的に上手く統合出来る。そしてそれは確固とした論理性であり論理力である。
と言うことが言いたかった訳です。
0674デフォルトの名無しさん
2007/01/27(土) 22:01:54>670によると、プライドが高いか低いかどっちかだそうだけどどっち?
0675デフォルトの名無しさん
2007/01/27(土) 22:05:17これは非常に簡単な理由がある。
こららは所謂理数系と言われる教科なのに、教えている教師は
所謂文系と言う数学物理を苦手とする傾向のある連中が多いからだ。
0676デフォルトの名無しさん
2007/01/27(土) 22:17:17×ナルシスト
○ナルシシスト
0677デフォルトの名無しさん
2007/01/27(土) 22:52:52大学以降は抽象化されたものを扱う機会が増えたが、どちらかというとそういう考え方は苦手だった。
そんな自分でも高校のときは数学や物理には苦労しなかったので、
受験レベルの内容に抽象的な思考はそれほど必要ないと考えている。
>>629には数学や物理ができない原因が本当に抽象的な思考の部分にあるのかを考えてみてほしいな。
あと元々はプログラミングが役に立つかという話だったけど、
プログラミングに手を出したコストに対して得られる対価が小さいかなと思う。
もっと楽にスキルを上げる方法はありそう。
少しでも興味を持ったのなら、受験は抜きにして趣味でやってみるのもいいかもしれない。
嫌だったら辞めればいいだけの話だし。
0678デフォルトの名無しさん
2007/01/27(土) 22:57:28それが終わったら赤に行け
0679661
2007/01/27(土) 23:40:15そうか、すまん。
お前が論理性があるというので、学校の数学から論理性を取り除いて残るものと言ったら計算しか思いつかなかった。
それだけのことだ。
0680デフォルトの名無しさん
2007/01/28(日) 11:10:15→ナルシシズム、ナルシシスト
ナルス
→ナルシズム、ナルシスト
http://ja.wikipedia.org/wiki/%E3%83%8A%E3%83%AB%E3%82%B9
お前同郷かw 頑張れよw
抽象化じゃなくて一般化という用語を使えばもっと楽だろうに。
〜ist の用法を一般化したらナルシシストになりそうなことは想像できるだろう。
ナルシスを知らなければしょうがないが。
0681デフォルトの名無しさん
2007/01/28(日) 14:02:550682デフォルトの名無しさん
2007/01/28(日) 15:52:17○machiavellian
×machiavellist
0683デフォルトの名無しさん
2007/01/28(日) 15:56:560684デフォルトの名無しさん
2007/01/30(火) 05:08:18かなり役に立った。特に微積に関してはかなり深く勉強できたと思う。
でも今はCOBOLのプログラマだ・・・orz
全く数値計算法なんて役に立たねえっす。
0685629
2007/01/31(水) 05:51:03良かったら教えてくださいまし。
0686デフォルトの名無しさん
2007/01/31(水) 10:30:540687デフォルトの名無しさん
2007/01/31(水) 11:03:43ttp://www.netlib.org/lapack/index.html
0688686
2007/01/31(水) 11:13:39LAPACKでググったら、直線のサンプルもすぐ見つかりました。
ttp://www.sip.eee.yamaguchi-u.ac.jp/kou/lapack.html
実際にやりたいのは真円のフィッティングだから、もうちょっと調査しないと。
0689デフォルトの名無しさん
2007/01/31(水) 11:33:59数値解を繰り返しで求めるか
ttp://www.tensyo.com/urame/prog/linealgo.htm
中心からの距離の2乗の誤差の2乗で代用すれば解けるそうだ
0690デフォルトの名無しさん
2007/01/31(水) 11:40:36ニューメリカルレシピ・イン・シー
0691686
2007/01/31(水) 11:42:04やりたいのは4点(真円になるか”?”)から円の中心と半径を決めたいです。
>中心からの距離の2乗の誤差の2乗で代用
ちょっと難?
関係ないけど今回の場合、
ttp://okwave.jp/qa1812202.html
は参考にならないような気がしました。
0692デフォルトの名無しさん
2007/01/31(水) 12:02:124点なら最小2乗を使うほどにはデータ量が多くない。
4点から3点を取り出しては公式で中心を求めて
4つ求まった中心の平均を出したらどうだろ?
その中心から4点への平均距離を求めて、それを半径の代用としたら?
0693デフォルトの名無しさん
2007/01/31(水) 12:09:42>ちょっと難?
これはデータの特性にもよるだろな。 使った感じでは十分だよ。
もっとも画像からのデータでデータ量は十分あったが
何より式一発で求まるのは非常にありがたい
0694691
2007/01/31(水) 12:19:57ごめんなさい、書き方が悪かったです。
そのやり方が悪いという意味じゃなくて、文章難しくて、
その文章どおりにするには何をどうすれば良いのか考えなきゃという意味です。
>>692
その通りやってみたんですが、トンデモない円が出来て(3点だと上手く行く)。。。
どういう算数でやろう?
最小二乗法は?
「中心からの距離の2乗の誤差の2乗で代用」ってどういう意味? ← 今ここ
といった感じです。
0695デフォルトの名無しさん
2007/01/31(水) 12:22:31最小2乗を使ってもとんでもない結果しか出ないよ
これがデータ100点あって1個だけ飛び離れてるような状態ならマシだけどさ
0697デフォルトの名無しさん
2007/01/31(水) 13:38:16もし論文なんかで出したら一笑されるぞ。
少数観測点からデータを復元する話は GPS などの位置計測で
よくある話で、そういう場合にまともな推定をしようとすれば
観測点の状態が分かってることが必要になる。
たとえば、真の値からのずれがある確率分布に従うと仮定を
置けば、最尤推定などによって最もそれっぽい場所が出る。
正規分布に従うと仮定してよい場合は最尤推定は最小自乗法に
一致するため、多くの GPS ではこの仮定を(知ってか知らずか)置いている。
多量のデータで最小自乗法が有効なのは、中心極限定理によって
分布が正規分布に漸近するため。
0698691
2007/01/31(水) 13:40:21>★最小2乗法による円弧推定
ってそのまんまのタイトルがあった。
3点から円を求めるのもこのページから見つけたのに気付かなかったOTL
0699691
2007/01/31(水) 13:43:33>何も考えずにたった4点で最小二乗法をやるのは馬鹿げてる。
>もし論文なんかで出したら一笑されるぞ。
了解、納得でつ。
必ずしも円になるのか分からないデータの解析なので、
逆にデータ解析結果から観測点が円にならないという結論を出しますが。
ある部分は円、ある部分は円にならず従って応力受けてる、みたいな。
0700デフォルトの名無しさん
2007/01/31(水) 13:48:36これはいえない。最小自乗法では円にならなくても、ほかの方法では
円になる場合があるかもしれない。
0701699
2007/01/31(水) 13:57:36じゃあ、ほかの方法教えて。
ところで、
>★最小2乗法による円弧推定
を読んだら、
>この3つの方程式を解けば良い.
で終わってる。
解かないといけないのか。
0702デフォルトの名無しさん
2007/01/31(水) 14:12:46上で述べられていた任意の三点の中心の平均や重心は有力。
また、697 で述べたように、最尤推定はひとつの別の方法。
ほかにも最小自乗や最尤推定を部分に含むが、適当な
ペナルティ関数を設定して最小化問題として解く手もある。
最尤推定やペナルティ法は、円になるかどうか分からないものの
判定なら、おそらく最小自乗を用いるよりも良い結果が得られる。
(最小自乗では、一直線上のデータは半径が(ほぼ)無限大の円と
判断されるが、通常の応用ではそれを弾きたいことが多い)
0703701
2007/01/31(水) 14:33:09これらはライブラリにもなってなさそうだね。
実装難しそう...
0704701
2007/01/31(水) 18:42:23>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)
>この3つの方程式を解けば良い.
やっぱこれ何度読んでも分からないんだけど、
4点をどうやって代入するわけでしょうか?
解答が目の前にありながら使いこなせないorz
0705デフォルトの名無しさん
2007/02/01(木) 12:41:01X → 肺i XX → 肺i*xi XXX → 肺i*xi*xi
Y → 輩i YY → 輩i*yi YYY → 輩i*yi*yi
以下略
の意じゃね? xi, yi が配列要素で直線上にない3点が最低限必要と。
0706デフォルトの名無しさん
2007/02/05(月) 17:23:56後は、1) 2)式から a,b2変数の連立一次方程式を解けば中心が求まるでしょ?
それを3)に代入すれば半径も求まる
0707デフォルトの名無しさん
2007/02/14(水) 10:09:264点中 3点で求まる円から、 残った1点までの円からの距離2乗(|中心からの距離-半径|)を出して
それを4つ加算した結果、
またはそれを半径で割って正規化したもの
でいいんじゃないのか?
0708デフォルトの名無しさん
2007/02/28(水) 21:55:010709デフォルトの名無しさん
2007/03/25(日) 15:02:16for(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 の円の時にしか点が等間隔になりません。
どうしたら良いでしょうか?
知っていなくてはいけない公式などありましたら教えてくださいませ。
0710デフォルトの名無しさん
2007/03/25(日) 15:24:580711デフォルトの名無しさん
2007/03/25(日) 15:26:040712デフォルトの名無しさん
2007/03/25(日) 15:44:51それを等間隔に分割というのは、その厄介な問題に挑戦するという事になる
数式では基本的に解けないから、数値積分で頑張るしかないよ
0713デフォルトの名無しさん
2007/03/25(日) 15:46:43int x = cx + W * cos(i * 2 * PI / N);
int y = cy + H * sin(i * 2 * PI / N);
でどうでっしゃろ?
0714713
2007/03/25(日) 15:47:200715デフォルトの名無しさん
2007/03/25(日) 16:06:21それつかって円周を求めた(ここではLとおく)として
始点(仮に中心からx軸方向に幅/2移動した点としよう)から
円周にそってl ( = n×L, nは適当な実数)分移動したときの座標(x,y)ってのは楽に求まるの?
精度気にしないのならこれでいけると思うんだけど、どうでしょうか?
0716デフォルトの名無しさん
2007/03/25(日) 16:20:25Nの個数を30倍くらいに増やし一度テーブルを作成して
補間で、近似値を求めるのが一番簡単だろう
0717デフォルトの名無しさん
2007/03/25(日) 16:40:25そのテーブルの中味と、
何のデータを元にして、
何の式で補間するか書いてください
お願いします(-∧-)
0718709
2007/03/25(日) 16:51:16>>715 円周に沿って l の座標を求める方法がよくわからないです。すみません。
>>716 で言う方法かはわかりませんが、自分ならこうやりますがどうでしょう?
1. 30N の全ての点の座標を求たテーブルを作る
2. 隣り合う点との距離を三平方の定理で求めることを全ての点について行い、円周の近時値を計算する
3. 円周を N等分した長さを求め、これを当初の N個の頂点間の距離 d とする
4. 30N のテーブル中の頂点の中から、それぞれの点の円周上の位置(始点からの距離)に近いものを選んで点を打つ
0719デフォルトの名無しさん
2007/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
を角度にするという感じ
0720デフォルトの名無しさん
2007/03/25(日) 17:04:31>円周に沿って l の座標を求める方法
楕円積分の逆関数だから、楕円関数使って求まる気がする。
もちろん、楕円関数の値は数値的に計算するものだけど。
0721719
2007/03/25(日) 17:13:07double 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++;
0722719
2007/03/25(日) 17:16:37ax = 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);
}
0723709
2007/03/25(日) 17:30:02丁寧にどうもありがとうございます。
よくわかりました。
>>720
まずは楕円積分というのを勉強しなくてはいけなそうです。
おいおいスキルアップしたいと思います。
>>710
よろしければそのスレを教えてもらえないでしょうか。
放物線を速さ一定で進むってのも興味あります。
0724デフォルトの名無しさん
2007/03/25(日) 17:44:59ttp://www.tensyo.com/urame/prog/ALGO.HTM
の後ろの方に計算方法が書いてあったけど
コレは今回の問題には応用できないな
0725しょう
2007/03/25(日) 20:12:390726デフォルトの名無しさん
2007/03/25(日) 20:19:03宿題は他をあたりなさい。
0727デフォルトの名無しさん
2007/03/25(日) 20:22:522) 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
0728デフォルトの名無しさん
2007/03/25(日) 20:23:340729しょう
2007/03/25(日) 20:27:120730デフォルトの名無しさん
2007/03/25(日) 20:28:17せめて何学年まであるか、一学年何クラスかぐらいはわからないと・・・
0731しょう
2007/03/25(日) 20:31:330732デフォルトの名無しさん
2007/03/25(日) 20:34:090733デフォルトの名無しさん
2007/03/25(日) 20:41:40中学生じゃ解けないよ
大学入って掃き出し法を習うまでとっておきなさい
0734デフォルトの名無しさん
2007/03/25(日) 21:04:47>727間違ってる
1)a*2=b+c
だな。
で、2)はb:c=5:6から出した式ですね
あとは連立方程式
0735デフォルトの名無しさん
2007/03/25(日) 21:31:43俺数学選択で早稲田受かったよ
0736デフォルトの名無しさん
2007/03/25(日) 21:38:140737デフォルトの名無しさん
2007/03/25(日) 21:42:19まぁ政経だけどねwショボいけど
0738デフォルトの名無しさん
2007/03/25(日) 22:12:110739デフォルトの名無しさん
2007/03/26(月) 00:05:29これだった
http://pc11.2ch.net/test/read.cgi/gamedev/1170387853/
0740デフォルトの名無しさん
2007/03/26(月) 00:10:22それを言うなら、 a = (a + b + c) * 1/3 だろ。実際のところ、3学年しかないという前提はどこにもないわけだが。
0741デフォルトの名無しさん
2007/03/26(月) 00:26:09楕円の周長計算、
「使えるだけでいい」という姿勢なら
単に数値計算ライブラリからソースコピればいける気が。
0742デフォルトの名無しさん
2007/03/27(火) 01:25:31http://www.youtube.suppa.jp/?uk5LgzILKQU59%+0kBtqQp0RJf30%+8rTFpttjtSU89%+eKWwj9JoAln29%+OSjDO5IUnWJ14%+amKCPh1YrDR14%+KmClRwiAMOy29%+4gIU8IV2hvR14%+CxM26WyV1_d29%+mIBUzAsDgyE132%+@1@_AutomaticPlay
0743デフォルトの名無しさん
2007/03/27(火) 02:01:51★☆YouTubeのCMを連続動画に☆★
//ame.x0.com/main/070327015730.html
http://ame.x0.com/main/070327015730.html
0744デフォルトの名無しさん
2007/03/28(水) 11:26:380745デフォルトの名無しさん
2007/03/28(水) 11:31:260746デフォルトの名無しさん
2007/03/28(水) 11:37:360747デフォルトの名無しさん
2007/03/28(水) 11:43:15>>689のリンク先のページは見た?
0748デフォルトの名無しさん
2007/03/28(水) 13:40:09その離散点使って、数値積分的なことすれば?
Σ(点から直線への距離 × 傾き)
みたいなのを。
0749デフォルトの名無しさん
2007/03/28(水) 19:00:370750デフォルトの名無しさん
2007/03/28(水) 23:08:36>あ、曲線は自由曲線(ドットというか。。。)です。
ドットを数えろ。
0751デフォルトの名無しさん
2007/03/29(木) 08:38:460752デフォルトの名無しさん
2007/03/29(木) 08:57:37案1、 (点と中心との距離-半径)^2 の平均 ÷ 半径^2
案2、 ( 点と中心との距離^2 - 半径^2 )^2 の平均 ÷ 半径^4
あたりだろう。
0753デフォルトの名無しさん
2007/03/29(木) 14:03:24基準円よりも広いところで円をなす点たちと
基準円上のある点にのみ集中する点たちで
どっちが「円からズレている」かは一概には言えない。
0754デフォルトの名無しさん
2007/03/29(木) 14:22:100755デフォルトの名無しさん
2007/03/29(木) 15:48:26残差自乗みたいな簡単な尺度では「(概形が)円からズレてない」
みたいな位相的な構造はとても捕まえられないはずだけども。
0756デフォルトの名無しさん
2007/03/29(木) 16:07:10指定された円とのズレという量があればそれを最小にする半径、中心も求められるわけで・・・
それとも何か素晴らしいアイデアをお持ちで?
0757デフォルトの名無しさん
2007/03/29(木) 16:33:07だから「ズレている」を定義してくれと言ってるんだけどな。
一点に集中してても「ズレてない」とするなら残差自乗で十分だし、
そうでないならより輪郭線抽出などの手法が要るかもしれない。
0758デフォルトの名無しさん
2007/03/29(木) 16:37:32もう少しややこしくするなら、点同士がどれだけ中心からの角度で分散しているかの数値を入れるかい?
0759デフォルトの名無しさん
2007/03/29(木) 16:44:25だとすれば>750で充分だと思うのだけど。
#目的も判らずに数学的な意味を見出そうとしても虚しいばかりだ。
0760デフォルトの名無しさん
2007/03/29(木) 16:52:340761ラジアンの比較で躓いてます
2007/03/29(木) 16:57:43質問の回答として、以下のコードを教えてもらいました。
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; //外積
とスカラー値を返しているのですが、今一つやっている意味が判りません。
http://yuki.to/math_cgi/prybbs.html?log=3&mode=res&no=36824
ぐぐったらこんな掲示板見つけたけど、回答者の答えがイマイチ判りません。
コードを通して、ベクトルを理解したいのですが、誰か教えてもらえませんか?
0762デフォルトの名無しさん
2007/03/29(木) 17:08:28v1,v2が平面上なのでZ成分のみとなります。 だから省略したのでしょう
0763デフォルトの名無しさん
2007/03/29(木) 17:24:431つは、>>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本の貼る平行四辺形の符号付面積を外積と呼ぶのかも。
0764ラジアンの比較で躓いてます
2007/03/29(木) 17:36:19なるほど・・・時計回りだとZ値が下方向
反時計回りだと上方向になる性質を利用して、
Z値だけで判断すればいいという事ですね!!
0765ラジアンの比較で躓いてます
2007/03/29(木) 17:39:31double norm1 = v1.x * v1.x + v1.y * v1.y; // |v1|^2
double norm2 = v2.x * v2.x + v2.y * v2.y; // |v2|^2
0766デフォルトの名無しさん
2007/03/29(木) 17:47:02え、えっと、変数名もnormだし、コメントもnormだし、normを計算しているんじゃないかなあ。
0767デフォルトの名無しさん
2007/03/29(木) 17:53:54(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)
を変形してくと・・・・
って内積計算した方が計算量少ないかもしれないが
まあ、折角 外積計算したからって所じゃないのかな
0768ラジアンの比較で躓いてます
2007/03/29(木) 18:03:57一瞬正規化?とか思っちゃいましたが、「ノルム」でしたか・・・orz
0769デフォルトの名無しさん
2007/03/29(木) 18:19:33って日本語の数学では何だっけ?
0770デフォルトの名無しさん
2007/03/29(木) 18:27:370771デフォルトの名無しさん
2007/03/29(木) 18:37:01参考までに、中国語では「范数」と書くそうな。
0772デフォルトの名無しさん
2007/03/29(木) 18:41:31ノルムはノルムじゃない?
数学用語としてじゃなくて、一般には基準とか模範って訳すけど。
ノルムに似たので(というか、絶対値の一般化)付値ってのがあるけど、
それは英語でも valuation。
0773デフォルトの名無しさん
2007/03/30(金) 08:56:16片方を90度回転したベクトルで内積をとっても、時計回り判定はできる。
90度回転操作を (x, y) → (-y, x) とすると、外積とコメントされた式と同じになる。
好きな方で解釈するといい。
0774デフォルトの名無しさん
2007/03/30(金) 09:20:27片方のベクトルが水平(y成分が0)になるように回転変換し、
もう一方のベクトルのx成分が0なら垂直で y成分の符号を見ればいい
でも、それが内積と外積になっちゃうんだけどね
0775デフォルトの名無しさん
2007/03/30(金) 11:56:04その都度、数学の教養を身に着ける順序でも遅くなくねえ?
問題集をひたすら解くだけの抽象的な数学の本ばかり読んで
いると、生きることの意味がわからなくなってくるよ。
0776デフォルトの名無しさん
2007/03/30(金) 11:57:48半径10センチの球表面の座標(XYZ)をファイルに出力したいと考えております
点の間隔は0.1センチぐらいでいいかな、と
ファイルに落とす部分は、わかっているんですが
座標を算出するアルゴリズムが、さっぱりわからなくて
お分かりになる方、御教授いただけると助かります
0777デフォルトの名無しさん
2007/03/30(金) 12:20:47それが許されればな。
>>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)
解説: http://hp.vector.co.jp/authors/VA030421/fdd06.htm
で角度(th1, th2)でループまわすことも考えられる。
角度で回すのなら、球面三角法
http://ja.wikipedia.org/wiki/%E7%90%83%E9%9D%A2%E4%B8%89%E8%A7%92%E6%B3%95
も見おくべし。
0779デフォルトの名無しさん
2007/03/30(金) 20:41:06できてもさしつかえないならいいけど、もしさしつかえあるなら
ユニバーサルメルカトル図法みたいな雰囲気で局所座標系を
とったりするとよさそうな気が
0780デフォルトの名無しさん
2007/03/31(土) 01:53:290781デフォルトの名無しさん
2007/03/31(土) 02:39:25http://www.cubido.at/Blog/tabid/176/EntryID/86/Default.aspx
ジオデシックドームの頂点求めるんだたら↓これ?
http://www2.tokai.or.jp/tomo-kun/readmeJ.html
0782デフォルトの名無しさん
2007/03/31(土) 04:54:500783デフォルトの名無しさん
2007/04/11(水) 11:07:50ヒントはないでしょうか?
0784デフォルトの名無しさん
2007/04/11(水) 11:16:08単純に円弧で近似したいのがどういう状況か判らないのだけど
曲線の場合は、微分が一致するように接続してゆけばいいのだけど
円弧の場合は2点と半径で求まってしまうので、
どうやっても接続点が尖がってしまう
0785デフォルトの名無しさん
2007/04/12(木) 00:12:150786デフォルトの名無しさん
2007/04/12(木) 08:40:070787デフォルトの名無しさん
2007/04/13(金) 04:39:530788784
2007/04/13(金) 06:06:45接続点で中心の方向が一致すればいいと解けば
接続点が尖らないように出来るか
0789デフォルトの名無しさん
2007/04/13(金) 06:16:25)
(
)
(
こういうこと?
0790デフォルトの名無しさん
2007/04/13(金) 13:54:27直線も半径を限りなく大きな値にすれば可能っぽいね。
適当に曲線から3点抽出して、3点を通る円を求めればいいんじゃね?
必要精度に達していなければ間隔を短くし、足りてれば長くして情報量を落とせばいけそう。
0791デフォルトの名無しさん
2007/04/13(金) 20:32:080792デフォルトの名無しさん
2007/04/14(土) 00:14:380|
―+―――→x
|
|
↓
y
を、
y
↑
|
|
―+―――→x
0 |
に変換する行列教えてください。
0793デフォルトの名無しさん
2007/04/14(土) 00:26:240794デフォルトの名無しさん
2007/04/14(土) 15:36:47>about:configで「network.standard-url.encode-utf8」を「true」にします。
上記の設定で、無事、日本語になりました。
気になるのは、IE7では、『"』⇒『"』でしたが、
http://forum.mozilla.gr.jp/?mode=new&no=0&X='国際化'
http://forum.mozilla.gr.jp/?mode=new&no=0&X="国際化"
FireFox2では、『"』⇒『%22』になっていました、少々オシイです。
http://forum.mozilla.gr.jp/?mode=new&no=0&X='国際化'
http://forum.mozilla.gr.jp/?mode=new&no=0&X=%22国際化%22
『%22』を『"』に戻す作業が残ってしまいます。
共通化としてOpera9の国際化URL設定も分かると良いと思います。
0795デフォルトの名無しさん
2007/04/18(水) 10:03:30ゲームプログラムは数学的要素満載だと思うのですが。
0796デフォルトの名無しさん
2007/04/18(水) 10:46:250797デフォルトの名無しさん
2007/04/18(水) 13:57:03どんな分野であっても、スレタイが読める程度の日本語力は必要。
0798デフォルトの名無しさん
2007/04/18(水) 21:51:26事務web系って何ですか?
0799デフォルトの名無しさん
2007/04/18(水) 22:46:08ネット通販の注文を受け付けて倉庫に発送を指示するとか
0800デフォルトの名無しさん
2007/04/19(木) 02:44:550801デフォルトの名無しさん
2007/04/19(木) 02:48:30出来るだけ長く働ける分野が良いんです
0802デフォルトの名無しさん
2007/04/19(木) 11:29:460803デフォルトの名無しさん
2007/04/19(木) 11:38:45http://pc11.2ch.net/test/read.cgi/tech/1171457185/
を汚した奴だろ?
あちこち汚してゆくような奴が奴隷の待遇の上下を聞いたって嫌われるだけだろうに
0804デフォルトの名無しさん
2007/04/19(木) 23:48:040805デフォルトの名無しさん
2007/04/21(土) 06:19:39あなたの息子を舐め舐めしますよ
0806デフォルトの名無しさん
2007/04/21(土) 09:01:30勉強方法が知りたい。大学受験レベルの数学問題をひたすら解くのがいいのか、
公式を理解しつつ大学レベルの数学を勉強するのがいいのか。
どんな風にすればいいか教えてください。。。。。
0807デフォルトの名無しさん
2007/04/21(土) 09:31:49プログラミング・数学の能力と、役立てようとしている
方面によって勉強すべきものは変わる。
目的に近いところを読んでから、不足してそうな部分を
適宜補うってのが順当なごく普通の勉強の仕方。
どんなのがやりたいか言ってくれれば本などは紹介するよ。
因みに受験レベルの数学は単純計算すら覚束ないなら
仕方が無いけれど、そうでなければやる必要なしと思う。
0808デフォルトの名無しさん
2007/04/21(土) 16:00:13これは例えるなら「スポーツに役立つ練習方法を知りたい」と言ってるようなもの。
スポーツの種目によって練習方法が違うように、プログラミングの分野によって役立つ
数学の分野も違う。
# プログラミングで飯を食っていきたいなら、数学より簿記でもやった方が潰しがきくような気がしないでもない。
0809806
2007/04/22(日) 13:35:06レスサンクス。
プログラミングはCの基本的な文法を理解している程度で、数学に関しては
高校レベルもサパーリな状況です。
簿記の知識って業務系だとそんなに役に立つのか。。。。
それなら一応日商の簿記一級を持ってるんで、この知識を生かして業務系に
進んだほうがいいのかなと思ってるんだが、
「生涯現役でプログラマなんだぜ?」
な漏れとしては組込み系や制御系のほうがいいのかなと考えたり、、、
今現在はどの分野にいこうか迷ってるところです。
どうやらプログラミングに必要な数学は各分野ごとにまちまちで、漏れの
やりたいこと自体もまだ定まってないので、とりあえず高校レベルの数学を
勉強する。そしてその時の勉強法としては大学受験を目指す感じのやり方
じゃなく、基本的な公式を理解する程度でよいってことでFA?
0810デフォルトの名無しさん
2007/04/22(日) 21:45:37全部頼るのは良くない。
自分で決める、たぶんそれが一番大事だ。
0811デフォルトの名無しさん
2007/04/22(日) 22:39:13日本語でOK。
0812デフォルトの名無しさん
2007/04/24(火) 00:45:08毎年毎年現れる上に、人によって状況が違うし、いちいち答えてたらキリが無い
0813全知全能者
2007/04/25(水) 05:10:58原始時代は「筋力」
江戸時代は「家柄」
そして現代は「金」
現代社会では金を持っている人間が強い。
革新的なパラダイムの転換が無い限りこの価値観は変わらない。
0814デフォルトの名無しさん
2007/04/25(水) 15:49:140815デフォルトの名無しさん
2007/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バイト数に拡張します
2) tmp = (short)(tmp * rate); //< 演算結果の下位2バイトを結果として受け取ります
3) tmp = (word)tmp >> 8; //< 無符号型としてシフト()
4) byte result = (byte)(dst + (byte)tmp); //< tmpの下位1バイトのみを足し込みます
これだけ見ると変に複雑に見えますが、
実は計算にはmmxを使っていて4要素まとめて演算します。
そこで、1要素につき2バイトの範囲内で
(src - dst){-255 〜 255} * rate{0 〜 255}
の符号付き乗算をしないといけないため、このようなことになっています。
0816815
2007/04/29(日) 00:24:22例えば、src=0, dst=255, rate=255の場合、
結果としては0が期待されますが、実際に計算してみると
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
とほぼ正しい結果が得られます。
上の式は有名なライブラリで使われている式でもあり、
正しいことはほとんど保証されているのですが、
これがなぜ正しいのか証明できる方はいませんか?
少なくとも自分には理解できないです。
文献でもいいです。よろしくお願いします。
0817デフォルトの名無しさん
2007/04/29(日) 00:48:54どうしてもなにも、定義そのままの計算式じゃん。
> 1) short tmp = (short)src - (short)dst; //< 符号付き2バイト数に拡張します
(src-dst) は、-255〜+255 なので、short に納まる。
> 2) tmp = (short)(tmp * rate); //< 演算結果の下位2バイトを結果として受け取ります
(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 で処理しても結果は同じ。
0818デフォルトの名無しさん
2007/04/29(日) 14:31:47おお!
8ビットまで精度を保証すればうまくいくってことね^^
どうもありがとう
0819デフォルトの名無しさん
2007/04/29(日) 18:19:52一番大事なのは物事を合理的に考える事が出来るか。
0820デフォルトの名無しさん
2007/04/29(日) 19:08:450821デフォルトの名無しさん
2007/04/29(日) 19:46:50ヒントはないでしょうか? と、質問した者です。
点列が1個の円や直線にフィットするかどうかは、最小二乗法などで、
解決できると思います。
したがって、どういう風なグループで円弧や直線にすると、さらに最小になるか、
という問題になるような気がするのですが、そういう問題には、どのような考え方で臨めばよいでしょうか。
0822デフォルトの名無しさん
2007/04/29(日) 20:18:42円は3点で決まるから、3点毎に円弧を描いても、近似と言い張れば近似になる。
0823デフォルトの名無しさん
2007/04/29(日) 20:40:030824デフォルトの名無しさん
2007/04/30(月) 01:30:110825デフォルトの名無しさん
2007/05/02(水) 18:06:14それを、精度をあたえることで、コンピュータに計算で円弧+直線に解かせられないでしょうかね。
0826デフォルトの名無しさん
2007/05/02(水) 18:27:41使う円弧の数を減らしたいという要求があるなら
単純に端から3点でフィッティングして、次の点が誤差の範囲内ならと処理してもいいし
最小2乗円を求めては、誤差の範囲が収まるならとやってもいいとおもう
ただし、誤差だけでは、つなぎ目がガクガクに見えるという事になる
だから、それが嫌なら、点との誤差が幾ら以内で、接続が滑らかであるというような条件を追加しないといけない
0827デフォルトの名無しさん
2007/05/02(水) 19:42:29用途を書いてくれればもう少し具体的なアドバイスが出てくると思われ
0828デフォルトの名無しさん
2007/05/02(水) 20:19:07滑らかじゃなくていいなら既に色々レスついてるからそれ参照。誤差ゼロ。
滑らかである必要があるなら、確かに精度(許容誤差)を与える必要はあるだろう。
滑らかって何とか言うならもっと勉強しなさい。
0829デフォルトの名無しさん
2007/05/04(金) 00:50:190830デフォルトの名無しさん
2007/05/15(火) 23:20:02の高速除算なのですが、
どうしてこれで正しく計算できるのか分かる方いませんか?
あと、割られる数が負数のときでも上手くいくような
つまり
X / D = Q (Qは整数)
のときに
-X / D = -Q
となるような高速除算の方法を知っている方はいますか?
0831デフォルトの名無しさん
2007/05/16(水) 01:13:49定数除算なら、誤差項を1ビット誤差以下に抑えればいいから、こういうことは出来るけど。
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
0832デフォルトの名無しさん
2007/05/16(水) 01:16:49nが負なら最後に1足してくんなまし。
0833デフォルトの名無しさん
2007/05/16(水) 07:36:56Web上にある情報を、どれも正しいと思うな
Q=trunc(X/D +0.5) と計算したい筈なのに
m = 2*n-------------------1)
R = (2^m - 1) / D ---------2)
Q = (R * X + 2^n) >> m ------3) から>>m を2のベキに変更して
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 は1ではない
2^(n-m) は 0.5 ではない
これは単に使いたかった範囲で巧くいっただけだろ
0834デフォルトの名無しさん
2007/05/16(水) 11:48:36ありがとうございます。
一応、プログラムを組んで実験してみたのですが、
違う結果しか得られませんでした。
自分の理解が足りていないのかもしれません。
チェック用プログラムをアップしますので、確認していただければ幸いです。
http://kissho.xii.jp/1/src/1jyou6388.lzh
>>833
整数における割り算なので、やりたいのは四捨五入ではなく切り捨て、
つまりQ=trunc(X/D) です。
あと、使える値の範囲も>>830のページ上に書いてありますよ。
0835デフォルトの名無しさん
2007/05/16(水) 16:28:10831に書いた方法は、
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 の決定がそんなに単純じゃないんだー。
だから定数除算って書いた。
っていうかハッカーのたのしみって本を買えばこういうことが書いてある。
興味あるならオススメ。
0836デフォルトの名無しさん
2007/05/16(水) 17:04:05http://swox.com/~tege/divcnst-pldi94.pdf
> Division by Invariant Integers using Multiplication
>
> 4 Unsigned division
> 5 Signed division, quotient rounded towards 0
> 6 Signed division, quotient rounded towards -∞
0837デフォルトの名無しさん
2007/05/16(水) 18:06:560838デフォルトの名無しさん
2007/05/16(水) 18:13:07実測してみないとなんともだけど割り算の結果を遅延評価するといいかもね。
0839デフォルトの名無しさん
2007/05/16(水) 23:59:10そうなんですか^^
ハッカーの教科書は持っていたのですがこんな話題があるとは知りませんでした。
ちょっと研究してみようと思います。
>>836
そのページを見てみたのですが、正確にやるのは結構大変なんですね^^
よく使う定数の除算のみを最適化して
あとはテーブルなりなんなりで処理することにしました。
>>837
http://pc11.2ch.net/test/read.cgi/tech/1168399966/99-
↑を見て欲しいのですが、x86のdiv命令は
おそらくもっとも効率化がめざましい命令の一種です。
最新の石ならばそのままが一番速いと思うのですが、
昔のCPUを考慮するとなかなかそういうわけにもいきません。
0840デフォルトの名無しさん
2007/05/17(木) 00:08:18この手の高速化は高級言語で記述してもコンパイラが糞なら無意味だし、そもそも可読性が落ちるので歓迎されない。
そうなるとアセンブラでの記述に限定されるわけだが、肝心のターゲットが書かれていないので議論するだけ無駄。
というかそもそもこのスレ向きの話題じゃない。
0841デフォルトの名無しさん
2007/05/17(木) 00:55:370842デフォルトの名無しさん
2007/05/17(木) 01:26:26その理屈で行くとハッカーのたのしみや>>836の内容は
まったく無意味ってことになりますね。
もちろんそんなことはなくて、
ビット演算やコンピュータ用数学(?)にはある程度の普遍性があります。
だから、ターゲットが明確に決まっていなくてもある程度は意味はあります。
(一応、特定アーキテクチャ専用の話題は避けたつもりなのですが、
div命令のレイテンシの話題はたしかにスレ違いですね^^)
>>841
ここでいう「昔のCPU」とは、Pentium4やAthlonXPクラスの
(多分)今一番普及しているCPUのことです。
さすがに考慮しないわけにはいかないです^^
0843デフォルトの名無しさん
2007/05/17(木) 06:19:35アセンブラ系スレなら有意義なんだろうが、それ以外の場所じゃ完全に無意味。
>ビット演算やコンピュータ用数学(?)にはある程度の普遍性があります。
これもない。現にシフト命令が異様に遅いコアが存在する。
5年10年前なら確かにその言い分は通用したが、複雑化したx86系CPUで
除算命令のような数クロックを稼ごうという時にターゲットも指定しないなんてありえないだろ?
0844デフォルトの名無しさん
2007/05/17(木) 07:11:01> アセンブラ系スレなら有意義なんだろうが、それ以外の場所じゃ完全に無意味。
んなことない。定数除算→乗算の最適化をやってくれない環境はあるし、自分
が仕事で使った RISC とかは加減シフトが実効 1 サイクル、乗算が 2 サイク
ルに対して除算命令は 33 サイクル(中で何やってるかわかるな。しかも除算
命令実行中に割り込みが入ると割り込み終了後に再度除算命令を実行し直す)。
> 現にシフト命令が異様に遅いコアが存在する。
後学のために教えてくれ。
0845デフォルトの名無しさん
2007/05/17(木) 07:46:16だから、それはアセンブラ系のスレでやってくれということでしょ
0846デフォルトの名無しさん
2007/05/17(木) 08:04:05速度だけの問題なら他でやって欲しい所。
0847デフォルトの名無しさん
2007/05/17(木) 09:09:48自分がもともと質問していたことは、除算→乗算+シフトの変換方法です。
(このスレのおかげで定数除算の最適化という分野があることを知りましたし、
特定の除算の最適化もできそうな感じです。
その点は本当に感謝しています)
確かに石によって演算速度の違いがあることは事実なのですが、
そういったこともふまえて、対象となるCPUをつらつら挙げて
CPU別の最適化方法をここで聞いた方がよかったのですか?
もちろん、自分的には大歓迎なのですが^^
>>846
収穫もあったことですし、自分はそろそろ去りたいと思います^^
0848デフォルトの名無しさん
2007/05/17(木) 10:32:12最適化スレもあればアセンブラスレもあるのに?
0849デフォルトの名無しさん
2007/05/17(木) 10:51:11そういう命令セットの定義とかも数学の範疇じゃないのか?
0850デフォルトの名無しさん
2007/05/17(木) 11:24:26>CPU別の最適化方法をここで聞いた方がよかったのですか?
は>>843の
>除算命令のような数クロックを稼ごうという時にターゲットも指定しないなんてありえないだろ?
に対するレスです。
当然、自分の考えは>>848さんと同じで「それはスレ違いでしょ?」です。
>>849
命令セットの定義が数学の範疇だというのははじめて聞きました。
そうなのかもしれないのですが、
自分は特定アーキテクチャの最適化の話はここでするべきでは無いと考えたのです。
0851デフォルトの名無しさん
2007/05/17(木) 11:45:33R>0 Xは正負のint 範囲で
Q=trunc(X/D) を
Q= (int)(((long)X*(long)R + k) >> (long)m )
で計算したいが 適切な R と m を求むって事だろ?
0852デフォルトの名無しさん
2007/05/17(木) 11:45:45求める問題ってのは、特に並列計算の分野で見ることがある研究だね。
実際のCPUくらいのコストモデルでやるのはかなり大変だと思うけれど
0853デフォルトの名無しさん
2007/05/17(木) 11:46:45少なくともここはアセンブラの話題振るスレじゃねえことくらいは理解してんだろ?
0854デフォルトの名無しさん
2007/05/17(木) 21:17:53ただ、その実装手段が大抵の場合アセンブラしかないだけ。
CPUの除算が遅いとか言い出したヤツが最初にスレ違い。
0855デフォルトの名無しさん
2007/05/17(木) 23:01:21除算は多倍長文字列でやった事あるけど、
別と比べて極端に難しかった。
0856デフォルトの名無しさん
2007/05/18(金) 01:50:28ここがム板である限り。
0857デフォルトの名無しさん
2007/05/19(土) 05:32:05除算は四則のうち、唯一「トライ&エラー」的要素が必要になってくるからな。
「除法の定理」という、それほど単純ではない性質に則った結果にしなければ
いけないという縛りがその原因だろう。
0858デフォルトの名無しさん
2007/05/19(土) 15:26:240859デフォルトの名無しさん
2007/05/21(月) 08:21:460860デフォルトの名無しさん
2007/05/21(月) 08:28:330861デフォルトの名無しさん
2007/05/21(月) 08:37:57> a/b
>ニュートン法でcを更新すれば良い結果が得られます
> ニュートン法x-f(x)/f'(x)を適用するのに
> c*b-2^n=0 f(x)=b*x-2^n f'(x)=b とやったのではダメです
> b-2^n/c=0 f(x)=b-2^n/x f'(x)= 2^n/x^2
> x:=x-b*x^2/2^n+x=(x+x-b*x^2/2^n)=x+(x*(1-b*x/2^n));
> c:=c+(c*(2^n-b*c))/2^n と更新します
0862デフォルトの名無しさん
2007/05/21(月) 08:44:08除算をわざわざ自作しなくても簡単に判ることですが?
それとも、r = p / d を r * d = p に変形することもできないんですか?
0863デフォルトの名無しさん
2007/05/21(月) 09:10:47それだと
0 / 0 が定義出来てしまいますが?
0864デフォルトの名無しさん
2007/05/21(月) 10:23:45あんた馬鹿? dが0のときにpが0になるrは、一意には決まらないんですが。
つまり、0 / 0は数学的には解があるってこと。1 / 0とは事情が違う。
0865デフォルトの名無しさん
2007/05/21(月) 11:06:44全くだなw
0866デフォルトの名無しさん
2007/05/21(月) 13:17:07アルキメデスの原理を学んでいれば理解出来る筈だ。
実際にグラスに氷を入れて、水をすりきれ一杯にして、水がこぼれない実験をするまでは判らないのだろう
0で割るとどうなるかも、そりゃ数式を弄くれば予想は出来るが
それと、実際に計算させてみての実感とは違うものなのだろう
0867デフォルトの名無しさん
2007/05/21(月) 14:43:02地面の上に乗ってる氷が溶けると海面上がるよな。
北極って、大陸ではないけど、
小さな島状の部分とか浅瀬もなく、完全に氷が海に浮いてる状態?
0868デフォルトの名無しさん
2007/05/21(月) 14:50:03>>867
どこまでを北極というかによって変わるが、基本的に北極海には問題になるほどの陸地はない。
勿論北米北欧ロシアの海岸線を含むわけだが、面積比で言えば問題にならない。
0869デフォルトの名無しさん
2007/05/21(月) 15:19:40映画でアルゴアが言ってた
>>866
グラスに氷とかのすぐできるような実験もしないで
「本に書いてあるからそうなんです」って胸を張るような
中世の暗黒時代さながらの学問態度は取りたくないものだな
馬の歯は何本あるんだっけ
0870デフォルトの名無しさん
2007/05/21(月) 15:47:10でも、氷河期が終わるんなら、氷床も無くなって当然では?
0871デフォルトの名無しさん
2007/05/21(月) 16:23:24>アルキメデスの原理を学んでいれば
うみのみずはしょっぱいです。でも、ほっきょくのこおりはしょっぱくないです。
0872デフォルトの名無しさん
2007/05/21(月) 16:41:13その上に降り積もった部分は違うだろうけど
0873デフォルトの名無しさん
2007/05/21(月) 16:43:46面白い考え方だな
0874デフォルトの名無しさん
2007/05/21(月) 17:09:45塩分濃度が35ppt で 温度15℃の時の比重が1.0263
つまり、塩が溶けると 1.035/1.0263 = 1.0085倍に膨張するといいたいわけ?
0875デフォルトの名無しさん
2007/05/21(月) 17:15:29つ[浮力]
0876デフォルトの名無しさん
2007/05/21(月) 17:15:36>海の水が凍ったものは 塩を含む
もしそうなら、半導体価格は今より何倍も高いはずだ。
0877デフォルトの名無しさん
2007/05/21(月) 17:16:370878デフォルトの名無しさん
2007/05/21(月) 17:17:150879デフォルトの名無しさん
2007/05/21(月) 17:18:24そんなことは聞いてねぇよ、北極の氷はプログラムと関係あんのかって聞いてんだよ
0880デフォルトの名無しさん
2007/05/21(月) 17:20:10あると思うの? んなことマジで聞いてるの?
0881デフォルトの名無しさん
2007/05/21(月) 17:23:48雑談でスレ埋めるな、って嫌味も通じない馬鹿か…
0882デフォルトの名無しさん
2007/05/21(月) 17:31:330883デフォルトの名無しさん
2007/05/21(月) 17:32:49さて
http://blog.livedoor.jp/yuge15978/archives/8126608.html
のように、真水が凍った氷が溶けて水面が上がのが正しそうに見えるけど
氷が溶けると、塩水の濃度は下がるのだから、濃度が下がる事により容積は下がるので
その効果が打ち消す可能性もあり、やっぱり実験しないと判らない
0884デフォルトの名無しさん
2007/05/21(月) 18:35:54海氷の塩分濃度は10〜12gだそうだ。 海水の1/3程度だとはいえ、まだ塩分はそれなりにあるみたい
0885デフォルトの名無しさん
2007/05/21(月) 18:42:39たまに、海水は真水より比重が大きいのだから、
水銀の上で氷が溶ければ水面が上がるように、水面が上がるという人がいる。
でも、これは海水の主成分が水である事を忘れている。
氷が溶ければ真水ではなく、塩水になり、塩水は薄まる。
塩分濃度と比重は変化が小さい時には比例するので、
トータル水面は変わらないのさ。
極端に言えば、海水中の塩分と、水を別けて計算すればいい。
0886デフォルトの名無しさん
2007/05/21(月) 20:59:20氷の比重0.917
水の比重1.000
海水の比重1.024
海水上に出る氷の容積は、1-0.917/1.024 = 0.104
氷が水になると1-1/1.024 = 0.023
一旦は水面は氷*2.3%だけ持ち上がるよね? その後塩水と混じると水面が下がるわけ?
0887デフォルトの名無しさん
2007/05/21(月) 22:48:440888デフォルトの名無しさん
2007/05/21(月) 22:52:390889デフォルトの名無しさん
2007/05/21(月) 22:54:080890デフォルトの名無しさん
2007/05/21(月) 22:56:040891デフォルトの名無しさん
2007/05/21(月) 22:59:280892デフォルトの名無しさん
2007/05/21(月) 23:06:040893デフォルトの名無しさん
2007/05/21(月) 23:16:00ヘンリーラウールの法則によると凍った海水には塩が含まれる
0894デフォルトの名無しさん
2007/05/22(火) 00:05:390895デフォルトの名無しさん
2007/05/22(火) 00:08:060896デフォルトの名無しさん
2007/05/22(火) 00:17:380897デフォルトの名無しさん
2007/05/22(火) 00:35:250898デフォルトの名無しさん
2007/05/22(火) 00:57:150899デフォルトの名無しさん
2007/05/22(火) 02:43:23はずかしいボケをしたお友達かな?
0900デフォルトの名無しさん
2007/05/22(火) 08:17:300901デフォルトの名無しさん
2007/05/23(水) 08:18:110902デフォルトの名無しさん
2007/06/02(土) 20:26:08/\ / /
\ \ / /______ _____ ___ ___ ___
.\// _______/ /__ __/ [][] _| |_| |__ _| |_
/\ /__/ ______ / / | _ | |_ レ'~ ̄|
\ \ / / | |___  ̄| | / / / /| |
.\/ / / ̄/ / ̄/ / \__| | |  ̄ /_ / | |_
/ ̄ ̄  ̄  ̄  ̄/ |_| |__| \/
/\ ̄/ / ̄/ / ̄/ / ̄
/ / /  ̄  ̄ ̄  ̄/
\/  ̄ ̄ ̄ ̄ ̄/ / ̄ 2007 Summer @ Shou-nan
/ ̄ / http://www.2ch.net/sea/
 ̄ ̄
0903デフォルトの名無しさん
2007/06/04(月) 19:59:530904デフォルトの名無しさん
2007/06/04(月) 21:16:030905デフォルトの名無しさん
2007/06/07(木) 14:35:14これは 引分、勝ち負け を0,1,2 の3値で表現して、14試合分を ランダム発券するというものです
売り上げの40%が一等原資、10%を2等以下で別けます。
300円でキャリーが無ければ上限3億円 キャリーがあれば6億円が1等賞金の上限です
でいつ買えば得なのか、還元率を求めたいと思いました。
2等以下はまず平均に出るとして 300円中 30円 とし、
1等原資は 売り上げX口として300*X*0.4+CY (CYはキャリー金額)
1等が1口以上出る確率は(1-(1-1/3^14)^X なので
30+(300*X*0.4 +CY)*(1-(1-1/3^14)^X)/X
と求まります。
が、上限が6億なので、155円以上の結果は正しくないと思えます。
どういうプログラムを作ればいいでしょう?
0906デフォルトの名無しさん
2007/06/07(木) 15:14:240907デフォルトの名無しさん
2007/06/07(木) 15:56:57P = 1/3^14 として
X 中 1個が出る確率* 当選金
X 中 2個が出る確率* 当選金
X 中 3個が出る確率* 当選金
を2項分布で確率を求めて累積したらいいんじゃないの?
もしかして、nCr とかオーバーフローするのかな?
ならポアソン分布で近似すればいいって、これも オーバーフローするか・・・・困ったね
0908デフォルトの名無しさん
2007/06/07(木) 16:06:19やっぱりポアソン分布の方がオーバフローし易いみたいだな。
n個中k個の確率は2項分布なら
n!/(k! * (n-k)!)*p^k*(1-p)^(n-k) で、階乗をそのまま計算したらオーバフローするから
kの大きい所まで計算する必要はないから、 for 文で計算すればオーバーフローしないんじゃないかな
と無責任に言ってみる
0909デフォルトの名無しさん
2007/06/07(木) 17:26:18いつ買っても損をする。
なぜなら賭博とはそういうふうに出来ているから。
totoがダメだったのは、賞金頭割りだから順当な結果の回は1等数千円とかザラ、荒れるとビックリするほど当らない。
BIGはどうなるかわからんけど、totoで客離れちゃったからキャリー期待するのは難しいんじゃないかなぁ。
計算すんのは楽しいけどねw
0910905
2007/06/07(木) 18:38:35ありがとう。その式で計算させました。 JavaScriptですが
function combination(n, k)
{
var i,c=1;
for(i=1;i<=k;i++){ c=c*(n+1-i)/i;}
return c;
}
function P(n,k)
{
var p0=1.0/Math.pow(3, 14);
if(k==0) return Math.pow( 1-p0 ,n );
return combination(n, k) * Math.pow(p0,k) * Math.pow( 1-p0 ,n-k );
}
function kangen(CY , X)
{
var i;
var gensi=X*300*0.4+CY;
var umax=6E8;
if(CY==0) umax=3E8;//キャリーが無い時は上限3億円
var sum=0;
for(i=1;i<20;i++) { //iが1当の数
var a=Math.min(gensi/i,umax); //1等の配当は原資を当選数で割った金と上限の小さい方
p=P(X,i);
sum += i*a*p;//当選金総額x確率を掛算
if(p<1E-6)break;//計算誤差が6桁以下になれば終わり
}
return 30+sum/X;
}
0911905
2007/06/07(木) 18:43:51キャリー 100万口 200万口 300万口 400万口 500万口 600万口 700万口 800万口 900万口
0億円 53円 71円 79円 82円 86円 86円 88円 89円 89円
1億円 72円 88円 101円 112円 113円 115円 117円 120円 123円
2億円 90円 105円 117円 119円 119円 121円 123円 126円 127円
3億円 109円 122円 126円 124円 125円 126円 129円 130円 130円
4億円 128円 134円 130円 129円 130円 132円 134円 133円 133円
5億円 144円 137円 135円 134円 136円 138円 136円 136円 136円
6億円 145円 140円 139円 140円 142円 140円 139円 139円 139円
7億円 147円 144円 143円 145円 143円 142円 142円 142円 143円
8億円 149円 147円 148円 147円 145円 144円 144円 145円 144円
9億円 151円 150円 150円 148円 147円 146円 147円 147円 146円
10億円 153円 153円 151円 149円 149円 149円 149円 148円 147円
11億円 155円 153円 152円 151円 150円 151円 150円 149円 148円
12億円 155円 154円 152円 152円 152円 151円 150円 150円 150円
キャリーが4億でも還元率は50%の 150円ありません。
そして、いくらキャリーが増えてもやっぱり155円から還元率は増えません。
9億くらいキャリーが溜まって、やっと50%の還元率です。
計算があってるとすればヒドイ設計のクジです。
0912905
2007/06/07(木) 18:50:06<HTML><HEAD>
</HEAD>
<BODY>
/////////////// この間に >>910 のコードを入れる。///////////
var x,y;
document.write('<TABLE BORDER><TR><TD>キャリー');
for(x=100;x<1000;x+=100)document.write('<TD>'+x+'万口' );
for(y=0;y<20;y=y+1)
{
document.write('<TR>');
document.write('<TD ALIGN=RIGHT>'+y+'億円</TD>');
for(x=100;x<1000;x=x+100){
document.write('<TD ALIGN=RIGHT>'+Math.round(kangen(y*1E8,x*1E4))+'</TD>');
}
document.write('</TR>');
}
</SCRIPT>
</BODY>
</HTML>
0913デフォルトの名無しさん
2007/06/07(木) 18:50:42宝くじって大体そんなもんだよ。
もっとひどいこともある。
競馬は還元率75%らしいけどね。
0914905
2007/06/07(木) 19:02:03宝くじやロトよりも悪いです。
これでは、キャリーが8億超えないと買う気になれず、
だから売れないのでしょう。
14試合ではなく13試合にしておけば、1等の平均賞金は2億程度になりますが
5億キャリーすれば還元率が100%を越える事になり、もっと面白いクジになったでしょうね。
0915デフォルトの名無しさん
2007/06/07(木) 20:30:170916デフォルトの名無しさん
2007/06/07(木) 20:47:34微塵でも還元率が100%超えちゃう可能性があるものはギャンブルとしてはちょっとどうかと。
0917905
2007/06/07(木) 21:08:52還元率が100%超えてる間は、当然みんなブッコムでしょ?
だから、心配は要らない。 というか売れてみんなハッピーだと思うけどな
13試合にすると平均は1億9千万だけど、この効果で売り上げが上がるだろうから
1等6億は逆に出易かったろう。
14試合にしたせいで、>>911のように全然面白くないクジになってしまったかと
だってキャリーが11億無いと、買うのは馬鹿ってクジだもの。
0918デフォルトの名無しさん
2007/06/07(木) 21:16:440919デフォルトの名無しさん
2007/06/07(木) 21:21:15そこでは買い手が馬鹿であることを仮定するからねえ。
0920デフォルトの名無しさん
2007/06/07(木) 21:28:57そもそもマジでもうけたい人は宝くじなんてしないさ。
0921905
2007/06/07(木) 21:37:32馬鹿でも欲は持ってるから、還元率が宝くじの45%より低いと買わないだろう。
回せる財布は限られてるからね。
で、買われないからキャリーが溜まらないの悪循環。
だから去年末からはじめて、1等がまだ2回しか出ていない惨状。
せめて、完全ランダムじゃなく、全部の組み合わせ 3^14 をシャッフルして
発行する方法にしておけば 480万口に1つは確実に1等が出るから
ずっと早くリミットの155円に達するのに、馬鹿だよね。
0922デフォルトの名無しさん
2007/06/07(木) 21:53:46具体的になんか変えたいならこんなとこに書いてないで適切なところに行ってね。
0923デフォルトの名無しさん
2007/06/07(木) 22:18:49まず、胴元が儲けを確保してから残りをどう分配するかっていう事を考えてる訳だから。
計算とかするまでもなく、ね。
BIGだろうがジャンボだろうが確率なんてあって無きがごとしなので、どっちが有利とかないな。
一枚も買わなければ確率0だから、とにかく一口買って、リターンがデカイ方に張るのが賢いやり方。
ジャンボを10枚セットで買う奴とかはどうかしてると思う。
もっとも、一口買う時点でノータリンの貧乏人と言われればそれまでだが。
0924905
2007/06/07(木) 22:37:02資産200万以下で、貯蓄中なら、収入の2%程度までを 還元率が悪くても
億くらいのリターンのある宝くじに入れるのは悪くないと思うんだ。
コツコツ貯めるだけだと老後の資金でせいっぱい、生きてるだけでせいっぱいだもんね。
でも、そういう目的だと、宝くじの3等とか4等っていらないんだよね。
トトBIGは1等に40%集中だから、いいかと思ったけど、やっぱりダメだね
0925デフォルトの名無しさん
2007/06/07(木) 22:42:1014試合中引き分けが5試合以上とか、強いチームが負けまくってるとか、
つまりそういう、「見た瞬間にはずれがわかるくじ」が氾濫しまくってることが原因だと思う。
0926905
2007/06/07(木) 22:47:07まあ、最近は役所仕事は信じられないからアレだけど
それより >>911 の結果を見てしまうと
キャリーが9億以下の時に買う気になれない方が大きいと思うのだが?
300円で5円くらいは損しても許せるが、それ以上は俺は許せんのだけど
0927デフォルトの名無しさん
2007/06/07(木) 22:56:57株でも貯蓄でもいいけど、還元率が100%を超える金融商品なんてそれこそたくさんある。
じゃあ何でわざわざ宝くじを買うかっていうと、
やっぱり一時であっても夢を見れるからだと思うんだよね。
「もし3億当たったら何買おう?」とか考えるだけでも楽しいしwww
ただ、totoBIGの場合はそれが難しいみたいだね。
そもそも1/3の確率で引き分けが出るなら14試合中4試合が引き分けになるけど、
その時点でほとんどはずれだって分かる。
0928905
2007/06/07(木) 23:10:56俺はこれから貯めるわけだからさ。
でもさ、株も確実に儲かるわけじゃないでしょ?
というか小金で株を景気の良い今買って景気が悪くなって金必要になって売ったら損するわけだしさ、
貯金は、サトウキビ燃料のせいで物価上がってインフレになりそうな今は目減りするだけのようだし
といって、土地は買えるだけ金ないしで、
まあ、ほんと、勉強でもして自分に投資でもするしかないけど
投資しても使い捨てられるだけのような感じだし。
0929デフォルトの名無しさん
2007/06/07(木) 23:44:120930デフォルトの名無しさん
2007/06/08(金) 01:05:59いや、貯める目的ならなら余計に宝くじはやめとけ。
貯蓄なら確実だし、株だって還元率50%は無いから安心しろw
0931デフォルトの名無しさん
2007/06/08(金) 03:02:59熊が人間に問いかけている
0932デフォルトの名無しさん
2007/06/08(金) 03:17:55プログラムが当たり前に書けるレベルの知能があるなら、完全ランダムの賭博はしないだろ。
そのへんが全く理解出来てない時点で905の程度が知れる。
0933デフォルトの名無しさん
2007/06/08(金) 12:49:48それで騙されて買う馬鹿がいるんだから
胴元にしてみればすばらしい設計ではないか
0934デフォルトの名無しさん
2007/06/11(月) 09:04:43(最大値-最小値)/(32-1)で32段階になると思うのですが、
これ、あってますか・・・?
0935デフォルトの名無しさん
2007/06/11(月) 09:17:250936デフォルトの名無しさん
2007/06/11(月) 09:22:460937デフォルトの名無しさん
2007/06/11(月) 09:37:510938デフォルトの名無しさん
2007/06/11(月) 14:42:470939デフォルトの名無しさん
2007/06/11(月) 15:28:54入力をXとして
trunc( (X-最小値) / (最大値-最小値)* (32-1)) とやるか
round( (X-最小値) / (最大値-最小値)* 32 ) とやるか という話?
0940デフォルトの名無しさん
2007/06/11(月) 15:36:420941デフォルトの名無しさん
2007/06/11(月) 17:15:56普通に考えると、こうなる。
(1) Y = trunc(X * 32 / 2.0) (X が 2.0 を含まない場合)
(2) Z = X * 31 / 2.0, Y = trunc(Z) (X が 2.0 を含む場合)
(3) Z = X * 31 / 2.0, Y = trunc(Z + 0.5) ((2)だとXがきっかり2.0のときしかYが31にならないので少し改良)
X が 2.0 を含まない場合は全く問題ない。
しかしそうで無い場合は、(3)を使っても
0.0 <= Z < 0.5 → Y = 0
0.5 <= Z < 1.5 → Y = 1
30.5 <= Z <= 31.0 → Y = 31 (Zの最大値は31のため)
つまりY=0, 31の場合だけ、とれるZの範囲が半分になってしまう。
この場合の量子化を正確に行うのは少し難しい。
改良案として、
Z = X * (32 - delta) / 2.0 (deltaは適度に小さい数。Xが整数値ならdelta = 1)
Y = trunc(Z)
としてYを計算する方法がある。
つまり、0.0 <= X <= 2.0 の値を 0.0 <= Z <= 31.99999... に写像してしまえばいい。
もう一つは、X = 2.0 や Y = 32 の値だけを特別扱いする。
つまり、
Y = trunc(X * 32 / 2.0)
if (Y == 32) Y = 31 (Y -= (Y >> 5))
としてしまう方法もある。
0942デフォルトの名無しさん
2007/06/13(水) 17:57:43[0, 32)
なんだから
量子化後は
[0, 31]
でいいよ
0943デフォルトの名無しさん
2007/06/13(水) 18:23:44それ[0.0 〜 2.0]と書かない
0944デフォルトの名無しさん
2007/06/14(木) 01:06:42って言ってるんじゃないの。
0945デフォルトの名無しさん
2007/06/14(木) 02:46:100946デフォルトの名無しさん
2007/06/24(日) 20:14:41http://ja.wikipedia.org/wiki/%E5%8C%BA%E9%96%93_%28%E6%95%B0%E5%AD%A6%29
0947デフォルトの名無しさん
2007/08/08(水) 07:42:250948デフォルトの名無しさん
2007/08/12(日) 03:44:24どうしてもわからないので質問したいのですが
DirectXで3Dのゲームを組んでるのですが、ある点からポリゴン面までの
最短距離をはかりたいのですが、数学的知識があまりなくうまくいきません
なんとか距離をはかる方法はないでしょうか?
ポリゴン面はとりあえず3点で構成されているただの一枚のポリゴンです。
0949デフォルトの名無しさん
2007/08/12(日) 07:01:48ある点=点P、ポリゴン面=�僊BCとする
次の3通りの可能性がある
1.面上の何処かが一番近い
>面の法線ベクトルで点Pを通る直線がポリゴンと交わる時、交わる点が一番近い
2.3辺の線上の何処かが一番近い
>1を満たさず、各辺と点Pからなる三角形の内角が90度を越えない場合、点Pと該当する辺に垂直に交わる点が一番近い
3.3頂点の何処かが一番近い
>1と2を満たさない場合、3頂点のうち一番距離が近い頂点が一番近い
0950デフォルトの名無しさん
2007/08/12(日) 07:12:19とかのほうがかえって早いとかは?
0951デフォルトの名無しさん
2007/08/12(日) 19:14:06両方試させて頂きます。
ただ、今のところ面の法線をだして、面上の交点をだそうと試みてるのですが、
なかなか上手くいきません・・
とにかくもう少しがんばってみます
0952デフォルトの名無しさん
2007/08/13(月) 05:37:180953デフォルトの名無しさん
2007/09/11(火) 11:29:540954デフォルトの名無しさん
2007/09/12(水) 11:37:52商業系の高校行って、大学進学したもんで普通科高校でやる数学I-III、A-Cを
まったくやってないんだけど、中学3年の証明でつまずいてたヤツが、新課程の
数学I-III、A-Cを独学するに辺り、最適な参考書、数学攻略本は何があるでしょうか。
0955デフォルトの名無しさん
2007/09/12(水) 11:55:15面白そうなの1冊、単なる小説1冊、ちょっと難しそうなの1冊として
毎週続くようにね。
0956デフォルトの名無しさん
2007/09/13(木) 01:03:060957デフォルトの名無しさん
2007/09/13(木) 01:27:06俺もちょうど今日会社で数セミの紹介記事見て面白そうって思ったところだ
アマゾンで買ってみるか
0958デフォルトの名無しさん
2007/09/29(土) 05:14:160959デフォルトの名無しさん
2007/10/10(水) 18:28:53離散データを微分したいのですが、
C/C++系で微分ライブラリって何がありますか?
オイラーみたいな?
0960959
2007/10/10(水) 18:35:40単純差分が1次、Eulerが2次、Taylorが3次、Runge-Kutteが4次、であってまつか?
0961デフォルトの名無しさん
2007/10/10(水) 18:58:15数値微分じゃなくて、微分方程式の解、要するに求積の方のアルゴリズムじゃない?
微分の方は、離散点のデータしかないとなると、
差分近似するか、適当な関数で補完して微分するかくらいしか思いつかない。
0962デフォルトの名無しさん
2007/10/10(水) 22:05:41http://ja.wikipedia.org/wiki/PID%E5%88%B6%E5%BE%A1
0963959
2007/10/11(木) 08:36:050964デフォルトの名無しさん
2007/10/11(木) 09:54:37fs/4以上含まないとしても、荒い近似しか出来ない。
0965デフォルトの名無しさん
2007/10/11(木) 14:04:49959は近似をしたいのか、差分を出したいのか、なんなのか。
0966959
2007/10/11(木) 14:25:082階微分って、もしかして、2回微分するってことですか?
それで微分とは等間隔の場合、隣の値と差分取るだけで良いのですか?
0967デフォルトの名無しさん
2007/10/11(木) 14:30:21>2階微分って、もしかして、2回微分するってことですか?
そう。
何回微分したかは rank、階っていう。
>それで微分とは等間隔の場合、隣の値と差分取るだけで良いのですか?
差分とるだけでいいというか、差分で近似するしか手がない。
サンプリング定理を満たすようなデータ列なら、
補完っていうかローパスフィルタで連続信号にできるから
それ使って数値微分(これも差分近似するだけだけど)で精度はあげれるけど。
0968959
2007/10/11(木) 14:32:040969デフォルトの名無しさん
2007/10/11(木) 19:51:430970デフォルトの名無しさん
2007/12/02(日) 00:53:440971デフォルトの名無しさん
2007/12/02(日) 00:57:380972デフォルトの名無しさん
2007/12/02(日) 01:06:180973デフォルトの名無しさん
2007/12/02(日) 01:22:290974デフォルトの名無しさん
2007/12/02(日) 01:23:00線形代数しか思い浮かばんぜ
0975デフォルトの名無しさん
2007/12/02(日) 01:25:310976デフォルトの名無しさん
2007/12/02(日) 01:26:26運動方程式を解くための微分方程式についての知識
0977デフォルトの名無しさん
2007/12/02(日) 01:41:08立てようと思ったが、関連スレが存在するのかすらわかんね。
0978デフォルトの名無しさん
2007/12/02(日) 01:43:30プログラミングの為の数学と算数 vol.2
プログラムに必要な数学、算数に関する話題について
語りましょう。TIPS/Q&Aスレです。
関連板:
http://science6.2ch.net/math/
これでいいんじゃ?
0979デフォルトの名無しさん
2007/12/03(月) 19:08:46それとも自分でつくっていますか?ガンマ関数とかいろいろ
ライブラリ作るの面倒くさいよ・・・
購入したほうがいいのかな・・・
0980デフォルトの名無しさん
2007/12/03(月) 19:33:100981デフォルトの名無しさん
2007/12/05(水) 00:37:190982デフォルトの名無しさん
2007/12/06(木) 00:29:060983デフォルトの名無しさん
2007/12/06(木) 01:03:120984デフォルトの名無しさん
2007/12/06(木) 18:10:28↑上位1ビットの値 というのは一番左端の値であり、1ですか?
0985デフォルトの名無しさん
2007/12/06(木) 18:20:510986デフォルトの名無しさん
2007/12/06(木) 21:48:310987デフォルトの名無しさん
2007/12/06(木) 23:51:460988デフォルトの名無しさん
2007/12/06(木) 23:57:44左をみろ!!→→→→→→→→→→→→→→→→→→→→→→→→→→
0989デフォルトの名無しさん
2007/12/07(金) 00:37:15十進数でも、
1234なら、1(千の位)が一番上の桁でしょうや。
0990デフォルトの名無しさん
2007/12/07(金) 04:23:48LSB
何の略だっけ
0991デフォルトの名無しさん
2007/12/07(金) 09:56:310992デフォルトの名無しさん
2007/12/07(金) 22:08:450993デフォルトの名無しさん
2007/12/07(金) 22:28:550994デフォルトの名無しさん
2007/12/07(金) 22:37:49ドイツ語で自分は孤独だという嘆き。
0995デフォルトの名無しさん
2007/12/08(土) 06:07:080996デフォルトの名無しさん
2007/12/08(土) 06:07:380997デフォルトの名無しさん
2007/12/08(土) 06:08:090998デフォルトの名無しさん
2007/12/08(土) 06:08:390999デフォルトの名無しさん
2007/12/08(土) 06:09:101000小倉優子 ◆YUKOH0W58Q
2007/12/08(土) 06:23:2010011001
Over 1000Threadもう書けないので、新しいスレッドを立ててくださいです。。。
レス数が1000を超えています。これ以上書き込みはできません。