結果
問題 |
No.1319 最強とんがりコーン
|
ユーザー |
![]() |
提出日時 | 2020-12-17 02:10:47 |
言語 | Rust (1.83.0 + proconio) |
結果 |
AC
|
実行時間 | 4 ms / 2,000 ms |
コード長 | 1,093 bytes |
コンパイル時間 | 13,397 ms |
コンパイル使用メモリ | 377,732 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-09-20 06:01:43 |
合計ジャッジ時間 | 13,431 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 70 |
ソースコード
fn calc(d: f64) -> f64 { let func = |r: f64| -> f64 { if 2.0 * r < d { return 0f64; } let p = (d / (2.0 * r)).acos() * 2.0; let pp = p * p; let mut v = 0f64; let m = 17; let mut de = -1.0 / (1..=m).fold(1f64, |s, a| s * a as f64); for i in 0..=7 { v = v * pp + de; de = -de * ((m - 2 * i) * (m - 2 * i - 1)) as f64; } v * r * r * pp * p }; let n = 1 << 17; let f = 1f64 / n as f64; let mut sum = 0f64; for i in (0..n).rev() { let r = (n - i) as f64 * f; let l = (n - 1 - i) as f64 * f; sum += func((l + r) * 0.5); } sum * f } fn run() { let mut s = String::new(); std::io::stdin().read_line(&mut s).unwrap(); let mut it = s.trim().split_whitespace(); let r: f64 = it.next().unwrap().parse().unwrap(); let h: f64 = it.next().unwrap().parse().unwrap(); let d: f64 = it.next().unwrap().parse().unwrap(); let ans = r * r * h * calc(d / r); println!("{:.7}", ans); } fn main() { run(); }