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

**デフォルトの名無しさん** · 2021/06/19(土) 00:02:57.84

プログラミングのお題スレです。

【出題と回答例】
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/

宿題は宿題スレがあるのでそちらへ。

※前スレ
プログラミングのお題スレ Part19
https://mevius.5ch.net/test/read.cgi/tech/1606662245/

**デフォルトの名無しさん** · 2021/11/12(金) 20:37:34.22

”2021”というフォルダがある
"2020"というフォルダもある
　：
"1950"というフォルダもある
そしてそれぞれに"01"～"12"のサブフォルダがあって
そしてそのそれぞれに"01"～"31"というサブフォルダがあって
更に"00"～"23"という・・・・

つまり、１時間ごとのフォルダがあって
そこには多くのファイルが詰まっている

そのファイルのタイムスタンプに変化があったら
「どこのファイルが書き換えられたぞ」と検出するプログラムを作りたい
片っ端からタイムスタンプを調べるしかないのか
これはファイルシステムの問題でしょうが
なんとか、前回調べてから変化があったか、短時間で検出できる方法はありませんか

ファイルに変化があったら、フォルダのタイムスタンプを変えてくれればいいんだけどね

**デフォルトの名無しさん** · 2021/11/12(金) 21:45:17.50

そういう検出機能はOSかハードウェアかについてなかったか　忘れたが

**デフォルトの名無しさん** · 2021/11/12(金) 21:49:29.72

NTFSの USN とは？
USN　ジャーナルとは
NTFSのUSN機能とはNTFSボリュームに対しての変更の記録です。
Cドライブ、Dドライブ、Eドライブごとなど、ドライブごとに加えられた変更を記録したものです。

USNとはUpdate Sequence Number Journalの略称です。
主に以下の目的で使用されます。
・FRS(File Replication Service)によるファイル複製時の差分洗い出し。
・スナップショットバックアップ時の差分洗い出し。
・ウィルス対策ソフトによるファイル検査。
・インデックスサービス向け。

http://tooljp.com/qa/what-is-USN-36BE.html

**デフォルトの名無しさん** · 2021/11/12(金) 22:31:42.38

Ruby では、文字列の中から、URL を抜き出す関数が非推奨になった。
頻繁に規格が変わるからかな？

正規表現でも多分、何万行も掛かる。
作るのに10年ぐらい掛かるので、誰も作らない

今あるものでも、正しく判定されるかどうか、分からない

**デフォルトの名無しさん** · 2021/11/12(金) 22:39:57.54

>>430
5ch のバグで、

そのURL をプログラム板全体で見た時は、&q とつながっているのに、

個別のスレッドで見ると、&;q と、&の直後に、;が入って、
リンクが&までになってしまう

**デフォルトの名無しさん** · 2021/11/13(土) 02:42:42.42

>>433
ドライブごとに加えられた変更がログに保存されているとは知りませんでした
データドライブに分けているので、ログから探すのも速いかもしれません
ありがとうございました

**デフォルトの名無しさん** · 2021/11/15(月) 00:16:09.47

お題: JSON形式のテキストをXML形式のテキストに変換しなさい

**デフォルトの名無しさん** · 2021/11/15(月) 21:06:19.54

>>437 C++
https://ideone.com/vkSCmj

**デフォルトの名無しさん** · 2021/11/23(火) 21:56:09.59

age

**デフォルトの名無しさん** · 2021/11/24(水) 11:45:59.40

ideone って boost は include 出来るのか
codepad はずっと落ちてるな

**ハノン** ◆QZaw55cn4c · 2021/11/24(水) 20:41:27.64

>>440
もう codepad はおしまい、codepad にアップした私のコードも永遠に失われました…

96 · 2021/12/02(木) 04:03:51.64

元のソースをローカルに保存しておかなかったのか

**デフォルトの名無しさん** · 2021/12/02(木) 12:54:02.17

短いもんだから直接打ち込んじゃったんじゃないかな。
そして消滅。

**デフォルトの名無しさん** · 2021/12/02(木) 16:48:13.78

>>442
デスクトップで書いたものはローカルのどこかに残っている（と思う）
ノートで書いたものはノートPCごとどこかにいってしまった…
いや、まあ、最初から書き捨てでしたけど、惜しいと思うものも多少ありますね

96 · 2021/12/02(木) 21:30:16.60

いま
http://codepad.org/
復活してない？

