[表示 : 全て 最新50 1-99 101- 201- 301- 401- 501- 601- 701- 801- 901- 1001- 2ch.scのread.cgiへ]
Update time : 05/07 12:57 / Filesize : 360 KB / Number-of Response : 1043
[このスレッドの書き込みを削除する]
[+板 最近立ったスレ&熱いスレ一覧 : +板 最近立ったスレ/記者別一覧] [類似スレッド一覧]


↑キャッシュ検索、類似スレ動作を修正しました、ご迷惑をお掛けしました

プログラミングのお題スレ Part12



1 名前:デフォルトの名無しさん mailto:sage [2018/09/28(金) 10:09:07.13 ID:phwOkayR.net]
プログラミングのお題スレです。

【出題と回答例】
1 名前:デフォルトの名無しさん
  お題:お題本文

2 名前:デフォルトの名無しさん
  >>1 使用言語
  回答本文

【ソースコードが長くなったら】 (オンラインでコードを実行できる)
https://ideone.com/
codepad.org/
compileonline.com/
rextester.com/runcode
https://runnable.com/
https://code.hackerearth.com/
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/

855 名前:デフォルトの名無しさん mailto:sage [2018/12/13(木) 09:40:48.46 ID:wBrOfppZ.net]
>>822
c
https://ideone.com/YkmUdB

c++
https://ideone.com/WXXvlU

856 名前:デフォルトの名無しさん mailto:sage [2018/12/15(土) 18:26:25.25 ID:RxHR1YVb.net]
お題:標準入力から読み取った行を出力せよ
ただし42を読んだ場合は出力せずに終了する
https://ideone.com/vdWddB

857 名前:デフォルトの名無しさん mailto:sage [2018/12/15(土) 18:51:24.63 ID:Hf86DN+J.net]
>>829
ruby -pe'~/^42$/&&exit'

858 名前:デフォルトの名無しさん [2018/12/15(土) 20:46:54.11 ID:Du8iBFv1.net]
>>829 Common Lisp
https://ideone.com/QQzEQh

read-lineの第三引数にnil以外を使ったの初めてだ

859 名前:デフォルトの名無しさん mailto:sage [2018/12/15(土) 23:34:07.15 ID:2jxK776v.net]
>>829
bash
https://paiza.io/projects/gIZdCMpBEoNvFnnjVvGGtA

860 名前:デフォルトの名無しさん mailto:sage [2018/12/15(土) 23:59:19.52 ID:QTag+jm2.net]
>>829 PowerShell
while (($var = (Read-Host)) -ne 42) { $var }

861 名前:デフォルトの名無しさん mailto:sage [2018/12/16(日) 03:36:20.47 ID:X+FaDx+z.net]
>>829
python 3
https://ideone.com/it2Axq

862 名前:デフォルトの名無しさん mailto:sage [2018/12/16(日) 13:54:15.20 ID:ob8ozoeg.net]
[お題] 来年と素数
今年も残りわずか、来年は2019年で平成31年。

1)"2019"の省略形の"19"について。
 素数の和で19を作る、すべての素数配列を列挙せよ(できれば辞書順で)。
 同じ素数を何個使ってもよいが、同じ素数同士は区別しない。
 ・例えば対象が"11"だと以下の6つ
 {2, 2, 2, 2, 3}
 {2, 2, 2, 5}
 {2, 2, 7}
 {2, 3, 3, 3}
 {3, 3, 5}
 {11}


 以下 2)3)4)は種類計のみ答える(明細は多いので略)。
2)来年の初まりは平成31年で、"31"について。
 素数の和で31を作る、その種類はいくつか。条件は1)と同様。

3)素数の和で2019を作る、その種類はいくつか。条件は1)と同様。

4)2019と31を続けた数 201931(=2019*100+31)について。
 素数の和で201931を作る、その種類はいくつか。
 但し、使用していい素数は31以下の素数かつ、
 同じ素数は最大2019個までしか使えない。同じ素数は区別しない。

 ※ 3)4)は64bit整数を超えるので、下10桁だけの回答も可。

863 名前: mailto:sage [2018/12/16(日) 14:12:11.90 ID:IjCemHtZ.net]
>>835
年忘れ課題の時期になったんですね…



864 名前:デフォルトの名無しさん mailto:sage [2018/12/16(日) 19:04:23.02 ID:P931WLXH.net]
>>835 Java
https://ideone.com/LgkDWF
あってんのかなー

865 名前:830 mailto:sage [2018/12/16(日) 19:54:33.15 ID:GrZg6kve.net]
>>837

2)の31からして自分のと違う。
自分のでは111だった。

javaのやついじって、112個の明細を出してみたら、[2,2,27] があった。
27は素数じゃないよね。いじり方が悪い?

866 名前:デフォルトの名無しさん mailto:sage [2018/12/16(日) 19:56:19.00 ID:pjubjjb0.net]
>>838
自分も31は111だったわ

867 名前:デフォルトの名無しさん mailto:sage [2018/12/16(日) 20:37:34.86 ID:VvrWecHB.net]
>>835 Ruby 1)は省略

require 'prime'

def fuge(n, cand = Prime.to_a(n), h = {}, succ = 0)
c0 = cand[0]
return h[[n, c0]] if h[[n, c0]]
return 0 if n == 0 || !c0 || n < c0
return 1 if n == c0
x = succ == 2019 ? 0 : fuge(n - c0, cand, h, succ + 1)
h[[n, c0]] = fuge(n, cand[1..-1], h) + x
end

p fuge(19) # => 23
p fuge(31) # => 111
p fuge(2019) # => 576202207044176168646563
p fuge(201931, Prime.to_a(31)) # => 4021686887140718864271667825968903

868 名前:デフォルトの名無しさん mailto:sage [2018/12/16(日) 20:50:18.98 ID:P931WLXH.net]
>>838
まじだ・・・修正した
犯人は118行目だ!orz

869 名前:830 mailto:sage [2018/12/17(月) 20:45:38.05 ID:4p2KDXiR.net]
>>837 836
 当方の回答と一致しました。

>>840
4) 201931 だけ異なる(自分のがあっているとは限らないが……)。
※再帰の動きがよくわからない

870 名前:830 mailto:sage [2018/12/20(木) 20:26:37.78 ID:W0v1JICZ.net]
>>835 python https://ideone.com/OFR7bn

 間が空いたので、出題者コメント

 1)問題文説明用に書いてみて、問題にもできると思った。

 2)3)は"オイラー 31 DP"でググってください。
 違いはコインの額面が素数に変わったくらいと、おおきさ。
 なかには一つ一つ数えて、本問の2)しか溶けない回答があるので見分けて。

 4)については、ソース上三種類書いて考察しています。
  ("201931"という変な数値は、3重ループ解を落とすのが目的だった)

871 名前:デフォルトの名無しさん mailto:sage [2018/12/21(金) 08:20:43.27 ID:hMEdBbLv.net]
お題:
真理値表から2入力NANDで最小ゲート数回路を作れ。実用上最も重要な問題群

872 名前:デフォルトの名無しさん mailto:sage [2018/12/21(金) 08:38:52.94 ID:choQhZIj.net]
遅延時間とかファンアウトとかは考慮不要ってことでいいのかな

873 名前:デフォルトの名無しさん [2018/12/22(土) 11:29:36.65 ID:7UfbjIEE.net]
参考にどうぞ
https://eng.kice.tokyo/logic/logic_ex3/



874 名前: mailto:sage [2018/12/22(土) 16:28:53.58 ID:vmp1HvU+.net]
>>846
昔カルノー図がかける範囲(3入力?4入力くらいだったかな?)で簡約化ツールを作ったことがあります
当時は書き捨て上等、でやっていたのですが、今は困ってしまっています、当時のプログラム記述能力は当時は今後もずっと保存されるものと当時は考えていたのですが、それは完全な誤りでした…

