結果

問題 No.89 どんどんドーナツどーんといこう！
ユーザー 👑 p-adicp-adic
提出日時 2022-08-12 12:28:25
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 1,933 bytes
コンパイル時間 857 ms
コンパイル使用メモリ 77,724 KB
実行使用メモリ 4,372 KB
最終ジャッジ日時 2023-10-23 22:02:01
合計ジャッジ時間 1,515 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,372 KB
testcase_01 AC 2 ms
4,372 KB
testcase_02 AC 2 ms
4,372 KB
testcase_03 AC 2 ms
4,372 KB
testcase_04 AC 2 ms
4,372 KB
testcase_05 AC 2 ms
4,372 KB
testcase_06 AC 2 ms
4,372 KB
testcase_07 AC 2 ms
4,372 KB
testcase_08 AC 2 ms
4,372 KB
testcase_09 AC 2 ms
4,372 KB
testcase_10 AC 2 ms
4,372 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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 0 

template <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 / 4
  DOUBLE( 6 , ( C * ( Rin + Rout ) * ( Rout - Rin ) * ( Rout - Rin ) ) * ( 3.1415926535898 * 3.1415926535898 ) / 4.0 );

}
0