サンプルURLとして前々ｽﾚの過去の回答

プログラミングのお題スレ Part18
http://mevius.5ch.net/test/read.cgi/tech/1594702426/174
>174 名前：デフォルトの名無しさん[] 投稿日：2020/08/04(火) 19:24:17.96 ID:lImhU2CT
> >>173
> あーアホだったw
> タブ消去の仕様だとPythonだと分けわからんことになるなw
> ↓貼り直した
> http://codepad.org/BWRzr8fW
> ※会員未登録の為、一定期間後に自動消去

も表示されたし。

保持しておきたいcodeがあれば今のうちにsaveするか、
Internet archiveに記録しとくのも手だよ。

**デフォルトの名無しさん** · 2021/12/08(水) 14:07:00.92

お題が出なくなったので全く実用性がないと思われる頭の体操的なやつを一つ。

お題: 与えられた数値の配列の 1 になっているビットをカウントする。

例えばC言語で unsigned char の配列とか使って、処理する関数はポインタで渡された所から何要素みたいに指定されるようにする。。

それで例えば 0xde, 0x96 (11011110, 10010110) なら答えは 10。

**デフォルトの名無しさん** · 2021/12/08(水) 17:14:42.64

.replace("0", "").length

**デフォルトの名無しさん** · 2021/12/08(水) 17:44:20.69

お題: 日付が与えられるので4月ならtrueを出力しなさい

< 2021/04/01
> true

< 04/01
> true

< 2021/04
> true

< 2021/05/01
> false

**デフォルトの名無しさん** · 2021/12/08(水) 19:09:00.93

>>446
最適なコードが環境に大きく依存するからなあ
シフトの速度とか専用命令有無とか

あまり良いお題じゃないし荒れそう

**デフォルトの名無しさん** · 2021/12/09(木) 05:41:43.04

英式　1/4/2021
米式　4/1/2021

**デフォルトの名無しさん** · 2021/12/09(木) 08:53:58.33

日 21/04/01
米 04/01/21
欧 01/04/21
年2桁で入り乱れたら更にわけわかめ

**デフォルトの名無しさん** · 2021/12/09(木) 13:38:41.08

逆に、年2桁の時は、
その数値で年月日が特定できるか問題があった。

02 12 29 --> 12年02月29日と 29年02月12日の2通りがある
02 13 29 --> 13年02月29日はないので、29年02月13日のみ

**デフォルトの名無しさん** · 2021/12/09(木) 13:44:59.59

>>452
特定できるかどうかだけならいいが……
反例を間違った　02 12 29 3通りある

**デフォルトの名無しさん** · 2021/12/09(木) 14:42:16.44

そういえば食品の賞味期限がどう読んだらいいか分からなかったことがあったなあ。

11/02/13

2013年2月11日なのか、それとも2013年11月2日なのか、まさか2011年2月13日だったりしたらどうしよう、みたいな。
普通に輸入されているものなら日本語ラベルがあるから分かるだろうけど、友達が現地のお土産屋で買ってきた
お菓子みたいなものだとそのまんまだったりする(まあ、アメリカかアメリカ以外かで考えればだいたいは大丈夫
なのかも知れないが)。
日本の元号みたいな独自年号まで使われているのが混在していたら余計に分からないよな。今の日本では
あまりないだろうけど昭和の頃はあったような気がするなあ。

**デフォルトの名無しさん** · 2021/12/09(木) 22:53:27.33

>>446 Ruby

p [0b1110, 0b1010111].reduce(0){|s, n| s + n.digits(2).count(1)}

# => 8

**デフォルトの名無しさん** · 2021/12/09(木) 23:18:40.40

>>446
https://ideone.com/qmIeDU

**デフォルトの名無しさん** · 2021/12/10(金) 00:33:54.80

>>446
C
https://paiza.io/projects/YQ8ORfELqUt1xkKqmysViQ

入力を配列に入れてからカウントさせている。(uint32_t 型の配列)
カウント部分は>>456に似てる。しかし最初に引くのは思いつかなかった。

**デフォルトの名無しさん** · 2021/12/10(金) 01:03:30.79

>>448
Perl
https://paiza.io/projects/o4VTdU_kK-S8T4G5HvkIRw

年月日順で年がある場合は必ず4桁でなければならない。
4だけでも4月に判定されるが、まあいいか。

**デフォルトの名無しさん** · 2021/12/10(金) 09:12:55.36