875 名前:デフォルトの名無しさん mailto:sage [2018/12/23(日) 16:48:28.85 ID:V9Fp8lZV.net]
>>846
俺全然わかってない人なんだけど、
そのURLの一番下の、入力部分のNOTゲートを無視しない場合、
9個のNANDゲートで構成できるって書いてあるけど、
NOTゲートの部分を2個減らして7個で作るとかはあかんのん?

876 名前:デフォルトの名無しさん mailto:sage [2018/12/23(日) 16:56:13.54 ID:jMuRFGAa.net]
半導体製造のコストがある
NANDゲートは作り易い、というか造りが簡素で単純
NOTゲートも作り易い

なので厳密には製造コストも考えないといけない

877 名前:デフォルトの名無しさん mailto:sage [2018/12/23(日) 17:29:24.88 ID:gPP7+XoH.net]
>>848
ホントだ、AとCのNOTを一つずつ減らせるね
というかこれが最小になる保証もないよね

878 名前:デフォルトの名無しさん mailto:sage [2018/12/23(日) 18:34:59.46 ID:m1FIlPHW.net]
>>846
カルノー図からの簡略結果が最終形ではない。
最終形は
Z=~ABC+A~C+~BC

879 名前:デフォルトの名無しさん mailto:sage [2018/12/23(日) 21:11:36.62 ID:9AXglOVX.net]
>>851
それA=0, B=0, C=0の結果から間違ってるんだが

880 名前:デフォルトの名無しさん mailto:sage [2018/12/23(日) 23:39:28.74 ID:DR7zt9Cz.net]
>>852
うむ

最終形は
Z=~ABC+A~C+~B~C
ですね

2入力NAND回路は
Step 9, Solution 1

0 : -1, -1: 0, 0, 0, 0, 1, 1, 1, 1

1 : -1, -1: 0, 0, 1, 1, 0, 0, 1, 1

2 : -1, -1: 0, 1, 0, 1, 0, 1, 0, 1

3 : 0, 0: 1, 1, 1, 1, 0, 0, 0, 0

4 : 3, 1: 1, 1, 0, 0, 1, 1, 1, 1

5 : 4, 2: 1, 0, 1, 1, 1, 0, 1, 0

6 : 5, 2: 1, 1, 1, 0, 1, 1, 1, 1

7 : 5, 4: 0, 1, 1, 1, 0, 1, 0, 1

8 : 7, 6: 1, 0, 0, 1, 1, 0, 1, 0

ともう一つでした。

881 名前:デフォルトの名無しさん mailto:sage [2018/12/24(月) 08:22:13.10 ID:bL5ZoHRn.net]
>>853
> Z=~ABC+A~C+~B~C
それA=0, B=1, C=0の結果が違うぞ

882 名前:デフォルトの名無しさん mailto:sage [2018/12/24(月) 11:32:09.70 ID:7aEQUcJ6.net]
>>854
~ABC → C=0なので0
A~C → A=0なので0
~B~C → B=1なので0
結果0であってるんじゃないのん?

883 名前:デフォルトの名無しさん mailto:sage [2018/12/24(月) 13:42:11.05 ID:bL5ZoHRn.net]
>>855
すまん、合ってるわ
ちょっとボケてたな



884 名前:デフォルトの名無しさん mailto:sage [2018/12/26(水) 16:33:56.35 ID:GYS6kcBm.net]
>>853
wxMaximaでお絵描きした。
/* [wxMaxima: input start ] */
load (graphs)$
gg() := draw_graph(
net,
show_weight=true,
vertex_size=3,
show_id=true,
show_vertices=[0,1,2,8],
show_vertex_type=filled_square,
head_length=0.2,
head_angle=5,
edge_color="dark-green",
text_color=blue
);
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
Solution[1]:[[[ 0, 3],2],[[ 3, 4],1],[[ 1, 4],1],
[[ 4, 5],1],[[ 2, 5],1],[[ 5, 6],1],[[ 2, 6],1],
[[ 5, 7],1],[[ 4, 7],1],[[ 7, 8],1],[[ 6, 8],1]]$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
aa:1$ for g:aa thru 1 do (
net : create_graph(9,Solution[g],directed=true),
print("Solution",aa,Solution[g]),aa:aa+1,
gg());
/* [wxMaxima: input end ] */

885 名前:デフォルトの名無しさん mailto:sage [2019/01/06(日) 20:37:20.68 ID:uX30aAuK.net]
エモい歌うまバーチャルユーチューバー
https://script.google.com/macros/s/AKfycbxuAt1ejx43_zMRfK7Xplwgwn7cXDQAdZTJkDdLrVCLo8SIdhI/exec

今、見られている話題のVtuberを知りたいなら、ハルカス式Vtuberランキングを見ましょう
【ハルカス式Vtuberランキング】※毎週・毎月自動更新
https://script.google.com/macros/s/AKfycbzZO93mJcZ5M_dPmIZOzuVQNr4mcGwQQdiT7GiLld_1Uews7uE/exec

年末年始の一週間12/27〜01/02までのハルカス式Vtuberランキング
https://docs.google.com/spreadsheets/d/1jdVQ2r5yTrAykMRF1nuX6_5DXpK0167D4JcVCQZTJdM/edit#gid=157464142

年末の一ヶ月12/1〜12/31までのハルカス式Vtuberランキング
https://docs.google.com/spreadsheets/d/1Yc6HD_h4lfurme9b3T0ZXBf2-iBpFSmMGg-1bdgSe0E/edit#gid=1971024984

第4回バーチャルYouTuber人気投票募集中(全304名、2018年11月23日〜)

あなたの好きなVtuberは?5問(配信内容、トーク、歌うま、ゲーム実況、キャラデザ)
https://script.google.com/macros/s/AKfycbwSwNBm8qYD4_kZN2uJLeqRIP8Mwpbo3YDTUEpaSMU02BDAR3jh/exec

・Googleフォームは標準仕様では不正投票が可能な為、改造して対策を行いました。
・2018/11/23時点でチャンネル登録数10,000人以上のuserlocalに登録されたVtuberを対象としています。
・一つのチャンネルでVtuberが複数人いる場合は、それぞれ分けています。
・このフォームに投票するにはグーグルアカウントでログインする必要はありません。
・並び順は前回の得票数の高い順と新人は登録者数の多い順に並んでいます。
・今回も1,000人の方が投票するまで継続します。

なお投票結果はスプレッドシートにリアルタイムで表示されるようにしました。
また【概要】と【詳細】でシートを分けております。
https://docs.google.com/spreadsheets/d/1udB81Vnia9CoIgPaan3a5XFeaCfcnBXkRoPTuxvRpAc/

grげrげ

886 名前:デフォルトの名無しさん [2019/01/18(金) 22:48:43.15 ID:sD1souQd.net]
あげ

887 名前:デフォルトの名無しさん [2019/01/19(土) 05:22:33.88 ID:M8Q3zGyy.net]
過疎ってるので別スレのこれをお題にする。
なお、クラスではなく関数でもサブルーチンでも良い。

https://mevius.5ch.net/test/read.cgi/tech/1544839627/332

332 デフォルトの名無しさん (ワッチョイ f11f-t+p0) sage 2019/01/18(金) 22:57:58.37 ID:1ZiD/8V00
文字列から数字をtmpWとtmpHに抜き出すクラスを作りたいのですがどうしたらいいですか?
演算子記号は+−もある、数字の桁数は不明で小数点もあり、数字が入るときは小文字一字始まり、引数はw or hで抜き出す戻り値を示す。


対象文字
@w900×HH
AWW×h1000
BWW×HH
Cw900/h100
結果 引数がwなら左  hなら右
@tmpW=900 tmpH=WW
AtmpW=WW tmpH=1000
BtmpW=WW tmpH=WW
CtmpW=900 tmpH=1000

