#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( ANSWER ) cerr << ANSWER << 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 ) ) #define LIBRARY_SEARCH if( LibrarySearch() != 0 ){ QUIT; }; #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( ANSWER ) #define COUT( ANSWER ) cout << ANSWER << "\n" #define ASSERT( A , MIN , MAX ) assert( ( MIN ) <= A && A <= ( MAX ) ) #define LIBRARY_SEARCH #endif #include using namespace std; using uint = unsigned int; using ll = long long; using ull = unsigned long long; #define ATT __attribute__( ( target( "sse4.2,fma,avx2,popcnt,lzcnt,bmi2" ) ) ) #define TYPE_OF( VAR ) decay_t #define CEXPR( LL , BOUND , VALUE ) constexpr LL BOUND = VALUE #define CIN( LL , A ) LL A; cin >> A #define CIN_ASSERT( A , MIN , MAX ) CIN( TYPE_OF( MAX ) , A ); ASSERT( A , MIN , MAX ) #define SET_ASSERT( A , MIN , MAX ) cin >> A; ASSERT( A , MIN , MAX ) #define GETLINE( A ) string A; getline( cin , A ) #define GETLINE_SEPARATE( A , SEPARATOR ) string A; getline( cin , A , SEPARATOR ) #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 FOREQINV( VAR , INITIAL , FINAL ) for( TYPE_OF( INITIAL ) VAR = INITIAL ; VAR >= FINAL ; VAR -- ) #define AUTO_ITR( ARRAY ) auto itr_ ## ARRAY = ARRAY .begin() , end_ ## ARRAY = ARRAY .end() #define FOR_ITR( ARRAY ) for( AUTO_ITR( ARRAY ) , itr = itr_ ## ARRAY ; itr_ ## ARRAY != end_ ## ARRAY ; itr_ ## ARRAY ++ , itr++ ) #define REPEAT( HOW_MANY_TIMES ) FOR( VARIABLE_FOR_REPEAT_ ## HOW_MANY_TIMES , 0 , HOW_MANY_TIMES ) #define QUIT return 0 #define SET_PRECISION( DECIMAL_DIGITS ) cout << fixed << setprecision( DECIMAL_DIGITS_ ) #define RETURN( ANSWER ) COUT( ( ANSWER ) ); QUIT inline void AlertAbort( int n ) { cerr << "abort関数が呼ばれました。assertマクロのメッセージが出力されていない場合はオーバーフローの有無を確認をしてください。" << endl; } template inline T Absolute( const T& a ){ return a > 0 ? a : -a; } template inline T Residue( const T& a , const T& p ){ return a >= 0 ? a % p : ( a % p ) + p; } #define POWER( ANSWER , ARGUMENT , EXPONENT ) \ static_assert( ! is_same::value && ! is_same::value ); \ TYPE_OF( ARGUMENT ) ANSWER{ 1 }; \ { \ TYPE_OF( ARGUMENT ) ARGUMENT_FOR_SQUARE_FOR_POWER = ( ARGUMENT ); \ TYPE_OF( EXPONENT ) EXPONENT_FOR_SQUARE_FOR_POWER = ( EXPONENT ); \ while( EXPONENT_FOR_SQUARE_FOR_POWER != 0 ){ \ if( EXPONENT_FOR_SQUARE_FOR_POWER % 2 == 1 ){ \ ANSWER *= ARGUMENT_FOR_SQUARE_FOR_POWER; \ } \ ARGUMENT_FOR_SQUARE_FOR_POWER *= ARGUMENT_FOR_SQUARE_FOR_POWER; \ EXPONENT_FOR_SQUARE_FOR_POWER /= 2; \ } \ } \ #define POWER_MOD( ANSWER , ARGUMENT , EXPONENT , MODULO ) \ ll ANSWER{ 1 }; \ { \ ll ARGUMENT_FOR_SQUARE_FOR_POWER = ( MODULO + ( ( ARGUMENT ) % MODULO ) ) % MODULO; \ TYPE_OF( EXPONENT ) EXPONENT_FOR_SQUARE_FOR_POWER = ( EXPONENT ); \ while( EXPONENT_FOR_SQUARE_FOR_POWER != 0 ){ \ if( EXPONENT_FOR_SQUARE_FOR_POWER % 2 == 1 ){ \ ANSWER = ( ANSWER * ARGUMENT_FOR_SQUARE_FOR_POWER ) % MODULO; \ } \ ARGUMENT_FOR_SQUARE_FOR_POWER = ( ARGUMENT_FOR_SQUARE_FOR_POWER * ARGUMENT_FOR_SQUARE_FOR_POWER ) % MODULO; \ EXPONENT_FOR_SQUARE_FOR_POWER /= 2; \ } \ } \ #define FACTORIAL_MOD( ANSWER , ANSWER_INV , INVERSE , MAX_INDEX , CONSTEXPR_LENGTH , MODULO ) \ static ll ANSWER[CONSTEXPR_LENGTH]; \ static ll ANSWER_INV[CONSTEXPR_LENGTH]; \ static ll INVERSE[CONSTEXPR_LENGTH]; \ { \ ll VARIABLE_FOR_PRODUCT_FOR_FACTORIAL = 1; \ ANSWER[0] = VARIABLE_FOR_PRODUCT_FOR_FACTORIAL; \ FOREQ( i , 1 , MAX_INDEX ){ \ ANSWER[i] = ( VARIABLE_FOR_PRODUCT_FOR_FACTORIAL *= i ) %= MODULO; \ } \ ANSWER_INV[0] = ANSWER_INV[1] = INVERSE[1] = VARIABLE_FOR_PRODUCT_FOR_FACTORIAL = 1; \ FOREQ( i , 2 , MAX_INDEX ){ \ ANSWER_INV[i] = ( VARIABLE_FOR_PRODUCT_FOR_FACTORIAL *= INVERSE[i] = MODULO - ( ( ( MODULO / i ) * INVERSE[MODULO % i] ) % MODULO ) ) %= MODULO; \ } \ } \ // 二分探索テンプレート // EXPRESSIONがANSWERの広義単調関数の時、EXPRESSION >= TARGETの整数解を格納。 #define BS( ANSWER , MINIMUM , MAXIMUM , EXPRESSION , TARGET , INEQUALITY , UPDATE_U , UPDATE_L , UPDATE_ANSWER ) \ static_assert( ! is_same::value && ! is_same::value ); \ ll ANSWER = MINIMUM; \ if( MINIMUM <= MAXIMUM ){ \ ll VARIABLE_FOR_BINARY_SEARCH_L = MINIMUM; \ ll VARIABLE_FOR_BINARY_SEARCH_U = MAXIMUM; \ ANSWER = ( VARIABLE_FOR_BINARY_SEARCH_L + VARIABLE_FOR_BINARY_SEARCH_U ) / 2; \ ll VARIABLE_FOR_DIFFERENCE_FOR_BINARY_SEARCH; \ while( VARIABLE_FOR_BINARY_SEARCH_L != VARIABLE_FOR_BINARY_SEARCH_U ){ \ VARIABLE_FOR_DIFFERENCE_FOR_BINARY_SEARCH = ( EXPRESSION ) - ( TARGET ); \ CERR( "二分探索中: " << VARIABLE_FOR_BINARY_SEARCH_L << "<=" << ANSWER << "<=" << VARIABLE_FOR_BINARY_SEARCH_U << ":" << EXPRESSION << "-" << TARGET << "=" << VARIABLE_FOR_DIFFERENCE_FOR_BINARY_SEARCH ); \ if( VARIABLE_FOR_DIFFERENCE_FOR_BINARY_SEARCH INEQUALITY 0 ){ \ VARIABLE_FOR_BINARY_SEARCH_U = UPDATE_U; \ } else { \ VARIABLE_FOR_BINARY_SEARCH_L = UPDATE_L; \ } \ ANSWER = UPDATE_ANSWER; \ } \ CERR( "二分探索終了: " << VARIABLE_FOR_BINARY_SEARCH_L << "<=" << ANSWER << "<=" << VARIABLE_FOR_BINARY_SEARCH_U << ":" << EXPRESSION << "-" << TARGET << ( EXPRESSION > TARGET ? ">0" : EXPRESSION < TARGET ? "<0" : "0" ) ); \ } else { \ CERR( "二分探索失敗: " << MINIMUM << ">" << MAXIMUM << "です。" ); \ } \ // 単調増加の時にEXPRESSION >= TARGETの最小解を格納。 #define BS1( ANSWER , MINIMUM , MAXIMUM , EXPRESSION , TARGET ) \ BS( ANSWER , MINIMUM , MAXIMUM , EXPRESSION , TARGET , >= , ANSWER , ANSWER + 1 , ( VARIABLE_FOR_BINARY_SEARCH_L + VARIABLE_FOR_BINARY_SEARCH_U ) / 2 ) \ // 単調増加の時にEXPRESSION <= TARGETの最大解を格納。 #define BS2( ANSWER , MINIMUM , MAXIMUM , EXPRESSION , TARGET ) \ BS( ANSWER , MINIMUM , MAXIMUM , EXPRESSION , TARGET , > , ANSWER - 1 , ANSWER , ( VARIABLE_FOR_BINARY_SEARCH_L + 1 + VARIABLE_FOR_BINARY_SEARCH_U ) / 2 ) \ // 単調減少の時にEXPRESSION >= TARGETの最大解を格納。 #define BS3( ANSWER , MINIMUM , MAXIMUM , EXPRESSION , TARGET ) \ BS( ANSWER , MINIMUM , MAXIMUM , EXPRESSION , TARGET , < , ANSWER - 1 , ANSWER , ( VARIABLE_FOR_BINARY_SEARCH_L + 1 + VARIABLE_FOR_BINARY_SEARCH_U ) / 2 ) \ // 単調減少の時にEXPRESSION <= TARGETの最小解を格納。 #define BS4( ANSWER , MINIMUM , MAXIMUM , EXPRESSION , TARGET ) \ BS( ANSWER , MINIMUM , MAXIMUM , EXPRESSION , TARGET , <= , ANSWER , ANSWER + 1 , ( VARIABLE_FOR_BINARY_SEARCH_L + VARIABLE_FOR_BINARY_SEARCH_U ) / 2 ) \ // 圧縮用 #define TE template #define TY typename #define US using #define ST static #define IN inline #define CL class #define PU public #define OP operator #define CE constexpr #define CO const #define NE noexcept #define RE return #define WH while #define VO void #define VE vector #define LI list #define BE begin #define EN end #define SZ size #define MO move #define TH this #define CRI CO int& #define CRUI CO uint& #define CRL CO ll& #define ASK_DETAILS( ... ) \ CERR( "問題の区分は以下の中で何番に該当しますか?" ); \ problems = { __VA_ARGS__ }; \ problems_size = problems.size(); \ FOR( i , 0 , problems_size ){ \ CERR( i << ": " << problems[i] ); \ } \ cin >> num; \ CERR( "" ); \ num_temp = 0; \ if( num < 0 || num >= problems_size ){ \ CERR( "返答は" << problems_size - 1 << "以下の非負整数にしてください。" ); \ CERR( "終了します。" ); \ CERR( "" ); \ return -1; \ } \ int LibrarySearch( int num = -1 ) { vector problems{}; int problems_size = 13; int num_temp = 0; string reply{}; if( num == -1 ){ CERR( "ライブラリーを探索しますか?[y/n]" ); cin >> reply; if( reply == "n" ){ CERR( "ライブラリーを探索せずに続行します。" ); CERR( "" ); return 0; } else if( reply != "y" ){ CERR( "y/nのいずれかで答えてください。" ); CERR( "終了します。" ); CERR( "" ); return -1; } CERR( "" ); CERR( "ライブラリーを探索します。" ); ASK_DETAILS( "数に関する問題。" , "配列に関する問題。" , "文字列に関する問題。" , "順列に関する問題。" , "矩形領域に関する問題。" , "グラフに関する問題。" , "部分和問題。" , "確率/期待値に関する問題。" , "ゲームに関する問題。" , "論理に関する問題。" , "半順序集合に関する問題。" , "関数適用に関する問題。" , "構築問題。" ); } else { CERR( "" ); } CEXPR( int , num_graph , 5 ); CEXPR( int , num_subsequence_sum , 6 ); CEXPR( int , num_game , 8 ); if( num == num_temp++ ){ CERR( "入力は1つの数か、1つの数と法を表す数ですか?[y/n/c]" ); cin >> reply; CERR( "" ); if( reply == "y" ){ CERR( "まずは小さい入力の場合を愚直に計算し、OEISで検索しましょう。" ); CERR( "https://oeis.org/?language=japanese" ); CERR( "" ); CERR( "次に出力の定義と等価な式を考察しましょう。" ); CERR( "- 単調ならば、冪乗や階乗" ); CERR( "- 定義にp進法が使われていれば、各種探索アルゴリズム" ); CERR( "- 入力が素数に近い場合に規則性があれば、p進付値、p進法、" ); CERR( " オイラー関数、約数の個数など" ); CERR( "を検討しましょう。" ); } else if( reply == "n" ){ ASK_DETAILS( "求解問題。" , "最大化問題。" , "数え上げ問題。" , "その他。" ); if( num == num_temp++ ){ CERR( "まず" ); CERR( "- 単調関数は二分探索" ); CERR( "- 可微分関数はニュートン法" ); CERR( "- 一次関数は掃き出し法" ); CERR( "を検討しましょう。" ); CERR( "" ); CERR( "それらが間に合わない時は" ); CERR( "- 探索範囲を絞れそうな場合は全探策" ); CERR( "- 多変数の合成関数で表せる場合は半分全列挙" ); CERR( "を検討しましょう。" ); } else if( num == num_temp++ ){ CERR( "- 凸関数は三分探索" ); CERR( "- 可微分関数は微分のニュートン法" ); CERR( "を検討しましょう。" ); } else if( num == num_temp++ ){ CERR( "- 変数の対称性があれば大小関係を制限した全探策" ); CERR( "- 何らかの約数となるなど動く範囲が狭い変数があればそれらを決め打った全探策" ); CERR( "- 多変数の合成関数で表せる場合は半分全列挙" ); CERR( "を検討しましょう。" ); } else if( num == num_temp++ ){ CERR( "このケースのライブラリー探索は不完全です。" ); } } else { CERR( "y/nのいずれかで答えてください。" ); CERR( "終了します。" ); CERR( "" ); return -1; } CERR( "" ); CERR( "マルチテストケースの場合は以下の前計算を検討しましょう;" ); CERR( "素数列挙、約数列挙、サブゴールとなる関係式を満たす解列挙。" ); } else if( num == num_temp++ ){ ASK_DETAILS( "区間処理問題。" , "最大化問題。" , "最長部分列問題。" , "数え上げ問題。" , "部分和問題。" ); if( num == num_temp++ ){ ASK_DETAILS( "可換群構造+を使う問題。" , "可換羃等モノイド構造∨を使う問題。" , "モノイド構造*を使う問題。" , "非結合的マグマ構造*を使う問題。" , "集合へのマグマ作用(*,\\cdot)を使う問題。" , "モノイドへのマグマ作用(+,\\cdot)を使う問題。" , "定数とのmaxを取った値の区間和取得を使う問題。" ); if( num == num_temp++ ){ CERR( "- 区間加算/区間取得が必要ならば可換群BIT" ); CERR( " \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\BIT\\Template" ); CERR( "- 一点代入/一点加算/区間取得が必要ならば可換群平方分割" ); CERR( " \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\SqrtDecomposition\\Template" ); CERR( "- 区間以外の領域で加算/全更新後の一点取得が必要ならば階差数列" ); CERR( " \\Mathematics\\SetTheory\\DirectProduct\\Tree\\DifferenceSeqeuence" ); CERR( "を検討しましょう。" ); } else if( num == num_temp++ ){ CERR( "- 一点代入/区間加算/一点取得/区間取得が必要ならば可換羃等モノイドBIT" ); CERR( " \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\BIT\\IntervalMax\\Template" ); CERR( "を検討しましょう。" ); } else if( num == num_temp++ ){ CERR( "- 一点代入/区間取得が必要ならばモノイドBIT" ); CERR( " \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\BIT\\Template\\Monoid" ); CERR( "- 一点加算/区間加算/一点取得/区間取得が必要ならばモノイド平方分割" ); CERR( " \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\SqrtDecomposition\\Template\\Monoid" ); CERR( "- 一点代入/一点取得/区間取得が必要ならばモノイドセグメント木" ); CERR( " \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\SegmentTree" ); CERR( "を検討しましょう。" ); } else if( num == num_temp++ ){ CERR( "- 写像のコード化" ); CERR( " \\Mathematics\\Function\\Encoder" ); CERR( "によりモノイドに帰着させることを検討しましょう。" ); } else if( num == num_temp++ ){ CERR( "- 一点作用/区間作用/一点取得が必要ならば双対平方分割" ); CERR( " \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\SqrtDecomposition\\Template\\Dual" ); CERR( "を検討しましょう。" ); } else if( num == num_temp++ ){ CERR( "- 区間代入/区間作用/区間加算/一点取得/区間取得が必要な場合は遅延評価平方分割" ); CERR( " \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\SqrtDecomposition\\Template\\LazyEvaluation" ); CERR( "を検討しましょう。" ); } else if( num == num_temp++ ){ CERR( "maxで全体更新でなく区間更新をする場合の汎用的な解法は分かりません。" ); CERR( "例えば区間が包含について単調でmaxを取る値も単調であれば全体更新と" ); CERR( "同様の処理ができます。状況に応じた考察をしましょう。" ); CERR( "" ); CERR( "maxで全体更新をする場合、maxを取る値は単調である場合に帰着できます。" ); CERR( "maxで全体更新をしない場合、つまりただmaxの区間和を処理するだけの場合、" ); CERR( "クエリの順番を入れ替えることができるので、単調な全体更新に帰着できます。" ); CERR( "従って以下では単調な全体更新の問題を考えます。" ); CERR( "" ); CERR( "maxを取る定数を変数化し、元の値との大小を表す{0,1}値の係数を考えます。" ); CERR( "すると区間作用前後の値は統一的にその係数と変数を使って表せます。" ); CERR( "配列の各成分の係数の値が変化するイベントとクエリをソートして管理し、" ); CERR( "クエリがイベントを跨ぐたびに係数を更新することを検討しましょう。" ); CERR( "" ); CERR( "例えばクエリB_qに対するmax(A_i,B_q)の区間和は、" ); CERR( "- 優先度つきキューA'={(A_i,i)|i}(構築O(N log N))" ); CERR( "- (B_q,q)_qをソートしたB'(構築O(Q log Q))" ); CERR( "- 長さNの数列C=(0,...,0)(構築O(N))" ); CERR( "を用意し、B'を前から探索して順に各クエリ(B_q,q)を処理します。" ); CERR( "具体的にはA'を前から探索して順にA_i \ class Interval ## MAX ## BIT \ { \ private: \ T m_init; \ T m_a[N]; \ T m_fenwick_0[N + 1]; \ T m_fenwick_1[N + 1]; \ \ public: \ inline Interval ## MAX ## BIT( const T& n ); \ inline Interval ## MAX ## BIT( const T& n , const T ( &a )[N] ); \ inline Interval ## MAX ## BIT( const T& n , T ( &&a )[N] ); \ \ inline const T& operator[]( const int& i ) const; \ inline const T& Get( const int& i ) const; \ T Interval ## MAX( const int& i_start , const int& i_final ) const; \ \ void Set( const int& i , const T& n ); \ void Set ## MAX( const int& i , const T& n ); \ void IntervalSet ## MAX( const int& i_start , const int& i_final , const T& n ); \ \ int BinarySearch( const T& n ) const; \ \ }; \ #define DEFINITION_OF_INTERVAL_MAX_BIT( MAX , INEQUALITY , OP ) \ template inline Interval ## MAX ## BIT::Interval ## MAX ## BIT( const T& n ) \ : m_init( n ) , m_a() , m_fenwick_0() , m_fenwick_1() \ { \ \ if( m_a[0] != m_init ){ \ \ for( int i = 0 ; i < N ; i++ ){ \ \ m_a[i] = m_init; \ \ } \ \ for( int j = 1 ; j <= N ; j++ ){ \ \ m_fenwick_0[j] = m_fenwick_1[j] = m_init; \ } \ \ } \ \ } \ \ template inline Interval ## MAX ## BIT::Interval ## MAX ## BIT( const T& n , const T ( &a )[N] ) : m_init( n ) , m_a() , m_fenwick_0() , m_fenwick_1() \ { \ \ for( int i = 0 ; i < N ; i++ ){ \ \ m_a[i] = a[i]; \ \ } \ \ for( int i = 0 ; i < N ; i++ ){ \ \ int j = i + 1; \ T& fenwick_0i = m_fenwick_0[j]; \ fenwick_0i = m_a[i]; \ const int j_llim = j - ( j & -j ); \ j--; \ \ while( j > j_llim ){ \ \ const T& tj = m_fenwick_0[j]; \ fenwick_0i INEQUALITY tj ? fenwick_0i = tj : fenwick_0i; \ j -= ( j & -j ); \ \ } \ \ } \ \ for( int i = N - 1 ; i >= 0 ; i-- ){ \ \ int j = i + 1; \ T& fenwick_1i = m_fenwick_1[j]; \ fenwick_1i = m_a[i]; \ const int j_ulim = min( j + ( j & -j ) , N + 1 ); \ j++; \ \ while( j < j_ulim ){ \ \ const T& tj = m_fenwick_1[j]; \ fenwick_1i INEQUALITY tj ? fenwick_1i = tj : fenwick_1i; \ j += ( j & -j ); \ \ } \ \ } \ \ } \ \ template inline Interval ## MAX ## BIT::Interval ## MAX ## BIT( const T& n , T ( &&a )[N] ) : m_init( n ) , m_a() , m_fenwick_0() , m_fenwick_1() \ { \ \ swap( m_a , a ); \ \ for( int i = 0 ; i < N ; i++ ){ \ \ int j = i + 1; \ T& fenwick_0i = m_fenwick_0[j]; \ fenwick_0i = m_a[i]; \ const int j_llim = j - ( j & -j ); \ j--; \ \ while( j > j_llim ){ \ \ const T& tj = m_fenwick_0[j]; \ fenwick_0i INEQUALITY tj ? fenwick_0i = tj : fenwick_0i; \ j -= ( j & -j ); \ \ } \ \ } \ \ for( int i = N - 1 ; i >= 0 ; i-- ){ \ \ int j = i + 1; \ T& fenwick_1i = m_fenwick_1[j]; \ fenwick_1i = m_a[i]; \ const int j_ulim = min( j + ( j & -j ) , N + 1 ); \ j++; \ \ while( j < j_ulim ){ \ \ const T& tj = m_fenwick_1[j]; \ fenwick_1i INEQUALITY tj ? fenwick_1i = tj : fenwick_1i; \ j += ( j & -j ); \ \ } \ \ } \ \ } \ \ template inline const T& Interval ## MAX ## BIT::operator[]( const int& i ) const { return m_a[i]; } \ template inline const T& Interval ## MAX ## BIT::Get( const int& i ) const { return m_a[i]; } \ \ template \ T Interval ## MAX ## BIT::Interval ## MAX( const int& i_start , const int& i_final ) const \ { \ \ T answer = m_init; \ const int j_min = i_start < 0 ? 1 : i_start + 1; \ const int j_max = i_final < N ? i_final + 1 : N; \ int j = j_min; \ int j_next = j + ( j & - j ); \ \ while( j_next <= j_max ){ \ \ const T& tj = m_fenwick_1[j]; \ answer INEQUALITY tj ? answer = tj : answer; \ j = j_next; \ j_next += ( j & -j ); \ \ } \ \ const T& a_centre = m_a[j-1]; \ ( j_min <= j_max && answer < a_centre ) ? answer = a_centre : answer; \ j = j_max; \ j_next = j - ( j & - j ); \ \ while( j_next >= j_min ){ \ \ const T& tj = m_fenwick_0[j]; \ answer INEQUALITY tj ? answer = tj : answer; \ j = j_next; \ j_next -= ( j & -j ); \ \ } \ \ return answer; \ \ } \ \ template \ void Interval ## MAX ## BIT::Set( const int& i , const T& n ) \ { \ \ T& ai = m_a[i]; \ \ if( n INEQUALITY ai ){ \ \ int j = i + 1; \ \ while( j <= N ){ \ \ const int lsb = ( j & -j ); \ m_fenwick_0[j] = OP( OP( Interval ## MAX( j - lsb + 1 , i - 1 ) , n ) , Interval ## MAX( i + 1 , j ) ); \ j += lsb; \ \ } \ \ j = i + 1; \ \ while( j > 0 ){ \ \ const int lsb = ( j & -j ); \ m_fenwick_0[j] = OP( OP( Interval ## MAX( j , i - 1 ) , n ) , Interval ## MAX( i + 1 , j + lsb - 1 ) ); \ j -= lsb; \ \ } \ \ ai = n; \ \ } else { \ \ Set ## MAX( i , n ); \ } \ \ return; \ \ } \ \ template \ void Interval ## MAX ## BIT::Set ## MAX( const int& i , const T& n ) \ { \ \ T& ai = m_a[i]; \ ai INEQUALITY n ? ai = n : ai; \ int j = i + 1; \ \ while( j <= N ){ \ \ T& tj = m_fenwick_0[j]; \ tj INEQUALITY n ? tj = n : tj; \ j += ( j & -j ); \ \ } \ \ j = i + 1; \ \ while( j > 0 ){ \ \ T& tj = m_fenwick_1[j]; \ tj INEQUALITY n ? tj = n : tj; \ j -= ( j & -j ); \ \ } \ \ return; \ \ } \ // 最大(最小)元による初期化O(N) // 配列による初期化O(N) // 一点取得O(1) // 区間max(min)取得O(log_2 N) // 一点更新O((log_2 N)^2) // max(min)による一点更新O(log_2 N) // max(min)による区間更新O(i_final-i_start+log_2 N) // t以上(以下)となる要素の添字の最小値の二分探索O(log_2 N) // そのうちの区間min取得と一点更新は // M. Dima, R. Ceterchi, Efficient Range Minimum Queries using Binary Indexed Trees, Olympiads in Informatics, 2015, Vol. 9, 39--44 // の手法をもとに実装 DECRALATION_OF_INTERVAL_MAX_BIT( Max ); DECRALATION_OF_INTERVAL_MAX_BIT( Min ); DEFINITION_OF_INTERVAL_MAX_BIT( Max , < , max ); DEFINITION_OF_INTERVAL_MAX_BIT( Min , > , min ); #define SFINAE_FOR_DOUBLING_BODY( DEFAULT ) enable_if_t >* DEFAULT template class DoublingBody { private: int m_length; map m_memory; vector m_memory_inv; protected: int m_size; int m_doubling[digit][size_max]; public: template inline DoublingBody( const int& size ); // n < 2のdigit乗 の場合のみサポート template T IteratedComposition( T t , INT n ); private: virtual T e( const int& i ); virtual int e_inv( const T& t ); }; template class Doubling : public DoublingBody { public: inline Doubling( const int& size ); private: inline int e( const int& i ); inline int e_inv( const int& t ); }; template class MemorisationDoubling : public DoublingBody { public: inline MemorisationDoubling( const int& size ); }; template class EnumerationDoubling : public DoublingBody { public: inline EnumerationDoubling( const int& size ); private: inline T e( const int& i ); inline int e_inv( const T& t ); }; template template inline DoublingBody::DoublingBody( const int& size ) : m_length() , m_memory() , m_memory_inv() , m_size( size ) , m_doubling() {} template inline Doubling::Doubling( const int& size ) : DoublingBody( size ) { using base = DoublingBody; for( int d = 0 ; d < digit ; d++ ){ int ( &doubling_d )[size_max] = base::m_doubling[d]; for( int i = 0 ; i < base::m_size ; i++ ){ doubling_d[i] = -1; } } } template inline MemorisationDoubling::MemorisationDoubling( const int& size ) : DoublingBody( size ) { using base = DoublingBody; for( int d = 0 ; d < digit ; d++ ){ int ( &doubling_d )[size_max] = base::m_doubling[d]; for( int i = 0 ; i < base::m_size ; i++ ){ doubling_d[i] = -1; } } } template inline EnumerationDoubling::EnumerationDoubling( const int& size ) : DoublingBody( size ) { using base = DoublingBody; for( int d = 0 ; d < digit ; d++ ){ int ( &doubling_d )[size_max] = base::m_doubling[d]; // enum_Tがm_size未満への全単射とは限らないのでsize_maxまで初期化する。 for( int i = 0 ; i < size_max ; i++ ){ doubling_d[i] = -1; } } } template template T DoublingBody::IteratedComposition( T t , INT n ) { int i = e_inv( t ); int d = 0; while( n != 0 ){ assert( d < digit ); int ( &doubling_d )[size_max] = m_doubling[d]; const int& doubling_d_i = doubling_d[i]; if( doubling_d_i == -1 ){ int j = i; if( d == 0 ){ while( doubling_d[j] == -1 ){ j = doubling_d[j] = e_inv( t = f( t ) ); } } else { int ( &doubling_d_minus )[size_max] = m_doubling[d - 1]; while( doubling_d[j] == -1 ){ j = doubling_d[j] = doubling_d_minus[doubling_d_minus[j]]; } } } ( n & 1 ) == 1 ? i = doubling_d_i : i; n >>= 1; d++; } return e( i ); } template T DoublingBody::e( const int& i ) { assert( i < m_length ); return m_memory_inv[i]; } template inline int Doubling::e( const int& i ) { return i; } template inline T EnumerationDoubling::e( const int& i ) { return enum_T( i ); } template int DoublingBody::e_inv( const T& t ) { if( m_memory.count( t ) == 0 ){ assert( m_length < m_size ); m_memory_inv.push_back( t ); return m_memory[t] = m_length++; } return m_memory[t]; } template inline int Doubling::e_inv( const int& t ) { return t; } template inline int EnumerationDoubling::e_inv( const T& t ) { return enum_T_inv( t ); } inline CEXPR( int , bound_N , 200000 ); // 0が5個 IntervalMaxBIT* pT; inline int f( const int& Ti ) { return max( Ti , pT->IntervalMax( 0 , Ti ) ); } int main() { UNTIE; LIBRARY_SEARCH; // CEXPR( int , bound_T , 100000 ); // CIN_ASSERT( T , 1 , bound_T ); // CEXPR( int , bound_N , 100000 ); // 0が5個 // CEXPR( ll , bound_N , 1000000000 ); // 0が9個 // CEXPR( ll , bound_N , 1000000000000000000 ); // 0が18個 CIN_ASSERT( N , 2 , bound_N ); CEXPR( int , bound_HT , 1000000000 ); // 0が9個 static tuple HTI[bound_N + 1]; FOREQ( i , 1 , N ){ CIN_ASSERT( Hi , 1 , bound_HT ); HTI[i] = { Hi , 0 , i }; } FOREQ( i , 1 , N ){ CIN_ASSERT( Ti , 1 , bound_HT ); get<1>( HTI[i] ) = Ti; } sort( HTI + 1 , HTI + N + 1 ); HTI[0] = { 0 , 0 , 0 }; static int T_prep[bound_N + 1]; T_prep[0] = 0; FOREQ( i , 1 , N ){ int& Ti = get<1>( HTI[i] ); BS2( j , 0 , N , get<0>( HTI[j] ) , Ti ); T_prep[i] = j; } static IntervalMaxBIT T{ 0 , move( T_prep ) }; pT = &T; FOREQ( i , 1 , N ){ get<0>( HTI[get<2>( HTI[i] )] ) = i; } static Doubling d{ N + 1 }; CEXPR( int , bound_Q , 200000 ); CIN_ASSERT( Q , 1 , bound_Q ); // CEXPR( int , bound_M , 100000 ); // 0が5個 // // CEXPR( ll , bound_M , 1000000000 ); // 0が9個 // // CEXPR( ll , bound_M , 1000000000000000000 ); // 0が18個 // CIN_ASSERT( M , 0 , bound_M ); REPEAT( Q ){ CIN_ASSERT( Aq , 1 , N ); CIN_ASSERT( Bq , 1 , N ); int& iq = get<0>( HTI[Aq] ); int& jq = get<0>( HTI[Bq] ); const int& Tiq = T[iq]; CERR( iq << "," << jq << "," << Tiq ); if( Tiq >= jq ){ CERR( "accessible at once: " << Tiq << ">=" << jq ); COUT( 1 ); } else { if( d.IteratedComposition( Tiq , N ) < jq ){ CERR( "inaccessible: " << d.IteratedComposition( Tiq , N ) << "<" << jq ); COUT( -1 ); } else { BS1( n , 0 , N , d.IteratedComposition( Tiq , n ) , jq ); CERR( "accessible: " << d.IteratedComposition( Tiq , n ) << ">=" << jq ); COUT( ++n ); } } } // REPEAT( T ){ // COUT( N ); // } QUIT; }