結果
問題 | No.2443 特殊線形群の標準表現 |
ユーザー | 👑 p-adic |
提出日時 | 2023-08-05 21:42:08 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 125 ms / 3,000 ms |
コード長 | 3,651 bytes |
コンパイル時間 | 3,333 ms |
コンパイル使用メモリ | 220,984 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-04-30 13:48:48 |
合計ジャッジ時間 | 6,848 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 3 ms
6,816 KB |
testcase_01 | AC | 4 ms
6,940 KB |
testcase_02 | AC | 4 ms
6,940 KB |
testcase_03 | AC | 4 ms
6,940 KB |
testcase_04 | AC | 4 ms
6,940 KB |
testcase_05 | AC | 5 ms
6,940 KB |
testcase_06 | AC | 4 ms
6,940 KB |
testcase_07 | AC | 4 ms
6,944 KB |
testcase_08 | AC | 5 ms
6,940 KB |
testcase_09 | AC | 4 ms
6,944 KB |
testcase_10 | AC | 4 ms
6,944 KB |
testcase_11 | AC | 4 ms
6,944 KB |
testcase_12 | AC | 5 ms
6,940 KB |
testcase_13 | AC | 7 ms
6,940 KB |
testcase_14 | AC | 20 ms
6,944 KB |
testcase_15 | AC | 125 ms
6,940 KB |
testcase_16 | AC | 122 ms
6,944 KB |
testcase_17 | AC | 121 ms
6,944 KB |
testcase_18 | AC | 122 ms
6,940 KB |
testcase_19 | AC | 122 ms
6,940 KB |
testcase_20 | AC | 113 ms
6,944 KB |
ソースコード
#ifdef DEBUG #define _GLIBCXX_DEBUG #define UNTIE ios_base::sync_with_stdio( false ); cin.tie( nullptr ); signal( SIGABRT , &AlertAbort ) #define DEXPR( LL , BOUND , VALUE , DEBUG_VALUE ) CEXPR( LL , BOUND , DEBUG_VALUE ) #define CERR( MESSAGE ) cerr << MESSAGE << endl; #define COUT( ANSWER ) cout << ANSWER << endl #define ASSERT( A , MIN , MAX ) CERR( "ASSERTチェック: " << ( MIN ) << ( ( MIN ) <= A ? "<=" : ">" ) << A << ( A <= ( MAX ) ? "<=" : ">" ) << ( MAX ) ); assert( ( MIN ) <= A && A <= ( MAX ) ) #else #pragma GCC optimize ( "O3" ) #pragma GCC optimize( "unroll-loops" ) #pragma GCC target ( "sse4.2,fma,avx2,popcnt,lzcnt,bmi2" ) #define UNTIE ios_base::sync_with_stdio( false ); cin.tie( nullptr ) #define DEXPR( LL , BOUND , VALUE , DEBUG_VALUE ) CEXPR( LL , BOUND , VALUE ) #define CERR( MESSAGE ) #define COUT( ANSWER ) cout << ANSWER << "\n" #define ASSERT( A , MIN , MAX ) assert( ( MIN ) <= A && A <= ( MAX ) ) #endif #include <bits/stdc++.h> using namespace std; using ll = long long; #define MAIN main #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 CIN_ASSERT( A , MIN , MAX ) TYPE_OF( MAX ) A; SET_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 FOREQ( VAR , INITIAL , FINAL ) for( TYPE_OF( FINAL ) VAR = INITIAL ; VAR <= FINAL ; VAR ++ ) #define REPEAT( HOW_MANY_TIMES ) FOR( VARIABLE_FOR_REPEAT_ ## HOW_MANY_TIMES , 0 , HOW_MANY_TIMES ) #define QUIT return 0 #define SET_ASSERT( A , MIN , MAX ) cin >> A; ASSERT( A , MIN , MAX ) #ifdef DEBUG inline void AlertAbort( int n ) { CERR( "abort関数が呼ばれました。assertマクロのメッセージが出力されていない場合はオーバーフローの有無を確認をしてください。" ); } #endif template <typename T> inline T Residue( const T& a , const T& p ){ return a >= 0 ? a % p : p - 1 - ( ( - ( a + 1 ) ) % p ); } #define OO first.first #define OI first.second #define IO second.first #define II second.second ll B; using Matrix = pair<pair<ll,ll>,pair<ll,ll> >; inline Matrix operator*( const Matrix& M , const Matrix& N ) { return { { ( M.OO * N.OO + M.OI * N.IO ) % B , ( M.OO * N.OI + M.OI * N.II ) % B } , { ( M.IO * N.OO + M.II * N.IO ) % B , ( M.IO * N.OI + M.II * N.II ) % B } }; } int MAIN() { UNTIE; DEXPR( int , bound_N , 100000 , 100 ); // 0が5個 CIN_ASSERT( N , 1 , bound_N ); CEXPR( ll , bound_ABx , 1000000000 ); // 0が9個 SET_ASSERT( B , 1 , bound_ABx ); DEXPR( int , bound_Q , 100000 , 100 ); // 0が5個 CIN_ASSERT( Q , 1 , bound_Q ); Matrix temp{ { 1 , 0 } , { 0 , 1 } }; Matrix prod[bound_N + 1] = { temp }; FOREQ( n , 1 , N ){ CIN_ASSERT( AOO , -bound_ABx , bound_ABx ); CIN_ASSERT( AOI , -bound_ABx , bound_ABx ); CIN_ASSERT( AIO , -bound_ABx , bound_ABx ); CIN_ASSERT( AII , -bound_ABx , bound_ABx ); assert( AOO * AII - AOI * AIO == 1 ); Matrix A = { { AOO , AOI } , { AIO , AII } }; prod[n] = temp = A * temp; } REPEAT( Q ){ CIN_ASSERT( Lq , 0 , N ); CIN_ASSERT( Rq , Lq , N ); CIN_ASSERT( x , -bound_ABx , bound_ABx ); CIN_ASSERT( y , -bound_ABx , bound_ABx ); Matrix& prod_Lq = prod[Lq]; Matrix& prod_Rq = prod[Rq]; temp = { { prod_Lq.II , -prod_Lq.OI } , { -prod_Lq.IO , prod_Lq.OO } }; temp = prod_Rq * temp; ll z = Residue( temp.OO * x + temp.OI * y , B ); ll w = Residue( temp.IO * x + temp.II * y , B ); COUT( z << " " << w ); } QUIT; }