888 名前:デフォルトの名無しさん mailto:sage [2019/01/19(土) 08:17:19.40 ID:9NMu+pMT.net]
お題スレじゃなくて、あなたの宿題無料でやりますよになってる

889 名前:デフォルトの名無しさん [2019/01/19(土) 09:54:02.71 ID:9IabhBF0.net]
>>860 Ruby
f = -> str {str.scan(/([A-Z]{2}|[a-z])((?:\d+\.)?\d+)?/).each_with_object({}){|(m, num), h| h["tmp#{m[0].upcase}"] = !num ? ' ' : num.match?(/\./) ? num.to_f : num.to_i}}

%w[
  w900×HH WW×h1000 WW×HH w900/h100 w3.14/h2.72
].each{|e| h = f[e]; puts "#{e} -> " + h.map{|v| '%s = %p' % v}.join(', ')}
# =>
w900×HH -> tmpW = 900, tmpH = " "
WW×h1000 -> tmpW = " ", tmpH = 1000
WW×HH -> tmpW = " ", tmpH = " "
w900/h100 -> tmpW = 900, tmpH = 100
w3.14/h2.72 -> tmpW = 3.14, tmpH = 2.72

890 名前:デフォルトの名無しさん mailto:sage [2019/01/19(土) 12:58:12.91 ID:XYN5JTgF.net]
> CtmpW=900 tmpH=1000

891 名前:デフォルトの名無しさん [2019/01/19(土) 15:20:58.20 ID:M8Q3zGyy.net]
>>861
お題をやるかどうかは自由だし使う言語も自由なのでここは宿題向きスレではない。それに宿題用のスレは別にある。

892 名前:デフォルトの名無しさん [2019/01/20(日) 17:00:10.01 ID:XepLmY3z.net]
お題
チェス盤をかく

□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□

893 名前:デフォルトの名無しさん mailto:sage [2019/01/20(日) 17:28:18.15 ID:bv4HmLR/.net]
for i in range(4):
print("□■□■□■□■")
print("■□■□■□■□")



894 名前:デフォルトの名無しさん mailto:sage [2019/01/20(日) 18:45:01.40 ID:vla32kuj.net]
console.log(`
□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□
`);

895 名前:デフォルトの名無しさん mailto:sage [2019/01/20(日) 19:21:05.18 ID:QLQt0ffh.net]
>>865
85kg 170cmの俺がJavascriptで回答。
[85, 170]
.map(d => [...d.toString(2).padStart(8,0)]
.map(b => +b ? '□' : '■')
.join`` + '\n')
.join``
.repeat(4)

896 名前:デフォルトの名無しさん mailto:sage [2019/01/20(日) 19:39:43.28 ID:CUY2SNg1.net]
>>865 Ruby

puts 72.times.map{|i| i % 9 / 8 * 10 + 15 / (i % 9 + 8) * (i & 1 ^ 9633)}.pack('U*')
# =>
□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□

897 名前:デフォルトの名無しさん mailto:sage [2019/01/20(日) 19:45:27.08 ID:z2xvgkTe.net]
>>868
>>869
お前ら天才だな
普段何やってんの

898 名前:デフォルトの名無しさん mailto:sage [2019/01/20(日) 19:58:44.02 ID:vla32kuj.net]
>>870
ヒマ人どもの公開オナニーに騙されるな。
ベンチ取ってみれば分かるけど
>>868より>>867のほうが処理速度も速く、CPU、メモリ消費も勝ってる。

899 名前:デフォルトの名無しさん mailto:sage [2019/01/20(日) 22:22:45.08 ID:n75q6Hj+.net]
解法や工夫は色々あっていい

900 名前:デフォルトの名無しさん mailto:sage [2019/01/20(日) 23:40:01.19 ID:3oCDYmJz.net]
>>871
可読性も最強だしなw

901 名前:デフォルトの名無しさん [2019/01/21(月) 00:16:11.01 ID:evkqklOK.net]
前々からつくろうとしていまだ決定版ができてないのがある。
コマンド fileA fileBで、重複行を除いて連結と、
A-Bで重複行を除くやつ。入力が数百Mだときつい。

902 名前:デフォルトの名無しさん [2019/01/21(月) 00:20:48.31 ID:q5TMCwho.net]
きついのは口臭だけにしとけよ

903 名前:デフォルトの名無しさん mailto:sage [2019/01/21(月) 04:33:49.89 ID:ZCrWhalc.net]
>>874
Linux とかでよく使われるシェル(bash とか)ならこんな感じだろうな。

( uniq fileA ; uniq fileB ) > fileC

シェルスクリプトにするならこんな内容のファイルを作ってそれの実行パーミッション立てておくだけ。

#!/bin/sh
uniq $1
uniq $2

ようするに uniq コマンド使ってるだけ。

Windows でも似たようなものはあるかも知れないが WSL 入れて Ubuntu とか動かしちゃった方が何も考えずに Linux 風にできるので楽かも知れない。
/mnt/c/ にCドライブマウントされてるからファイルのやり取りを考える必要もない。
改行が LF なのとデフォルトのエンコーディングが UTF-8 なのは考慮する必要があるが。



904 名前:デフォルトの名無しさん [2019/01/21(月) 07:09:30.84 ID:evkqklOK.net]
そんな機能あったっけ、と確認したみたけど意図する動作しない、連続しない重複も対象



Linuxコマンド集 【 uniq 】 ソート済みのファイルから重複した行を削除する
ファイルで重複している行を削除する。ただし,ファイルはソートしてある必要がある
https://tech.nikkeibp.co.j

905 名前:p/it/article/COLUMN/20060227/230918/


連続していない離れた重複行も削除したければ、sortコマンドコマンドで予めソートする必要がある。その代わり順番は保存されない。
https://hydrocul.github.io/wiki/commands/uniq.html
[]
[ここ壊れてます]

906 名前:デフォルトの名無しさん mailto:sage [2019/01/21(月) 08:10:55.04 ID:rd0KMO5T.net]
要するにAの末尾とBの先頭で一致するとこだけを削除したいってこと?

907 名前:デフォルトの名無しさん [2019/01/21(月) 08:19:59.40 ID:evkqklOK.net]
2つあって。
一つは連結した単一ファイルで、後方にある重複行を取り除く。
もう一つはBに存在する行を取り除いたAの重複行を取り除く。
それぞれuniq拡張(A+B) 、uniq拡張(A-B)のような。

908 名前:デフォルトの名無しさん [2019/01/21(月) 09:47:36.96 ID:evkqklOK.net]
スクリプトやめてC++でやってみたら
ロードと書き込みの速度とほぼ同程度でできるようだ。

