結果

問題 No.3127 Multiple of Twin Prime
コンテスト
ユーザー srtry
提出日時 2025-05-04 16:17:22
言語 Rust
(1.83.0 + proconio)
結果
AC  
実行時間 113 ms / 2,500 ms
コード長 2,159 bytes
コンパイル時間 12,030 ms
コンパイル使用メモリ 401,312 KB
実行使用メモリ 14,488 KB
最終ジャッジ日時 2025-05-04 16:17:38
合計ジャッジ時間 14,957 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 12
権限があれば一括ダウンロードができます
コンパイルメッセージ
warning: unused import: `proconio::input`
 --> src/main.rs:6:5
  |
6 | use proconio::input;
  |     ^^^^^^^^^^^^^^^
  |
  = note: `#[warn(unused_imports)]` on by default

warning: value assigned to `mid` is never read
  --> src/main.rs:12:13
   |
12 |     let mut mid:usize = (left+right+1)/2;
   |             ^^^
   |
   = help: maybe it is overwritten before being read?
   = note: `#[warn(unused_assignments)]` on by default

ソースコード

diff # 

/* 
 *      Author:    srtry
 *      Created:   2025-04-30T04:38:41+09:00
 */

use proconio::input;
use std::io::{stdout,Write,BufWriter};

fn solve(&n:&usize, twin_primes_prod:&Vec<usize>) -> Option<usize> {
    let mut left:usize = 0;
    let mut right:usize = twin_primes_prod.len()-1;
    let mut mid:usize = (left+right+1)/2;
    loop {
        if right-left <= 1 {
            break;
        }
        mid = (left+right+1)/2;
        if twin_primes_prod[mid] <= n {
            left = mid;
        } else {
            right = mid-1;
        }
    }

    if twin_primes_prod[right] <= n {
        return Some(twin_primes_prod[right])
    } else if twin_primes_prod[left] <= n {
        return Some(twin_primes_prod[left]);
    } else {
        return None;
    }
}

fn main() {
    let stdin = std::io::read_to_string(std::io::stdin()).unwrap();
    let mut stdin = stdin.split_whitespace();
    
    let mut max_max_p:usize = 0;
    let t:usize = stdin.next().unwrap().parse::<usize>().unwrap();
    let case = {
        let mut tmp:Vec<usize> = vec![0;t];
        for i in 0..t {
            tmp[i] = stdin.next().unwrap().parse::<usize>().unwrap();
            if tmp[i] > max_max_p {
                max_max_p = tmp[i];
            }
        }
        tmp
    };
    max_max_p = (max_max_p as f64).sqrt() as usize; 

    // エラトステネスの篩
    let mut is_prime:Vec<bool> = vec![true;max_max_p+3];
    is_prime[0] = false;
    is_prime[1] = false;
    for i in 2..is_prime.len() {
        if is_prime[i] {
            for j in i..=(max_max_p/i) {
                is_prime[j*i] = false;
            }
        }
    }

    // 双子素数の積
    let mut twin_primes_prod:Vec<usize> = Vec::new();
    for i in 0..=max_max_p {
        if is_prime[i] && is_prime[i+2] {
            twin_primes_prod.push(i*(i+2));
        }
    }

    
    let out = stdout();
    let mut out = BufWriter::new(out.lock());

    let mut ans;
    for &c in case.iter() {
        ans = solve(&c, &twin_primes_prod);
        match ans {
            Some(x) => write!(out, "{}\n", x).unwrap(),
            None => write!(out, "-1\n").unwrap(),
        }
    }
}
0