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

プログラミングのお題スレ Part22 (788ﾚｽ)
上下前次 1-新
通常表示 512ﾊﾞｲﾄ分割ﾚｽ栞

ﾘﾛｰﾄﾞ規制です｡10分ほどで解除するので､他のﾌﾞﾗｳｻﾞへ避難してください｡

683: 9 [sage] 2025/03/18(火) 16:07:07.28 ID:GYPHuJM6

>>682
実行結果 (改行数を減らすため適度につなげてます)
$ perl 22_680_hex_substr_1.pl
38d7ea4c67fff 999999999999999
1 .. 0x3897994967999
Tue Mar 18 09:15:31 2025
1, 0x1　　 2, 0x2　　 3, 0x3　　 4, 0x4　　 5, 0x5　　 6, 0x6　　 7, 0x7　　 8, 0x8　　 9, 0x9
357440, 0x57440　　 357441, 0x57441　　 357442, 0x57442　　 357443, 0x57443　　 357444, 0x57444
357445, 0x57445　　 357446, 0x57446　　 357447, 0x57447　　 357448, 0x57448　　 357449, 0x57449
1079653, 0x107965　　 1081713, 0x108171　　 1122966, 0x112296　　 1123079, 0x112307　　 1123080, 0x112308
2246166, 0x224616　　 3369253, 0x336925　　 3371313, 0x337131　　 3412566, 0x341256
4494393, 0x449439　　 4494400, 0x449440　　 4535653, 0x453565
5658739, 0x565873　　 5658740, 0x565874　　 5660793, 0x566079　　 5660800, 0x566080　　 5702166, 0x570216
6783879, 0x678387　　 6783880, 0x678388　　 6784000, 0x678400
6825253, 0x682525　　 7948339, 0x794833　　 7948340, 0x794834　　 7950393, 0x795039　　 7950400, 0x795040
2182104640, 0x82104640　　 2182104641, 0x82104641　　 2182104642, 0x82104642　　 2182104643, 0x82104643　　 2182104644, 0x82104644
2182104645, 0x82104645　　 2182104646, 0x82104646　　 2182104647, 0x82104647　　 2182104648, 0x82104648　　 2182104649, 0x82104649
1263629042727, 0x12636290427　　 1307655353654, 0x13076553536　　 2573583194436, 0x25735831944　　 2617616245848, 0x26176162458
3330782168640, 0x30782168640　　 3330782168641, 0x30782168641　　 3330782168642, 0x30782168642　　 3330782168643, 0x30782168643
3330782168644, 0x30782168644　　 3330782168645, 0x30782168645　　 3330782168646, 0x30782168646　　 3330782168647, 0x30782168647
3330782168648, 0x30782168648　　 3330782168649, 0x30782168649　　 3883544086630, 0x38835440866　　 3927569962533, 0x39275699625
3927570397557, 0x39275703975

Core i7-8559U で6時間ほど実行してここまで高々13桁。
やはり想定通り気の利いた高速解放が要りますﾃﾍﾍﾟﾛ。

http://mevius.5ch.net/test/read.cgi/tech/1691038333/683

697: デフォルトの名無しさん [sage] 2025/03/21(金) 12:56:29.73 ID:CgJbZEAu

>>680 　Rust
fn odai_680() -> Vec<i128> {
　let mut answer = vec![0];
　let n_max = (0..).find(|&n| pow16(n + 1) > pow10(36)).unwrap();
　for s in (0..).take_while(|&s| pow16(n_max) >= pow10(n_max + s + 1)) {
　　let c = (0..=n_max).map(|i| pow16(i) - pow10(i + s)).collect::<Vec<_>>();
　　let rmax = c.iter().scan(0, |s, &c| { *s += if c > 0 { c * 9 } else { 0 }; Some(*s) }).collect::<Vec<_>>();
　　let rmin = c.iter().scan(0, |s, &c| { *s += if c < 0 { c * 9 } else { 0 }; Some(*s) }).collect::<Vec<_>>();
　　let (mut i, mut n, mut d, mut ct) = (0, 1, vec![0; c.len()], vec![0; c.len() + 1]);
　　loop {
　　　d[i] += 1;
　　　if d[i] < 10 {
　　　　let m = pow10(n as u32 + s); ct[i] = c[i] * d[i] + ct[i+1];
　　　　if i == 0 {
　　　　　if ct[0] >= 0 && ct[0] % m < pow10(s) { answer.push(d.iter().take(n).rev().fold(0, |sum, &d| sum * 16 + d)) }
　　　　} else {
　　　　　let (max, min) = (ct[i] + rmax[i-1], ct[i] + rmin[i-1]);
　　　　　if max >= 0 && (max - min > m || pow10(s) > min % m || min % m > max % m) { i -= 1; }
　　　　}
　　　} else { d[i] = -1; i += 1; if i == n { if n == d.len() { break; } n += 1; } }
　　}
　}
　answer.sort(); answer
}

http://mevius.5ch.net/test/read.cgi/tech/1691038333/697

750: デフォルトの名無しさん [sage] 2025/04/11(金) 07:38:20.09 ID:oaeJuxMT

>>738 に手を加えて10倍速くしてみた
fn solve(n: usize, limit: usize) -> Vec<usize> {
　let mut answer = Vec::new();
　let mut pnt = generate_primes(n).windows(2).rev().map(|s| (s[1], s[0], 0, 0)).collect::<Vec<_>>();
　let (mut ci, mut cn, mut ct) = (0, n, 1_usize);
　'advance: loop {
　　pnt[ci..].iter_mut().for_each(|(_p, _q, n, t)| (*n, *t) = (cn, ct));
　　if cn & 1 == 0 && ct.leading_zeros() >= (cn >> 1) as u32 {
　　　ct <<= cn >> 1; if ct <= limit { answer.push(ct); }
　　}
　　'back: for (i, (p, q, n, t)) in pnt.iter_mut().enumerate().rev() {
　　　'again: loop {
　　　　if *n < *p { continue 'back; }
　　　　*n -= *p; *t *= *p;
　　　　if *n ==1 || *t > limit { continue 'back; }
　　　　if *n == 0 { answer.push(*t); continue 'back; }
　　　　if *q > 3 {
　　　　　let mut tt = *t * (*n % *q);
　　　　　for _ in 0..(*n / *q) { tt *= *q; if tt > limit { continue 'again; } }
　　　　}; break 'again;
　　　}; (ci, cn, ct) = (i, *n, *t); continue 'advance;
　　}; break 'advance;
　}; answer.sort(); answer
}

http://mevius.5ch.net/test/read.cgi/tech/1691038333/750

上下前次 1-新書関写板覧索設栞歴

ｽﾚ情報赤ﾚｽ抽出画像ﾚｽ抽出歴の未読ｽﾚ AAｻﾑﾈｲﾙ

ぬこの手ぬこTOP 0.026s