909 名前:デフォルトの名無しさん [2019/01/21(月) 11:56:42.83 ID:ub5sAEWO.net]
>>865 J
}:("1) 8 9 $ '_#'

_#_#_#_#
#_#_#_#_
_#_#_#_#
#_#_#_#_
_#_#_#_#
#_#_#_#_
_#_#_#_#
#_#_#_#_

910 名前:デフォルトの名無しさん mailto:sage [2019/01/21(月) 13:07:07.09 ID:jV3tC4jb.net]
JってUnicode使えないの?

911 名前:デフォルトの名無しさん [2019/01/21(月) 13:14:23.06 ID:ub5sAEWO.net]
使えるんだけど私の方が勉強不足で。

912 名前:デフォルトの名無しさん [2019/01/21(月) 14:03:05.80 ID:chW+Guj1.net]
>>865
Perl
https://paiza.io/projects/06iuVCc-HeWgTb7j2TrTMQ

913 名前:デフォルトの名無しさん mailto:sage [2019/01/21(月) 14:14:22.79 ID:P73TaUU+.net]
>>879
A-Bの部分は
grep -x -v -f B A
でできるけど順序を保存したuniqはコマンドだけではできないのかな



914 名前:デフォルトの名無しさん [2019/01/21(月) 18:14:39.49 ID:chW+Guj1.net]
そういうコマンドを別途作るか、全部作っちゃうかかな。

915 名前:デフォルトの名無しさん mailto:sage [2019/01/21(月) 21:03:38.83 ID:mEBdwP2z.net]
Ruby なら、ハッシュで重複行を排除できる。
ARGF は、a.txt b.txt (ARGV)を連結した、仮想ファイル

ruby script.rb a.txt b.txt

以下は、script.rb の内容

hash = { }

ARGF.each_line { |line| hash[ line ] = true } # true には、特に意味がない

puts hash.keys

-- a.txt

a

b

-- b.txt
x

a


-- 出力

a
b
x


916 名前:デフォルトの名無しさん [2019/01/21(月) 21:49:42.62 ID:QH90tNLy.net]
速度出るのできた。CRCと文字列長のハッシュしかみてなく重複らしいのはduplicate.txtへ。


標準出力
deldup.exe fileA fileB ・・・ はマージして重複削除
deldup.exe /D fileA fileB ・・・ はA -B -・・・の重複削除
https://ideone.com/ABJNAy

917 名前:デフォルトの名無しさん mailto:sage [2019/01/21(月) 22:56:33.96 ID:uToUSKRf.net]
>>865 octave
https://ideone.com/KyX1rR

918 名前:デフォルトの名無しさん mailto:sage [2019/01/21(月) 23:04:52.30 ID:PFiT+1XJ.net]
>>226
Javaはunsignedの概念がない

919 名前:デフォルトの名無しさん mailto:sage [2019/01/22(火) 02:53:38.24 ID:t1xtxA/a.net]
>>879
後方にある重複業削除なら awk で連想配列(要するにハッシュ)使ってこんな風にすればできるな。

awk '{if(!n[$0]){n[$0]=1;print}}' fileA fileB

ただし全行をメモリ上に置くことになるのでファイルがでかいとまともに動かない環境があるかも。

超大きいファイルの場合は行ごとに MD5 や SHA1 等のハッシュ値計算してそれだけ保管しておいて比較した方が良いかもね。

920 名前:デフォルトの名無しさん [2019/01/22(火) 10:22:28.48 ID:f5q/almD.net]
>>865 Lua
for i=1,72 do
if i%9==0 then
print("")
else
io.write(({"■","□"})[i%2+1])
end
end

921 名前:デフォルトの名無しさん mailto:sage [2019/01/22(火) 11:22:34.61 ID:ez/DcBYq.net]
普通に行ごとにソートして重複を除去すればいいだけやん

922 名前:デフォルトの名無しさん mailto:sage [2019/01/22(火) 12:39:35.18 ID:bOf9t ]
[ここ壊れてます]

923 名前:fZi.net mailto: >>865 javascript
[...function*(){for(var i=0;i<64;i++)yield (i&7)==7?'\n':'■□'[i&8?i%2:+!(i%2)]}()].join``
[]
[ここ壊れてます]



924 名前:デフォルトの名無しさん [2019/01/22(火) 12:48:37.38 ID:J1mDHp9Y.net]
>>893
順序変えなくないんだって。

925 名前:デフォルトの名無しさん [2019/01/22(火) 18:32:37.01 ID:f5q/almD.net]
>>865 Ruby

puts (["□■"*4,"■□"*4]*4)

926 名前:デフォルトの名無しさん mailto:sage [2019/01/22(火) 18:48:04.82 ID:bOf9tfZi.net]
>>881
J分からない人向け解説キボンヌ

927 名前:デフォルトの名無しさん mailto:sage [2019/01/22(火) 22:15:16.89 ID:4TOwjU0o.net]
>>865 javascript
['□■', '■□'].map(s => s.repeat(4) + '\n').join``.repeat(4)

928 名前:デフォルトの名無しさん [2019/01/24(木) 21:04:14.42 ID:HE3o+0uB.net]
>>865 Lua
print (((("_#"):rep(36)):gsub("(........).","%1\n")))

929 名前:デフォルトの名無しさん mailto:sage [2019/01/24(木) 22:20:43.56 ID:HLNzgCZV.net]
>>865 javascript
Array(72).fill().map((v, i)=>i%9?'■□'[i%2]:'\n').join``

>>894はボード横幅が7になってる間違いでしたごめんなさい。

930 名前:デフォルトの名無しさん [2019/01/25(金) 20:09:30.47 ID:s1700d1K.net]
お題
n個からr個選ぶ組み合わせの個数を求める

931 名前:デフォルトの名無しさん mailto:sage [2019/01/25(金) 22:06:56.43 ID:7NAKvkXa.net]
数学の問題が叩かれるのはここのアホどもには難しいからなんだ!
決してスレ違いだからじゃないんだ!
だから高校1年の教科書の定義引き写すだけの簡単な問題なら叩かれないんだ!

932 名前:デフォルトの名無しさん mailto:sage [2019/01/25(金) 22:30:26.19 ID:oIz/D2lY.net]
教科書捨てたから定義ググりました。
const p = (n, r) => (n < 2 || r < 1) ? 1 : n * p(n - 1, r - 1);
const c = (n, r) => p(n, r) / p(r);
c(1200, 100);
//=> 1.2734481790871909e+148

chromeでBigInt実装されてるのでせっかくだからと思ってpを
const p = (n, r) => (n < 2n || r < 1n) ? 1n : n * p(n - 1n, r - 1n);
に改良しようとしたらエラーでる…
なんでか分かる方いたら教えてください。

933 名前:デフォルトの名無しさん mailto:sage [2019/01/26(土) 03:39:53.55 ID:mo79YSSI.net]
>>901 Ruby

c = -> (n, k, m = []) {
  m[n] = [1] unless m[n]
  k = n - k if n < k * 2
  return m[n][k] if m[n][k]
  m[n][k] = c[n - 1, k - 1, m] + c[n - 1, k, m]
}

p c[5, 2] # => 10
p c[100, 50] # => 100891344545564193334812497256

>>903
undefined - 1n でエラーが起きているので
const c = (n, r) => p(n, r) / p(r, r);



934 名前:デフォルトの名無しさん [2019/01/26(土) 03:52:04.18 ID:NFcmjJe3.net]
>>901
効率はともかく、関数プログラミングって本でこれ見た時は感動したな。(言語は数学とSMLを合わせた仮想言語だったが)

Haskell

c n 0 = 1
c n m |n == m = 1
c n m = c (n - 1) m + c (n - 1) (m - 1)

実行結果
5 `c` 3
>10

935 名前:デフォルトの名無しさん [2019/01/26(土) 04:12:24.35 ID:NFcmjJe3.net]
高校の教科書版

Haskell

c' n r = product [1..n] `div` (product [1..(n - r)] * product [1..r])

解説
n C r
= n P r / r!
= n! / ((n - r) ! * r!) (n P r = n ! / (n - r)!に分解)

936 名前:デフォルトの名無しさん mailto:sage [2019/01/26(土) 04:32:41.56 ID:Ve68vOks.net]
>>904
ありがとう!


937 名前:ずかしい… []
[ここ壊れてます]

938 名前:デフォルトの名無しさん mailto:sage [2019/01/26(土) 05:11:41.67 ID:iYrqGgfh.net]
>>901 C++
inline uint64_t gcd(uint64_t a, uint64_t b) {
uint64_t r;
while ((r = a % b) != 0) {
a = b; b = r;
}
return b;
}

uint64_t binomial(uint64_t n, uint64_t m) {
if (n < m)
return 0;
m = (m > n / 2) ? n - m : m;
uint64_t result = 1;
for (uint64_t i = 1; i <= m; i++, n--) {
uint64_t d = gcd(result, i);
result /= d;
result *= n / (i / d);
}
return result;
}

939 名前:デフォルトの名無しさん [2019/01/26(土) 10:34:54.19 ID:ls0y1m5K.net]
>>901 Lua
function cc(n,r)
local a = 1
for i = 1, r do
a=a*(n-i+1)/i
end
return a
end
print(cc(20,10))

184756.0

940 名前:デフォルトの名無しさん mailto:sage [2019/01/26(土) 11:22:21.12 ID:tq46h9Ux.net]
お題
複数の線分が平面上で重なりがあるかないかしらべる。

941 名前:デフォルトの名無しさん mailto:sage [2019/01/26(土) 12:02:35.05 ID:Owe/QLUY.net]
>>865
@Mathematica

Table[If[EvenQ[i + j], 0, 1], {i, 8}, {j, 8}] //
 ArrayPlot[#, Mesh -> True] &

942 名前:デフォルトの名無しさん [2019/01/26(土) 13:54:20.56 ID:+oP10x+/.net]
>>910
適当に連立方程式とくだけじゃなかったです?
数学出来ないんですけど。

943 名前:デフォルトの名無しさん [2019/01/26(土) 14:10:31.09 ID:t2YAFdPa.net]
人狼AI と 人狼ゲーム



944 名前:デフォルトの名無しさん mailto:sage [2019/01/26(土) 14:39:24.25 ID:YRte2Dpq.net]
>>865 Perl5
print ((qw{□■}x4, "\n", qw{■□}x4, "\n")x4);

実行結果
$ perl 12_860.pl
□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□

945 名前:デフォルトの名無しさん mailto:sage [2019/01/26(土) 15:04:32.49 ID:vn81VQ6e.net]
>>865
Pharo Smalltalk
String cr join: ((1 to: 4) collect: [:i | {'□■' repeat: 4 . '■□' repeat: 4} ]) flattened
実行結果

□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□
□■□■□■□■
■□■□■□■□

946 名前:デフォルトの名無しさん mailto:sage [2019/01/26(土) 17:16:36.10 ID:5cGQ6lSj.net]
お題: マークダウン記法の見出しをHTMLに変換しなさい

947 名前:デフォルトの名無しさん [2019/01/26(土) 18:50:28.91 ID:OZN9QunV.net]
>>865
rfo-basic

For i=1 To 44
m$=m$+Chr$((9622+Mod(i,2))*(Mod(i,9)<>0)+10)
Next
Print m$

End

プロは奇をてらわない。
bai 砂坊主

948 名前:デフォルトの名無しさん mailto:sage [2019/01/26(土) 21:39:22.15 ID:Ve68vOks.net]
byでは?

949 名前:デフォルトの名無しさん mailto:sage [2019/01/26(土) 21:56:29.11 ID:lXnCadeF.net]
レス古事記に構うなよ

950 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 00:46:03.01 ID:EZ97VlcS.net]
>>865
@Mathematica

セルオートマトン(ルール32)を使って、、

CellularAutomaton[32, {0, 1, 0, 1, 0, 1, 0, 1}, 8] //
 ArrayPlot[#, Mesh -> True] &

951 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 01:01:04.00 ID:u288xi28.net]
>>865
python3
[print("".join("■" if (x+y) % 2 else "□" for y in range(8))) for x in range(8)]

952 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 04:52:31.63 ID:rwuNPA7J.net]
>>865 C
for (int r = 0; r < 8; r++) {
  for (int c = 0; c < 8; c++) {
    printf((r + c) % 2 ? "■" : "□");
  }
  printf("\n");
}

