プログラミングの為の数学と算数 vol.2
■ このスレッドは過去ログ倉庫に格納されています
0001デフォルトの名無しさん
NGNG語りましょう。TIPS/Q&Aスレです。
0339デフォルトの名無しさん
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 ヽ_. ' `'ﾞ`'‐'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‐‐‐'--､_..... ｌ_,.．-'''""'- "
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大学（慶応、上智、東京、明治、立教、早稲田）のどっかに入らないと
イカンわけだろ？
いやー回答者の資格に制限がある問題は難しいわ、ホント難しい、
全然答えられそうにないよｗｗｗ
最近のセンター試験ってこんな答えるのが難しい問題出すの？（・∀・）
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近似的に解くべきだなｗ
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);
だろが！
■ このスレッドは過去ログ倉庫に格納されています