>>446

Haskell

sumCntBits = id
. length
. filter odd
. ( >>= ( takeWhile ( /= 0 ) . iterate ( flip div 2 ) ) )

main = do
print $ sumCntBits [ 0xde, 0x96 ]
print $ sumCntBits [ 0x12345, 0x6789a,0xbcdef ]

----
10
32

**デフォルトの名無しさん** · 2021/12/11(土) 11:47:37.10

>>446 rust
https://ideone.com/9yFeU3
fn main() {
let f = |a: &[u32]| a.iter().map(|n| n.count_ones()).sum::<u32>();
println!("{}", f(&[0xde, 0x96]));
}

**デフォルトの名無しさん** · 2021/12/11(土) 20:38:05.35

>>446
Kotlin
https://paiza.io/projects/7hQyycMNf5nuV5HcBg6iQg

普通に作るのは出尽くした感があるのでちょっと変わったやり方にした。
入力から Int の List を作り、それを 1 ビットづつの Boolean のリスト(というか Iterator) にしてから true のみをカウントしている。

**デフォルトの名無しさん** · 2021/12/11(土) 22:30:09.92

>>446 octave
https://ideone.com/BNjv3I
f = @(a) sum(dec2bin(a)(:) - '0');
f([0xde 0x96])

**デフォルトの名無しさん** · 2021/12/14(火) 17:40:23.45

もうすぐ、2022年
[お題] 2022は"x3y1数"(造語)?

以下の二つを満たす正の整数を"x3y1数"と呼ぶ
・各桁の数値が、二種類のみの数字からなる
・上の二数の個数比は 3:1

　該当例：1112, 2212, 2022, 32222223, 999999999888
　ダメな例：2213(種類), 4444(種類), 33232(個数比), 0222(先頭ゼロ)

　整数A,Bが与えられる。A以上B以下の"x3y1数"はいくつあるか？
　制約: 0 < A <= B <= 10^18

1) 2923 3311 --> 8
　　該当は [2999, 3000, 3033, 3111, 3133, 3222, 3233, 3303]
2) 2021 77977796 --> 2022
3) 3000 6666566566566555 --> ?
4) 999999999889 1000000000000110 --> ?
5) 1 1000000000000000000 --> ?

※実行時間：上の5問をクエリーとし全部を2秒程度

**デフォルトの名無しさん** · 2021/12/16(木) 03:59:00.84

お題：自分用double-double演算ライブラリ
最低限、通常のdoubleとの相互変換は可能であること。それに加えて、
1)加減算
2)加減算 + 乗算
3)四則演算
数字が大きいもの程上級者向けです。

演算子のオーバーロードなどは任意とします。

**デフォルトの名無しさん** · 2021/12/16(木) 07:13:07.04

>>464
多倍長演算ライブラリ、のことですか？

**デフォルトの名無しさん** · 2021/12/16(木) 07:37:15.94

>>465
アホ

**デフォルトの名無しさん** · 2021/12/16(木) 20:36:30.04

>>463 c
https://ideone.com/SjWUYZ
・数字を数えて判定
・範囲内の全ての整数をチェック
・想像以上に遅くてダメだった

>>463 ruby
https://ideone.com/JYQPIF

**デフォルトの名無しさん** · 2021/12/16(木) 20:39:48.12

>>466
double の演算を自分で実装するという意味ですか？
sum(double, double)
diff(double, double)
mul(double, double)
div(double, double)
を自分で実装する、という話でいいですか？
あと double のフォーマットは IEEE754　でいいですか？

**デフォルトの名無しさん** · 2021/12/16(木) 20:57:05.50

問題が説明不足では？

**デフォルトの名無しさん** · 2021/12/16(木) 21:53:18.15

お題: テキストを読み込みそれをクリスマスツリーにして出力しなさい
クリスマスツリーに見えれば形は自由とする

入力

本日は良いお日柄ですね

出力

___本
__日は
_良いお
日柄です
___ね

**デフォルトの名無しさん** · 2021/12/16(木) 22:32:14.13

>>470
文字コードは何を仮定すればいいのですか？

**デフォルトの名無しさん** · 2021/12/16(木) 22:34:05.93

>>471
UTF-8
日本語の扱いが難しい言語では英語のみの対応も良しとする

**デフォルトの名無しさん** · 2021/12/17(金) 00:19:35.20

>>470 octave
https://ideone.com/RseGCJ