953 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 06:36:29.52 ID:rwuNPA7J.net]
>>910 C
double abac = (bx - ax) * (cy - ay) - (by - ay) * (cx - ax);
double abad = (bx - ax) * (dy - ay) - (by - ay) * (dx - ax);
double cdca = (dx - cx) * (ay - cy) - (dy - cy) * (ax - cx);
double cdcb = (dx - cx) * (by - cy) - (dy - cy) * (bx - cx);
return (abac * abad < 0 && cdca * cdcb < 0);



954 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 06:44:59.11 ID:mVMeZ4rl.net]
>>923
直線じゃなくて線分だぞ

955 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 06:46:02.03 ID:z0hosMYm.net]
>>916 javascript(es2018 regexネームドグル

956 名前:ープのため)
const atxMatcher = line => line.trim().match(/^(?<directive>\S*)[ \t]+(?<contents>.*)/)
const atxRules = [{
name: 'HnTag',
regex: /^#{1,6}/,
process: (hnDirective, hnContents) => {
const hnEl = document.createElement('h'+ hnDirective.length)
hnEl.textContent = hnContents
return hnEl.outerHTML
}
}]
const makeHtml = line => {
const atxMatch = atxMatcher(line)
if (!atxMatch) return line
const {directive, contents} = atxMatch.groups
let tag
for (const rule of atxRules) {
const directiveMatch = directive.match(rule.regex)
if (directiveMatch) {
tag = rule.process(directiveMatch[0], contents)
break
}
}
return tag ? tag : line
}
`# aaa
##\tbbb
###### \t ccc
################ ddd
`.trim().split`\n`.map(makeHtml).join``
結果:
<h1>aaa</h1><h2>bbb</h2><h6>ccc</h6><h6>ddd</h6>
[]
[ここ壊れてます]

957 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 08:20:11.42 ID:rwuNPA7J.net]
>>924
線分ですよ

958 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 11:17:31.03 ID:/cT7Kgd3.net]
>>923
> return (abac * abad < 0 && cdca * cdcb < 0);
return (abac * abad <= 0 && cdca * cdcb <= 0);
端点も含まんとあかんのとちゃう?

959 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 12:22:09.16 ID:rwuNPA7J.net]
重なりの定義によるかも
接しているのは重なりなしとしました

960 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 12:42:45.94 ID:rwuNPA7J.net]
と思ったけど=ありのが正解でした

961 名前:デフォルトの名無しさん [2019/01/27(日) 13:15:46.24 ID:4tSdie01.net]
お題
半径が整数の3つの円ある。
どの円も他の2つの円に外接している。
3つの円の中心点を結んでできる
三角形がピタゴラスの三角形に
なる場合の3つの円の半径を
100未満についてすべて求める。
例えば3つの円の半径が1,2,3の時
辺の長さが3,4,5のピタゴラスの
三角形ができる。

962 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 13:54:53.07 ID:wgzaqL4l.net]
>>930 C++
https://ideone.com/TvwEx3

963 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 14:02:09.42 ID:/cT7Kgd3.net]
>>930 Java
https://ideone.com/vHvOJu



964 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 14:08:22.33 ID:u288xi28.net]
>>930
python3
https://ideone.com/PPyoFC

965 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 14:21:11.01 ID:/cT7Kgd3.net]
>>932
しまった…半径100未満なのに100まで調べてるw
> r(20, 30, 100) = p(50, 120, 130)
が余分やww

>>930 Java 修正しました
https://ideone.com/BkN6Ls

966 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 17:44:27.67 ID:rwuNPA7J.net]
>>927
=の追加で要件をクリアしたと思ってましたが、対象が線分から直線に変わってました
端点については、以下のコードの呼び出しが4つ必要になります

int abap = (b.x - a.x) * (p.y - a.y) - (b.y - a.y) * (p.x - a.x);
if (abap == 0 &&
    (p.x >= min(a.x, b.x) && p.x <= max(a.x, b.x)) &&
    (p.y >= min(a.y, b.y) && p.y <= max(a.y, b.y))) {
  return 1;
}
return 0;

967 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 18:36:00.36 ID:rE8NR22a.net]
>>930 Ruby 2.5.0

def kotae(max)
  (1..max).each_with_object([]) do |ra, ary|
    x, y = max - ra, 2 * ra**2
    1.step(Integer.sqrt(x), 2) do |n|
      next if x * n**2 < y
      rx = ra + n**2
      m, r = (ra + rx).divmod(n)
      ry = ra * m / n
      next if m.even? || r != 0 || m.gcd(n) != 1
      pcpy = [ra, rx, ry].sort
      ary.push(*(1..max / pcpy[2]).map{|i| pcpy.map{|e| e * i}})
    end
  end
