結果
問題 | No.89 どんどんドーナツどーんといこう! |
ユーザー |
👑 |
提出日時 | 2022-08-12 12:28:25 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 5,000 ms |
コード長 | 1,933 bytes |
コンパイル時間 | 542 ms |
コンパイル使用メモリ | 76,644 KB |
最終ジャッジ日時 | 2025-01-30 20:24:48 |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 9 |
ソースコード
#include <iostream>#include <list>#include <vector>#include <string>#include <stdio.h>#include <stdint.h>#include <iomanip>using namespace std;using uint = unsigned int;using ll = long long;#define CIN( LL , A ) LL A; cin >> A#define GETSEGMENT( A , SEPARATOR ) string A; getline( cin , A , SEPARATOR )#define GETLINE( A ) string A; getline( cin , A )#define FOR_LL( VAR , INITIAL , FINAL_PLUS_ONE ) for( ll VAR = INITIAL ; VAR < FINAL_PLUS_ONE ; VAR ++ )#define FOR_ITR( ARRAY , ITR , END ) for( auto ITR = ARRAY .begin() , END = ARRAY .end() ; ITR != END ; ITR ++ )#define RETURN( LL , ANSWER ) const LL answer_for_contest = ( ANSWER ); cout << answer_for_contest << endl; return 0#define DOUBLE( DIGITS , ANSWER ) const double answer_for_contest = ( ANSWER ); cout << fixed << setprecision( DIGITS ) << answer_for_contest << endl;return 0template <typename T> inline T Distance( const T& a , const ll& b ){ return a < b ? b - a : a - b; }int main(){CIN( ll , C );CIN( ll , Rin );CIN( ll , Rout );// integral( real r , w , theta ; ( r - ( Rin + Rout ) / 2 ) ^ 2 + w ^ 2 <= ( ( Rout - Rin ) / 2 ) ^ 2 && 0 <= theta && theta < 2 * pi ; r++ , w++, theta++ ) r dr dw dtheta// = 2 * pi * integral( real r ; Rin <= r && r <= Rout ; r++ ) 2 * r * sqrt( ( ( Rout - Rin ) / 2 ) ^ 2 - ( r - ( Rin + Rout ) / 2 ) ^ 2 ) dr// = 2 * pi *// (// difference( - ( 2 / 3 ) * ( ( ( Rout - Rin ) / 2 ) ^ 2 - ( X - ( Rin + Rout ) / 2 ) ^ 2 ) ^ ( 3 / 2 ) , Rin , Rout ]// - integral( real r ; Rin <= r && r <= Rout ; r++ ) - ( Rin + Rout ) * sqrt( ( ( Rout - Rin ) / 2 ) ^ 2 - ( r - ( Rin + Rout ) / 2 ) ^ 2 ) dr// )// = 2 * pi * ( 0 + ( Rin + Rout ) * pi * ( ( Rout - Rin ) / 2 ) ^ 2 / 2 )// = ( Rin + Rout ) * ( Rout - Rin ) ^ 2 * pi ^ 2 / 4DOUBLE( 6 , ( C * ( Rin + Rout ) * ( Rout - Rin ) * ( Rout - Rin ) ) * ( 3.1415926535898 * 3.1415926535898 ) / 4.0 );}