**デフォルトの名無しさん** · 2021/12/17(金) 04:20:32.14

二種類のみの数字からなり個数比は 3:1
引数の範囲は 1-10^18 = 1001-9999999999998888(16桁)
以下の範囲に限られる
1000-9998
1000 0001-9999 9988
1000 0000 0011-9999 9999 9888
1000 0000 0000 0111-9999 9999 9999 8888

「二種類のみの数字からなる」を計算式で判定する方法ある？

**デフォルトの名無しさん** · 2021/12/17(金) 05:36:58.41

1([\d&&[^1]])\1{2} 最上位桁が比1
111[\d&&[^1]],11[\d&&[^1]]1,1[\d&&[^1]]11 最上位桁が比3
一般化 (\d)(?!\1)(\d)\2{2}|(\d)\1{2}(?!\1)\d|(\d)\1(?!\1)\d\1|(\d)(?!\1)\d\1{2}
4桁ならこれでもいいけど8桁以上になると複雑化するし
地道に数えるより 4の倍数桁,数字2種,比率1:3 のルールで生成する方が速そう

**463** · 2021/12/17(金) 16:22:06.93

>>463
　https://ideone.com/xTDtME

　想定解は、事前に4,8,12,16桁の"x3y1"数を全列挙して作っておく。
　プログラミング的には、各言語の順列や組合せを使って、作れるだろう。
　(想定解例では組合せは2ベキとpopcountから作っている)

　「それは、全列挙数が小さいとわかっているからでは..？」に対して
　プログラムで出すのなら、雑に最も大きい16桁が4つあるとして計算
　　10P2 * 16C4 * 4 < 70万なので、全列挙可能
　まじめに計算すると 10P2 * (16C4 + 12C3 + 8C4 + 4C1) * 9 /10 = 167,832

　列挙済みならば、クエリー5件程度なら、16.7万*5 チェックで間に合う。
　ちゃんとやるなら、ソートして二分探索すれば、数千単位のクエリーに対応できる。
　(想定解例では後者でやっている)

**デフォルトの名無しさん** · 2021/12/17(金) 20:17:58.00

>>468
>>466

**467** · 2021/12/17(金) 20:31:12.28

>>463 c
https://ideone.com/bmYThw
・Ruby版の移植
・組み合わせの列挙方法は丸パクリ
・Ralph William Gosper Jr. 氏に感謝

**デフォルトの名無しさん** · 2021/12/17(金) 23:34:29.99

>>463 c
https://ideone.com/oPUphG
・>>478から若干の整理
・組み合わせ列挙用バッファ廃止

**デフォルトの名無しさん** · 2021/12/18(土) 16:20:56.22

>>464 C++
https://ideone.com/tkSsy4

**464** · 2021/12/18(土) 16:29:44.66

>>465
惜しい

>>468
IEEE754の倍精度(binary64)を整数演算で実装するのではありません。
binary64を二つ使って、上位53ビットと下位53ビットとで106ビットの浮動小数に見立てたものが
double-double演算です。

Wikipediaの「四倍精度浮動小数点数」の項に少しだけ載ってますです。

**デフォルトの名無しさん** · 2021/12/18(土) 16:50:56.06

> Wikipediaの「四倍精度浮動小数点数」の項に少しだけ載ってますです。
一般的な用語じゃないんだから初めからこれ書いとけよ

**465,468** · 2021/12/19(日) 00:01:36.35

>>481
I see.

**デフォルトの名無しさん** · 2021/12/19(日) 21:10:50.32

>>448 octave
https://ideone.com/2NglYm

**デフォルトの名無しさん** · 2021/12/21(火) 19:32:25.39

128ビットあるのに106ビットしか使わんの？
もったいなくね？

**デフォルトの名無しさん** · 2021/12/21(火) 19:47:17.61

素人はだまってろ

**デフォルトの名無しさん** · 2021/12/21(火) 21:17:48.02

>>485
もったいなよ
ただ、既存のdouble計算リソースが使えるという利点がある

**デフォルトの名無しさん** · 2021/12/22(水) 04:27:38.39

double-doubleはFMAがFMAとして役立つ数少ない用途だな
積和じゃなくて3個の和のfused命令も欲しくなる

**デフォルトの名無しさん** · 2021/12/23(木) 07:32:47.64

お題: 1つの整数から規則性のある複数の整数を生成せよ
生成される整数は再現性がなければならない

