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

プログラミングのお題 ..

658:デフォルトの名無しさん
17/11/16 21:25:39.77 /xRbPsNU.net
数学ムズイ。。。
PGも算数で解いてるからな。

659:デフォルトの名無しさん
17/11/17 00:16:40.63 5DUWZGJy.net
紙とペ

660:ンで考えてみたところ 0以外の任意の整数なら3,5,7で表わせるから問題として不適なのでは？

661:デフォルトの名無しさん
17/11/17 00:20:03.76 M2EFWWXH.net
17

662:デフォルトの名無しさん
17/11/17 10:03:52.61 oe8UBfUe.net
>>637
3,5,7は１回までって回数制限があるから
表せる数は限られるよ。

663:デフォルトの名無しさん
17/11/17 10:44:12.02 a6b9gyRQ.net
17が出来ない

664:デフォルトの名無しさん
17/11/17 19:35:00.15 5DUWZGJy.net
ああ、本当だ。17はどうやっても作れないね
しかしこれをどうやってコードで計算すんだろう
^2があるから全探査はできないし
自分は「+または-」をいくつ使うかで場合分けして一個一個可能性を消していったんだけれども

665:デフォルトの名無しさん
17/11/17 20:04:40.56 M2EFWWXH.net
コードはアップしないけど出来たよ

666:デフォルトの名無しさん
17/11/17 20:07:16.14 M2EFWWXH.net
独自有理数クラス
演算回数を1回ずつ増やしていって、
出来た値に対応するフラグをセット

667:デフォルトの名無しさん
17/11/17 20:09:06.68 M2EFWWXH.net
数値をstd::multisetで保持
演算n回目のmultisetをstd::setで保持

668:デフォルトの名無しさん
17/11/18 17:51:34.79 6foiYhRZ.net
ABC4D
-E3FG
-----
77777
A～G は、1～9 の異なる数字。
ただし、3, 4 ではない

669:デフォルトの名無しさん
17/11/18 17:56:36.72 R4dFDjUs.net
はい、そうですか。

670:デフォルトの名無しさん
17/11/18 19:08:28.21 8fhXEikQ.net
>>645 Java
2年前の問題と俺の回答
ｽﾚﾘﾝｸ(tech板:451番)
を改造したもの (50-54行目と標準入力)
URLﾘﾝｸ(ideone.com)

671:デフォルトの名無しさん
17/11/18 19:44:36.57 oFg54zrO.net
>>645 Ruby
f = ->a, b, c, d, e, f, g{10000*a + 1000*b + 100*c + d - (1000*e + 10*f + g) == 78037}
[1, 2, 5, 6, 7, 8, 9].permutation{|a| puts "%d%d%d4%d - %d3%d%d == 77777" % a if f[*a]}
#=>87142 - 9365 == 77777

672:デフォルトの名無しさん
17/11/18 20:40:01.77 6foiYhRZ.net
>>647
何も、そこまで作り込まなくても良いだろw
色々な覆面算に対応するため、汎用的に書いたのか

673:デフォルトの名無しさん
17/11/19 22:39:02.74 oda4btU4.net
500, 100, 50, 10, 5, 1円のすべての種類の硬貨を、1枚以上使って、
合計15枚で750円にする時、10円硬貨は何枚になるか？

A～E の5人のランナーが走った結果、
完走したのは、1着とべべの2人で、残りの3人は、途中で棄権した
ここで、完走した2人は、必ず真実を言い、
棄権した3人は、必ず嘘をつくものとする
(つまり、事実に対して、真偽値を取る)
A: D は棄権した
B: A は、べべだった
C: E は棄権した
D: C は、べべだった
E: B は完走した
A～Eがこのように答えた時、1着は誰か？

674:デフォルトの名無しさん
17/11/20 01:25:03.39 Z32/GYkn.net
先に答えやそれに至る式がわかっててコードに書き直すだけになっちゃうから
数学的に道筋立てて答えが出せるものはあんまりおもしろくないんだよな

675:デフォルトの名無しさん
17/11/20 03:34:46.27 GkhyFhEh.net
アルゴリズムとは、数式の完全コピー
最初に、数式を考えて、その数式が間違っていれば、
撃墜モードでは、そこを突かれて撃墜される
結局、数式の証明が大事。
証明に、勘違いが無いかどうか

676:デフォルトの名無しさん
17/11/20 11:24:27.08 7i/OQPcC.net
べべってなんぞ？

677:デフォルトの名無しさん
17/11/20 11:38:36.82 OJcNabXy.net
>>653
たぶんこの場合は大阪の方言

678:デフォルトの名無しさん
17/11/20 11:45:30.45 7i/OQPcC.net
俺地方の人間だからわかんない。

679:デフォルトの名無しさん
17/11/20 18:24:16.80 Slkhafwt.net
べべって最下位のことじゃないか
どべ

680:デフォルトの名無しさん
17/11/21 23:50:05.58 zUV8sDjk.net
>>645
こんなのどうかな。Kotlin で作った。
URLﾘﾝｸ(paiza.io)

681:デフォルトの名無しさん
17/11/23 10:05:44.85 zWeuVerg.net
お題
1から99を表示する
お題：1から999を出力する
ただし0を含む数は除く

682:デフォルトの名無しさん
17/11/23 10:11:40.40 TrZHjzbP.net
>>658
1000.times{|i|p i unless i.to_s[?0]}

683:デフォルトの名無しさん
17/11/23 12:15:31.91 6AL/1aep.net
>>658 GNU Smalltalk
1 to: 999 do: [:n | (n asString includes: $0) ifFalse: [n displayNl]]

684:デフォルトの名無しさん
17/11/23 12:38:42.26 KUvGqrz8.net
>>658 F#
let () = seq { 1..999 } |> Seq.iter (printfn "%d")

685:デフォルトの名無しさん
17/11/23 12:40:24.86 KUvGqrz8.net
>>661
すまん、問題よく読んでなかった･･･

686:デフォルトの名無しさん
17/11/23 13:30:23.43 jBvfUrCY.net
>>658
URLﾘﾝｸ(ideone.com)
C++。ほかの言語だと一行で書けるんだけどなぁ。
まぁ過去に比べれば大分短くなったけど。

687:デフォルトの名無しさん
17/11/23 15:45:02.09 ys+VuKpG.net
>>658
文字も数もその場に合わせて適当に解釈してくれる言語だと楽だね。
perl だとこれでできる。
for(1..999){print"$_\n"unless(/0/)}

