プログラミングのお題スレです。
【出題と回答例】
1 名前：デフォルトの名無しさん
お題：お題本文
2 名前：デフォルトの名無しさん
>>1 使用言語
回答本文
【ソースコードが長くなったら】 (オンラインでコードを実行できる)
https://ideone.com/
http://codepad.org/
http://compileonline.com/
http://rextester.com/runcode
https://runnable.com/
https://code.hackerearth.com/
http://melpon.org/wandbox
https://paiza.io/
宿題は宿題スレがあるのでそちらへ。
前スレ
プログラミングのお題スレ Part10
https://mevius.5ch.net/test/read.cgi/tech/1514772904/
プログラミングのお題スレ Part11
https://mevius.5ch.net/test/read.cgi/tech/1524570314/
プログラミングのお題スレ Part12
レス数が1000を超えています。これ以上書き込みはできません。
2018/09/28(金) 10:09:07.13ID:phwOkayR
952デフォルトの名無しさん
2019/01/31(木) 22:25:01.31ID:xJsSt9Re Python では、機械学習などの周辺ライブラリ一式が揃っているから、使われる
Ruby では、ベクトル演算のNArray があって、
処理速度はOctave にも匹敵するけど、一式揃っていない
北大の湊教授が作った、ZDD はあるけど
Ruby では、ベクトル演算のNArray があって、
処理速度はOctave にも匹敵するけど、一式揃っていない
北大の湊教授が作った、ZDD はあるけど
953デフォルトの名無しさん
2019/01/31(木) 22:46:32.84ID:2IY59Hh/ scipy/scipy/special/cephes at master ・ scipy/scipy
ttps://github.com/scipy/scipy/tree/master/scipy/special/cephes
この辺り見てみるとフォートランとCがそのまま使われてる
scipy/_comb.pyx at master ・ scipy/scipy
ttps://github.com/scipy/scipy/blob/master/scipy/special/_comb.pyx
このあたりをそのまま書ければいい
ttps://github.com/scipy/scipy/tree/master/scipy/special/cephes
この辺り見てみるとフォートランとCがそのまま使われてる
scipy/_comb.pyx at master ・ scipy/scipy
ttps://github.com/scipy/scipy/blob/master/scipy/special/_comb.pyx
このあたりをそのまま書ければいい
954デフォルトの名無しさん
2019/02/01(金) 04:58:43.05ID:IMBPiIxm お題
平方数かどうか判定する
平方数かどうか判定する
955デフォルトの名無しさん
2019/02/01(金) 10:08:25.26ID:VdCPb4pG C++
#include <iostream>
#include <limits>
template <typename T>
bool is_safely_multiplicable(T a, T b) {
if (b == 0) return true;
return a <= (std::numeric_limits<T>::max() / b);
}
uint64_t sqrt_int(uint64_t x) {
if (x == 0) return 0;
uint64_t a = 1, b = x, c;
while (b - a > 1) {
c = (b - a) / 2 + a;
if (is_safely_multiplicable(c, c) && x >= c * c) a = c;
else b = c;
}
return a;
}
bool is_square_number(uint64_t x) {
uint64_t rt_x = sqrt_int(x);
return rt_x * rt_x == x;
}
#include <iostream>
#include <limits>
template <typename T>
bool is_safely_multiplicable(T a, T b) {
if (b == 0) return true;
return a <= (std::numeric_limits<T>::max() / b);
}
uint64_t sqrt_int(uint64_t x) {
if (x == 0) return 0;
uint64_t a = 1, b = x, c;
while (b - a > 1) {
c = (b - a) / 2 + a;
if (is_safely_multiplicable(c, c) && x >= c * c) a = c;
else b = c;
}
return a;
}
bool is_square_number(uint64_t x) {
uint64_t rt_x = sqrt_int(x);
return rt_x * rt_x == x;
}
956デフォルトの名無しさん
2019/02/01(金) 18:32:58.66ID:iYm26gLc >>954
Haskell
isSqr a = isSqr' a [n * n | n <- [0..a]]
where
isSqr' x (l:_) |(x == l) = True
isSqr' x (l:_) |(x < l) = False
isSqr' x (_:ls) = isSqr' x ls
Haskell
isSqr a = isSqr' a [n * n | n <- [0..a]]
where
isSqr' x (l:_) |(x == l) = True
isSqr' x (l:_) |(x < l) = False
isSqr' x (_:ls) = isSqr' x ls
957デフォルトの名無しさん
2019/02/01(金) 20:32:42.61ID:iYm26gLc Python
def isSqr(x):
n = 0
a = 0
while x <= a:
if x == a:
return True
n += 1
a = n * n
return False
リスト内包表記だとすぐ落ちた。。。
Haskellより持たない。。。
def isSqr(x):
n = 0
a = 0
while x <= a:
if x == a:
return True
n += 1
a = n * n
return False
リスト内包表記だとすぐ落ちた。。。
Haskellより持たない。。。
958デフォルトの名無しさん
2019/02/01(金) 21:10:55.66ID:qFuZrTt3 これはいかに巨大数を高速判定するかだろ
959デフォルトの名無しさん
2019/02/01(金) 21:19:45.40ID:qFuZrTt3 目星はこれを荒く使えばつきそう
逆数と平方根を求める高次収束アルゴリズム
http://www.finetune.co.jp/~lyuka/technote/fract/sqrt.html
高速根号計算
http://takashiijiri.com/study/miscs/fastsqrt.html
逆数と平方根を求める高次収束アルゴリズム
http://www.finetune.co.jp/~lyuka/technote/fract/sqrt.html
高速根号計算
http://takashiijiri.com/study/miscs/fastsqrt.html
960デフォルトの名無しさん
2019/02/01(金) 22:13:46.18ID:kNvVsHFY961デフォルトの名無しさん
2019/02/01(金) 23:10:20.49ID:iYm26gLc >>958
m = ？
n = m * m
m=10の時nの1/10で、m=100だとnの1/100
探すべきmが高速で小さくなって行く。
大きな数になる程、0<=mなループの始まり求めるの無理ゲー。。。
m = ？
n = m * m
m=10の時nの1/10で、m=100だとnの1/100
探すべきmが高速で小さくなって行く。
大きな数になる程、0<=mなループの始まり求めるの無理ゲー。。。
962デフォルトの名無しさん
2019/02/01(金) 23:12:44.04ID:VdCPb4pG ttps://ideone.com/xr1vcE
Newton-Raphson法ほど速くはないがせめてO(log(N))くらいで実装せにゃお題として出てくる意味がなかろう
Newton-Raphson法ほど速くはないがせめてO(log(N))くらいで実装せにゃお題として出てくる意味がなかろう
963デフォルトの名無しさん
2019/02/02(土) 00:09:40.70ID:/6KX0oFw >>954 Lua
function isSquare(n)
local a = 0
local i = 1
while a < n do
a = a + i
i = i + 2
end
return a == n
end
print(isSquare(100000))
print(isSquare(10000))
function isSquare(n)
local a = 0
local i = 1
while a < n do
a = a + i
i = i + 2
end
return a == n
end
print(isSquare(100000))
print(isSquare(10000))
964デフォルトの名無しさん
2019/02/02(土) 03:47:32.54ID:i2SNxKFt いっそsqrt 関数実装。
頭悪いんで数学的じゃない。
256以上入れるとフリーズするけど、その範囲なら精度抜群。
精度を1つ落とすと(range(15)の値を1つ減らすと)2桁くらい大きな数を入れてもフリーズしなくなる。
def sqrt(x):
i = 1
if x == 0:
return 0
while x >= (i * i):
if x / i == i:
return i
i += 1
i -= 1
a = 0.1
for j in range(15):
while x >= (i * i):
i += a
i -= a
a *= 0.1
return (i - a)
頭悪いんで数学的じゃない。
256以上入れるとフリーズするけど、その範囲なら精度抜群。
精度を1つ落とすと(range(15)の値を1つ減らすと)2桁くらい大きな数を入れてもフリーズしなくなる。
def sqrt(x):
i = 1
if x == 0:
return 0
while x >= (i * i):
if x / i == i:
return i
i += 1
i -= 1
a = 0.1
for j in range(15):
while x >= (i * i):
i += a
i -= a
a *= 0.1
return (i - a)
965デフォルトの名無しさん
2019/02/02(土) 09:51:18.58ID:/6KX0oFw >>954 J
f =: = *: @ <. @ %:
f =: = *: @ <. @ %:
966デフォルトの名無しさん
2019/02/02(土) 13:18:27.64ID:CwD+xRo8 平方数かどうかを高速に判定する方法 - hnwの日記
https://hnw.hatenablog.com/entry/20140503
GNU MPのmpz_perfect_square_p関数の実装
GNU MPのソースコードを確認してみたところ、次のような処理だとわかります。
引数nのmod 256を計算し、平方数ではない数を判定する
平方数のmod 256は44種類の値しか出現しないので、入力の82.8%は平方数でないと判定できる
同様に入力nのmod 9, 5, 7, 13, 17（64-bitシステムではmod 97も）を計算し、平方数ではない数を判定する
これにより入力の99.25%（64-bitシステムでは99.62%）について平方数でないと判定できる
最後に、平方根を計算して平方数かどうか確認する
平方数ではない数の多くを事前にふるい落とし、判定できなかった数だけ真面目に平方根を求める、という方針だとわかります。
もちろん、GNU MPの場合のmod pの選び方は多倍長整数演算ならではだと言えます。
mod 256はサイズNにかかわらずO(1)で計算できるので、最初に行うことで全体の高速化に貢献できます。
また、2ステップ目の計算も2^24-1 = 9 * 5 * 7 * 13 * 17 * ...であることを利用し、多倍長整数であっても比較的高速に計算できるような実装になっています。
10進整数でmod 9を求める場合に全部の桁を足し合わせてからmod 9しても同じ結果になるのと同様、
下位から24bit区切りの数を足し合わせてからmod 2^24-1を計算することで、元の数のmod 2^24-1が計算できるのです。
https://hnw.hatenablog.com/entry/20140503
GNU MPのmpz_perfect_square_p関数の実装
GNU MPのソースコードを確認してみたところ、次のような処理だとわかります。
引数nのmod 256を計算し、平方数ではない数を判定する
平方数のmod 256は44種類の値しか出現しないので、入力の82.8%は平方数でないと判定できる
同様に入力nのmod 9, 5, 7, 13, 17（64-bitシステムではmod 97も）を計算し、平方数ではない数を判定する
これにより入力の99.25%（64-bitシステムでは99.62%）について平方数でないと判定できる
最後に、平方根を計算して平方数かどうか確認する
平方数ではない数の多くを事前にふるい落とし、判定できなかった数だけ真面目に平方根を求める、という方針だとわかります。
もちろん、GNU MPの場合のmod pの選び方は多倍長整数演算ならではだと言えます。
mod 256はサイズNにかかわらずO(1)で計算できるので、最初に行うことで全体の高速化に貢献できます。
また、2ステップ目の計算も2^24-1 = 9 * 5 * 7 * 13 * 17 * ...であることを利用し、多倍長整数であっても比較的高速に計算できるような実装になっています。
10進整数でmod 9を求める場合に全部の桁を足し合わせてからmod 9しても同じ結果になるのと同様、
下位から24bit区切りの数を足し合わせてからmod 2^24-1を計算することで、元の数のmod 2^24-1が計算できるのです。
967デフォルトの名無しさん
2019/02/02(土) 13:43:47.68ID:OgiywF+Q >>944
えぇ…
えぇ…
968デフォルトの名無しさん
2019/02/02(土) 16:31:05.91ID:9W2pTWu+ >>967 その心は？
969デフォルトの名無しさん
2019/02/02(土) 16:42:03.23ID:OgiywF+Q970デフォルトの名無しさん
2019/02/02(土) 17:31:46.69ID:rEiZ26fd 本質的に高度な数学的知識が要求され
る問題は板違いかと
言語によっては、問題にジャストフィット
するライブラリが標準で添付されているとか
ライブラリ管理コミュニティを通じて
容易に入手できるとかあるかもしれないが
その場合は紹介程度に。
言語処理系の言語処理の為の機能提供
である場合を除いて多くは本質的に多言
語に対応しており〜言語用のライブラリと
いう表現は兎も角〜言語のライブラリと表
現すると曖昧で誤解を招きやすい表現
として嫌われる場合もありえることに注意。
（別言語で記述される場合もある）
る問題は板違いかと
言語によっては、問題にジャストフィット
するライブラリが標準で添付されているとか
ライブラリ管理コミュニティを通じて
容易に入手できるとかあるかもしれないが
その場合は紹介程度に。
言語処理系の言語処理の為の機能提供
である場合を除いて多くは本質的に多言
語に対応しており〜言語用のライブラリと
いう表現は兎も角〜言語のライブラリと表
現すると曖昧で誤解を招きやすい表現
として嫌われる場合もありえることに注意。
（別言語で記述される場合もある）
971デフォルトの名無しさん
2019/02/02(土) 17:36:46.77ID:XgXX/tZQ なんしたのきゅうに
972デフォルトの名無しさん
2019/02/02(土) 17:57:26.00ID:rEiZ26fd セルフホスティング対応な言語でライブラリを
自前で構成しているもので、他言語に移植されて
いない独自の機能を持つものもあるかもしれない
し、もっと極端に言えば独自性の高いライブラリに
最適化された言語を自前で作ってそれで記述され
ているものもあるかもしれないけど一般的には入手
は容易ではないかも
自前で構成しているもので、他言語に移植されて
いない独自の機能を持つものもあるかもしれない
し、もっと極端に言えば独自性の高いライブラリに
最適化された言語を自前で作ってそれで記述され
ているものもあるかもしれないけど一般的には入手
は容易ではないかも
973デフォルトの名無しさん
2019/02/02(土) 19:52:10.51ID:hDNgHqpo じゃあさ、各々複数行の
a.txt と、b.txtを並べて表示して。
a.txt と、b.txtを並べて表示して。
974デフォルトの名無しさん
2019/02/02(土) 23:25:24.42ID:g8xy/J6N 辺に沿って動くとき、AからBまでの最短経路はいくつあるか
┏┳┳┳┓B
┣╋╋╋┫
┏┳╋╋╋╋┫
┏┳┳┳╋╋╋╋╋┻┛
┣╋╋╋╋╋╋╋┫
┣╋╋╋╋╋╋╋┫
┣╋╋╋╋╋╋╋┛
┗┻┻┻┻┻┻┛
A
┏┳┳┳┓B
┣╋╋╋┫
┏┳╋╋╋╋┫
┏┳┳┳╋╋╋╋╋┻┛
┣╋╋╋╋╋╋╋┫
┣╋╋╋╋╋╋╋┫
┣╋╋╋╋╋╋╋┛
┗┻┻┻┻┻┻┛
A
975デフォルトの名無しさん
2019/02/03(日) 00:45:38.22ID:UGH880J+ おねえさんがロボットになるやつかｗ
976デフォルトの名無しさん
2019/02/03(日) 00:55:19.58ID:UTpNqxd4 中学受験の算数の問題に経路を数える問題があるよな
つまり小学生でも手計算で解ける
つまり小学生でも手計算で解ける
977デフォルトの名無しさん
2019/02/03(日) 00:57:22.00ID:72eosYJ+ >>970 アホちゃうの？
プログラムというのは、問題解決のための言語であって、問題が解決できなければどんなに美しい言語でも存在価値はない。
最速で美しく問題解決にたどり着ける言語が一番。 言語そのものの美しさなんて二の次。
英語でも日本語でも何でも良い。 利用できるものは何でも利用すれば良い。 好き嫌い言ってる奴はアホ。
どんなに綺麗な国の言葉でも高等教育が自国語でできない国がほとんど、それは利用できる教材が自国語で書いたものがないから。
日本語だったら英語を知らなくてもノーベル賞が取れるだけの教材(ライブラリ)が揃ってる。
プログラムというのは、問題解決のための言語であって、問題が解決できなければどんなに美しい言語でも存在価値はない。
最速で美しく問題解決にたどり着ける言語が一番。 言語そのものの美しさなんて二の次。
英語でも日本語でも何でも良い。 利用できるものは何でも利用すれば良い。 好き嫌い言ってる奴はアホ。
どんなに綺麗な国の言葉でも高等教育が自国語でできない国がほとんど、それは利用できる教材が自国語で書いたものがないから。
日本語だったら英語を知らなくてもノーベル賞が取れるだけの教材(ライブラリ)が揃ってる。
978デフォルトの名無しさん
2019/02/03(日) 01:53:34.29ID:wxxHWwaf 別に高度な数学が必要になることは構わないけど、数学だけで話が閉じちゃうような問題はつまらないとは思う。数学的に解を求める解法自体が主題で、コーディングはただその解法を(日本語、英語などの自然言語と同様に)とあるプログラミング言語で記述しただけのような。
もちろん数学に価値がないといっているのではなく、話のウェイトとしてプログラムである部分がほぼないのなら、その話、ここじゃなくてもいいじゃんと思うよ。
もちろん数学に価値がないといっているのではなく、話のウェイトとしてプログラムである部分がほぼないのなら、その話、ここじゃなくてもいいじゃんと思うよ。
979デフォルトの名無しさん
2019/02/03(日) 04:02:35.76ID:JyP+XfGy つまらない問題だと思ったら無視すればいいだけじゃないかな
その問題に対して回答する人がいるなら、回答者(と出題者)にとっては
つまらない問題ではないってことなんだろうし
その問題に対して回答する人がいるなら、回答者(と出題者)にとっては
つまらない問題ではないってことなんだろうし
980デフォルトの名無しさん
2019/02/03(日) 04:46:34.72ID:1lu6X4vo >>977
ＯＳ記述の言語は言語じゃないのか？
コンパイラなどの言語処理系記述の言語
も言語じゃないってわけだな
問題解決手段の一つではあるかもしれないが
通信系ソフトウェアの存在はどう説明？
既存の資産を維持しより効率的に活用する研究が
新しいプログラミング言語が次々と生成されてくる原
動力じゃないかと
ＯＳ記述の言語は言語じゃないのか？
コンパイラなどの言語処理系記述の言語
も言語じゃないってわけだな
問題解決手段の一つではあるかもしれないが
通信系ソフトウェアの存在はどう説明？
既存の資産を維持しより効率的に活用する研究が
新しいプログラミング言語が次々と生成されてくる原
動力じゃないかと
981デフォルトの名無しさん
2019/02/03(日) 05:08:39.33ID:l3Qt7IvN >>977
>日本語だったら英語を知らなくてもノーベル賞が取れるだけの教材(ライブラリ)が揃ってる。
ダウト。ノーベル賞を受賞した研究者はみんな英語で論文書いている。
あと、今の日本語での教育教材の蓄積があるのは、大雑把に言って、
ここ10年でノーベル賞受賞した面々の世代の研究者や技術者が書いた。
もちろん彼らは英語の文献を読み、勉強し、日本語の文献を書いた。
だから、英語の文献なんて要らないみたいな話はアホ。
ソファでポテチ食いながら「ジャガイモなんて買えばいいんだから、育ててる奴はアホ」と言ってるようなもの。
>日本語だったら英語を知らなくてもノーベル賞が取れるだけの教材(ライブラリ)が揃ってる。
ダウト。ノーベル賞を受賞した研究者はみんな英語で論文書いている。
あと、今の日本語での教育教材の蓄積があるのは、大雑把に言って、
ここ10年でノーベル賞受賞した面々の世代の研究者や技術者が書いた。
もちろん彼らは英語の文献を読み、勉強し、日本語の文献を書いた。
だから、英語の文献なんて要らないみたいな話はアホ。
ソファでポテチ食いながら「ジャガイモなんて買えばいいんだから、育ててる奴はアホ」と言ってるようなもの。
982981
2019/02/03(日) 05:15:20.75ID:l3Qt7IvN プログラミング言語も一緒。
Cは泥臭い実用言語で、その前にはAlgolなど「綺麗な」言語が下敷きとしてあった。
Cは現在でも多くの問題を解決するための道具として活用されているが、
現在では使われていないAlgolなどの「綺麗な」言語なしにCは生まれなかった。
Cは泥臭い実用言語で、その前にはAlgolなど「綺麗な」言語が下敷きとしてあった。
Cは現在でも多くの問題を解決するための道具として活用されているが、
現在では使われていないAlgolなどの「綺麗な」言語なしにCは生まれなかった。
983デフォルトの名無しさん
2019/02/03(日) 05:39:02.55ID:1lu6X4vo 個々の言語の歴史観の講釈はスレ違い
確かに問題解決の道具であることが
中心かもしれないが、扱ってきた問
題の種類によって文法やライブラリ・
その取り扱い方に差異が生じている
が、同じ問題を別言語で「解く」と
優劣の違いがわかって面白いかもし
れない。
が、あんまし長くやってると不毛な
言語比較論、文化比較論になったり
して色々ヤバいからそろそろ一旦
お開きにしたら？
続きはそれぞれの言語別スレッドで
ということで
確かに問題解決の道具であることが
中心かもしれないが、扱ってきた問
題の種類によって文法やライブラリ・
その取り扱い方に差異が生じている
が、同じ問題を別言語で「解く」と
優劣の違いがわかって面白いかもし
れない。
が、あんまし長くやってると不毛な
言語比較論、文化比較論になったり
して色々ヤバいからそろそろ一旦
お開きにしたら？
続きはそれぞれの言語別スレッドで
ということで
984デフォルトの名無しさん
2019/02/03(日) 07:04:30.87ID:LaZtKDWq お題：プログラムの実行時刻が午前なら「おはようございます、ご主人様！」、午後なら「お疲れ様です、ご主人様！」と表示させる
985デフォルトの名無しさん
2019/02/03(日) 07:48:34.73ID:AEg+fU/i >>984 C
time_t now = time(NULL);
struct tm *p = localtime(&now);
if (p->tm_hour * 60 + p->tm_min < 12 * 60) {
printf("おはようございます、ご主人様！\n");
} else if (p->tm_hour * 60 + p->tm_min > 12 * 60) {
printf("お疲れ様です、ご主人様！\n");
}
time_t now = time(NULL);
struct tm *p = localtime(&now);
if (p->tm_hour * 60 + p->tm_min < 12 * 60) {
printf("おはようございます、ご主人様！\n");
} else if (p->tm_hour * 60 + p->tm_min > 12 * 60) {
printf("お疲れ様です、ご主人様！\n");
}
986デフォルトの名無しさん
2019/02/03(日) 08:56:17.85ID:l3Qt7IvN >>984
Pharo Smalltalk
Smalltalk ui inform: (Time now meridianAbbreviation = 'AM' ifTrue: [ 'おはようございます、ご主人様！' ] ifFalse: [ 'お疲れ様です、ご主人様！' ])
Pharo Smalltalk
Smalltalk ui inform: (Time now meridianAbbreviation = 'AM' ifTrue: [ 'おはようございます、ご主人様！' ] ifFalse: [ 'お疲れ様です、ご主人様！' ])
987デフォルトの名無しさん
2019/02/03(日) 09:19:47.58ID:72eosYJ+ python
from datetime import datetime
if datetime.now().hour < 12:
print('おはようございますご主人様')
else:
print('お疲れ様です、ご主人様')
from datetime import datetime
if datetime.now().hour < 12:
print('おはようございますご主人様')
else:
print('お疲れ様です、ご主人様')
988デフォルトの名無しさん
2019/02/03(日) 09:23:27.86ID:1lu6X4vo 午前12時=00:00
午後12時=12:00
23:59の後は00:00
午前12時=12:00
午後12時=24:00
24:00の後は00:01
午後12時=12:00
23:59の後は00:00
午前12時=12:00
午後12時=24:00
24:00の後は00:01
989デフォルトの名無しさん
2019/02/03(日) 09:40:26.58ID:cfde/ig7990デフォルトの名無しさん
2019/02/03(日) 09:58:12.21ID:I0qputsI >>964 のHaskell版。
負の数の場合の処理、負の数含め、絶対値が256以上だった場合エラー吐く様に処理を追加。
※Haskellは整数と少数を明確に分ける為、渡す数に小数点が無いとエラーになる。
mysqrt x = mysqrt' x 0
where
mysqrt' x m |x < 0 = - mysqrt (abs x)
mysqrt' x m |x == m * m = m
mysqrt' x m |x < m * m = fsqrt x (m - 1) 0.1 15
mysqrt' x m = mysqrt' x (m + 1)
fsqrt _ a _ 0 = a
fsqrt v _ _ _ | v > 256 = error "\"fsqr\":out of range 0..256"
fsqrt v a f n | v <= a * a = fsqrt v (a - f) (f * 0.1) (n - 1)
fsqrt v a f n = fsqrt v (a + f) f n
使用例
main = print (mysqrt x) >> print (mysqrt x * mysqrt x)
結果
1.41421356237309
2.0
負の数の場合の処理、負の数含め、絶対値が256以上だった場合エラー吐く様に処理を追加。
※Haskellは整数と少数を明確に分ける為、渡す数に小数点が無いとエラーになる。
mysqrt x = mysqrt' x 0
where
mysqrt' x m |x < 0 = - mysqrt (abs x)
mysqrt' x m |x == m * m = m
mysqrt' x m |x < m * m = fsqrt x (m - 1) 0.1 15
mysqrt' x m = mysqrt' x (m + 1)
fsqrt _ a _ 0 = a
fsqrt v _ _ _ | v > 256 = error "\"fsqr\":out of range 0..256"
fsqrt v a f n | v <= a * a = fsqrt v (a - f) (f * 0.1) (n - 1)
fsqrt v a f n = fsqrt v (a + f) f n
使用例
main = print (mysqrt x) >> print (mysqrt x * mysqrt x)
結果
1.41421356237309
2.0
991デフォルトの名無しさん
2019/02/03(日) 10:00:13.11ID:I0qputsI fsqrt v _ _ _ | v > 256 = error "\"fsqr\":out of range 0..256"
fsqrt v a f n | v <= a * a = fsqrt v (a - f) (f * 0.1) (n - 1)
fsqrt v a f n = fsqrt v (a + f) f n
fsqrt v a f n | v <= a * a = fsqrt v (a - f) (f * 0.1) (n - 1)
fsqrt v a f n = fsqrt v (a + f) f n
992デフォルトの名無しさん
2019/02/03(日) 10:13:31.71ID:xEPkQ4sk この改行長文おじさんはどこからきたの。笑
>>977
＞日本語だったら英語を知らなくてもノーベル賞が取れるだけの教材(ライブラリ)が揃ってる。
最近はそうでもないようですよ…ペーパーは英語だし、教科書＝テキストレベルでも英語でしか発刊されない状況といいます
haskell をやろうとして圏論の教科書を探しましたが、欧米の本の和訳ばかりで日本人が書いた圏論の教科書はありませんでした
＞日本語だったら英語を知らなくてもノーベル賞が取れるだけの教材(ライブラリ)が揃ってる。
最近はそうでもないようですよ…ペーパーは英語だし、教科書＝テキストレベルでも英語でしか発刊されない状況といいます
haskell をやろうとして圏論の教科書を探しましたが、欧米の本の和訳ばかりで日本人が書いた圏論の教科書はありませんでした
994デフォルトの名無しさん
2019/02/03(日) 10:35:41.62ID:/jO+7TC8 誰か次スレ頼む
995デフォルトの名無しさん
2019/02/03(日) 11:07:15.54ID:72eosYJ+ お題1: 現在地の緯度、経度を出せ
緯度:、、、、
経度:、、、、
お題2: 東京都新宿区西新宿2丁目8-1 の緯度、経度を出せ
緯度:、、、
経度:、、、
お題3: お題2で求めた緯度経度から住所を出せ
郵便番号:、、、
住所:東京都、、、、
緯度:、、、、
経度:、、、、
お題2: 東京都新宿区西新宿2丁目8-1 の緯度、経度を出せ
緯度:、、、
経度:、、、
お題3: お題2で求めた緯度経度から住所を出せ
郵便番号:、、、
住所:東京都、、、、
996デフォルトの名無しさん
2019/02/03(日) 11:22:46.76ID:72eosYJ+997デフォルトの名無しさん
2019/02/03(日) 11:23:36.56ID:72eosYJ+ >>995 は、次スレに移動させるね。
998デフォルトの名無しさん
2019/02/03(日) 17:37:12.27ID:csrqlAvs うめ
999デフォルトの名無しさん
2019/02/03(日) 17:37:50.59ID:oUppVF8S1000デフォルトの名無しさん
2019/02/03(日) 17:38:30.96ID:csrqlAvs10011001
Over 1000Thread このスレッドは１０００を超えました。
新しいスレッドを立ててください。
life time: 128日 7時間 29分 24秒
新しいスレッドを立ててください。
life time: 128日 7時間 29分 24秒
10021002
Over 1000Thread 5ちゃんねるの運営はプレミアム会員の皆さまに支えられています。
運営にご協力お願いいたします。
───────────────────
《プレミアム会員の主な特典》
★ 5ちゃんねる専用ブラウザからの広告除去
★ 5ちゃんねるの過去ログを取得
★ 書き込み規制の緩和
───────────────────
会員登録には個人情報は一切必要ありません。
月300円から匿名でご購入いただけます。
▼ プレミアム会員登録はこちら ▼
https://premium.5ch.net/
▼ 浪人ログインはこちら ▼
https://login.5ch.net/login.php
運営にご協力お願いいたします。
───────────────────
《プレミアム会員の主な特典》
★ 5ちゃんねる専用ブラウザからの広告除去
★ 5ちゃんねるの過去ログを取得
★ 書き込み規制の緩和
───────────────────
会員登録には個人情報は一切必要ありません。
月300円から匿名でご購入いただけます。
▼ プレミアム会員登録はこちら ▼
https://premium.5ch.net/
▼ 浪人ログインはこちら ▼
https://login.5ch.net/login.php
レス数が1000を超えています。これ以上書き込みはできません。
ニュース
- 橋下徹氏 外務省幹部の訪中受け「口だけ番長」へ痛烈指摘 「喧嘩は日本の完敗…なんとかっこ悪い日本か」 [冬月記者★]
- 【外国人問題】小野田紀美担当相「不法就労や不法滞在は許さない」 [シャチ★]
- 【野球】井端監督 大谷翔平、山本由伸らのWBCへの参加 「1日も早く返事ほしい」「待っててといっても、国内組が遅くなってしまう」★3 [冬月記者★]
- 経団連会長、日中は建設的対話を 経済3団体が高市首相と初会談も日中関係は話題に登らず… [BFU★]
- 中国で「クレしん」公開延期 対日報復、エンタメに波及 [蚤の市★]
- 東京株式市場 インバウンド関連株が下落 中国政府の渡航自粛要請で [バイト歴50年★]
- 中国「高市が頭を下げて謝罪しない限り、絶対に許さない」 [329329848]
- 🏡
- 有識者「高市総理が発言を撤回したり、辞職するしかないと言っている人は、それで日中関係が今まで通りになると思ってる？」 [834922174]
- 減税は低所得者差別
- さわやかハンバーグって何が凄いの？？ [134367759]
- 【動画】陽キャ、パルクール失敗