**デフォルトの名無しさん** · 2021/12/23(木) 07:37:01.70

>>489

function f: Integer -> Integer{
return 0;
}

**デフォルトの名無しさん** · 2021/12/23(木) 09:01:28.76

>>489 Ruby
def sequence( seed, number )
srand( seed )
Array.new( number ){ rand(100) }
end

p sequence( 123, 10 ) #=> [66, 92, 98, 17, 83, 57, 86, 97, 96, 47]
p sequence( 123, 10 ) #=> [66, 92, 98, 17, 83, 57, 86, 97, 96, 47]

**デフォルトの名無しさん** · 2021/12/23(木) 19:26:32.49

>>486 ocaml
https://ideone.com/NzF5f2
let f =
let rec fib = function
0 -> 0 | 1 -> 1 | n -> fib (n - 1) + fib (n - 2)
and aux acc = function
-1 -> acc | m -> aux (fib m :: acc) (m - 1)
in aux []

**デフォルトの名無しさん** · 2021/12/23(木) 19:27:19.05

>>489 ocaml
https://ideone.com/NzF5f2
let f =
let rec fib = function
0 -> 0 | 1 -> 1 | n -> fib (n - 1) + fib (n - 2)
and aux acc = function
-1 -> acc | m -> aux (fib m :: acc) (m - 1)
in aux []

**デフォルトの名無しさん** · 2021/12/24(金) 17:16:46.67

>>488
Juliaのhypot()でもFMA使ってますです

**デフォルトの名無しさん** · 2021/12/24(金) 22:51:57.83

>>494
ただ積と和を1命令にして高速化しただけの積和
の効果だけじゃなくて
融合(fused)の効果が効く用途の話

**デフォルトの名無しさん** · 2021/12/25(土) 03:52:36.33

>>489
Kotlin
https://paiza.io/projects/xmMY6y8BGb8zlhn5QEmKvQ

何も考えずにただ Random 使っただけ。

**デフォルトの名無しさん** · 2021/12/27(月) 20:25:00.12

[お題] 平均が2022な素数数列

　5000以下のあい異なる素数で、加算平均がぴったり 2022 の数列を作る。
　数列の項数(要素数)を最大化する、最大はいくつか。
　最大数と数列を表示する。
　
※解答例(もちろん4以上がある)
　4
　[1747, 2099, 2113, 2129]
※実行時間は素数生成を含めて、4秒以内
　最大な数列は複数通りあると思うので、一例のみで

**デフォルトの名無しさん** · 2021/12/29(水) 15:58:08.20

答えは595個？
計算機+理詰めで595個っぽいけど

**デフォルトの名無しさん** · 2021/12/29(水) 16:38:34.79

「Log4j」2.17.0にもリモートコード実行の脆弱性

**497** · 2021/12/29(水) 20:37:12.14

>>497 c++
　https://ideone.com/UBbtWd

　"素数-2022"で適当に最大化DPすれば、合計0相当の所に答えが……。
　個人的には他の復元方法を見たかった。
　pythonは遅いのであきらめた(高速化方法を知らない)

>>498
595個でした。手作業でできるレベルなら……

**デフォルトの名無しさん** · 2021/12/29(水) 22:36:34.62

オレがやったのは
p[n]をn番目の素数、
P[n]をn番目までの素数の集合、
s[n]をP[n]の和
a=2022として
まず素数のn元集合で平均が1番小さくなるのはP[n]でその平均値s[n]/nは単調増加だからs[n]/n>aの時元数n以上の解はない
s[596]/596>aは計算機で確認
なので595元集合で解があればそれが最大
s[595]/595<aなのでP[595]はダメ
s[596]-595a = 1205525 - 1203090 = 2435は素数ではないのでP[597]から一個消すのはダメ
s[597]-595a = 1209898 - 1203090 = 6808 はp597=4373以下の素数2459,4349の和で表すことができる
よってP[597]＼{2459,4349}の和は595aとなる
完全に全自動で探索するプログラムも作れそうだけど答え出たらもういいかなと手が止まってしまった

**デフォルトの名無しさん** · 2021/12/30(木) 12:24:49.52

何個取り除いたら平均2022にできるか考えると
確か74個だったけな？
JavaScriptで1秒もかかんなかったか

**502** · 2021/12/30(木) 12:30:19.06

表示時間除いたら1秒かかってないな
探索はほぼ一回でボトムまで到達して
発見された

