結果

問題 No.1262 グラフを作ろう！
ユーザー akakimidoriakakimidori
提出日時 2020-10-16 22:55:00
言語 Rust
(1.77.0)
結果
AC  
実行時間 1,371 ms / 3,000 ms
コード長 3,378 bytes
コンパイル時間 1,989 ms
コンパイル使用メモリ 155,332 KB
実行使用メモリ 33,724 KB
最終ジャッジ日時 2023-09-28 04:48:52
合計ジャッジ時間 48,564 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 591 ms
24,320 KB
testcase_01 AC 413 ms
16,756 KB
testcase_02 AC 692 ms
27,140 KB
testcase_03 AC 290 ms
21,064 KB
testcase_04 AC 91 ms
15,628 KB
testcase_05 AC 430 ms
18,284 KB
testcase_06 AC 59 ms
17,440 KB
testcase_07 AC 89 ms
6,396 KB
testcase_08 AC 593 ms
23,408 KB
testcase_09 AC 172 ms
13,456 KB
testcase_10 AC 320 ms
14,748 KB
testcase_11 AC 624 ms
24,864 KB
testcase_12 AC 486 ms
23,040 KB
testcase_13 AC 339 ms
21,588 KB
testcase_14 AC 333 ms
15,184 KB
testcase_15 AC 110 ms
11,724 KB
testcase_16 AC 610 ms
22,756 KB
testcase_17 AC 533 ms
20,040 KB
testcase_18 AC 773 ms
28,844 KB
testcase_19 AC 134 ms
14,828 KB
testcase_20 AC 372 ms
21,692 KB
testcase_21 AC 731 ms
27,756 KB
testcase_22 AC 51 ms
12,568 KB
testcase_23 AC 254 ms
14,072 KB
testcase_24 AC 618 ms
23,412 KB
testcase_25 AC 93 ms
9,448 KB
testcase_26 AC 187 ms
12,776 KB
testcase_27 AC 408 ms
21,096 KB
testcase_28 AC 350 ms
16,840 KB
testcase_29 AC 247 ms
17,348 KB
testcase_30 AC 290 ms
13,396 KB
testcase_31 AC 570 ms
25,072 KB
testcase_32 AC 443 ms
17,872 KB
testcase_33 AC 107 ms
7,872 KB
testcase_34 AC 547 ms
24,244 KB
testcase_35 AC 398 ms
21,520 KB
testcase_36 AC 83 ms
7,920 KB
testcase_37 AC 817 ms
29,476 KB
testcase_38 AC 500 ms
20,288 KB
testcase_39 AC 215 ms
13,240 KB
testcase_40 AC 322 ms
14,556 KB
testcase_41 AC 300 ms
18,744 KB
testcase_42 AC 719 ms
28,736 KB
testcase_43 AC 404 ms
18,268 KB
testcase_44 AC 509 ms
17,612 KB
testcase_45 AC 701 ms
29,452 KB
testcase_46 AC 223 ms
20,644 KB
testcase_47 AC 118 ms
17,492 KB
testcase_48 AC 630 ms
26,984 KB
testcase_49 AC 754 ms
28,700 KB
testcase_50 AC 684 ms
27,808 KB
testcase_51 AC 152 ms
18,372 KB
testcase_52 AC 278 ms
19,244 KB
testcase_53 AC 133 ms
7,748 KB
testcase_54 AC 390 ms
22,900 KB
testcase_55 AC 405 ms
21,460 KB
testcase_56 AC 289 ms
17,444 KB
testcase_57 AC 874 ms
31,052 KB
testcase_58 AC 319 ms
21,216 KB
testcase_59 AC 289 ms
18,152 KB
testcase_60 AC 813 ms
29,752 KB
testcase_61 AC 282 ms
13,948 KB
testcase_62 AC 643 ms
26,100 KB
testcase_63 AC 321 ms
22,148 KB
testcase_64 AC 397 ms
16,616 KB
testcase_65 AC 156 ms
10,728 KB
testcase_66 AC 301 ms
14,584 KB
testcase_67 AC 269 ms
13,676 KB
testcase_68 AC 190 ms
18,596 KB
testcase_69 AC 304 ms
13,248 KB
testcase_70 AC 528 ms
25,736 KB
testcase_71 AC 267 ms
17,412 KB
testcase_72 AC 547 ms
26,340 KB
testcase_73 AC 327 ms
24,944 KB
testcase_74 AC 60 ms
21,412 KB
testcase_75 AC 294 ms
22,648 KB
testcase_76 AC 237 ms
16,916 KB
testcase_77 AC 172 ms
22,148 KB
testcase_78 AC 348 ms
22,972 KB
testcase_79 AC 589 ms
27,516 KB
testcase_80 AC 527 ms
21,740 KB
testcase_81 AC 135 ms
21,720 KB
testcase_82 AC 208 ms
23,792 KB
testcase_83 AC 487 ms
20,804 KB
testcase_84 AC 251 ms
17,252 KB
testcase_85 AC 612 ms
27,988 KB
testcase_86 AC 281 ms
22,972 KB
testcase_87 AC 470 ms
26,084 KB
testcase_88 AC 129 ms
21,588 KB
testcase_89 AC 432 ms
19,684 KB
testcase_90 AC 45 ms
25,688 KB
testcase_91 AC 44 ms
25,780 KB
testcase_92 AC 44 ms
25,692 KB
testcase_93 AC 47 ms
25,784 KB
testcase_94 AC 46 ms
25,788 KB
testcase_95 AC 1,371 ms
33,724 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// ---------- begin input macro ----------
// reference: https://qiita.com/tanakh/items/0ba42c7ca36cd29d0ac8
macro_rules! input {
    (source = $s:expr, $($r:tt)*) => {
        let mut iter = $s.split_whitespace();
        input_inner!{iter, $($r)*}
    };
    ($($r:tt)*) => {
        let s = {
            use std::io::Read;
            let mut s = String::new();
            std::io::stdin().read_to_string(&mut s).unwrap();
            s
        };
        let mut iter = s.split_whitespace();
        input_inner!{iter, $($r)*}
    };
}