end

pys = kotae(100 - 1)
p pys.size # => 82
p pys # => [[1, 2, 3],
 [2, 3, 10],
 [2, 4, 6],
 # 中

968 名前:略 #
 [33, 66, 99],
 [35, 84, 85],
 [36, 84, 90]]
[]
[ここ壊れてます]

969 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 18:48:49.63 ID:z0hosMYm.net]
>>930 javascript
https://ideone.com/S1PkhD

970 名前:デフォルトの名無しさん mailto:sage [2019/01/27(日) 21:02:27.16 ID:u288xi28.net]
>>910
python3
ax+by=c (d<x<e)
px+qy=r (s<x<t)
とする

import numpy as np
co = np.array([[a,b],[p,q]])
in = np.array([c,r])
ans = np.linalg.solve(co,in)
print(d<ans[0]<e and s<ans[0]<t)

本当は直線が平行なときも調べないといけないが略

971 名前:デフォルトの名無しさん mailto:sage [2019/01/28(月) 02:42:20.01 ID:Ht+H3Zdt.net]
こちらで検証したパターンです
line ab = {{0, 0}, {100, 0}};
line cd1 = {{0, 0}, {100, 0}};    // 重なりあり
line cd2 = {{1, 0}, {100, 0}};    // 重なりあり
line cd3 = {{0, 0}, {99, 0}};    // 重なりあり
line cd4 = {{-1, 0}, {100, 0}};    // 重なりあり
line cd5 = {{0, 0}, {101, 0}};    // 重なりあり
line cd6 = {{-1, 0}, {99, 0}};    // 重なりあり
line cd7 = {{1, 0}, {101, 0}};    // 重なりあり
line cd8 = {{-100, 0}, {0, 0}};    // 重なりあり
line cd9 = {{100, 0}, {200, 0}};    // 重なりあり
line cd10 = {{-100, 0}, {-1, 0}};    // 重なりなし
line cd11 = {{101, 0}, {200, 0}};    // 重なりなし
line cd12 = {{0, -100}, {0, -1}};    // 重なりなし
line cd13 = {{0, -100}, {0, 0}};    // 重なりあり
line cd14 = {{0, -100}, {0, 1}};    // 重なりあり

972 名前:デフォルトの名無しさん mailto:sage [2019/01/28(月) 19:06:17.59 ID:ZiBPaW1B.net]
>>910
2直線を
(x,y)=(x1,y1)+t(x2*x1,y2*y1)
(x,y)=(x3,y3)+s(x4*x3,y4*y3)
としたときt、sが0以上1以下の条件で解を持つことであろう。

973 名前:デフォルトの名無しさん [2019/01/28(月) 20:30:57.09 ID:BBM43iQf.net]
>>901
>>906の改良版。
product [(n - (r - 1))..n]はそのままn P rにもなる。

p n r = product [(n - (r - 1))..n]

Haskell

c'' n r = product [(n - (r - 1))..n] `div` product [1..r]

Python

def fact(n):
a = 1
i = 1
while i <= n:
a *= i
i += 1
return a

def p(n,r):
a = 1
i = n - (r - 1)
while i <= n:
a *= i
i += 1
return a

def c(n,r):
return (p(n,r) // fact(r))



974 名前:デフォルトの名無しさん mailto:sage [2019/01/30(水) 16:02:43.18 ID:WqNbA033.net]
お題
古のスタートレックゲームを作る

975 名前:デフォルトの名無しさん [2019/01/30(水) 16:19:17.11 ID:gOIe9eSU.net]
>>942
そのゲームがどんなものだったのか詳細が書いてある所はある?
昔々 Apple ][ でテキストのやつを人がやってる所は見たことあるのだが細かいルールがどんなんだったか忘れた。
確かワープして移動してってクリンゴンが撃ってくるとかクリンゴンに向けて撃つだとかそんなのを繰り返すんだよな?

976 名前:デフォルトの名無しさん [2019/01/30(水) 19:13:41.13 ID:32wJF+eS.net]
>>865
RFOーBASIC

TEXT.INPUT S$,"","チェス番を入力してください。"
PRINT S$
END

977 名前:デフォルトの名無しさん mailto:sage [2019/01/30(水) 21:23:13.66 ID:zHXVmep5.net]
>>943
この辺りを参考にすればいいかと
https://ja.m.wikipedia.org/wiki/スタートレック_(マイコンゲーム)

978 名前:デフォルトの名無しさん [2019/01/31(木) 17:40:00.64 ID:uIyAv6WO.net]
>>865 Lua
print("░")
横から見た

979 名前:デフォルトの名無しさん mailto:sage [2019/01/31(木) 17:57:03.10 ID:6YZEfuaO.net]
>>941 python の場合、これじゃいかんのか? まだ全くpythonを知らん人間だが。

scipy.special.comb(n, k, exact = False)

python はライブラリを使ってなんぼでしょ。 生の言語で比較するのは時代遅れ。 なんで膨大な処理が必要なAI でpython が主流になったか。

980 名前:デフォルトの名無しさん mailto:sage [2019/01/31(木) 18:21:11.67 ID:JEItRzDd.net]
ライブラリを紹介するスレじゃないんだが。
mathjs
math.combinations(7, 5)

こういうの延々と続けるの?
バカじゃないのwww

981 名前:デフォルトの名無しさん mailto:sage [2019/01/31(木) 18:41:17.06 ID:6YZEfuaO.net]
>>948 別に構わんだろ。 そもそも言語なんてやりたいことを解決するのに便利かどうかだけの話なんだから。

そういう意味では、普通のことをやるにはどの言語でも大差はない。

FORTRAN なんて大昔の言語だがいまだに科学技術計算で使われ続けてきているのは、言語仕様ではなく高速なライブラリが整ってるからに他ならない。

なんでAI や科学技術計算にスクリプト言語のpython がもてはやされてるのかもう一度考えてみる必要があるのでは?

982 名前:デフォルトの名無しさん mailto:sage [2019/01/31(木) 18:50:20.38 ID:6YZEfuaO.net]
いくら言語が優れていても、ある特定の環境だけでしか実用性がない言語と言うのは、別環境用のライブラリが整っていないからに他ならない。

例えばSwift これは、Mac環境下ではObjective-C の全てが使えるが、
言語的にはWindows などでコンパイルできても、ライブラリがないから実用として使うにはまだ程遠い。
API が整っていることが重要。

python がもてはやされてるのは、ライブラリが整ってるからだろ。 言語の良さがどうかなんてあまり関係ない。
ま、個人的な好き嫌いが出てくるのは仕方ないが、実用性があるかどうかはまた別の問題。

983 名前:デフォルトの名無しさん mailto:sage [2019/01/31(木) 18:58:10.77 ID:6YZEfuaO.net]
そもそも、言語なんてアセンブラの時代からマクロで自分好みの文をかけるようになってたんだから、今更何を議論してるのかという気もする。

標準化された文法というのは評価はするが、それはどっちかと言うと頭でっかちな学者が考えること。
現場はいかにバグが少なく早く問題を解決できるかにかかってる。 つまり総合環境だろ。



984 名前:デフォルトの名無しさん mailto:sage [2019/01/31(木) 18:59:37.75 ID:aRQ4w2UO.net]
溜まってたものを一気に語ったな

985 名前:デフォルトの名無しさん [2019/01/31(木) 19:46:27.30 ID:REgWEVPY.net]
うむ。これからの時代はファミリーBASICだな。

986 名前:デフォルトの名無しさん mailto:sage [2019/01/31(木) 19:59:18.32 ID:6YZEfuaO.net]
>>953 最初の頃のBASICは、10数個位しか命令は無かった。
ま、電卓がわりみたいなものでは有るが、基本的なプログラミング言語としては成り立ってた。

987 名前:デフォルトの名無しさん [2019/01/31(木) 20:11:31.45 ID:4JULsXj7.net]
>>947
たまたま、ライブラリで解決出来る問題だからってマウント取られてもね。。。
それに、C/C++みたいに多倍長整数サポートしてない言語と、サポートしてる言語でのアルゴリズムの違いや、
式変形して引き算を足し算に変えられる所は変えた方が速度出るとか、色々勉強にはなってるよ。
>>941の例じゃ(n - (r - 1))を(n - r + 1)にした方が速い。
(Haskellでは速いだけでなく、2の補数表現を(元の数値と別に)生成してるらしく、メモリ効率も上がる)

988 名前:デフォルトの名無しさん mailto:sage [2019/01/31(木) 21:19:24.42 ID:Te48fFDx.net]
>(n - (r - 1))を(n - r + 1)にした方が速い。

>>955 あのさ、そんなのはコンパイラの最適化次第だろ。 勿論ど素人のプログラマじゃ無ければそれなりの工夫はするはずだが、そんな事は些細な事。 くだらんことを言うな。 むしろそんなコンパイラがクソ。 言語なんてのは言語道断。

989 名前:デフォルトの名無しさん mailto:sage [2019/01/31(木) 22:25:01.31 ID:xJsSt9Re.net]
Python では、機械学習などの周辺ライブラリ一式が揃っているから、使われる

Ruby では、ベクトル演算のNArray があって、
処理速度はOctave にも匹敵するけど、一式揃っていない

北大の湊教授が作った、ZDD はあるけど

990 名前:デフォルトの名無しさん mailto:sage [2019/01/31(木) 22:46:32.84 ID:2IY59Hh/.net]
scipy/scipy/special/cephes at master ・ scipy/scipy

991 名前:
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

このあたりをそのまま書ければいい
[]
[ここ壊れてます]

992 名前:デフォルトの名無しさん [2019/02/01(金) 04:58:43.05 ID:IMBPiIxm.net]
お題
平方数かどうか判定する

993 名前:デフォルトの名無しさん mailto:sage [2019/02/01(金) 10:08:25.26 ID:VdCPb4pG.net]
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;
}