**502** · 2021/12/30(木) 12:38:50.05

可能と不可能が極端に分かれている問題なので
事前のブランチカットだけが重要なヘンなお題ｗ
>>352はまだ未解決

**デフォルトの名無しさん** · 2021/12/30(木) 20:10:02.90

>>502
最大の個数を求めよやろ？

**デフォルトの名無しさん** · 2021/12/30(木) 20:24:56.63

千葉興業銀行、4月から副業解禁　県内地銀初

南都銀行、4月から行員の副業制度導入　ウェブ制作など

荘内銀、行員の副業・兼業解禁

フィデアHD、副業・兼業制度を導入

横浜銀行、10月から従業員の副業・兼業解禁

鹿児島銀、副業解禁を検討　九州FGと肥後銀は10月導入

肥後銀行が副業制度導入へ　多様な働き方認める　10月から

**デフォルトの名無しさん** · 2021/12/31(金) 15:04:12.13

>>497
haskell
https://ideone.com/GLMXRV
>>501のアルゴリズムを自動化してみた
すげー簡単なところでどハマりして半日かかった
まだ解なしの場合とかの動作チェックとかしてないけどもうどうでもいい

**デフォルトの名無しさん** · 2021/12/31(金) 15:09:29.33

ちなみに出力形式は
(最大個数、最大素数の通し番号、最大素数までの間での素数で除外する素数のリスト)

try 67%5
→(5,8,[2,3,5])
は最初の個数8個[2,3,5,7,11,13,17,19]から[2,3,5]を抜いた[7,11,13,17,19]の5個が平均が67/5になる素数のリストの一つ
長さ
6以上はない

**デフォルトの名無しさん** · 2022/01/08(土) 11:42:09.36

お題：
xをゼロ以上の浮動小数点数として
2^floor(log2(x))
の計算。ただし、x == 0 の場合はゼロとする。

**デフォルトの名無しさん** · 2022/01/08(土) 18:45:49.69

>>481
＞binary64を二つ使って、上位53ビットと下位53ビットとで106ビットの浮動小数に見立てたものが
＞double-double演算です。
現在検討中ですが、binary64 中には仮数部に使用できるビット幅は 52 bits しかありません。つまりケチビット表現です
53+53 とのことですが、実際には 53 + 52 = 105 しか格納できないのではないでしょうか？

**464** · 2022/01/08(土) 21:07:22.52

>>510
上位の±1/2ulp相当が下位になります

**デフォルトの名無しさん** · 2022/01/09(日) 01:28:46.61

>>509
perlでワンライナー。入力は標準入力からする。

perl -MPOSIX -ne 'chomp;$n=$_?2**floor(log($_)/log(2)):0;print "$n\n"'

でも、こんなので良いの？自分ではほとんど何も考えてないんだが。
(log2()がないからlog(n)/log(2)でやるって所ぐらいしか工夫がない)

**デフォルトの名無しさん** · 2022/01/09(日) 07:56:15.63

>>509 C++
環境+コンパイルオプション依存 little endian, double = 64bit, long double = 128bit

double fl2( double x )
{
*( (uint64_t*) &x ) &= 0xFFF0000000000000LLU;
return x;
}

long double fl2( long double x )
{
*( (__uint128_t*) &x ) &= *( (__uint128_t*)"\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\xFF\xFF" );
return x;
}

**デフォルトの名無しさん** · 2022/01/09(日) 07:59:19.16

2進表記した時に先頭のビット以外を0にすればいいだけだがワンライナーで書ける気がしない

**デフォルトの名無しさん** · 2022/01/09(日) 08:00:08.09

先頭のビットというより「(存在するなら)0じゃない一番位の大きいビット」だな

**デフォルトの名無しさん** · 2022/01/09(日) 08:27:45.00

ビットいじる方法だと、非正規化数が0に、NaNが無限大になる
そういえば、無限大やNaNはどうすればいいのだろうか

**509** · 2022/01/09(日) 09:16:00.82

Python
でかい値だとうまくいかないやつ

def foo(x):
_c = x * float(2 ** 52 + 1)
_xh = c - (c - x)
_if xh > x: return(xh * 0.5)
_return(xh)

>>516
深く考えてませんでしたｗ

**509** · 2022/01/09(日) 09:53:47.91

