プログラミングのお題スレです。
【出題と回答例】
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/
宿題は宿題スレがあるのでそちらへ。
※前スレ
https://mevius.2ch.net/test/read.cgi/tech/1538096947/
探検
プログラミングのお題スレ Part13
■ このスレッドは過去ログ倉庫に格納されています
1デフォルトの名無しさん
2019/02/03(日) 11:21:53.20ID:72eosYJ+ >>36
C++ 多桁長演算(加減乗除)の一環として、mpz_class との互換を目指していました
https://mevius.5ch.net/test/read.cgi/tech/1434079972/51
C++ 多桁長演算(加減乗除)の一環として、mpz_class との互換を目指していました
https://mevius.5ch.net/test/read.cgi/tech/1434079972/51
2019/02/04(月) 20:50:12.06ID:Ly7rB5Pz
39デフォルトの名無しさん
2019/02/04(月) 21:00:40.83ID:mK0I6q3a modをみる
nが平方数なら、n=x^2 だが、n=x^2 (mod m)でもある
逆にmod で平方数でなければ、元々も平方数ではない
mod 3だと0 1 は平方数だが、2はちがう。3i + 2 は平方数にはならない
nが平方数なら、n=x^2 だが、n=x^2 (mod m)でもある
逆にmod で平方数でなければ、元々も平方数ではない
mod 3だと0 1 は平方数だが、2はちがう。3i + 2 は平方数にはならない
40デフォルトの名無しさん
2019/02/04(月) 22:02:57.78ID:eX/1kX5o >>38
元は小学生にプログラミングを通じて、割り算への理解を深めてもらえないかと考えたんで、単純に演算子を置き換えて欲しくないかも。。。
元は小学生にプログラミングを通じて、割り算への理解を深めてもらえないかと考えたんで、単純に演算子を置き換えて欲しくないかも。。。
41デフォルトの名無しさん
2019/02/05(火) 10:18:54.39ID:ij/1zyvC 小学生の時テストで0点をもらった間違った割り算の
やり方をプログラムにしてみた。
--Lua
function myDivMod(a, b)
local r = 0
while a > 0 do
a = a - b
r = r + 1
end
return r, -a
end
print(myDivMod(10,3))
実行結果
4 2
やり方をプログラムにしてみた。
--Lua
function myDivMod(a, b)
local r = 0
while a > 0 do
a = a - b
r = r + 1
end
return r, -a
end
print(myDivMod(10,3))
実行結果
4 2
2019/02/05(火) 12:28:40.98ID:aE6b0ZPr
>>35 結局 整数のsqrt を作って、範囲の中に納まる最小最大の整数のsqrt を取り出し、その差(+1)がその範囲の中にある平方数の個数と言う作りにした。
ポイントとなった整数のsqrt が秀逸だったのでここに書いておく。 どんなに巨大な数字でも数回のシフト操作だけで終わるから極端にスピードが速い。
ソースは、gist.github.com/bnlucas/5879594
# integer square root
def isqrt_2(n):
if n < 0:
raise ValueError('Square root is not defined for negative numbers.')
x = int(n)
if x == 0:
return 0
a, b = divmod(x.bit_length(), 2)
#divmod(a, b)は(a // b, a % b)のタプルを返す。
#平方数は半分のビット数以下だからそれを最大値で計算開始
n = 2 ** (a + b)
while True:
y = (n + x // n) >> 1 #1bit右にシフト
if y >= n:
return n
n = y
-----------------
#結果 0.0 秒 count= 1000000000
#範囲 999999999000000000000000000000000000 1000000001000000000000000000000000000
#入力bit_length()=120 入力bit_length()=120
平方数範囲 999999999500000000 1000000000499999999
上の二乗 999999999000000000250000000000000000 1000000000999999998249999999000000001
ポイントとなった整数のsqrt が秀逸だったのでここに書いておく。 どんなに巨大な数字でも数回のシフト操作だけで終わるから極端にスピードが速い。
ソースは、gist.github.com/bnlucas/5879594
# integer square root
def isqrt_2(n):
if n < 0:
raise ValueError('Square root is not defined for negative numbers.')
x = int(n)
if x == 0:
return 0
a, b = divmod(x.bit_length(), 2)
#divmod(a, b)は(a // b, a % b)のタプルを返す。
#平方数は半分のビット数以下だからそれを最大値で計算開始
n = 2 ** (a + b)
while True:
y = (n + x // n) >> 1 #1bit右にシフト
if y >= n:
return n
n = y
-----------------
#結果 0.0 秒 count= 1000000000
#範囲 999999999000000000000000000000000000 1000000001000000000000000000000000000
#入力bit_length()=120 入力bit_length()=120
平方数範囲 999999999500000000 1000000000499999999
上の二乗 999999999000000000250000000000000000 1000000000999999998249999999000000001
43デフォルトの名無しさん
2019/02/05(火) 12:33:09.36ID:NCwCR2JI >>41
おしいなw
おしいなw
2019/02/05(火) 13:45:18.58ID:jFne2R1T
掛け算は簡単だけど除算は結構面倒
ttps://ideone.com/EheeYM
ttps://ideone.com/EheeYM
2019/02/05(火) 18:45:08.63ID:63VtM8MC
>>42 の isqrt_2 を使ったパフォーマンステスト。
次のようなのを継ぎ足してテストした。 例によってインデント部は全角空白に変換してるから、逆変換しないと動かない。
def isSqrt(n):
return n == isqrt_2(n)**2
v0 = 12345678901234567890
v = v0**2 # 整数平方される対象の数値
loopc = 100000 # をこの回数繰り返す。
isqr=0
start =time.process_time()
for i in range(loopc): isqr=isqrt_2(v)
end =time.process_time()
print('#整数平方(v)の結果',end-start,'秒')
print(' 繰返し数の回数',loopc),print(),print('#v0 ',v0)
print('#v=v0**2=',v),
print('#isqrt(v)',isqr)
print('#上の**2',isqr**2)
print('対象数vのビット数',v.bit_length(),'bit')
print('vが平方数かどうかの判定',isSqrt(v))
-----
#整数平方(v)の結果 0.22398700000002236 秒
繰返し数の回数 100000
#v0 12345678901234567890
#v=v0**2= 152415787532388367501905199875019052100
#isqrt(v) 12345678901234567890
#上の**2 152415787532388367501905199875019052100
対象数vのビット数 127 bit
vが平方数かどうかの判定 True
次のようなのを継ぎ足してテストした。 例によってインデント部は全角空白に変換してるから、逆変換しないと動かない。
def isSqrt(n):
return n == isqrt_2(n)**2
v0 = 12345678901234567890
v = v0**2 # 整数平方される対象の数値
loopc = 100000 # をこの回数繰り返す。
isqr=0
start =time.process_time()
for i in range(loopc): isqr=isqrt_2(v)
end =time.process_time()
print('#整数平方(v)の結果',end-start,'秒')
print(' 繰返し数の回数',loopc),print(),print('#v0 ',v0)
print('#v=v0**2=',v),
print('#isqrt(v)',isqr)
print('#上の**2',isqr**2)
print('対象数vのビット数',v.bit_length(),'bit')
print('vが平方数かどうかの判定',isSqrt(v))
-----
#整数平方(v)の結果 0.22398700000002236 秒
繰返し数の回数 100000
#v0 12345678901234567890
#v=v0**2= 152415787532388367501905199875019052100
#isqrt(v) 12345678901234567890
#上の**2 152415787532388367501905199875019052100
対象数vのビット数 127 bit
vが平方数かどうかの判定 True
2019/02/06(水) 11:48:35.22ID:Cmz9AyOj
>>45 同じ条件で2分探索法でやると 5.5秒かかった。
2019/02/06(水) 16:53:02.43ID:XOfzhWu4
Wikipediaに Integer square root
https://en.wikipedia.org/wiki/Integer_square_root
があり、その中の 2.1 Using bitwise operations
の二つを試してみたが、
最初のrecursive call を使った方が 1.65秒
次の方が 2.05秒
早いことは早いが、>>42 >>45 のビットシフト法の方がかなり早い。
0.22秒
gmpのisqrt は早そうだが Pythonistaでは使えないので試していない。
https://en.wikipedia.org/wiki/Integer_square_root
があり、その中の 2.1 Using bitwise operations
の二つを試してみたが、
最初のrecursive call を使った方が 1.65秒
次の方が 2.05秒
早いことは早いが、>>42 >>45 のビットシフト法の方がかなり早い。
0.22秒
gmpのisqrt は早そうだが Pythonistaでは使えないので試していない。
2019/02/09(土) 01:16:57.63ID:VrkeVQvn
>>4 そう言う記述を度々見るんだが、具体的なプログラムを提示してくれない?
2019/02/09(土) 05:46:19.68ID:ZuP9aKcb
>>48
前スレでもGMPだって紹介されてたし
https://gmplib.org/manual/Perfect-Square-Algorithm.html
https://gmplib.org/repo/gmp/file/tip/mpn/generic/perfsqr.c
実際に使われる法はビルド時生成(https://gmplib.org/repo/gmp/file/tip/gen-psqr.c)だけど大抵の32/64bitシステムの場合は書いてある通り
前スレでもGMPだって紹介されてたし
https://gmplib.org/manual/Perfect-Square-Algorithm.html
https://gmplib.org/repo/gmp/file/tip/mpn/generic/perfsqr.c
実際に使われる法はビルド時生成(https://gmplib.org/repo/gmp/file/tip/gen-psqr.c)だけど大抵の32/64bitシステムの場合は書いてある通り
50デフォルトの名無しさん
2019/02/09(土) 06:52:32.52ID:y7fm8J5o なんでも放題にすればいい
お題
平方数は 256で割るとあまりの
パターンが44種類になるという。
この44種類を求める。
お題
平方数は 256で割るとあまりの
パターンが44種類になるという。
この44種類を求める。
51デフォルトの名無しさん
2019/02/09(土) 08:05:29.21ID:y7fm8J5o 放題は、お題の間違い
2019/02/09(土) 08:30:22.91ID:IKLi4q/e
python
{(x**2) % 256 for x in range(0,256)}
{(x**2) % 256 for x in range(0,256)}
2019/02/09(土) 10:10:19.89ID:DhY4ZB+P
Ruby
p 0x100.times.map{|i| i*i&0xFF}.uniq.sort
# => [0, 1, 4, [...], 233, 241, 249]
p 0x100.times.map{|i| i*i&0xFF}.uniq.sort
# => [0, 1, 4, [...], 233, 241, 249]
2019/02/09(土) 10:33:40.76ID:O4yJeWlE
javascript
[...new Set(function*(){for(let i=0;i<256;i++)yield i*i%256}())].sort((a,b)=>a-b)
[...new Set(function*(){for(let i=0;i<256;i++)yield i*i%256}())].sort((a,b)=>a-b)
55デフォルトの名無しさん
2019/02/09(土) 14:20:17.67ID:BaccQTUO お題:ハノイの塔の最少手数は一種類しかないのか
2019/02/09(土) 14:22:34.94ID:VrkeVQvn
>>50 その44種(mod256)の中に x%256 が一致するか
python
mod256=sorted({(i**2)%256 for i in range(256)})
x=123*123
print(x%256 in mod256)
# True
これでmod の話が理解できた。 要は完全な平方数じゃないのは平方根の計算をしないと言う話ね。
python
mod256=sorted({(i**2)%256 for i in range(256)})
x=123*123
print(x%256 in mod256)
# True
これでmod の話が理解できた。 要は完全な平方数じゃないのは平方根の計算をしないと言う話ね。
2019/02/09(土) 15:28:05.40ID:1XyVHoA8
>>50 Perl5
%a = map{$_=>1} map{$_*$_%256} 0..256;
@a = keys %a;
print "@a\n";
実行結果
~ $ perl 13_50.pl | wc -w
44
$ perl 13_50.pl
137 89 161 57 1 4 100 17 36 49 121 64 68 144 201 177 65 185 16 9 193 169 129
105 196 132 25 73 249 209 33 233 225 97 41 81 241 164 145 228 217 0 153 113
%a = map{$_=>1} map{$_*$_%256} 0..256;
@a = keys %a;
print "@a\n";
実行結果
~ $ perl 13_50.pl | wc -w
44
$ perl 13_50.pl
137 89 161 57 1 4 100 17 36 49 121 64 68 144 201 177 65 185 16 9 193 169 129
105 196 132 25 73 249 209 33 233 225 97 41 81 241 164 145 228 217 0 153 113
2019/02/09(土) 17:42:33.15ID:1XyVHoA8
2019/02/09(土) 18:17:43.26ID:HLblHrgV
2019/02/09(土) 18:57:59.38ID:luPnpF49
61デフォルトの名無しさん
2019/02/10(日) 06:54:30.87ID:qszHu1wC2019/02/10(日) 12:16:48.56ID:8pY6FeJB
お題
ある数 n とする。
下位から24bit区切りの数を足し合わせてからmod 2^24-1 した数が、元の数nのmod 2^24-1 と一致することを確認しなさい。
ある数 n とする。
下位から24bit区切りの数を足し合わせてからmod 2^24-1 した数が、元の数nのmod 2^24-1 と一致することを確認しなさい。
2019/02/10(日) 12:27:09.74ID:Mq5me4ef
>>62
任意の n に対してある自然数 N, 自然数列{a_k} が存在して
n = Σ_{k = 0}^{N} a_k * 2^(24 * k) と書けるので mod 2^24 - 1 として
n = Σ_{k = 0}^{N} a_k * 1^k
= Σ_{k = 0}^{N} a_k ■
任意の n に対してある自然数 N, 自然数列{a_k} が存在して
n = Σ_{k = 0}^{N} a_k * 2^(24 * k) と書けるので mod 2^24 - 1 として
n = Σ_{k = 0}^{N} a_k * 1^k
= Σ_{k = 0}^{N} a_k ■
2019/02/10(日) 14:27:03.54ID:H2rtpzeI
>>59,60
256はmodじゃなくて&255を取る
確率的には大きい順じゃなくて9,97,17,13,7,5が良いのでは?
9,97,17,13,7,5でmodを取る場合、大きい数からそのままmodを取るのではなく2^48-1でmodを取った数値に対してmod
これで速度どうなる?
256はmodじゃなくて&255を取る
確率的には大きい順じゃなくて9,97,17,13,7,5が良いのでは?
9,97,17,13,7,5でmodを取る場合、大きい数からそのままmodを取るのではなく2^48-1でmodを取った数値に対してmod
これで速度どうなる?
2019/02/10(日) 14:33:21.22ID:JrJlQQ/Q
>>63 それなんと言うプログラム?
2019/02/10(日) 14:41:11.23ID:fWGYOSi7
またこの流れ?
2019/02/10(日) 14:49:22.42ID:NRo2aHHT
mod 255にしたら遅くなるんじゃねーの
0 < n mod 255 < 254
だぞ
0 < n mod 255 < 254
だぞ
2019/02/10(日) 14:49:43.12ID:NRo2aHHT
0 <= n mod 255 <= 254 だった
2019/02/10(日) 15:12:41.36ID:95x0uvij
>>62 python
# n%(2**24-1) を求める
def mod224get(n):
bn=(n.bit_length()+7)//8 #byte長
bb=n.to_bytes(bn,'little')
s= sum([int.from_bytes(bb[i:i+3],'little')
for i in range(0,bn,3) ]) #24bit毎の合計
return s%(2**24-1)
v0=12345678901234567890
#v0=0
n=v0**2
loop = 100000
print('テスト範囲は ',n,'〜',n+loop-1,'loop回数=',loop)
start =time.process_time()
for i in range(n,n+loop):
if mod224get(n) != n%(2**24-1) :print('間違い見っけ',n)
end =time.process_time()
print('全問正解 かかった時間は、',end-start,'秒')
#
―― 結果
テスト範囲は 152415787532388367501905199875019052100 〜 152415787532388367501905199875019152099 loop回数= 100000
全問正解 かかった時間は、 0.2963889999999765 秒
# n%(2**24-1) を求める
def mod224get(n):
bn=(n.bit_length()+7)//8 #byte長
bb=n.to_bytes(bn,'little')
s= sum([int.from_bytes(bb[i:i+3],'little')
for i in range(0,bn,3) ]) #24bit毎の合計
return s%(2**24-1)
v0=12345678901234567890
#v0=0
n=v0**2
loop = 100000
print('テスト範囲は ',n,'〜',n+loop-1,'loop回数=',loop)
start =time.process_time()
for i in range(n,n+loop):
if mod224get(n) != n%(2**24-1) :print('間違い見っけ',n)
end =time.process_time()
print('全問正解 かかった時間は、',end-start,'秒')
#
―― 結果
テスト範囲は 152415787532388367501905199875019052100 〜 152415787532388367501905199875019152099 loop回数= 100000
全問正解 かかった時間は、 0.2963889999999765 秒
2019/02/10(日) 16:02:45.96ID:8pY6FeJB
>>64 アドバイスありがとう。 それは思ったんだけど、現代の言語がそんなところで手抜きはしていないだろうと信じてテストしていなかった。
今、&255 に変えてテストしてみたけど、スピードの差はなかった。 そう言う発想は昔は非常に重要だったけど、今は言語の中で吸収してるみたいだね。
その下のアドバイスに対しては、何故ご自分では試されないのですか?
あまりやるつもりはないのは、mod 2**24-1 と言うのが理解できていないからです。 これで早くなるのなら色々試してみたいんですが、このリストを作るだけでもかなりの時間がかかりめげてます。
今、&255 に変えてテストしてみたけど、スピードの差はなかった。 そう言う発想は昔は非常に重要だったけど、今は言語の中で吸収してるみたいだね。
その下のアドバイスに対しては、何故ご自分では試されないのですか?
あまりやるつもりはないのは、mod 2**24-1 と言うのが理解できていないからです。 これで早くなるのなら色々試してみたいんですが、このリストを作るだけでもかなりの時間がかかりめげてます。
2019/02/10(日) 16:33:26.08ID:H2rtpzeI
>>70
剰余の順番に関しては確率がこんなんだからやで
3 / 5 = 0.600000
4 / 7 = 0.571429
4 / 9 = 0.444444
7 / 13 = 0.538462
9 / 17 = 0.529412
49 / 97 = 0.505155
テーブルは9, 97, 17, 13, 7, 5の物で良いんやで?
多倍長整数の剰余より32bit整数/64bit整数の剰余のほうが計算量が少ないから、
(32bitの場合) 2^24-1で剰余を取ったものに対して9, 17, 13, 7, 5の剰余で平方数かどうかを調べる
(64bitの場合) 2^48-1で剰余を取ったものに対して9, 97, 17, 13, 7, 5の剰余で平方数かどうかを調べる
なしてこんなことができるかってーと、
2^24-1(=16777215)の因数に5, 7, 9, 13, 17が、2^48-1(=281474976710655)の因数に5, 7, 9, 13, 17, 97含まれているからやで
剰余の順番に関しては確率がこんなんだからやで
3 / 5 = 0.600000
4 / 7 = 0.571429
4 / 9 = 0.444444
7 / 13 = 0.538462
9 / 17 = 0.529412
49 / 97 = 0.505155
テーブルは9, 97, 17, 13, 7, 5の物で良いんやで?
多倍長整数の剰余より32bit整数/64bit整数の剰余のほうが計算量が少ないから、
(32bitの場合) 2^24-1で剰余を取ったものに対して9, 17, 13, 7, 5の剰余で平方数かどうかを調べる
(64bitの場合) 2^48-1で剰余を取ったものに対して9, 97, 17, 13, 7, 5の剰余で平方数かどうかを調べる
なしてこんなことができるかってーと、
2^24-1(=16777215)の因数に5, 7, 9, 13, 17が、2^48-1(=281474976710655)の因数に5, 7, 9, 13, 17, 97含まれているからやで
2019/02/10(日) 16:55:58.44ID:8pY6FeJB
>>71 あまり深入りするつもりはないけど、mod 2**24-1 でチェックしたら、
mod 9, 97, 17, 13, 7, 5 でチェックする必要はないと言う事?
ま、数学を解いてるつもりは全くなく、プログラムの練習だからいかに沢山の人が素晴らしいプログラムを見せてくれるかにしか興味はない。
プログラムを書かない人は自分にとってはなんの意味もない。
mod 9, 97, 17, 13, 7, 5 でチェックする必要はないと言う事?
ま、数学を解いてるつもりは全くなく、プログラムの練習だからいかに沢山の人が素晴らしいプログラムを見せてくれるかにしか興味はない。
プログラムを書かない人は自分にとってはなんの意味もない。
2019/02/10(日) 17:11:01.79ID:H2rtpzeI
>>72
ちゃうねん
mod 2**24-1をした数値に対してmod 9, 17, 13, 7, 5 でチェックするねん
もしくはmod 2**48-1をした数値に対してmod 9, 97, 17, 13, 7, 5 でチェックするねん
ちゃうねん
mod 2**24-1をした数値に対してmod 9, 17, 13, 7, 5 でチェックするねん
もしくはmod 2**48-1をした数値に対してmod 9, 97, 17, 13, 7, 5 でチェックするねん
2019/02/10(日) 19:32:43.48ID:/XsfFvRM
>> 73
2重に剰余を取るのはGMPみたく多倍長なら意味があるけど, 32/64bit固定長ならあまり意味はない
複数回剰余を確認する必要があるから多倍長から固定長(32/64bit)にしていて, 更に因数を使えば剰余を求めるための除算の代わりに乗算が使えるから因数の多い2^24 - 1や2^48 - 1を採用してる
2重に剰余を取るのはGMPみたく多倍長なら意味があるけど, 32/64bit固定長ならあまり意味はない
複数回剰余を確認する必要があるから多倍長から固定長(32/64bit)にしていて, 更に因数を使えば剰余を求めるための除算の代わりに乗算が使えるから因数の多い2^24 - 1や2^48 - 1を採用してる
2019/02/11(月) 00:35:41.13ID:8Hdd2FlG
>>62
ガウス少年が見出したように
Σ1,2,…,n-2,n-1=n *(n +1) /2
なので、
n の mod 2^24-1
と
Σ1,2,…,n-2,n-1 =n *(n +1) /2 の mod 2^24-1
が等しいのは自明だと思うけど、
そういう、ちょっとした数学を使わず
Σ1,2,…,n-2,n-1
をloopで和を算出し mod 2^24-1 して比較する
n の mod 2^24-1 と比較する
プログラムを作れという題なんだろうか…
ガウス少年が見出したように
Σ1,2,…,n-2,n-1=n *(n +1) /2
なので、
n の mod 2^24-1
と
Σ1,2,…,n-2,n-1 =n *(n +1) /2 の mod 2^24-1
が等しいのは自明だと思うけど、
そういう、ちょっとした数学を使わず
Σ1,2,…,n-2,n-1
をloopで和を算出し mod 2^24-1 して比較する
n の mod 2^24-1 と比較する
プログラムを作れという題なんだろうか…
2019/02/11(月) 00:37:45.16ID:8Hdd2FlG
2019/02/11(月) 00:57:54.78ID:HnU/OI7o
2019/02/11(月) 01:10:25.11ID:8Hdd2FlG
2019/02/11(月) 01:12:22.73ID:8Hdd2FlG
>>77
ゴメンなんか誤解したかも、よく読む
ゴメンなんか誤解したかも、よく読む
2019/02/11(月) 01:15:22.84ID:f+GXhEiR
ある数nのビット表記方法によって一致する/しないを答えればいいのかな
2019/02/11(月) 01:42:54.29ID:8Hdd2FlG
>>62 Perl5
use bignum (l=>GMP);
use feature say;
sub sum24 {
my $v = $_[0];
if ($v > 0) {
my $d = int($v / 2**24);
my $m = $v % 2**24; # $v - $d * $f6;
$m + sum24($d);
} else {
0;
}
}
$n = 12345678901234567890;
say $n % (2**24 -1);
say sum24($n) % (2**24 -1);
実行結果
~ $ perl 13_62.pl
13189905
13189905
use bignum (l=>GMP);
use feature say;
sub sum24 {
my $v = $_[0];
if ($v > 0) {
my $d = int($v / 2**24);
my $m = $v % 2**24; # $v - $d * $f6;
$m + sum24($d);
} else {
0;
}
}
$n = 12345678901234567890;
say $n % (2**24 -1);
say sum24($n) % (2**24 -1);
実行結果
~ $ perl 13_62.pl
13189905
13189905
2019/02/11(月) 01:47:30.91ID:8Hdd2FlG
2019/02/11(月) 01:53:44.12ID:C0KPLnD/
どうみても自明なんだから確認も糞もないけどな
2019/02/11(月) 01:57:28.89ID:8Hdd2FlG
お題を作ることの難しさだよな…
2019/02/11(月) 02:15:34.48ID:HnU/OI7o
2019/02/11(月) 02:26:56.38ID:8Hdd2FlG
そんな怒るなよ。
暖かくしてぐっすりお休みよ
暖かくしてぐっすりお休みよ
2019/02/11(月) 02:31:44.77ID:HnU/OI7o
しかしここまで複雑な処理をして本当に早くなるのかどうか疑問だけどな。 mod 2**24-1 って結構時間がかかりそうな気がする。
2019/02/11(月) 02:35:42.19ID:ucqIUq+7
>>85
一番能書き垂れてんのお前だろ
クソみたいな御託並べる前に自分のことを考えろっつったろうが
GMPが一体どこで
> n**2%(2**24-1) のリスト
なんか使ってんだ?91で割った場合のテーブルでさえ12byte必要だってのにどうやってそんな巨大なテーブル用意するんだ?
GMPの中身なんか数学の成果の塊だぞ?お前が数学したくないだか出来ないだか知らんがアルゴリズム考えるようなスレでクソみたいなこと喋ってんじゃねぇよ
お前はコードを書いても価値がない
一番能書き垂れてんのお前だろ
クソみたいな御託並べる前に自分のことを考えろっつったろうが
GMPが一体どこで
> n**2%(2**24-1) のリスト
なんか使ってんだ?91で割った場合のテーブルでさえ12byte必要だってのにどうやってそんな巨大なテーブル用意するんだ?
GMPの中身なんか数学の成果の塊だぞ?お前が数学したくないだか出来ないだか知らんがアルゴリズム考えるようなスレでクソみたいなこと喋ってんじゃねぇよ
お前はコードを書いても価値がない
2019/02/11(月) 02:35:55.83ID:8Hdd2FlG
単なるbitmaskで済まない様な場合
あるいは除算して剰余を求めるなら
さんざ研究されていると思うから自力で1から考える前に
先人の業績を知れってことだろ
アバヨ ノシ
あるいは除算して剰余を求めるなら
さんざ研究されていると思うから自力で1から考える前に
先人の業績を知れってことだろ
アバヨ ノシ
2019/02/11(月) 02:36:00.30ID:IhaR3BEX
お題:ポーカーダイス
通常のサイコロを5個振って出た目をポーカーの役になぞってそれぞれの出現確率を求める。
役は、5カード、4カード、ストレート、フルハウス、3カード、2ペア、1ペア、ブタ(ノーペア)
例えば出た目が 1,1,3,1,4 ならスリーカード。2,5,4,6,3 ならストレート。
5カードは6/(6^5)、4カードは(5*5*6)/(6^5)というように数学的に
求めても難しくはないのですが、ここはプログラミングのスレなので
全通り力技でチェックして求めてみてください。
解答例:C言語 https://ideone.com/4X62Am
通常のサイコロを5個振って出た目をポーカーの役になぞってそれぞれの出現確率を求める。
役は、5カード、4カード、ストレート、フルハウス、3カード、2ペア、1ペア、ブタ(ノーペア)
例えば出た目が 1,1,3,1,4 ならスリーカード。2,5,4,6,3 ならストレート。
5カードは6/(6^5)、4カードは(5*5*6)/(6^5)というように数学的に
求めても難しくはないのですが、ここはプログラミングのスレなので
全通り力技でチェックして求めてみてください。
解答例:C言語 https://ideone.com/4X62Am
2019/02/11(月) 03:04:24.50ID:8Hdd2FlG
6^5総当りせよってか…
native compiler系言語で力技か
native compiler系言語で力技か
2019/02/11(月) 03:20:03.89ID:K/18SmCD
Jニキはよ
2019/02/11(月) 03:29:32.98ID:ucqIUq+7
大した数じゃないからズルいことが出来る
https://ideone.com/yEcdPV
https://ideone.com/yEcdPV
2019/02/11(月) 04:00:05.29ID:8Hdd2FlG
お なかなか
2019/02/11(月) 08:16:46.17ID:b3B7Bg4u
python3
https://ideone.com/k6Ea4j
最後の出力部分はpython 3.6以降だと
for k,v in hand.items(): print("{} :\n {} / 7776 ({} %)".format(k,v, round(100*v/7776,2)))
でいけるけど実行環境が3.5なのでやむなく
https://ideone.com/k6Ea4j
最後の出力部分はpython 3.6以降だと
for k,v in hand.items(): print("{} :\n {} / 7776 ({} %)".format(k,v, round(100*v/7776,2)))
でいけるけど実行環境が3.5なのでやむなく
2019/02/11(月) 16:44:30.42ID:xTuBWJbc
なんか数学でもできる力技お題増えてきたな
もっとプログラミングじゃないとできないような良いお題無いんだろうか
もっとプログラミングじゃないとできないような良いお題無いんだろうか
2019/02/11(月) 17:22:02.16ID:7gZS39yo
>>96
そんなの存在しないんじゃない?
そんなの存在しないんじゃない?
2019/02/11(月) 17:28:00.80ID:6aFdKLEP
確率の問題でも特定の疑似乱数と種を使った偏りを求めるとかは数学では難しい
99さまよえる蟻人間 ◆T6xkBnTXz7B0
2019/02/11(月) 17:37:44.75ID:adj8EvAq お題: 日本語文字とカッコ { } とスラッシュ(/)で構成された入力文字列Sが与えられる。{ }で囲まれ、かつ
スラッシュで区切られた部分文字列について、それぞれ場合分けを行って、複数の文字列のリストに展開して改行区切りで出力せよ。
カッコの対応が間違っている場合はERRORを出力せよ。
(例1) {ひまわり/あさがお}は{植物/花}です。
(出力結果)
ひまわりは植物です。
あさがおは植物です。
ひまわりは花です。
あさがおは花です。
スラッシュで区切られた部分文字列について、それぞれ場合分けを行って、複数の文字列のリストに展開して改行区切りで出力せよ。
カッコの対応が間違っている場合はERRORを出力せよ。
(例1) {ひまわり/あさがお}は{植物/花}です。
(出力結果)
ひまわりは植物です。
あさがおは植物です。
ひまわりは花です。
あさがおは花です。
100さまよえる蟻人間 ◆T6xkBnTXz7B0
2019/02/11(月) 17:59:47.38ID:BEdrdhIs なお、展開の順序については問わない。カッコがなくなるまで繰り返し展開せよ。
101さまよえる蟻人間 ◆T6xkBnTXz7B0
2019/02/11(月) 18:20:25.29ID:BEdrdhIs (1) {あ{いう/え}/お{か/き}/く}け{こ}
(2) さ{し/す}せそ{{た/ち}つ/て}と
(2) さ{し/す}せそ{{た/ち}つ/て}と
102デフォルトの名無しさん
2019/02/11(月) 19:00:31.50ID:MkFOBvt9 ネストありかよ、ちょっと面倒だな
103さまよえる蟻人間 ◆T6xkBnTXz7B0
2019/02/11(月) 19:03:55.20ID:BEdrdhIs ヒント: まず、適当な場所でブツ切りにしてノードに分ける。
104デフォルトの名無しさん
2019/02/11(月) 19:20:26.84ID:Q78+FEDq >>101 で、どう言う結果を正解とするの?
105さまよえる蟻人間 ◆T6xkBnTXz7B0
2019/02/11(月) 19:32:34.31ID:BEdrdhIs (1)の答え(※ソート済み)
あいうけこ
あいうけこ
あえけこ
あえけこ
おかけこ
おかけこ
おきけこ
おきけこ
くけこ
くけこ
くけこ
くけこ
あいうけこ
あいうけこ
あえけこ
あえけこ
おかけこ
おかけこ
おきけこ
おきけこ
くけこ
くけこ
くけこ
くけこ
106さまよえる蟻人間 ◆T6xkBnTXz7B0
2019/02/11(月) 19:33:13.25ID:BEdrdhIs (2)の答え(※ソート済み)
さしせそたつと
さしせそちつと
さしせそてと
さしせそてと
さすせそたつと
さすせそちつと
さすせそてと
さすせそてと
さしせそたつと
さしせそちつと
さしせそてと
さしせそてと
さすせそたつと
さすせそちつと
さすせそてと
さすせそてと
107デフォルトの名無しさん
2019/02/11(月) 19:36:37.18ID:MkFOBvt9 これでいいのか?
> (1) {あ{いう/え}/お{か/き}/く}け{こ}
あ いう け こ
あ え け こ
お か け こ
お き け こ
く け こ
> (2) さ{し/す}せそ{{た/ち}つ/て}と
さ し せそ た つ と
さ す せそ た つ と
さ し せそ ち つ と
さ す せそ ち つ と
さ し せそ て と
さ す せそ て と
> (1) {あ{いう/え}/お{か/き}/く}け{こ}
あ いう け こ
あ え け こ
お か け こ
お き け こ
く け こ
> (2) さ{し/す}せそ{{た/ち}つ/て}と
さ し せそ た つ と
さ す せそ た つ と
さ し せそ ち つ と
さ す せそ ち つ と
さ し せそ て と
さ す せそ て と
108デフォルトの名無しさん
2019/02/11(月) 19:37:56.93ID:MkFOBvt9 あれ?
変化しないケースも出力するの?
変化しないケースも出力するの?
110さまよえる蟻人間 ◆T6xkBnTXz7B0
2019/02/11(月) 19:46:19.44ID:adj8EvAq ごめんなさい。間違えました。重複は単一化して下さい。
112さまよえる蟻人間 ◆T6xkBnTXz7B0
2019/02/11(月) 20:21:24.27ID:BEdrdhIs 単純に場所分けを樹木図で考えると重複ができるようだ。すみません。
113デフォルトの名無しさん
2019/02/11(月) 20:48:49.82ID:uHNor3GB お題:Aが真であるならばBが真である ことをプログラムしなさい。
114さまよえる蟻人間 ◆T6xkBnTXz7B0
2019/02/11(月) 21:04:41.39ID:BEdrdhIs アホちゃいまんねん、パーでんねん。
115デフォルトの名無しさん
2019/02/12(火) 00:09:52.31ID:VqanzRzk バカなのか?AとBに因果関係があるわけじゃないし、この世の全てがプログラム言語でマッピングできるわけじゃない、数学徒は帰れ
116デフォルトの名無しさん
2019/02/12(火) 00:29:24.58ID:xM7yD0R2 const A = true;
const B = A === true ? true : false;
console.log(B);
const B = A === true ? true : false;
console.log(B);
117デフォルトの名無しさん
2019/02/12(火) 01:58:51.98ID:ww6cDCbZ118デフォルトの名無しさん
2019/02/12(火) 01:59:52.12ID:/350tEey >>113
!A&&B
!A&&B
119デフォルトの名無しさん
2019/02/12(火) 02:31:32.61ID:qK/pLy4w >>118 python
B = A
B = A
120デフォルトの名無しさん
2019/02/12(火) 02:31:55.89ID:qK/pLy4w >>113 の間違い
121デフォルトの名無しさん
2019/02/12(火) 02:52:28.28ID:jwrsqhME {あ{いう/え}/お{か/き}/く}けこ
あいうけこ
あえけこ
おかけこ
おきけこ
くけこ
さ{し/す}せそ{{た/ち}つ/て}と
さしせそたつと
さしせそちつと
さしせそてと
さすせそたつと
さすせそちつと
さすせそてと
あいうけこ
あえけこ
おかけこ
おきけこ
くけこ
さ{し/す}せそ{{た/ち}つ/て}と
さしせそたつと
さしせそちつと
さしせそてと
さすせそたつと
さすせそちつと
さすせそてと
122デフォルトの名無しさん
2019/02/12(火) 07:13:30.28ID:WW36R8Qd123デフォルトの名無しさん
2019/02/12(火) 09:01:25.53ID:eC1lEXzI >>117も間違い。偽のときは未定義なんだからエラー吐かなきゃ
124デフォルトの名無しさん
2019/02/12(火) 10:15:35.66ID:EWuoyvxz 未定義じゃねえだろアホ
125デフォルトの名無しさん
2019/02/12(火) 10:18:36.24ID:eC1lEXzI126デフォルトの名無しさん
2019/02/12(火) 10:38:38.54ID:/lUdPPCt Aが偽の時はエラー吐かなきゃいけないとかBを偽にしてはいけない
とかいうのは勝手な拡大解釈でしかない
とかいうのは勝手な拡大解釈でしかない
127デフォルトの名無しさん
2019/02/12(火) 10:45:59.39ID:dUnMTtNo 真面目に考えるだけ時間の無駄
128デフォルトの名無しさん
2019/02/12(火) 11:03:54.75ID:8L309PqZ129デフォルトの名無しさん
2019/02/12(火) 11:18:15.32ID:eC1lEXzI >>128
> AならばBでAが偽ならばそれは真だっつーの
えっ、どういうことなの?
それは
AならばB
のとき
AでないならばB
ということ?
BはAに関わらず真ということ?
> AならばBでAが偽ならばそれは真だっつーの
の意味がよくわからん…
> AならばBでAが偽ならばそれは真だっつーの
えっ、どういうことなの?
それは
AならばB
のとき
AでないならばB
ということ?
BはAに関わらず真ということ?
> AならばBでAが偽ならばそれは真だっつーの
の意味がよくわからん…
130デフォルトの名無しさん
2019/02/12(火) 11:29:33.12ID:dUnMTtNo >>129
論理としては A => B (AならばB)は対偶論理 ¬B => ¬A (BでないならばAでない)を成り立たせるために通常 ¬A∨B (AでないかまたはBである) で定義される
つまり A => B という論理式は A が偽であれば B の真偽に依らず真になる
だから何だという話ではある
論理としては A => B (AならばB)は対偶論理 ¬B => ¬A (BでないならばAでない)を成り立たせるために通常 ¬A∨B (AでないかまたはBである) で定義される
つまり A => B という論理式は A が偽であれば B の真偽に依らず真になる
だから何だという話ではある
131デフォルトの名無しさん
2019/02/12(火) 11:30:20.12ID:7Ldk0kbC132デフォルトの名無しさん
2019/02/12(火) 11:36:13.30ID:puzbyhsI AならばBと
Aが真ならばBが真
とは違うだろ
Aが真ならばBが真
とは違うだろ
133デフォルトの名無しさん
2019/02/12(火) 11:43:29.03ID:puzbyhsI 「AならばB」
と言う命題は
「Aが真でBが真である
Aが偽であればBは真である」という命題の
上の文の3行目のはじめの部分をプログラムしろということだぞ
と言う命題は
「Aが真でBが真である
Aが偽であればBは真である」という命題の
上の文の3行目のはじめの部分をプログラムしろということだぞ
134デフォルトの名無しさん
2019/02/12(火) 11:45:33.01ID:puzbyhsI 4行目間違えた
135デフォルトの名無しさん
2019/02/12(火) 12:10:48.64ID:2r3VUiS2 A: 自然数 : 1,2,3,・・・・・
B: 整数 : ・・・・・ , -2,-1,0,1,2,3,・・・・・
AならばBである
AでなければBでもない
BでなければAでもない
B: 整数 : ・・・・・ , -2,-1,0,1,2,3,・・・・・
AならばBである
AでなければBでもない
BでなければAでもない
136デフォルトの名無しさん
2019/02/12(火) 12:11:50.75ID:/o8EBvgR■ このスレッドは過去ログ倉庫に格納されています
ニュース
- 外務省局長は無言で厳しい表情…日中の高官協議終了か 高市首相“台湾”発言で中国が強硬対応 発言撤回求めたか… [BFU★]
- 中国国営メディア「沖縄は日本ではない」… ★6 [BFU★]
- 政府、株式の配当など金融所得を高齢者の医療保険料や窓口負担に反映する方針を固めた [バイト歴50年★]
- ナイツ塙が指摘のローソンコーヒーカップ、ロゴ「L」で誤解生みデザイン変更へ 在庫使い切る3か月後にリニューアル [muffin★]
- バービー、 台湾有事の発言の波紋で「たまったもんじゃない」「高市さんに真意は聞きたい」「国民に向けて説明してほしい」 [muffin★]
- 【速報】 高市政権、「日本版DOGE」を立ち上げ 米国で歳出削減をした「政府効率化省(DOGE)」になぞらえたもの [お断り★]
- 高市早苗、岸田政権(当時)に「台湾有事は日本の有事か」という質問をしていた [175344491]
- 日銀植田総裁「高市早苗と為替について議論、政府と連携して注視することで一致した」 [888298477]
- 【悲報】中国→日本行きの航空チケット、高市有事の影響で50万人分がキャンセルされる [834922174]
- 岸田文雄「国民は不安よな 岸田、動きます」 非核三原則堅持の立場示す [175344491]
- 【悲報】早速高市首相のせいで全国の民泊でキャンセルラッシュwwwwwwwwwwww 経営者も嘆き「こんな事は初めてだ…」😲 [871926377]
- 【悲報】日銀植田「利上げはデータ次第」と高市にぴしゃり [733893279]
