結果

問題 No.2756 GCD Teleporter
コンテスト
ユーザー magurofly
提出日時 2024-05-10 23:06:10
言語 Rust
(1.83.0 + proconio)
結果
WA  
実行時間 -
コード長 4,072 bytes
コンパイル時間 13,900 ms
コンパイル使用メモリ 387,792 KB
実行使用メモリ 18,080 KB
最終ジャッジ日時 2024-12-20 07:16:48
合計ジャッジ時間 16,942 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 29 WA * 7
権限があれば一括ダウンロードができます
コンパイルメッセージ
warning: variable `N` should have a snake case name
 --> src/main.rs:9:9
  |
9 |         N: usize,
  |         ^ help: convert the identifier to snake case: `n`
  |
  = note: `#[warn(non_snake_case)]` on by default

warning: variable `A` should have a snake case name
  --> src/main.rs:10:9
   |
10 |         A: [usize; N],
   |         ^ help: convert the identifier to snake case: `a`

ソースコード

diff # 

pub use __cargo_equip::prelude::*;

use proconio::input;
use std::collections::*;
use acl_dsu::*;

fn main() {
    input! {
        N: usize,
        A: [usize; N],
    }
    
    let max = *A.iter().max().unwrap();
    let (_, sieve) = linear_sieve(max);
    let mut by_prime_factor = HashMap::new();
    for i in 0 .. N {
        let mut n = A[i];
        while n > 1 {
            let p = sieve[n];
            while n % p == 0 {
                n /= p;
            }
            by_prime_factor.entry(p).or_insert_with(Vec::new).push(i);
        }
    }
    
    let mut dsu = Dsu::new(N);
    for (_, vertices) in &by_prime_factor {
        for i in 1 .. vertices.len() {
            dsu.merge(vertices[0], vertices[i]);
        }
    }
    
    let mut min_prime_factor = vec![200000; N];
    let mut leaders = vec![];
    for v in 0 .. N {
        let r = dsu.leader(v);
        leaders.push(r);
        let p = &mut min_prime_factor[r];
        let p0 = *p;
        let q = sieve[A[v]];
        if q > 1 {
            *p = p0.min(q);
        }
    }
    leaders.sort();
    leaders.dedup();
    
    // eprintln!("{:?}", leaders);
    // eprintln!("{:?}", min_prime_factor);
    
    let mut ans = 0;
    for i in 1 .. leaders.len() {
        let p = min_prime_factor[leaders[0]];
        let q = min_prime_factor[leaders[i]];
        let pq = p.min(q);
        if 4 < pq {
            ans += 4;
            min_prime_factor[leaders[0]] = 2;
        } else {
            ans += pq;
            min_prime_factor[leaders[0]] = pq;
        }
    }
    
    println!("{}", ans);
}

fn linear_sieve(limit: usize) -> (Vec<usize>, Vec<usize>) {
    let mut primes = vec![];
    let mut table = vec![1; limit + 1];
    for d in 2 ..= limit {
        if table[d] == 1 {
            table[d] = d;
            primes.push(d);
        }
        for &p in &primes {
            if p * d > limit || p > table[d] { break; }
            table[p * d] = p;
        }
    }
    (primes, table)
}
// The following code was expanded by `cargo-equip`.

///  # Bundled libraries
/// 
///  - `ac-library-rs-parted-dsu 0.1.0 (path+██████████████████████████████████████████████████████████████████████)` published in https://github.com/qryxip/ac-library-rs-parted licensed under `CC0-1.0` as `crate::__cargo_equip::crates::acl_dsu`
#[cfg_attr(any(), rustfmt::skip)]
#[allow(unused)]
mod __cargo_equip {
    pub(crate) mod crates {
        pub mod acl_dsu {pub use self::dsu::*;mod dsu{pub struct Dsu{n:usize,parent_or_size:Vec<i32>,}impl Dsu{pub fn new(size:usize)->Self{Self{n:size,parent_or_size:vec![-1;size],}}pub fn merge(&mut self,a:usize,b:usize)->usize{assert!(a<self.n);assert!(b<self.n);let(mut x,mut y)=(self.leader(a),self.leader(b));if x==y{return x;}if-self.parent_or_size[x]< -self.parent_or_size[y]{std::mem::swap(&mut x,&mut y);}self.parent_or_size[x]+=self.parent_or_size[y];self.parent_or_size[y]=x as i32;x}pub fn same(&mut self,a:usize,b:usize)->bool{assert!(a<self.n);assert!(b<self.n);self.leader(a)==self.leader(b)}pub fn leader(&mut self,a:usize)->usize{assert!(a<self.n);if self.parent_or_size[a]<0{return a;}self.parent_or_size[a]=self.leader(self.parent_or_size[a]as usize)as i32;self.parent_or_size[a]as usize}pub fn size(&mut self,a:usize)->usize{assert!(a<self.n);let x=self.leader(a);-self.parent_or_size[x]as usize}pub fn groups(&mut self)->Vec<Vec<usize>>{let mut leader_buf=vec![0;self.n];let mut group_size=vec![0;self.n];for i in 0..self.n{leader_buf[i]=self.leader(i);group_size[leader_buf[i]]+=1;}let mut result=vec![Vec::new();self.n];for i in 0..self.n{result[i].reserve(group_size[i]);}for i in 0..self.n{result[leader_buf[i]].push(i);}result.into_iter().filter(|x|!x.is_empty()).collect::<Vec<Vec<usize>>>()}}}}
    }

    pub(crate) mod macros {
        pub mod acl_dsu {}
    }

    pub(crate) mod prelude {pub use crate::__cargo_equip::crates::*;}

    mod preludes {
        pub mod acl_dsu {}
    }
}
0