結果

問題 No.2973 シュニレルマン積分入門
コンテスト
ユーザー 👑 p-adic
提出日時 2023-08-20 18:29:11
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,022 bytes
コンパイル時間 2,811 ms
コンパイル使用メモリ 217,764 KB
最終ジャッジ日時 2025-02-16 11:52:34
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other WA * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

// 誤解法（単純にcomplex<double>を用いた誤差の大きい近似解）チェック
#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( PRECISION ) cout << fixed << setprecision( PRECISION )
#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;
  }
  ll m = 1;
  ll m_factorial = 1;
  while( m < 11 ){
    m++;
    m_factorial *= m;
  }
  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