macro_rules! input_inner {
    ($iter:expr) => {};
    ($iter:expr, ) => {};
    ($iter:expr, $var:ident : $t:tt $($r:tt)*) => {
        let $var = read_value!($iter, $t);
        input_inner!{$iter $($r)*}
    };
}

macro_rules! read_value {
    ($iter:expr, ( $($t:tt),* )) => {
        ( $(read_value!($iter, $t)),* )
    };
    ($iter:expr, [ $t:tt ; $len:expr ]) => {
        (0..$len).map(|_| read_value!($iter, $t)).collect::<Vec<_>>()
    };
    ($iter:expr, chars) => {
        read_value!($iter, String).chars().collect::<Vec<char>>()
    };
    ($iter:expr, bytes) => {
        read_value!($iter, String).bytes().collect::<Vec<u8>>()
    };
    ($iter:expr, usize1) => {
        read_value!($iter, usize) - 1
    };
    ($iter:expr, $t:ty) => {
        $iter.next().unwrap().parse::<$t>().expect("Parse error")
    };
}
// ---------- end input macro ----------

fn run() {
    input! {
        n: usize,
        m: usize,
        a: [usize; m],
    }
    let w = n;
    let mut f = (0..=w).collect::<Vec<_>>();
    for i in 2..=w {
        if i * i > w {
            break;
        }
        if i == f[i] {
            for j in i..=(w / i) {
                f[j * i] = i;
            }
        }
    }
    let mut memo = vec![0; n + 1];
    let mut ans = 0;
    for a in a {
        if memo[a] == 0 {
            let mut d = vec![1];
            let mut b = a;
            while b > 1 {
                let p = f[b];
                let mut c = 0;
                while b % p == 0 {
                    b /= p;
                    c += 1;
                }
                let len = d.len();
                let mut mul = 1;
                for _ in 0..c {
                    mul *= p;
                    for i in 0..len {
                        let v = d[i] * mul;
                        d.push(v);
                    }
                }
            }
            d.sort();
            let mut sum = vec![0; d.len()];
            for (s, d) in sum.iter_mut().zip(d.iter()) {
                *s = a / *d;
            }
            let mut b = a;
            while b > 1 {
                let p = f[b];
                while b % p == 0 {
                    b /= p;
                }
                let mut j = 0;
                for i in 0..d.len() {
                    if d[i] * p > a {
                        break;
                    }
                    if a / d[i] % p != 0 {
                        continue;
                    }
                    while d[i] * p != d[j] {
                        j += 1;
                    }
                    sum[i] -= sum[j];
                }
            }
            for (d, s) in d.iter().zip(sum.iter()) {
                memo[a] += ((a - 1) % *d) * *s;
            }
        }
        ans += memo[a];
    }
    println!("{}", ans);
}

fn main() {
    run();
}
0