688:デフォルトの名無しさん
17/11/23 16:16:10.88 JcpJJmmU.net
>>658
@Mathematica
nListWithoutZero[n_]:=n//
　Range[1,#]&//
　Map[ToString,#]&//
　StringCases[#,RegularExpression["^(?!.*0).*$"]]&//
　Flatten;
In[1] := nListWithoutZero[999]
Out[1] = (略)

689:デフォルトの名無しさん
17/11/23 16:34:54.15 2sNCCDGP.net
>>658
ダメだお題の意味がわからん
降参

690:デフォルトの名無しさん
17/11/23 16:42:00.54 QTAUjuBR.net
>>658 rust
URLﾘﾝｸ(ideone.com)
fn main() {
println!("{:?}", (1..1000).filter(|i| !i.to_string().contains("0")).collect::<Vec<_>>())
}

691:デフォルトの名無しさん
17/11/23 16:46:25.38 ys+VuKpG.net
>>658
Kotlin で文字列変換してやる場合
fun main(args: Array<String>) {
　for (i in 1..999)
　　if (! i.toString().contains('0', false))
　　　println(i)
}
数値のままやる場合
fun main(args: Array<String>) {
　for (i in 1..999)
　　if (i % 10 != 0
　　　　 && (i < 10 || i / 10 % 10 != 0)
　　　　 && (i < 100 || i / 100 % 10 != 0))
　　　　println(i)
}

692:デフォルトの名無しさん
17/11/23 17:05:13.53 fGVRHt7J.net
>>658
>>668 同じく数値のままやる場合
URLﾘﾝｸ(ideone.com)

693:デフォルトの名無しさん
17/11/23 17:08:37.18 fGVRHt7J.net
ｽﾚﾘﾝｸ(math板:722番)
お題：
n^2-1 = m^5
を満たす自然数 n, m は存在するか？
存在するという人と存在しないという人の両方が存在します

694:デフォルトの名無しさん
17/11/23 17:18:38.32 TrZHjzbP.net
>>670
(n, m) = (1, 0)
揚げ足取りはおいておいて、プログラミングで説く問題じゃないよね

695:デフォルトの名無しさん
17/11/23 18:35:09.23 zveldNvP.net
>>658
python
今回は必要ないかもだけど桁数増えた場合を考え再帰で
URLﾘﾝｸ(ideone.com)

696:デフォルトの名無しさん
17/11/23 19:09:48.15 fGVRHt7J.net
>>671
まあ自明な解はさておき、その他は見つからないのが不思議です

697:

698:デフォルトの名無しさん
17/11/23 20:21:54.46 /mQ4CZGQ.net
>>673
カタラン予想ですでに存在しないことが証明されているのに何が不思議なのかね

699:デフォルトの名無しさん
17/11/23 20:24:10.50 hjkeK8jf.net
>>673
その問題は数学的に解くものではないかな?
まあ、コンピュータなら力業でかなりの値を n, m に入れて計算して確認できるけどさ。

700:デフォルトの名無しさん
17/11/23 20:56:22.87 fGVRHt7J.net
>>674
カタラン予想というのがあるんですね、ありがとうございます！
URLﾘﾝｸ(en.wikipedia.org)

701:デフォルトの名無しさん
17/11/23 21:05:28.66 uF7hi9HH.net
>>658
#!/bin/sh
seq 999|grep -v 0

702:668
17/11/24 06:21:36.98 8wyGH9pr.net
>>658
Kotlin数値判定版。こんな風にも書けるなと後で気づいた。
fun f(n: Int): Boolean {
　var m = n;
　while (m != 0) {
　　if (m % 10 == 0)
　　　return false
　　m = m / 10
　}
　return true
}
fun main(args: Array<String>) {
　(1..999).filter(::f).forEach(::println)
}

703:デフォルトの名無しさん
17/11/24 07:42:25.55 MEHEP0+e.net
存在するしないをプログラミングで証明するのはお題として良くない

704:デフォルトの名無しさん
17/11/24 20:42:02.54 G34PGfZh.net
log 2 を2進数表記した時の小数点第 n 位から n + 9 位までを求めよ． (1 ≦ n ≦ 10^10)
cf. log 2 = 0.10110001...
*Sample input*
1
11
10000
31415926
314159265
*Sample output*
1011000101
1100100001
0010110110
1001010110
0111101001

705:デフォルトの名無しさん
17/11/24 23:31:00.22 r53+zpq0.net
>>680
c++で書いたけど小数第100億位を計算するのに5時間くらいかかりそうorz

706:
17/11/25 00:04:03.20 ROI3Hzdd.net
>>681
もう初手に届くとは劇速ですね

707:デフォルトの名無しさん
17/11/25 06:53:37.62 Uo3oYb2P.net
無条件でlogって書いたら普通自然対数だろ

708:デフォルトの名無しさん
17/11/25 07:01:13.88 Uo3oYb2P.net
ライブラリを使えばほとんど何も書かなくて良いけど
どこから書くことを求められてるの？

709:デフォルトの名無しさん
17/11/25 07:03:54.53 Uo3oYb2P.net
>>679
「良くない」じゃなくて「出来ない」でしょ

710:デフォルトの名無しさん
17/11/25 07:16:55.64 Uo3oYb2P.net
>>684
と思ったけど、普通に全桁計算したら終わらないな

711:デフォルトの名無しさん
17/11/25 07:34:12.57 Uo3oYb2P.net
Σ { 1 / (2^i × i) }
を使って10^10項位までを42bitくらいだけ計算すれば出来るかな？
1/nの周期性を考えないと計算量的に無理？
10^10が微妙に32bitを越えてるのがイヤだねえ

712:デフォルトの名無しさん
17/11/25 08:21:38.07 Uo3oYb2P.net
>>687
ダメだ
ざっと計算量を見積もったらとても5時間じゃ終わらない

713:デフォルトの名無しさん
17/11/25 13:10:34.24 l6j6CjYT.net
>>375
xxx@xxx-VirtualBox:~/casl$ casl -s -e -i stdlib.casl -i bigint.casl fact.casl
1
1
2
6
24
120
720
5040
40320
362880
　　　　　　　　　　途　中　省　略
1405006117752879898543142606244511569936384000000000
60415263063373835637355132068513997507264512000000000
2658271574788448768043625811014615890319638528000000000
119622220865480194561963161495657715064383733760000000000
5502622159812088949850305428800254892961651752960000000000

714:258623241511168180642964355153611979969197632389120000000000 12413915592536072670862289047373375038521486354677760000000000 608281864034267560872252163321295376887552831379210240000000000 30414093201713378043612608166064768844377641568960512000000000000 暇つぶしに書いてみたけど足算掛算割算しかできない引算は難しすぎるんで諦めた

715:デフォルトの名無しさん
17/11/25 16:35:37.76 J1zvm3XW.net
バイナリ法で最適化した結果なんとか1時間あれば10^10位は計算できるようになったがまだ縮められるかな

716:
17/11/25 16:41:37.19 ROI3Hzdd.net
備忘メモ
ｽﾚﾘﾝｸ(math板:890番)
ちゃんとお題にして出題しなおします

717:
17/11/25 21:49:51.82 ROI3Hzdd.net
>>680
指数関数のマクローリン展開で試してみたのですが、これは収束が遅すぎますね、それに収束半径を超えてるし…
なにか収束の早いよい方法はないものか…

718:
17/11/25 21:54:37.24 ROI3Hzdd.net
>>692
×指数関数
○対数関数

719:デフォルトの名無しさん
17/11/25 22:58:19.91 yyAYDlfh.net
対数関数のマクローリン展開？
そりゃ無理だ
log 0 が定義されてない

720:デフォルトの名無しさん
17/11/25 23:12:41.85 FRsJtlII.net
>>689
CASLで書いたの?
ソースコードは?

721:デフォルトの名無しさん
17/11/25 23:38:41.06 yyAYDlfh.net
>>680
log 2 = Σ_[i=1, 2, ...] { 1 / (2^i × i) }
冪剰余
でいける気がしてきた
しばらく暇がない
時間が空いたら
アセンブラ & C++ & OpenMP
でやってみる

722:デフォルトの名無しさん
17/11/26 02:33:02.82 T275kIwU.net
>>650
(setq aaa '(1 5 10 50 100 500))
(setq ddd 750)
(setq jjj 15)
(defun bbb (ccc iii)
(if (= iii 0)
ccc
(let (eee)
(dolist (fff ccc)
(dolist (ggg aaa)
(when (<= (+ (apply #'+ fff) ggg) ddd)
(push (cons ggg fff) eee))))
(bbb (remove-duplicates (mapcar (lambda (x) (sort (copy-seq x) #'<)) eee) :test 'equal) (1- iii)))))
(let* ((kkk (bbb '((0)) jjj))
(lll (mapcar (lambda (x) (remove 0 x)) kkk)))
(remove-if-not (lambda (x) (and
(= (apply #'+ x) ddd)
(= (length x) jjj)
(= (length (remove-duplicates x)) (length aaa))
)) lll))
((1 1 1 1 1 5 5 5 10 10 10 50 50 100 500))

723:デフォルトの名無しさん
17/11/26 02:34:12.86 T275kIwU.net
>>650
(setq aaa '(A B C D E))
(defun fff (ddd)
(if (null (cdar ddd))
ddd
(let (eee)
(dolist (jjj ddd)
(let ((bbb (car jjj))
(ccc (cdr jjj)))
(setq eee (append (mapcar (lambda (x) (cons (cons x bbb) (remove x ccc))) ccc) eee))))
(fff eee))))
(defun iii (kkk)
(if (< kkk 2) #'identity #'not))
(let* ((ggg (fff (list (cons nil aaa))))
(hhh (mapcar (lambda (x) (car x)) ggg)))
(remove-if-not (lambda (x) (and (funcall (iii (position 'A x)) (> (position 'D x) 1))
(funcall (iii (position 'B x)) (= (position 'A x) 1))
(funcall (iii (position 'C x)) (> (position 'E x) 1))
(funcall (iii (position 'D x)) (= (position 'C x) 1))
(funcall (iii (position 'E x)) (< (position 'B x) 2)))) hhh))
((D C B E A) (D C E B A) (D C A E B) (D C E A B) (D C A B E) (D C B A E))

724:
17/11/26 11:18:11.00 rNgJnhxq.net
>>694
いやいや
URLﾘﾝｸ(ja.wikipedia.org)
log(1-x) = - Σ((1/n)x^n) に x = -1 を機械的に代入しました、収束半径外ですが、この値は正しいらしい。

725:デフォルトの名無しさん
17/11/26 12:12:20.72 ffy1o2uq.net
お題
ASCIIコード表が載っている�

726:{をあげよ

727:デフォルトの名無しさん
17/11/26 12:35:31.87 tJzac9f2.net
>>695
うんCASL
全部で1200行かあ
ｘｘｘ@ｘｘｘ-VirtualBox:~/casl$ wc -l stdlib.casl bigint.casl fact.casl
274 stdlib.casl
851 bigint.casl
76 fact.casl
1201 合計
ソースはこういうのが延々続いててずっと眺めてるとゲシュタルト崩壊起こして
何が何だか分からなくなるよ
ld gr5,0,gr1
ld gr6,1,gr1
lad gr4,4,gr1
addl gr4,gr0
st gr4,0,gr1
st gr6,1,gr1
ld gr4,=1
st gr4,2,gr1
st gr0,3,gr1
ld gr6,gr1
ld gr1,0,gr1
st gr5,0,gr1
st gr6,1,gr1
xor gr4,gr4
st gr4,2,gr1
lad gr2,-4,gr2
subl gr2,gr0
st gr2,3,gr1
ld gr0,gr3

728:デフォルトの名無しさん
17/11/26 12:35:43.37 WgExDItE.net
>>699
他だ単に対数って言えば
log x のこと
これをマクローリン展開は無理
log (1-x) のマクローリン展開ならそう書かないと通じない
収束半径の外じゃなくて、収束半径丁度でしょ

729:デフォルトの名無しさん
17/11/26 12:49:31.90 WgExDItE.net
-log (1-x) のマクローリン展開に、
x = 1/2 を入れると
>>696 になる

730:デフォルトの名無しさん
17/11/26 13:21:30.43 jiBTwXK4.net
理解はしてないが、出てきたので貼っとく。

指数対数関数等の超越関数の多倍精度計算
本論文では、指数対数関数の高精度計算として Taylor 展開に BSA 法を使って高速化する方法提案する。
約 1000 桁以下の精度の計算では、 Taylor 展開を使った計算が Sasaki and Kanada[5] によって、様々な計算
法を比較して最も高速であることが示されているので、計算時間が問題となるのは、 1000 桁以上の精度の
計算である。ここで提案した Taylor 展開に BSA 法を適用して高速化した方法と Sasaki and Kanda によっ
て提案された方法を 1000 桁を超えた精度で比較し、その高速性を示した。
211 階乗計算例
10000! の計算を行う。この計算では、 BSA 法を使うだけでなく、 1600 桁以上の数値に対しては FFT を利用して乗算を行っている。
計算方法　計算時間(msec)
BSA 47
従来の方法　3578
このほか、三角関数、逆三角関数、双曲線関数など簡単な規則で各項の係数が表現でき、多くの関数がこの
行列の乗算形式に変形できます。Taylor 展開の係数が簡単な規則で表現できない $\tan x$ が例外的に表現できないだけである。
3 まとめ
指数関数や対数関数の Taylor 展開に BSA 法を適用することによって、 BSA を使わない従来の方法に比べ40 %程度の高速化ができた。
対数関数に対しては、 5000 桁程度の精度で最も高速な計算方法として知られた Sasaki and Kanada の方法を超えることを示した。
URLﾘﾝｸ(www.kurims.kyoto-u.ac.jp)

731:
17/11/26 13:37:36.96 rNgJnhxq.net
>>702
たしかに
収束半径は |x < 1| なのでは？

732:デフォルトの名無しさん
17/11/26 13:44:48.59 jiBTwXK4.net
とりあえず理解はできた計算方法として、logxの近似値などをaとおいたとき、
logx = a + log(x/e^a)という変形を用いる方法だ。
aが近似値だと、x≒e^aなので良いらしい。

733:デフォルトの名無しさん
17/11/26 13:58:53.69 WgExDItE.net
>>705
収束半径は1
収束半径は、収束するエリアと収束しないエリアの境目となる円の半径
収束半径丁度の時は収束する場合もあるししない場合もある

734:デフォルトの名無しさん
17/11/26 14:20:52.91 WgExDItE.net
>>706
計算する項が少なければ良いわけではなく、
各項の計算時間も重要

735:デフォルトの名無しさん
17/11/26 15:04:05.33 jiBTwXK4.net
exp(x)は、(exp(x/k))^k (kは2ベキ)、とするといいらしい。
k=2なら、括弧内を計算したやつ同士

736:の掛け算。

737:695
17/11/26 18:33:11.28 XHzckn9a.net
>>701
なるほどありがとう
怖いもの見たさがあった

738:
17/11/26 20:54:13.61 rNgJnhxq.net
>>689
えへへ、調べさせてもらったよｗ
後半は 42! から 50! までの値だね
この範囲なら、多数桁×ワンレジスタの計算で済みますね
多数桁×多数桁を実装すれば思いっきり褒めてあげるよ、えへへ：－）

739:デフォルトの名無しさん
17/11/26 21:11:54.45 tJzac9f2.net
>>711
ほい
xxx@xxx-VirtualBox:~/casl$ casl -s -e -i stdlib.casl -i bigint.casl bimul.casl
350306543997676425792
153864088327713953064
53899597027434699691252340823058767026688

740:
17/11/26 21:47:36.23 rNgJnhxq.net
>>712
おお、すごい、gmp で確認したがあってるぞ
私はまだ 10進文字列から多数桁型への変換は実装できていない、また仕事ができてしまったなあ

741:デフォルトの名無しさん
17/11/26 23:01:31.48 tJzac9f2.net
>>713
gmpで確かめてるのかあ
俺は確認はclispかrubyでやってるよ
10進文字列からの変換は10倍しながら足しこんでいくだけだから
そんなに難しくないでしょ
掛算なしでも(n<<3)+(n<<1)でできるし
逆の10進文字列への変換は割算が必要だから実装するの大変だったなあ
これができあがるまではメモリを16進ダンプして計算が合ってるか確かめてた

742:デフォルトの名無しさん
17/11/27 06:49:43.60 qqP20rnw.net
>>696がよさげだな
こういうことだろ?
Σ(1/n) >> n
小数にもシフトを適用するとして。
和の誤差がかわかれば、どこまで計算したらいいかわかるがどうやるんだ?

743:デフォルトの名無しさん
17/11/27 07:17:24.36 qqP20rnw.net
誤差しらべたら、テイラー展開は平均値の定理一般化だったか

数値計算とテイラー展開
ある区間において,関数 f(x)がn 回微分可能であるとし,定数aはこの区間に含まれるものとする.x もこの区間内に含まれるとき,
URLﾘﾝｸ(math-lab.main.jp)
をみたすa とx の間の実数c (a <c <x または x <c <a)が存在する
URLﾘﾝｸ(math-lab.main.jp)

744:デフォルトの名無しさん
17/11/27 07:40:30.25 qqP20rnw.net
log(1+x)の誤差項、剰余項は、(-1)^(n-1)/n * (x/(1+c))^n らしいので、
-log(1-x)では、1/n * (x/(1+c))^n か。
x=1/2で考えると、この項をなるべく大きくするならc=0で、誤差は(1/n) >> n以下か。ふたたびシフト使用。

745:デフォルトの名無しさん
17/11/27 07:44:14.48 qqP20rnw.net
まとめると、
log2 - (Σ(1/k) >> k) < (1/n) >> n 　(級数はn-1までの和)

746:デフォルトの名無しさん
17/11/27 08:54:49.81 qqP20rnw.net
どこか間違えてる　次数か?

747:デフォルトの名無しさん
17/11/27 08:59:59.63 qqP20rnw.net
いやあってるか。
A(k) = (1/k)>>kと置くと、
log2 - ΣA(k) < A(n) (級数はn-1までの和)　で
ΣA(k)　(n+1以上の和) < A(n) が成立するのか。

748:デフォルトの名無しさん
17/11/27 10:11:23.24 qqP20rnw.net
やはり、どこか間違ってるな。
上のとおりだと、log2 - ΣA(k)　(級数はn-1までの和)は、
A(n) を含むので、A(n)より小さいはずがない。

749:
17/11/27 23:07:50.22 b0u4s5jJ.net
>>701
＞ソースはこういうのが延々続いててずっと眺めてるとゲシュタルト崩壊起こして何が何だか分からなくなるよ
対象のコード書いてるときの「感覚」をハードディスクかどこかに貯めておいて、
必要に応じてまた脳みそに搭載できるようにならないものか…
そのマシン語の一行一行にも、もともとはなんらかの意味的構造があったのに、それが消えてしまうなんて損失以外のなにものでもないよね

750:デフォルトの名無しさん
17/11/27 23:36:53.65 bwTh4Bk9.net
>>680
>>696の方法で作ってみました
n=10^10の時に48.3秒くらいです

751:デフォルトの名無しさん
17/11/28 13:09:45.40 IU/PYwM1.net
>>723 はHaswell (4770) 3.4GHz固定での結果で
Skylake (6700K) 定格だと38.4秒でした
ちゃんとCPUも進化してるんですね

752:デフォルトの名無しさん
17/11/28 13:17:18.48 IU/PYwM1.net
>>681
>>690
C++だとどうやって計算してるかが非常に気になります
32bitを越える値同士の乗算(結果が64bitを越える)部分
アセンブラだと
64bit x 64bit ===> 128bit
128bit / 64bit ===> 64bit
等があるのでそれを使っちゃってますが

753:デフォルトの名無しさん
17/11/28 13:20:30.03 IU/PYwM1.net
冪剰余を求めるのに
(a * b) % c
みたいなのがたくさん出てきませんか？
aもbもcも32bitの範囲を微妙に越えてて

754:デフォルトの名無しさん
17/11/28 14:44:32.26 jzUFRHpN.net
誤差部分の間違いが判った。これでよさげだ。
ただし誤差評価を荒くやってはダメそうだが。一番最後の行のところ。

誤差項ありのマクローリン展開は、0<=c<=xが存在して
f(x) = Σ x^k * f(k)(0)/k! (kは0からn-1まで)　+　 x^n * f(n)(c)/n!
f(x) = -log(1-x)のn次導関数は、(n-1)!/(1-x)^n。
このときマクローリン展開は誤差項は x^n / (n*(1-c)^n)
x=1/2ならば、c=1/2のとき最大で、1/n

755:デフォルトの名無しさん
17/11/28 19:45:30.22 jzUFRHpN.net
これが収束速いようだ。

log(2) = 3log(81/80) + 5log(25/24) + 7log(16/15)
log((x+1)/(x-1))
= log((1+1/x)/(1-1/x))
= 2 Σ 1/((2n+1)*x^(2n+1))

756:デフォルトの名無しさん
17/11/28 22:18:22.38 7WoPw74F.net
>>728
1/log(2) ≒ 3.32
1/2log(161)+1/2log(49)+1/2log(31) ≒ 0.85
なので、計算に必要な項数は1/4程度
でも、1つの項の計算には時間がかかる
log(1-x)のマクローリン展開に0.5を入れた物は
分母が i * 2^i だから速く計算できるのだ

757:デフォルトの名無しさん
17/11/28 22:20:46.48 7WoPw74F.net
>>727
残りの項を等比数列と見なせば
簡単に誤差の上限が出ます

758:デフォルトの名無しさん
17/11/28 22:47:09.33 7WoPw74F.net
>>724
Haswellで33.96秒に縮まりました
シングルスレッドだと182.54秒で5.3倍
HTTが効くということは、
まだ多少改善の余地がありそう
一番内側のループは
vmulpd
vmulpd
vroundpd
vfmsub213pd
vfmsub132pd
vsubpd
なんと浮動小数点で計算してます

759:デフォルトの名無しさん
17/11/28 22:53:54.93 7WoPw74F.net
n=10000000000の時は
0000010101 でした
出題者さま、合ってます？

また、たまたまですが
n=10000000004では
0101010101
n=10000000005では
1010101010
になります

760:デフォルトの名無しさん
17/11/28 23:30:15.70 9HoDrqB3.net
一番内側のループのコード
URLﾘﾝｸ(fast-uploader.com)
PORT5がガラ空きで、処理のほとんどがPORT0,PORT1
こんなんでもHTTが効く
やっぱり浮動小数点はレイテンシがデカい
AVX512になれば
レジスタの数が倍になるので
8パラにしてレイテンシを隠蔽出来るんだけど
もちろんレジスタ長が倍になる方が大きい

761:デフォルトの名無しさん
17/11/29 13:17:33.09 mHyZby47.net
>>728は後半部分が間違ってるか。log((x+1)/x) = log(1+1/x) の展開を用いるのが正解で。
log(・)の中身を1に近づけた方が収束が早くなるが、
こういった分解　log(2) = 3log(81/80) + 5log(25/24) + 7log(16/15)はどうみつけるのか。
これは数値が(x+1)/x の形だけど、(x+1)/(x-1)の分解もあるのか。こっちだと計算するベキ項が一つ飛ばしにできる。>>728のように。

762:デフォルトの名無しさん
17/11/29 13:34:57.75 8/kTvoZy.net
2 = (81/80)^3 * (25/24)^5 * (16/15)^7
3 と 5 の指数の合計が0になる組み合わせを検索すれば良い

763:デフォルトの名無しさん
17/11/29 13:37:30.83 8/kTvoZy.net
log(81/80) = log(162/160) = log((161+1)/(161-1))
わかってて書いてるんだと思ったが
>>729のlogの中身はこの値

764:デフォルトの名無しさん
17/11/29 13:42:05.17 mHyZby47.net
>>728はそういうことか。みつけたやつのコピペで、そのとき考慮はしてなかった。

765:デフォルトの名無しさん
17/11/29 13:45:31.51 mHyZby47.net
指数も固定でなくていいはずで、
16/15よりかはたとえば1001/1000のほうが1に近いからそういうのはいくらでも見つけられるのかとおもった。

766:デフォルトの名無しさん
17/11/29 14:44:26.09 8/kTvoZy.net
分母分子の素因数の数と同じ項数が必要
例えば素因数が 2, 3, 5, 7 の4種類の場合、
1個差もしくは2個差のペアを4個探す
例えば
126/125
225/224
2401/2400
4375/4374
これらを適当に掛け算して2^nになるようにすると
項が4個の式がみつかる

767:デフォルトの名無しさん
17/11/30 00:31:43.08 H4qIjcIH.net
分母、分子とも 2, 3, 5, 7, 11, 13, 17 のみしか素因数を持たない形の場合、
以下が一番計算する項の数が少ないようです
log(2) = 72*log(126/125)+27*log(225/224)-19*log(2401/2400)+31*log(4375/4374)

768:デフォルトの名無しさん
17/11/30 05:07:54.11 fMs2N0Mh.net
>>740
その数値を検索すると、音楽のコンマというのが出てくる。関係あったり理論があったりするのか。

A.D. Fokker: Unison Vectors and Periodicity Blocks
URLﾘﾝｸ(www.huygens-fokker.org)
List of 7-prime limit accidentals - The?Sagittal?forum
URLﾘﾝｸ(forum.sagittal.org)

769:デフォルトの名無しさん
17/11/30 05:28:17.43 fMs2N0Mh.net
log(2)とは無関係で、単に一個差のやつで適当な素因数分解できるやつに名前がついてるだけ?

An Investigation into the Extraction of Melodic and Harmonic Features from Digital Audio
unit interval name
4375/4374 Ragisma
2401/2400 Breedsma
225/224 Septimal Kleisma
145/144 Difference between 29:16 and 9:5
126/125 Small Septimal Semicomma
121/120 Undecimal Seconds Comma
81/80 Syntonic Comma
URLﾘﾝｸ(scholar.sun.ac.za)

770:デフォルトの名無しさん
17/11/30 11:02:35.16 fMs2N0Mh.net
log2のほうは、分子・分母の素因数分解が似通ってないと成立しないってことで、
音楽のほうは小さい素数に限定して一個差ペアを求めたと理解。
log2のほうは、共通の素数で大きいやつを最初に固定すれば考えれば、よさげかと。

771:片山博文MZ
17/11/30 12:12:48.85 NsMGt5if.net
お題。横x[cm]、縦y[cm]の長方形のステンレスの１枚の板がある。この板からm枚の複数の長方形の部材を切り出す。
部材のサイズは配列で与えられる。
部材のサイズ(縦×横)はそれぞれだいたい決まっているが、1cm程度変わってもよい。
ただし、部材の縦または横が変わるとそれぞれ一点減点となる。
すべての部材を切り出すことができれば、減点がなるべく少ない方法の切り出し方法を出力せよ。
すべての部材を切り出すことができなければ、面積が広い順になるべくたくさん部材を切り出せ。

772:片山博文MZ
17/11/30 12:19:33.84 NsMGt5if.net
テストデータ。
x=10, y=10,
{
{5, 10}, {2, 2}, {2, 2}, {4, 3}, {6, 5}
}
x=5, y=12
{
{2, 5}, {3, 3}, {2,9}, {3, 2}, {4,3}
}

773:片山博文MZ
17/11/30 12:28:07.36 NsMGt5if.net
部材の縦と横は入れ替わってもよい。
可能ならば、切り出し方法をSVG形式で出力せよ。

774:片山博文MZ
17/11/30 12:51:59.04 NsMGt5if.net
切り出し方法は、
切り出す部材のx座標、y座標、幅、高さ
のリストとして出力せよ。

775:片山博文MZ
17/11/30 12:56:57.25 NsMGt5if.net
切り出しに余裕があるときは、なるべくx座標の大きい方、y座標の大きい方を残すようにせよ。

776:デフォルトの名無しさん
17/11/30 13:08:21.37 tkxPMdZc.net
斜めもOK？

777:デフォルトの名無しさん
17/11/30 13:08:55.87 tkxPMdZc.net
人間用のパズルで、斜めにしないと解けないのとかありそう

778:デフォルトの名無しさん
17/11/30 14:23:24.17 SHLZLl2M.net
問題文の条件が“だいたい”“なるべく”なんて
あいまい表現だらけ
これでプログラミングの問題かよ

779:片山博文MZ
17/11/30 16:18:01.24 NsMGt5if.net
斜めは考えなくてもよい。
訂正。
すべての部材を切り出すことができれば、減点が最小である切り出し方法を出力せよ。
すべての部材を切り出すことができなければ、面積が広い順に切り出せる部材の面積が最大になるよう部材を切り出せ。
切り出しに余裕があるときは、x座標の大きい方、y座標の大きい方を残すようにせよ。

780:680
17/11/30 17:10:37.24 8ZVWPbH7.net
>>732
そこにある3つとも正解です
当初は L = Σ1/(k*2^k) として
2^n * L の小数部分を愚直に求める方法を想定していました

781:デフォルトの名無しさん
17/11/30 17:26:52.63 r8WkgLop.net
普通に多倍長で計算したら計算量的に終わらないですよね？
n=314159265を求めるのに
冪剰余は使ってますよね？
おそらく私も同じような方法と思います
FMA3命令とOpenMPで高速化してるだけで

782:デフォルトの名無しさん
17/11/30 22:49:50.18 TklDiPhy.net
愚直という言い方は良く割りませんでしたね
仰る通り冪剰余は用います

783:片山博文MZ
17/12/01 14:26:32.50 fw1UFg83.net
再出題。横cx[cm]、縦cy[cm]の長方形または正方形のステンレスの１枚の板がある。この板からm枚の複数の長方形または正方形の部材を切り出す。
m枚の部材のサイズは(縦, 横)の配列で与えられる。
すべての部材を切り出すことができれば、切り出し方法を出力せよ。切り出しが不可能ならば「impossible」と出力せよ。
切り出し方法は、
(部材インデックス、部材の一番左のx座標、部材の一番上のy座標、幅、高さ)
のリストとして出力せよ。斜めの方向の切り出しは考えなくてもよい。
切り出しに余裕があるときは、x座標の大きい方、y座標の大きい方を残すようにせよ。ただし、x軸は右の向き、y軸は下向きとする。
テストデータ。
cx=10, cy=10, m=5, {5, 10}, {2, 2}, {2, 2}, {4, 3}, {6, 5}
cx=5, cy=12, m=5, {2, 5}, {3, 3}, {2, 8}, {3, 2}, {4, 3}

784:デフォルトの名無しさん
17/12/01 17:01:28.13 BqJQPTjH.net
頭の悪そうな文章だな
正方形⊂長方形
ステンレスの板の座標上の位置指定が無い
余裕がある場合の条件の意味が曖昧

785:片山博文MZ
17/12/01 18:41:36.20 fw1UFg83.net
ステンレスの板の左上座標は原点にあるものとする。
切り出しは、可能な限り、座標の小さい方を優先する(「余裕」の意味)。

786:デフォルトの名無しさん
17/12/01 20:32:10.86 9YSaSAW0.net
>>758
相変わらず曖昧な表現
全ての長方形の上辺と左辺がどこかに接していれば良いのか？
そうじゃないのか？
このような条件のを1個だけ見つければ良いのか
すべて見つけるのか

787:片山博文MZ
17/12/01 20:41:14.79 fw1UFg83.net
全ての部材の上辺と左辺が別の部材の辺、もしくは、元の板の端に接していること。
このような条件のをすべて見つけること。

788:片山博文MZ
17/12/01 20:43:57.24 fw1UFg83.net
なんか、工学関係でこのような問題があるらしいが、まだ解決策があるかどうかわからん。これが解ければ、実用化待ったなし。

789:デフォルトの名無しさん
17/12/01 21:15:23.73 9YSaSAW0.net
お題じゃなくて依頼かよ

790:デフォルトの名無しさん
17/12/01 21:19:54.38 9YSaSAW0.net
実用性なら
切り

791:やすさとか余り素材の形状とかそういうのが重要だろうに問題として中途半端過ぎる

792:デフォルトの名無しさん
17/12/01 21:20:21.36 qVeescqP.net
また片山博文MZ が乞食をやってるのか

793:デフォルトの名無しさん
17/12/01 21:23:07.33 9YSaSAW0.net
普通に考えて、NP問題だろう

794:デフォルトの名無しさん
17/12/01 21:48:16.44 8H4JUlF5.net
なにやら揉めてますね
そろそろうんざりなので次のお題どうぞ

795:デフォルトの名無しさん
17/12/01 23:36:29.58 N1IVcYDB.net
スポーツスケジューリング問題。

796:デフォルトの名無しさん
17/12/03 02:47:12.16 QazTjKaA.net
お題ってこういうのでもいいのかな
a[10] = { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 }
このように10個の整数を要素に持つ配列がある（整数の値は不定）
b = エントリーポイント
c = 移動する数量
d = 移動する距離
を与えることにより、配列の要素を移動させるプログラムを作りなさい
ただし配列の右端の要素を一つ移動させると、配列の左端に移動するものとする
例
a = 3, b = 1, c = 5, a = 0124567839
b = 1, c = 3, d = 1, a = 0412356789
b = 7, c = 1, d = 5, a = 1273456890
b = 0, c = 8, d = 1, a = 8012345679
b = 4, c = 5, d = 4, a = 6783901245
b = 9, c = 5, d = 4, a = 5679012384
b = 7, c = 3, d = 1, a = 9123456078

797:デフォルトの名無しさん
17/12/03 02:49:08.83 QazTjKaA.net
>>768
一番上の例は
b = 3, c = 1, d = 5, a = 0124567839
の間違いです

798:デフォルトの名無しさん
17/12/03 09:00:27.20 ucQfMVKf.net
>>768
URLﾘﾝｸ(ideone.com)
C++。なかなか素直な問題なのに邪悪な思念で見ていたため手間取った。
マニュアルを熟知していると簡単かもしれない。

799:デフォルトの名無しさん
17/12/03 09:04:48.71 b0AFqm0g.net
>>768 Ruby
URLﾘﾝｸ(ideone.com)

800:デフォルトの名無しさん
17/12/03 10:14:55.69 ZeeZWfzn.net
>>768
@Mathematica
URLﾘﾝｸ(ideone.com)

801:デフォルトの名無しさん
17/12/03 12:12:41.56 jbLAlLAh.net
>>768
b + c が10を超えることはありますか？

802:デフォルトの名無しさん
17/12/03 12:55:22.10 sC6phLBS.net
>>773訂正
b,c,dが10を越えることはありますか？

803:デフォルトの名無しさん
17/12/03 13:05:23.73 b0AFqm0g.net
>>771はb,c,dが任意の自然数でも大丈夫なようにしておいた

804:デフォルトの名無しさん
17/12/03 13:08:19.92 QazTjKaA.net
みなさんプログラム作るの早いですね
>>774
一応a、b、cの条件は、このようになると思います。
0 <= a <= 9
0 <= b <= 8
0 <= c <= 10 - b - 1
bとcは0では意味がないので、こっちのほうがいいのかな
0 <= a <= 9
1 <= b <= 8
1 <= c <= 10 - b - 1

805:デフォルトの名無しさん
17/12/03 13:12:58.80 QazTjKaA.net
>>776
すみません、間違えました、abcじゃなくてbcdですね
0 <= b <= 9
1 <= c <= 8
1 <= d <= 10 - b - 1

806:デフォルトの名無しさん
17/12/03 13:34:35.32 QazTjKaA.net
>>777
すみません、訂正の訂正です。dは10-c-1ですorz
0 <= b <= 9
1 <= c <= 8
1 <= d <= 10 - c - 1

807:デフォルトの名無しさん
17/12/03 17:47:25.64 Gq7SJlPX.net
Linked list 使うと楽そうな感じするな

808:デフォルトの名無しさん
17/12/03 18:26:21.27 QazTjKaA.net
>>770,771
せっかくなのでこのテストデータを作って検証してみましたが、全部合っていました、さすがですね
それにしても他人のプログラムって動かすことはできても、理解するのは困難ですね
>>772
Mathematicaは残念ながら持っていないので検証できませんでした
b = 8 c = 6 d = 2 a = 8901236745
b = 4 c = 3 d = 2 a = 0123784569
b = 4 c = 5 d = 2 a = 8123904567
b = 2 c = 5 d = 2 a = 0178234569
b = 0 c = 5 d = 3 a = 5670123489
b = 9 c = 3 d = 4 a = 3459016782
b = 5 c = 1 d = 7 a = 1253467890
b = 4 c = 6 d = 2 a = 8923014567
b = 6 c = 7 d = 2 a = 8901253467
b = 7 c = 5 d = 4 a = 5789016234
b = 4 c = 2 d = 2 a = 0123674589
b = 6 c = 4 d = 5 a = 4678950123
b = 7 c = 4 d = 4 a = 4789056123
b = 4 c = 5 d = 4 a = 6783901245
b = 8 c = 1 d = 5 a = 1238456790
b = 0 c = 2 d = 6 a = 2345670189
b = 3 c = 7 d = 1 a = 9120345678
b = 9 c = 4 d = 2 a = 4901256783
b = 8 c = 3 d = 6 a = 3456890712
b = 2 c = 5 d = 3 a = 0178923456

809:デフォルトの名無しさん
17/12/03 20:16:05.68 hV+xPFYR.net
>>768
ruby ワンライナー
URLﾘﾝｸ(ideone.com)

810:デフォルトの名無しさん
17/12/03 20:37:58.61 ucQfMVKf.net
>>780
お、あんた偉いな。
出題者でちゃんと答え合わせやってくれる人が稀でレスポンス薄いことが多いんだよ。
ちゃんと動いてよかったよ。

811:デフォルトの名無しさん
17/12/03 20:51:48.60 ucQfMVKf.net
>>774
自分は10超えてもらったら困るなぁ。
10っていうか、配列の要素数だな。
修正はそんなに難しくないけど。

812:デフォルトの名無しさん
17/12/03 23:00:23.32 QazTjKaA.net
>>781
1行プログラムすごいですね
さすがスクリプト言語
>>782
そうなんですか、でもいろいろな言語で回答されるから、出題者も大変かもしれないですね
私は問題の意味すらよく分からないことが多いので、もっぱら見る専ですが

813:デフォルトの名無しさん
17/12/04 02:36:11.06 iGjrIGoV.net
>>768
これでいいのかなあ？
Perl で書いた。
URLﾘﾝｸ(paiza.io)

814:デフォルトの名無しさん
17/12/04 06:47:18.44 Rc7ie/2s.net
>>785
b = 6、 c = 6、 d = 2
のときの動作がおかしいようです
8901452367
となるはずが
0167892345
となっており、移動する"678901"の並びが崩れてしまっています

815:785
17/12/05 02:32:35.27 LDxS5CId.net
>>786
問題勘違いしてました。素直にぐるぐる回すように修正しました。
URLﾘﾝｸ(paiza.io)
下の「入力」タブの方にスペース区切りで b, c, d の値を1行づつ並べて入れてから実行させると「出力」に結果が出ます。
とりあえず >>768 に書いてある値を入力にセットしたところ出力は同じになりました。

816:デフォルトの名無しさん
17/12/05 20:55:13.95 vy+ohhoY.net
>>787
入出力の確認をしただけですが、今回は問題ないようです
paiza使うと簡単に確認できて便利ですね

817:デフォルトの名無しさん
17/12/05 21:48:50.31 32LsMTj+.net
>>768
僕の頭だとどうしても右端から左端に行くときの動きがイメージできないので、解答締め切ったあとでもいいのでちょろっと教えてもらえると嬉しいです
上3つまではわかるけどそこから先がそもそも答えにたどり着けない……

818:デフォルトの名無しさん
17/12/05 23:47:14.24 ynbcQBXQ.net
>>789
ヒント。
(インデックス++)%配列の長さ
を繰り返すとどうなりますか？%は余剰デス。

819:デフォルトの名無しさん
17/12/06 00:25:49.49 gvvJf1Ph.net
>>789
上の３つまでは分かるということなので、３つめのcの値を一つずつ増やしてみると、こんな感じになります
"7"、"78"、"789"、"7890"と並ぶ数値が増えていっているのが分かると思います
c+d の合計は9が最高なので、これ以上cを増やすにはdを減らさなくてはなりません
b = 7 c = 1 d = 5 a = 1273456890
b = 7 c = 2 d = 5 a = 2378456901
b = 7 c = 3 d = 5 a = 3478956012
b = 7 c = 4 d = 5 a = 4578906123
今度はcを2に固定して、dの値を一つずつ増やすと、こんな感じになります
"78"の並びが一つずつ右へずれていっているのが分かると思います
c+d の合計は9が最高なので、これ以上dを増やすにはcを減らさなくてはなりません
b = 7 c = 2 d = 1 a = 0123456978
b = 7 c = 2 d = 2 a = 8123456907
b = 7 c = 2 d = 3 a = 7823456901
b = 7 c = 2 d = 4 a = 2783456901
b = 7 c = 2 d = 5 a = 2378456901
b = 7 c = 2 d = 6 a = 2347856901
b = 7 c = 2 d = 7 a = 2345786901
こんな感じで分かるでしょうか？私も自分の頭で考えると混乱しますw

820:デフォルトの名無しさん
17/12/06 00:31:57.65 +0bHqE6f.net
�

821:ｩ分の考え方は、配列の一部を切り取ってパッディングするんだけど、まず配列を回転させて切り取る第0インデックスを配列の最初に持って来る。すると、配列の後ろにはすでにパディングが終わった数列のができてる。で切り取って削除して尻尾にくっつける。で、さっき回した分を戻してやると完成。という方法で、>>770を解いた。

822:デフォルトの名無しさん
17/12/06 00:32:05.18 +0bHqE6f.net
自分の考え方は、配列の一部を切り取ってパッディングするんだけど、
まず配列を回転させて切り取る第0インデックスを配列の最初に持って来る。
すると、配列の後ろにはすでにパディングが終わった数列のができてる。
で切り取って削除して尻尾にくっつける。
で、さっき回した分を戻してやると完成。
という方法で、>>770を解いた。

823:デフォルトの名無しさん
17/12/06 00:32:32.42 +0bHqE6f.net
あちゃー
二重投稿になったある。
そんなに大事でもないのだけど。

824:デフォルトの名無しさん
17/12/06 00:38:22.54 +0bHqE6f.net
C++は比較的不自由な言語なので頭使うのは鍛えられるよ。
パーツが少ないのでほとんど自作しないといけない。

825:デフォルトの名無しさん
17/12/06 00:46:18.89 DokUEpLm.net
コンパイルが必要とかでスクリプト系の言語より
使うことに不便があるかもしれないが、
できるアプリが自由で高速・軽量！
ライブラリもＳＴＬやboost以外にも、探せば便利なのがいっぱい。

826:デフォルトの名無しさん
17/12/06 01:07:01.76 +0bHqE6f.net
よそのライブラリ使うとイデオンで動かないからなぁ。
まぁ、業務ではライセンスに合ったものを使えばいいよ。

827:デフォルトの名無しさん
17/12/06 02:01:03.49 QVT4XBLu.net
>>789
塊が移動すると考えれば良い。
例えば 7, 4, 2 だとすると、 789 と先頭の 0 が移動する塊だ。
で、ちょっと分かり易くするためにこの塊を伏せて * で書くとするとこうなる。
*123456***
それでこの塊を右に一つずらす。その時に1は食われて尻尾から吐き出される。
**234561**
移動量2なのでもう一回右にずらす。すると2が食われて尻尾から吐き出される。
***345612*
それでは * から 7890 に戻してみよう。
8903456127
できあがり。

828:デフォルトの名無しさん
17/12/06 02:58:25.26 +0bHqE6f.net
>>793
これ間違ってるな。
最後尻尾にくっつけるんじゃなくて、適所にインサートだった。
で戻す。
書いた日から時間たってるからすまん。

829:デフォルトの名無しさん
17/12/06 08:41:18.93 obBhCrma.net
>>790
>>791
>>798
ありがとうございます
やっと理解できました
がんばって解いてみます

830:デフォルトの名無しさん
17/12/06 17:48:37.81 T95E7suL.net
log2だけど、これ使うといいらしいよ。分割統治法のBSA法というやつらしい。
2個ずつの積を繰り返すことで計算回数が減らせる。

{{1, a0}, {0, r}}*{{1, a1}, {0, r}}を[ [○,P], [○,Q] ]とおくと
r*P/Qはa0 + a1/rであり、
{{1, a0}, {0, r}}*{{1, a1}, {0, r}}*{{1, a2}, {0, r}}を[ [○,P], [○,Q] ]とおくと
r*P/Qはa0 + a1/r + a2/r^2。
上の式を実際に計算してみる。
URLﾘﾝｸ(www.wolframalpha.com)

831:デフォルトの名無しさん
17/12/06 21:09:44.45 N+lgs3o9.net
>>801
普通のlog 2の計算には使えても
>>680 には使えないかと

832:デフォルトの名無しさん
17/12/07 00:07:53.25 YC8MWngg.net
>>731 から更に縮まりました
もう誰も興味が無いかも知れませんが
コードいる？

833:デフォルトの名無しさん
17/12/07 09:24:10.80 XIHsqoOR.net
>>801でできるはずだ。やってみる

834:デフォルトの名無しさん
17/12/07 10:09:03.65 DjN/Fbox.net
n=10^10で19.28秒にまで縮めた
n=10^13くらいまでなら�

次ページ