結果

問題 No.1319 最強とんがりコーン
ユーザー akakimidori
提出日時 2020-12-17 01:42:09
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 818 ms / 2,000 ms
コード長 849 bytes
コンパイル時間 400 ms
コンパイル使用メモリ 31,360 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-20 05:54:05
合計ジャッジ時間 28,672 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 70
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<stdio.h>
#include<stdlib.h>
#include<stdint.h>
#include<inttypes.h>
#include<math.h>

typedef int32_t i32;
typedef int64_t i64;

double PARAM = 0.0;

double func (double r) {
    if (2 * r < PARAM) {
        return 0;
    }
    double p = acos (PARAM / (2 * r)) * 2;
    return r * r * (p - sin(p));
}

double calc (void) {
    const i32 n = 1 << 23;
    double f = (double) 1 / n;
    double sum = 0;
    double i6 = (double)1 / 6;
    for (i32 i = n - 1; i >= 0; --i) {
        double r = (n - i) * f;
        double l = (n - 1 - i) * f;
        sum += (func(l) + 4 * func((l + r) * 0.5) + func(r)) * i6;
    }
    return sum * f;
}

void run (void) {
    double r, h, d;
    scanf ("%lf%lf%lf", &r, &h, &d);
    PARAM = d / r;
    double ans = r * r * h * calc();
    printf ("%7f", ans);
}

int main (void) {
    run();
    return 0;
}
0