Submission #2163433
Source Code Expand
/**
* _ _ __ _ _ _ _ _ _ _
* | | | | / / | | (_) | (_) | | (_) | |
* | |__ __ _| |_ ___ ___ / /__ ___ _ __ ___ _ __ ___| |_ _| |_ ___ _____ ______ _ __ _ _ ___| |_ ______ ___ _ __ _ _ __ _ __ ___| |_ ___
* | '_ \ / _` | __/ _ \ / _ \ / / __/ _ \| '_ ` _ \| '_ \ / _ \ __| | __| \ \ / / _ \______| '__| | | / __| __|______/ __| '_ \| | '_ \| '_ \ / _ \ __/ __|
* | | | | (_| | || (_) | (_) / / (_| (_) | | | | | | |_) | __/ |_| | |_| |\ V / __/ | | | |_| \__ \ |_ \__ \ | | | | |_) | |_) | __/ |_\__ \
* |_| |_|\__,_|\__\___/ \___/_/ \___\___/|_| |_| |_| .__/ \___|\__|_|\__|_| \_/ \___| |_| \__,_|___/\__| |___/_| |_|_| .__/| .__/ \___|\__|___/
* | | | | | |
* |_| |_| |_|
*
* https://github.com/hatoo/competitive-rust-snippets
*/
#[allow(unused_imports)]
use std::cmp::{max, min, Ordering};
#[allow(unused_imports)]
use std::collections::{BTreeMap, BTreeSet, BinaryHeap, HashMap, HashSet, VecDeque};
#[allow(unused_imports)]
use std::iter::FromIterator;
#[allow(unused_imports)]
use std::io::{stdin, stdout, BufWriter, Write};
mod util {
use std::io::{stdin, stdout, BufWriter, StdoutLock};
use std::str::FromStr;
use std::fmt::Debug;
#[allow(dead_code)]
pub fn line() -> String {
let mut line: String = String::new();
stdin().read_line(&mut line).unwrap();
line.trim().to_string()
}
#[allow(dead_code)]
pub fn chars() -> Vec<char> {
line().chars().collect()
}
#[allow(dead_code)]
pub fn gets<T: FromStr>() -> Vec<T>
where
<T as FromStr>::Err: Debug,
{
let mut line: String = String::new();
stdin().read_line(&mut line).unwrap();
line.split_whitespace()
.map(|t| t.parse().unwrap())
.collect()
}
#[allow(dead_code)]
pub fn with_bufwriter<F: FnOnce(BufWriter<StdoutLock>) -> ()>(f: F) {
let out = stdout();
let writer = BufWriter::new(out.lock());
f(writer)
}
}
#[allow(unused_macros)]
macro_rules ! get { ( $ t : ty ) => { { let mut line : String = String :: new ( ) ; stdin ( ) . read_line ( & mut line ) . unwrap ( ) ; line . trim ( ) . parse ::<$ t > ( ) . unwrap ( ) } } ; ( $ ( $ t : ty ) ,* ) => { { let mut line : String = String :: new ( ) ; stdin ( ) . read_line ( & mut line ) . unwrap ( ) ; let mut iter = line . split_whitespace ( ) ; ( $ ( iter . next ( ) . unwrap ( ) . parse ::<$ t > ( ) . unwrap ( ) , ) * ) } } ; ( $ t : ty ; $ n : expr ) => { ( 0 ..$ n ) . map ( | _ | get ! ( $ t ) ) . collect ::< Vec < _ >> ( ) } ; ( $ ( $ t : ty ) ,*; $ n : expr ) => { ( 0 ..$ n ) . map ( | _ | get ! ( $ ( $ t ) ,* ) ) . collect ::< Vec < _ >> ( ) } ; ( $ t : ty ;; ) => { { let mut line : String = String :: new ( ) ; stdin ( ) . read_line ( & mut line ) . unwrap ( ) ; line . split_whitespace ( ) . map ( | t | t . parse ::<$ t > ( ) . unwrap ( ) ) . collect ::< Vec < _ >> ( ) } } ; ( $ t : ty ;; $ n : expr ) => { ( 0 ..$ n ) . map ( | _ | get ! ( $ t ;; ) ) . collect ::< Vec < _ >> ( ) } ; }
#[allow(unused_macros)]
macro_rules ! debug { ( $ ( $ a : expr ) ,* ) => { println ! ( concat ! ( $ ( stringify ! ( $ a ) , " = {:?}, " ) ,* ) , $ ( $ a ) ,* ) ; } }
#[allow(dead_code)]
pub const M: u64 = 1_000_000_007;
#[allow(dead_code)]
pub fn gcd(a: u64, b: u64) -> u64 {
if b == 0 {
a
} else {
gcd(b, a % b)
}
}
#[allow(dead_code)]
pub fn lcm(a: u64, b: u64) -> u64 {
a / gcd(a, b) * b
}
#[allow(dead_code)]
/// (gcd, x, y)
pub fn extgcd(a: i64, b: i64) -> (i64, i64, i64) {
if b == 0 {
(a, 1, 0)
} else {
let (gcd, x, y) = extgcd(b, a % b);
(gcd, y, x - (a / b) * y)
}
}
#[allow(dead_code)]
pub fn mod_pow(x: u64, n: u64, m: u64) -> u64 {
let mut res = 1;
let mut x = x;
let mut n = n;
while n > 0 {
if n & 1 == 1 {
res = (res * x) % m;
}
x = (x * x) % m;
n >>= 1;
}
res
}
#[allow(dead_code)]
pub fn mod_inverse(a: u64, m: u64) -> u64 {
let (_, x, _) = extgcd(a as i64, m as i64);
((m as i64 + x) as u64 % m) % m
}
#[allow(dead_code)]
pub fn fact_table(len: usize, m: u64) -> Vec<u64> {
let mut res = vec![1; len + 1];
for i in 1..len + 1 {
res[i] = (i as u64 * res[i - 1]) % m;
}
res
}
#[allow(dead_code)]
/// Factorial and Inverse factorial table
pub fn fact_inv_table(size: usize, m: u64) -> (Vec<u64>, Vec<u64>) {
let mut fact = vec![1; size];
let mut fact_inv = vec![1; size];
for i in 2..size {
fact[i] = fact[i - 1] * i as u64 % m;
fact_inv[i] = m - ((m / i as u64) * fact_inv[(m % i as u64) as usize] % m);
}
for i in 1..size {
fact_inv[i] = fact_inv[i - 1] * fact_inv[i] % m;
}
(fact, fact_inv)
}
#[allow(dead_code)]
/// (a mod p, e when n! = a p\^e)
pub fn mod_fact(n: u64, p: u64, fact: &[u64]) -> (u64, u64) {
if n == 0 {
(1, 0)
} else {
let (a, b) = mod_fact(n / p, p, fact);
let pow = b + n / p;
if n / p % 2 != 0 {
(a * (p - fact[(n % p) as usize]) % p, pow)
} else {
(a * fact[(n % p) as usize] % p, pow)
}
}
}
#[allow(dead_code)]
/// C(n, k) % p
pub fn mod_comb(n: u64, k: u64, p: u64, fact: &[u64]) -> u64 {
if n < k {
0
} else {
let (a1, e1) = mod_fact(n, p, fact);
let (a2, e2) = mod_fact(k, p, fact);
let (a3, e3) = mod_fact(n - k, p, fact);
if e1 > e2 + e3 {
0
} else {
a1 * mod_inverse(a2 * a3 % p, p) % p
}
}
}
#[allow(dead_code)]
/// H(n, k) % p
pub fn mod_comb_repetition(n: u64, k: u64, p: u64, fact: &[u64]) -> u64 {
mod_comb(n - 1 + k, n - 1, p, fact)
}
#[derive(Eq, PartialEq, Clone, Debug)]
/// Equivalent to std::cmp::Reverse
pub struct Rev<T>(pub T);
impl<T: PartialOrd> PartialOrd for Rev<T> {
fn partial_cmp(&self, other: &Rev<T>) -> Option<Ordering> {
other.0.partial_cmp(&self.0)
}
}
impl<T: Ord> Ord for Rev<T> {
fn cmp(&self, other: &Rev<T>) -> Ordering {
other.0.cmp(&self.0)
}
}
#[allow(dead_code)]
fn main() {
let n = get!(usize);
let mut wh = get!(u64, u64; n);
wh.sort_by_key(|t| Rev(t.1));
let fact = fact_table(200000, M);
let mut ans = 0;
let inv = |x| mod_inverse(x, M);
let comb = |n, r| mod_comb(n, r, M, &fact);
let mut p = 0;
let mut q = 0;
for ((w, h), i) in wh.into_iter().zip(1..) {
if i >= 3 {
let l = (p + M - h * q % M) % M;
let r = 2 * fact[i] % M * comb(n as u64, i as u64) % M * fact[n - i] % M * w % M;
ans += l * r % M;
ans %= M;
}
if i > 1 {
let yy = inv(i as u64 * (i as u64 + 1));
p += h * yy;
p %= M;
q += yy;
q %= M;
}
}
println!("{}", ans);
}
Submission Info
Submission Time |
|
Task |
E - 天下一コップ |
User |
hatoo |
Language |
Rust (1.15.1) |
Score |
140 |
Code Size |
7588 Byte |
Status |
AC |
Exec Time |
103 ms |
Memory |
8444 KB |
Judge Result
Set Name |
Sample |
Subtask1 |
Subtask2 |
Subtask3 |
Score / Max Score |
0 / 0 |
15 / 15 |
55 / 55 |
70 / 70 |
Status |
|
|
|
|
Set Name |
Test Cases |
Sample |
sample_01.txt, sample_02.txt, sample_03.txt |
Subtask1 |
sample_01.txt, sample_02.txt, sample_03.txt, subtask1_01.txt, subtask1_02.txt, subtask1_03.txt, subtask1_04.txt, subtask1_05.txt, subtask1_06.txt, subtask1_07.txt, subtask1_08.txt, subtask1_09.txt, subtask1_10.txt, subtask1_11.txt, subtask1_12.txt, subtask1_13.txt, subtask1_14.txt |
Subtask2 |
sample_01.txt, sample_02.txt, sample_03.txt, subtask1_01.txt, subtask1_02.txt, subtask1_03.txt, subtask1_04.txt, subtask1_05.txt, subtask1_06.txt, subtask1_07.txt, subtask1_08.txt, subtask1_09.txt, subtask1_10.txt, subtask1_11.txt, subtask1_12.txt, subtask1_13.txt, subtask1_14.txt, subtask2_01.txt, subtask2_02.txt, subtask2_03.txt, subtask2_04.txt, subtask2_05.txt, subtask2_06.txt, subtask2_07.txt, subtask2_08.txt, subtask2_09.txt, subtask2_10.txt, subtask2_11.txt, subtask2_12.txt, subtask2_13.txt, subtask2_14.txt |
Subtask3 |
sample_01.txt, sample_02.txt, sample_03.txt, subtask1_01.txt, subtask1_02.txt, subtask1_03.txt, subtask1_04.txt, subtask1_05.txt, subtask1_06.txt, subtask1_07.txt, subtask1_08.txt, subtask1_09.txt, subtask1_10.txt, subtask1_11.txt, subtask1_12.txt, subtask1_13.txt, subtask1_14.txt, subtask2_01.txt, subtask2_02.txt, subtask2_03.txt, subtask2_04.txt, subtask2_05.txt, subtask2_06.txt, subtask2_07.txt, subtask2_08.txt, subtask2_09.txt, subtask2_10.txt, subtask2_11.txt, subtask2_12.txt, subtask2_13.txt, subtask2_14.txt, subtask3_01.txt, subtask3_02.txt, subtask3_03.txt, subtask3_04.txt, subtask3_05.txt, subtask3_06.txt, subtask3_07.txt, subtask3_08.txt, subtask3_09.txt, subtask3_10.txt, subtask3_11.txt, subtask3_12.txt, subtask3_13.txt, subtask3_14.txt |
Case Name |
Status |
Exec Time |
Memory |
sample_01.txt |
AC |
3 ms |
6396 KB |
sample_02.txt |
AC |
3 ms |
6396 KB |
sample_03.txt |
AC |
3 ms |
6396 KB |
subtask1_01.txt |
AC |
3 ms |
6396 KB |
subtask1_02.txt |
AC |
3 ms |
6396 KB |
subtask1_03.txt |
AC |
3 ms |
6396 KB |
subtask1_04.txt |
AC |
3 ms |
6396 KB |
subtask1_05.txt |
AC |
3 ms |
6396 KB |
subtask1_06.txt |
AC |
3 ms |
6396 KB |
subtask1_07.txt |
AC |
3 ms |
6396 KB |
subtask1_08.txt |
AC |
3 ms |
6396 KB |
subtask1_09.txt |
AC |
3 ms |
6396 KB |
subtask1_10.txt |
AC |
3 ms |
6396 KB |
subtask1_11.txt |
AC |
3 ms |
6396 KB |
subtask1_12.txt |
AC |
3 ms |
6396 KB |
subtask1_13.txt |
AC |
3 ms |
6396 KB |
subtask1_14.txt |
AC |
3 ms |
6396 KB |
subtask2_01.txt |
AC |
5 ms |
6396 KB |
subtask2_02.txt |
AC |
5 ms |
6396 KB |
subtask2_03.txt |
AC |
5 ms |
6396 KB |
subtask2_04.txt |
AC |
5 ms |
6396 KB |
subtask2_05.txt |
AC |
5 ms |
6396 KB |
subtask2_06.txt |
AC |
5 ms |
6396 KB |
subtask2_07.txt |
AC |
5 ms |
6396 KB |
subtask2_08.txt |
AC |
5 ms |
6396 KB |
subtask2_09.txt |
AC |
5 ms |
6396 KB |
subtask2_10.txt |
AC |
5 ms |
6396 KB |
subtask2_11.txt |
AC |
5 ms |
6396 KB |
subtask2_12.txt |
AC |
5 ms |
6396 KB |
subtask2_13.txt |
AC |
5 ms |
6396 KB |
subtask2_14.txt |
AC |
5 ms |
6396 KB |
subtask3_01.txt |
AC |
87 ms |
8444 KB |
subtask3_02.txt |
AC |
98 ms |
8444 KB |
subtask3_03.txt |
AC |
102 ms |
8444 KB |
subtask3_04.txt |
AC |
98 ms |
8444 KB |
subtask3_05.txt |
AC |
97 ms |
8444 KB |
subtask3_06.txt |
AC |
95 ms |
8444 KB |
subtask3_07.txt |
AC |
102 ms |
8444 KB |
subtask3_08.txt |
AC |
95 ms |
8444 KB |
subtask3_09.txt |
AC |
96 ms |
8444 KB |
subtask3_10.txt |
AC |
103 ms |
8444 KB |
subtask3_11.txt |
AC |
97 ms |
8444 KB |
subtask3_12.txt |
AC |
100 ms |
8444 KB |
subtask3_13.txt |
AC |
99 ms |
8444 KB |
subtask3_14.txt |
AC |
103 ms |
8444 KB |