結果

問題 No.2973 シュニレルマン積分入門
ユーザー 👑 p-adicp-adic
提出日時 2023-08-20 18:30:47
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,991 bytes
コンパイル時間 3,092 ms
コンパイル使用メモリ 217,664 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-11-29 20:50:41
合計ジャッジ時間 11,235 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 WA -
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

// 誤解法（m!の代わりに2^4*3^2*5*7*11*13*17 = 12252240を用いた近似解）チェック
#pragma GCC optimize ( "O3" )
#pragma GCC optimize( "unroll-loops" )
#pragma GCC target ( "sse4.2,fma,avx2,popcnt,lzcnt,bmi2" )
#include <bits/stdc++.h>
using namespace std;
using ll = long long;

#define UNTIE ios_base::sync_with_stdio( false ); cin.tie( nullptr )
#define TYPE_OF( VAR ) decay_t<decltype( VAR )>
#define CEXPR( LL , BOUND , VALUE ) constexpr LL BOUND = VALUE
#define CIN( LL , A ) LL A; cin >> A 
#define ASSERT( A , MIN , MAX ) assert( MIN <= A && A <= MAX ) 
#define CIN_ASSERT( A , MIN , MAX ) CIN( TYPE_OF( MAX ) , A ); ASSERT( A , MIN , MAX ) 
#define FOR( VAR , INITIAL , FINAL_PLUS_ONE ) for( TYPE_OF( FINAL_PLUS_ONE ) VAR = INITIAL ; VAR < FINAL_PLUS_ONE ; VAR ++ ) 
#define QUIT return 0 
#define COUT( ANSWER ) cout << ANSWER << "\n"
#define SET_PRECISION( DECIMAL_DIGITS ) cout << fixed << setprecision( DECIMAL_DIGITS )
#define RETURN( ANSWER ) COU( ( ANSWER ) ); QUIT 

int main()
{
  CEXPR( ll , bound , 20 );
  CIN_ASSERT( N , - bound , bound );
  bool negative;
  if( N < 0 ){
    negative = true;
    N *= -1;
  } else {
    negative = false;
  }
  CEXPR( ll , m_factorial , 12252240 );
  CEXPR( double , pi , 3.14159265358979323846264338 );
  double theta = 2 * pi / m_factorial;
  complex<double> answer{ 0.0 , 0.0 };
  complex<double> zeta{ cos( theta ) , sin( theta ) };
  complex<double> zeta_N{ 1.0 , 0.0 };
  complex<double> zeta_power = negative ? 1.0 / zeta : zeta;
  while( N > 0 ){
    ( N & 1 ) == 1 ? zeta_N *= zeta_power : zeta_N;
    zeta_power *= zeta_power;
    N >>= 1;
  }
  complex<double> zeta_N_power{ 1.0 , 0.0 };
  zeta_power = zeta_N_power;
  FOR( k , 0 , m_factorial ){
    answer += zeta_N_power / ( zeta_power * ( zeta_power + 1.0 ) + 10.0 );
    zeta_N_power *= zeta_N;
    zeta_power += zeta;
  }
  answer /= m_factorial;
  CEXPR( ll , d , 1000000000000000000 );
  ll n = answer.real() * d;
  COUT( n << "/" << d );
  QUIT;
}
0