プログラミングのお題スレ 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/ 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] 数学の問題が叩かれるのはここのアホどもには難しいからなんだ！ 決してスレ違いだからじゃないんだ！ だから高校１年の教科書の定義引き写すだけの簡単な問題なら叩かれないんだ！ 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 プロは奇をてらわない。 ｂａｉ　砂坊主 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まで調べてるｗ > r(20, 30, 100) = p(50, 120, 130) が余分やｗｗ >>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 ２直線を (x,y)=(x1,y1)+t(x2*x1,y2*y1) (x,y)=(x3,y3)+s(x4*x3,y4*y3) としたときｔ、ｓが０以上１以下の条件で解を持つことであろう。 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) こういうの延々と続けるの？ バカじゃないのｗｗｗ 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] おねえさんがロボットになるやつかｗ 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 ＯＳ記述の言語は言語じゃないのか？ コンパイラなどの言語処理系記述の言語 も言語じゃないってわけだな 問題解決手段の一つではあるかもしれないが 通信系ソフトウェアの存在はどう説明？ 既存の資産を維持しより効率的に活用する研究が 新しいプログラミング言語が次々と生成されてくる原 動力じゃないかと 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] このスレッドは１０００を超えました。 新しいスレッドを立ててください。 life time: 128日 7時間 29分 24秒 1042 名前：過去ログ ★ [[過去ログ]] ■ このスレッドは過去ログ倉庫に格納されています

---

---

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

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