994 名前:デフォルトの名無しさん [2019/02/01(金) 18:32:58.66 ID:iYm26gLc.net]
>>959
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

995 名前:デフォルトの名無しさん [2019/02/01(金) 20:32:42.61 ID:iYm26gLc.net]
Python

def isSqr(x):
n = 0
a = 0
while x <= a:
if x == a:
return True
n += 1
a = n * n
return False

リスト内包表記だとすぐ落ちた。。。
Haskellより持たない。。。

996 名前:デフォルトの名無しさん [2019/02/01(金) 21:10:55.66 ID:qFuZrTt3.net]
これはいかに巨大数を高速判定するかだろ

997 名前:デフォルトの名無しさん [2019/02/01(金) 21:19:45.40 ID:qFuZrTt3.net]
目星はこれを荒く使えばつきそう


逆数と平方根を求める高次収束アルゴリズム
www.finetune.co.jp/~lyuka/technote/fract/sqrt.html

高速根号計算
takashiijiri.com/study/miscs/fastsqrt.html

998 名前:デフォルトの名無しさん mailto:sage [2019/02/01(金) 22:13:46.18 ID:kNvVsHFY.net]
https://code.i-harness.com/ja-jp/q/4829b

999 名前:デフォルトの名無しさん [2019/02/01(金) 23:10:20.49 ID:iYm26gLc.net]
>>963
m = ?
n = m * m

m=10の時nの1/10で、m=100だとnの1/100

探すべきmが高速で小さくなって行く。
大きな数になる程、0<=mなループの始まり求めるの無理ゲー。。。

1000 名前:デフォルトの名無しさん mailto:sage [2019/02/01(金) 23:12:44.04 ID:VdCPb4pG.net]
ttps://ideone.com/xr1vcE

Newton-Raphson法ほど速くはないがせめてO(log(N))くらいで実装せにゃお題として出てくる意味がなかろう

1001 名前:デフォルトの名無しさん [2019/02/02(土) 00:09:40.70 ID:/6KX0oFw.net]
>>959 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))

1002 名前:デフォルトの名無しさん [2019/02/02(土) 03:47:32.54 ID:i2SNxKFt.net]
いっそ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)

1003 名前:デフォルトの名無しさん [2019/02/02(土) 09:51:18.58 ID:/6KX0oFw.net]
>>959 J
f =: = *: @ <. @ %:



1004 名前:デフォルトの名無しさん [2019/02/02(土) 13:18:27.64 ID:CwD+xRo8.net]
平方数かどうかを高速に判定する方法 - 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が計算できるのです。

1005 名前:デフォルトの名無しさん mailto:sage [2019/02/02(土) 13:43:47.68 ID:OgiywF+Q.net]
>>949
えぇ…

1006 名前:デフォルトの名無しさん mailto:sage [2019/02/02(土) 16:31:05.91 ID:9W2pTWu+.net]
>>972 その心は?

1007 名前:デフォルトの名無しさん mailto:sage [2019/02/02(土) 16:42:03.23 ID:OgiywF+Q.net]
>>970
えぇ…

の間違い。

1008 名前:デフォルトの名無しさん mailto:sage [2019/02/02(土) 17:31:46.69 ID:rEiZ26fd.net]
本質的に高度な数学的知識が要求され
る問題は板違いかと

言語によっては、問題にジャストフィット
するライブラリが標準で添付されているとか
ライブラリ管理コミュニティを通じて
容易に入手できるとかあるかもしれないが
その場合は紹介程度に。

言語処理系の言語処理の為の機能提供
である場合を除いて多くは本質的に多言
語に対応しており〜言語用のライブラリと
いう表現は兎も角〜言語のライブラリと表
現すると曖昧で誤解を招きやすい表現
として嫌われる場合もありえることに注意。
(別言語で記述される場合もある)

1009 名前:デフォルトの名無しさん mailto:sage [2019/02/02(土) 17:36:46.77 ID:XgXX/tZQ.net]
なんしたのきゅうに

1010 名前:デフォルトの名無しさん mailto:sage [2019/02/02(土) 17:57:26.00 ID:rEiZ26fd.net]
セルフホスティング対応な言語でライブラリを
自前で構成しているもので、他言語に移植されて
いない独自の機能を持つものもあるかもしれない
し、もっと極端に言えば独自性の高いライブラリに
最適化された言語を自前で作ってそれで記述され
ているものもあるかもしれないけど一般的には入手
は容易ではないかも

1011 名前:デフォルトの名無しさん mailto:sage [2019/02/02(土) 19:52:10.51 ID:hDNgHqpo.net]
じゃあさ、各々複数行の
a.txt と、b.txtを並べて表示して。

1012 名前:デフォルトの名無しさん [2019/02/02(土) 23:25:24.42 ID:g8xy/J6N.net]
辺に沿って動くとき、AからBまでの最短経路はいくつあるか
         ┏┳┳┳┓B
         ┣╋╋╋┫
      ┏┳╋╋╋╋┫
┏┳┳┳╋╋╋╋╋┻┛
┣╋╋╋╋╋╋╋┫
┣╋╋╋╋╋╋╋┫
┣╋╋╋╋╋╋╋┛
┗┻┻┻┻┻┻┛
A

1013 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 00:45:38.22 ID:UGH880J+.net]
おねえさんがロボットになるやつかw



1014 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 00:55:19.58 ID:UTpNqxd4.net]


1015 名前:学受験の算数の問題に経路を数える問題があるよな
つまり小学生でも手計算で解ける
[]
[ここ壊れてます]

