結果

問題 No.764 浮動点
ユーザー koba-e964
提出日時 2023-07-26 00:37:32
言語 Rust
(1.83.0 + proconio)
結果
AC  
実行時間 2 ms / 1,500 ms
コード長 3,167 bytes
コンパイル時間 14,220 ms
コンパイル使用メモリ 377,748 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-02 17:12:28
合計ジャッジ時間 13,709 ms
ジャッジサーバーID
(参考情報)
judge3 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#[allow(unused_imports)]
use std::cmp::*;
#[allow(unused_imports)]
use std::collections::*;
use std::io::{Write, BufWriter};
// https://qiita.com/tanakh/items/0ba42c7ca36cd29d0ac8
macro_rules! input {
($($r:tt)*) => {
let stdin = std::io::stdin();
let mut bytes = std::io::Read::bytes(std::io::BufReader::new(stdin.lock()));
let mut next = move || -> String{
bytes.by_ref().map(|r|r.unwrap() as char)
.skip_while(|c|c.is_whitespace())
.take_while(|c|!c.is_whitespace())
.collect()
};
input_inner!{next, $($r)*}
};
}
macro_rules! input_inner {
($next:expr) => {};
($next:expr,) => {};
($next:expr, $var:ident : $t:tt $($r:tt)*) => {
let $var = read_value!($next, $t);
input_inner!{$next $($r)*}
};
}
macro_rules! read_value {
($next:expr, ( $($t:tt),* )) => { ($(read_value!($next, $t)),*) };
($next:expr, [ $t:tt ; $len:expr ]) => {
(0..$len).map(|_| read_value!($next, $t)).collect::<Vec<_>>()
};
($next:expr, chars) => {
read_value!($next, String).chars().collect::<Vec<char>>()
};
($next:expr, usize1) => (read_value!($next, usize) - 1);
($next:expr, [ $t:tt ]) => {{
let len = read_value!($next, usize);
read_value!($next, [$t; len])
}};
($next:expr, $t:ty) => ($next().parse::<$t>().expect("Parse error"));
}
// Find the area of x^2 + y^2 <= a^2, (x-l)^2 + y^2 <= b^2
fn f(a: i64, b: i64, l: i64) -> f64 {
if a + b <= l {
return 0.0;
}
let pi = std::f64::consts::PI;
if a + l <= b {
return pi * (a * a) as f64;
}
if b + l <= a {
return pi * (b * b) as f64;
}
let x0 = (l * l + a * a - b * b) as f64 / (2 * l) as f64;
let height = (a as f64 * a as f64 - x0 * x0).sqrt();
let mut ans = l as f64 * height;
ans = -ans;
let adeg = height.atan2(x0);
ans += adeg * a as f64 * a as f64;
let bdeg = height.atan2(l as f64 - x0);
ans += bdeg * b as f64 * b as f64;
ans
}
// https://yukicoder.me/problems/no/764 (3.5)
// 2
// -> WA deg * r^2/2 pi * deg * r^2/2
fn main() {
let out = std::io::stdout();
let mut out = BufWriter::new(out.lock());
macro_rules! puts {($($format:tt)*) => (let _ = write!(out,$($format)*););}
input! {
n: usize,
suml: i64,
l: [i64; n + 1],
}
let mut lft = vec![(0, 0); n];
let mut ma = 0;
let mut sum = 0;
for i in 1..n + 1 {
sum += l[i - 1];
ma = max(ma, l[i - 1]);
lft[i - 1] = (max(0, 2 * ma - sum), sum);
}
let mut rgt = vec![(0, 0); n];
let mut ma = 0;
let mut sum = 0;
for i in (1..n + 1).rev() {
sum += l[i];
ma = max(ma, l[i]);
rgt[i - 1] = (max(0, 2 * ma - sum), sum);
}
for i in 0..n {
let (a, b) = lft[i];
let (c, d) = rgt[i];
puts!("{}\n", f(b, d, suml) - f(a, d, suml) - f(b, c, suml) + f(a, c, suml));
}
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0