>>512
log2(n)をlog(n)/log(2)で近似した際の「誤差」
(ぴったり整数値になって欲しいのにそれよりも数ulp小さい値になったとか、
数ulp大きくて、ぎりぎりfloor()で切り捨てられる筈の値が1大きくなったとか)
を補償するコードが欲しい。

**デフォルトの名無しさん** · 2022/01/09(日) 15:38:39.37

そこは性能とトレードオフになるしかない気はする
例えばfloor( log 32 / log2 ) = 5が使用上の動作だけど
log 32/log2 =4.9999999999978 × 10^0が返されて答えが4になってしまうのを避けるならおそらく大きすぎる場合は無視していいから
res' = floor( log x / log2 )
if (2^res'* 1.5) < x
// (1.5倍でも届かないなら不当な丸め誤差が出たと判断)
then res = res' + 1
else res = res'
とかちょっと汚い書き方するしかない希ガス

**デフォルトの名無しさん** · 2022/01/09(日) 21:50:51.38

なんか変なこと書いた

res' = floor( log x / log2 )
if 2^res' < x
then res = res' + 1
else res = res'

やな
res'の算出で丸め誤差は-1までと仮定して補正する
しかしもちろんlog(x)とか2^nハード的にFPUとかで高速にやってくれてしかも整数演算は誤差なしでやってくれる前提
この辺は高級言語プログラマレベルの話でなんとかなるもんではない
やるならアセンブリ言語レベルでやるしかない

**デフォルトの名無しさん** · 2022/01/10(月) 00:03:40.12

>>509 C
https://ideone.com/5kDuzA
非正規化数、NaN、無限大とかはそのまま返すようにした
やり方は>>513と変わらない

**509** · 2022/01/10(月) 00:53:57.15

いまどきのコンパイラなら、frexp()やldexp()をいい塩梅に最適化してくれるのだろうか？

from math import frexp, ldexp
def foo(x):
_m, e = frexp(x)
_if m == 0: return 0.0
_return ldexp(0.5, e)

**512** · 2022/01/10(月) 01:47:47.54

>>509
またPerlでワンライナー。

perl -ne 'chomp;if($_<=0){print"0\n"}else{for(my$i=0;;$i++){if((1<<$i)>$_){print 1<<($i-1),"\n";last}}}'

今度は計算している内容から考えて結果が同じになるようにした。浮動小数点演算をしていない。
また整数値が何ビットであるかも考慮しておらず、Perlの整数が32bitだった場合は2^32以上の値を入力されたら多分うまく動かない。
当然64bitだったら2^64以上の値の入力で多分うまく動かない。

**509** · 2022/01/10(月) 02:23:50.20

>>522の修正
from math import frexp, ldexp
def foo(x): return 0.0 if x == 0 else ldexp(0.5, frexp(x)[1])

**デフォルトの名無しさん** · 2022/01/10(月) 03:27:40.49

>>509
またまたPerlでワンライナー。

perl -MPOSIX -ne 'chomp;if($_==0.0){print"0\n"}else{print ldexp(0.5,(frexp($_))[1]),"\n"}'

これは>>524の真似(ていうかやってること同じ)。

**デフォルトの名無しさん** · 2022/01/10(月) 12:43:53.01

>>511
それは答えになっていないかと
質問を変えます。下位側の指数部も意味を持つようにインプリメントするべきでしょうか？

**464** · 2022/01/11(火) 02:38:04.94

>>526
先人の実装例だと、
上位 + 下位 = double doubleの数値
という事になってますね(上位側の指数が決まると、下位側の指数も決まる)。
tps://na-inet.jp/na/qd_ja.pdf

勿論、そう実装しないのもあり。

**デフォルトの名無しさん** · 2022/01/11(火) 03:06:24.10

>>527
となると、
>>510
＞binary64 中には仮数部に使用できるビット幅は 52 bits しかありません

よって下位側指数部無視なら 53bit + 52 bit = 105bit の実装となりますが？
下位側指数部有意ならば、下位側にもケチビットを適用できますが、今度は仮数部が 106 ビットとはいいきれなくなりますね（数によって変わる）

**デフォルトの名無しさん** · 2022/01/11(火) 03:08:13.70

>>527
失礼 pdf が紹介されていることを見落としていました、精査します、紹介ありがとうございます

**デフォルトの名無しさん** · 2022/01/30(日) 18:02:46.10

お題: 1より小さい実数を1以上2より下にせよ

< 0.123
> 1.23

< 0.0000123
> 1.23