1016 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 00:57:22.00 ID:72eosYJ+.net]
>>975 アホちゃうの?
プログラムというのは、問題解決のための言語であって、問題が解決できなければどんなに美しい言語でも存在価値はない。
最速で美しく問題解決にたどり着ける言語が一番。 言語そのものの美しさなんて二の次。

英語でも日本語でも何でも良い。 利用できるものは何でも利用すれば良い。 好き嫌い言ってる奴はアホ。
どんなに綺麗な国の言葉でも高等教育が自国語でできない国がほとんど、それは利用できる教材が自国語で書いたものがないから。
日本語だったら英語を知らなくてもノーベル賞が取れるだけの教材(ライブラリ)が揃ってる。

1017 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 01:53:34.29 ID:wxxHWwaf.net]
別に高度な数学が必要になることは構わないけど、数学だけで話が閉じちゃうような問題はつまらないとは思う。数学的に解を求める解法自体が主題で、コーディングはただその解法を(日本語、英語などの自然言語と同様に)とあるプログラミング言語で記述しただけのような。
もちろん数学に価値がないといっているのではなく、話のウェイトとしてプログラムである部分がほぼないのなら、その話、ここじゃなくてもいいじゃんと思うよ。

1018 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 04:02:35.76 ID:JyP+XfGy.net]
つまらない問題だと思ったら無視すればいいだけじゃないかな
その問題に対して回答する人がいるなら、回答者(と出題者)にとっては
つまらない問題ではないってことなんだろうし

1019 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 04:46:34.72 ID:1lu6X4vo.net]
>>982
OS記述の言語は言語じゃないのか?
コンパイラなどの言語処理系記述の言語
も言語じゃないってわけだな
問題解決手段の一つではあるかもしれないが
通信系ソフトウェアの存在はどう説明?

既存の資産を維持しより効率的に活用する研究が
新しいプログラミング言語が次々と生成されてくる原
動力じゃないかと

1020 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 05:08:39.33 ID:l3Qt7IvN.net]
>>982
>日本語だったら英語を知らなくてもノーベル賞が取れるだけの教材(ライブラリ)が揃ってる。

ダウト。ノーベル賞を受賞した研究者はみんな英語で論文書いている。
あと、今の日本語での教育教材の蓄積があるのは、大雑把に言って、
ここ10年でノーベル賞受賞した面々の世代の研究者や技術者が書いた。
もちろん彼らは英語の文献を読み、勉強し、日本語の文献を書いた。
だから、英語の文献なんて要らないみたいな話はアホ。
ソファでポテチ食いながら「ジャガイモなんて買えばいいんだから、育ててる奴はアホ」と言ってるようなもの。

1021 名前:981 mailto:sage [2019/02/03(日) 05:15:20.75 ID:l3Qt7IvN.net]
プログラミング言語も一緒。
Cは泥臭い実用言語で、その前にはAlgolなど「綺麗な」言語が下敷きとしてあった。
Cは現在でも多くの問題を解決するための道具として活用されているが、
現在では使われていないAlgolなどの「綺麗な」言語なしにCは生まれなかった。

1022 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 05:39:02.55 ID:1lu6X4vo.net]
個々の言語の歴史観の講釈はスレ違い

確かに問題解決の道具であることが
中心かもしれないが、扱ってきた問
題の種類によって文法やライブラリ・
その取り扱い方に差異が生じている
が、同じ問題を別言語で「解く」と
優劣の違いがわかって面白いかもし
れない。
が、あんまし長くやってると不毛な
言語比較論、文化比較論になったり
して色々ヤバいからそろそろ一旦
お開きにしたら?


1023 名前:きはそれぞれの言語別スレッドで
ということで
[]
[ここ壊れてます]



1024 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 07:04:30.87 ID:LaZtKDWq.net]
お題:プログラムの実行時刻が午前なら「おはようございます、ご主人様!」、午後なら「お疲れ様です、ご主人様!」と表示させる

1025 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 07:48:34.73 ID:AEg+fU/i.net]
>>989 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");
}

1026 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 08:56:17.85 ID:l3Qt7IvN.net]
>>989
Pharo Smalltalk

Smalltalk ui inform: (Time now meridianAbbreviation = 'AM' ifTrue: [ 'おはようございます、ご主人様!' ] ifFalse: [ 'お疲れ様です、ご主人様!' ])

1027 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 09:19:47.58 ID:72eosYJ+.net]
python

from datetime import datetime
if datetime.now().hour < 12:
print('おはようございますご主人様')
else:
print('お疲れ様です、ご主人様')

1028 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 09:23:27.86 ID:1lu6X4vo.net]
午前12時=00:00
午後12時=12:00
23:59の後は00:00

午前12時=12:00
午後12時=24:00
24:00の後は00:01

1029 名前:デフォルトの名無しさん [2019/02/03(日) 09:40:26.58 ID:cfde/ig7.net]
>>865 Common Lisp
https://pastebin.com/TLD3d9R1

実行結果
https://i.imgur.com/6Dne4jz.png

1030 名前:デフォルトの名無しさん [2019/02/03(日) 09:58:12.21 ID:I0qputsI.net]
>>969 の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

1031 名前:デフォルトの名無しさん [2019/02/03(日) 10:00:13.11 ID:I0qputsI.net]
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

1032 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 10:13:31.71 ID:xEPkQ4sk.net]
この改行長文おじさんはどこからきたの。笑

1033 名前: mailto:sage [2019/02/03(日) 10:31:38.09 ID:t4xt++Qj.net]
>>982
>日本語だったら英語を知らなくてもノーベル賞が取れるだけの教材(ライブラリ)が揃ってる。
最近はそうでもないようですよ…ペーパーは英語だし、教科書=テキストレベルでも英語でしか発刊されない状況といいます
haskell をやろうとして圏論の教科書を探しましたが、欧米の本の和訳ばかりで日本人が書いた圏論の教科書はありませんでした



1034 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 10:35:41.62 ID:/jO+7TC8.net]
誰か次スレ頼む

1035 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 11:07:15.54 ID:72eosYJ+.net]
お題1: 現在地の緯度、経度を出せ
緯度:、、、、
経度:、、、、
お題2: 東京都新宿区西新宿2丁目8-1 の緯度、経度を出せ
緯度:、、、
経度:、、、
お題3: お題2で求めた緯度経度から住所を出せ
郵便番号:、、、
住所:東京都、、、、

1036 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 11:22:46.76 ID:72eosYJ+.net]
立てたよ
プログラミングのお題スレ Part13
https://mevius.2ch.net/test/read.cgi/tech/1549160513/

1037 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 11:23:36.56 ID:72eosYJ+.net]
>>1000 は、次スレに移動させるね。

1038 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 17:37:12.27 ID:csrqlAvs.net]
うめ

1039 名前:デフォルトの名無しさん [2019/02/03(日) 17:37:50.59 ID:oUppVF8S.net]
>>969
今までの苦労は一体。。。

数学的な平方根の近似値は√x = x ^ (1/2)だった。。。

Haskell だとこんだけ。

mysqrt x = x ** 0.5

1040 名前:デフォルトの名無しさん mailto:sage [2019/02/03(日) 17:38:30.96 ID:csrqlAvs.net]
次スレ

プログラミングのお題スレ Part13
https://mevius.2ch.net/test/read.cgi/tech/1549160513/

1041 名前:1001 [Over 1000 Thread.net]
このスレッドは1000を超えました。
新しいスレッドを立ててください。
life time: 128日 7時間 29分 24秒

1042 名前:過去ログ ★ [[過去ログ]]
■ このスレッドは過去ログ倉庫に格納されています






[ 新着レスの取得/表示 (agate) ] / [ 携帯版 ]

前100 次100 最新50 [ このスレをブックマーク! 携帯に送る ] 2chのread.cgiへ
[+板 最近立ったスレ&熱いスレ一覧 : +板 最近立ったスレ/記者別一覧](;´∀`)<360KB

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef