結果
| 問題 | No.2242 Cities and Teleporters |
| ユーザー |
👑  p-adic
|
| 提出日時 | 2023-06-28 21:31:14 |
| 言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
| 結果 |
AC
|
| 実行時間 | 943 ms / 3,000 ms |
| コード長 | 45,856 bytes |
| コンパイル時間 | 6,604 ms |
| コンパイル使用メモリ | 256,492 KB |
| 実行使用メモリ | 24,576 KB |
| 最終ジャッジ日時 | 2024-07-05 15:00:10 |
| 合計ジャッジ時間 | 18,734 ms |
|
ジャッジサーバーID （参考情報） |
judge1 / judge3 |
（要ログイン）
テストケース
テストケース表示| 入力 | 結果 | 実行時間 実行使用メモリ |
|---|---|---|
| testcase_00 | AC | 17 ms
22,272 KB |
| testcase_01 | AC | 18 ms
22,144 KB |
| testcase_02 | AC | 19 ms
22,272 KB |
| testcase_03 | AC | 18 ms
22,272 KB |
| testcase_04 | AC | 18 ms
22,144 KB |
| testcase_05 | AC | 619 ms
24,448 KB |
| testcase_06 | AC | 614 ms
24,576 KB |
| testcase_07 | AC | 606 ms
24,448 KB |
| testcase_08 | AC | 943 ms
24,576 KB |
| testcase_09 | AC | 773 ms
24,448 KB |
| testcase_10 | AC | 337 ms
24,576 KB |
| testcase_11 | AC | 263 ms
24,576 KB |
| testcase_12 | AC | 267 ms
24,448 KB |
| testcase_13 | AC | 367 ms
24,576 KB |
| testcase_14 | AC | 529 ms
24,576 KB |
| testcase_15 | AC | 363 ms
24,576 KB |
| testcase_16 | AC | 509 ms
24,576 KB |
| testcase_17 | AC | 767 ms
24,576 KB |
| testcase_18 | AC | 214 ms
24,576 KB |
| testcase_19 | AC | 233 ms
24,576 KB |
| testcase_20 | AC | 356 ms
24,448 KB |
| testcase_21 | AC | 289 ms
24,320 KB |
| testcase_22 | AC | 214 ms
24,320 KB |
| testcase_23 | AC | 224 ms
24,576 KB |
| testcase_24 | AC | 216 ms
24,448 KB |
| testcase_25 | AC | 213 ms
24,576 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( 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 <bits/stdc++.h>
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<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 ) 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 <typename T> inline T Absolute( const T& a ){ return a > 0 ? a : -a; }
template <typename T> 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<TYPE_OF( ARGUMENT ),int>::value && ! is_same<TYPE_OF( ARGUMENT ),uint>::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<TYPE_OF( TARGET ),uint>::value && ! is_same<TYPE_OF( TARGET ),ull>::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<string> 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( "入力は１つの数か、１つの数と法を表す数ですか？[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<B_qとなる各iに対し" );
CERR( "- A'から(A_i,i)を削除（クエリ合計O(N)）" );
CERR( "- A_iを0に更新（クエリ合計O(N log N)）" );
CERR( "- C_iを1に更新（各クエリO(log N)）" );
CERR( "- A+C×B_qの区間和取得（各クエリO(log N)）" );
CERR( "とすれば合計O((N + Q)log N + Q log Q)で処理できます。" );
}
} else if( num == num_temp++ ){
CERR( "最大化すべき式のサブゴールfに表れる項xのうち決め打ちやすいものを探しましょう。" );
CERR( "- 配列の長さをiで打ち切った時のxの候補数をX(i)" );
CERR( "- 配列の長さをiで打ち切ってxを決め打った時の配列の長さi+1でのxの候補数をdX(i)" );
CERR( "と置きます。" );
CERR( "- O(sum_i X(i) dX(i))が通りそうでfがxからO(1)で計算できるならば、iとxに関する動的計画法" );
CERR( "- O(N log_2 X(N))が通りそうで" );
CERR( " - fがxからO(N)で計算できxに関して単調ならば、xの二分探索" );
CERR( " - fがxからO(N/x)で計算できるならば、xの全探索" );
CERR( "- O(N log_2 N)が通りそうでxを上手く並び替えるとfがxからO(log_2 N)で計算できるならば、" );
CERR( " 優先度つきキューなどでのxの管理" );
CERR( "を検討しましょう。" );
} else if( num == num_temp++ ){
CERR( "全順序か疎な半順序かで効率的な実装が違います。" );
CERR( "- 全順序ならば、条件を満たす部分列の長さの最大値をインデックスに持つ配列を用いて、" );
CERR( " それらの部分列の末尾である項を記録すること" );
CERR( "- 疎な半順序ならば、条件を満たす部分列の末尾をインデックスに持つ連想配列を用いて、" );
CERR( " それら部分列の長さの最大値を記録すること" );
CERR( "を検討しましょう。" );
} else if( num == num_temp++ ){
CERR( "gcdやmaxなどの羃等演算に関する等式を指定した数え上げは不等式の方が" );
CERR( "扱いやすいのでゼータ変換／メビウス変換を検討しましょう。" );
} else if( num == num_temp++ ){
return LibrarySearch( num = num_subsequence_sum );
}
} else if( num == num_temp++ ){
ASK_DETAILS(
"部分列マッチング問題。" ,
"回文探索問題。"
);
if( num == num_temp++ ){
CERR( "基本的には丁寧にループを回して解きましょう。" );
CERR( "- 比較対象が少ない場合、前または後ろから順に探索（貪欲法／動的計画法）" );
CERR( "- ワイルドカードを含む場合、" );
CERR( " - 前または後ろから順に場合分けをしてO(N)で処理できるか" );
CERR( " - 可能な代入方法を絞り込んでO(N)種類に落せるか" );
CERR( "- 比較回数が多い場合、ローリングハッシュ" );
CERR( " \\Utility\\String\\RollingHash" );
CERR( "- マッチングする文字列の最長化をする場合、Zアルゴリズム" );
CERR( " https://qiita.com/Pro_ktmr/items/16904c9570aa0953bf05" );
CERR( "を検討しましょう。" );
} else if( num == num_temp++ ){
CERR( "回文判定は長さに関して再帰的に計算できます。" );
CERR( "- O(N^2)が通る場合、愚直な再帰により前計算で全ての部分列の回文判定" );
CERR( "- O(N^2)が通らない場合、Manacherのアルゴリズムやローリングハッシュで前計算" );
CERR( " https://snuke.hatenablog.com/entry/2014/12/02/235837" );
CERR( "を検討しましょう。" );
}
} else if( num == num_temp++ ){
CERR( "符号の計算は転倒数の計算に帰着させましょう。" );
CERR( "符号と何かの積の和は行列式に帰着させましょう。" );
CERR( "- 行列式そのものなら行基本変形でO(N^3)" );
CERR( "- 余因子展開の途中の値はメモ化再帰でO(N 2^N)" );
CERR( "で計算できます。" );
CERR( "" );
CERR( "１つの順列の転倒数は、" );
CERR( "- O(N^2)が通りそうならば愚直な二重ループ" );
CERR( "- O(N log_2 N)が通りそうならば可換群BIT" );
CERR( " \\Mathematics\\Combinatorial\\Permutation" );
CERR( " \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\BIT" );
CERR( "で計算しましょう。" );
CERR( "" );
CERR( "条件を満たす順列全体をわたる転倒数の総和は、" );
CERR( "各i<jごとにそこで転倒が生じる順列の個数を計算し、その総和を取りましょう。" );
CERR( "条件が良ければ、転倒が生じる順列の個数は転倒が生じるとは限らない順列の個数の" );
CERR( "半分となります。" );
} else if( num == num_temp++ ){
ASK_DETAILS(
"探索問題。" ,
"絶対値の最小／最大化問題。" ,
"その他の最小／最大化問題。"
);
if( num == num_temp++ ){
CERR( "隣接関係の定める無向グラフの問題に帰着させましょう。" );
return LibrarySearch( num = num_graph );
} else if( num == num_temp++ ){
CERR( "符号を用いて絶対値を外しましょう。" );
CERR( "- 単調な式に帰着できる場合、二分探索を検討しましょう。" );
CERR( "- 最大化問題の場合、符号パターンの全探策を検討しましょう。" );
CERR( "- マンハッタン距離などは一次変換で簡単にすることを検討しましょう。" );
CERR( "複数のパラメータを決定すべき場合は、サブゴールの式の値を決め打ちましょう。" );
} else if( num == num_temp++ ){
CERR( "このケースのライブラリー探索は不完全です。" );
CERR( "- O(HW)が通りそうならば動的計画法" );
CERR( "- O(HW)が通らなさそうならば見方を変えて最短経路問題に帰着できないか" );
CERR( "を検討しましょう。" );
CERR( "" );
CERR( "例えば迷路の攻略可能性は" );
CERR( "- スタートとゴールが同一の弧状連結成分に属すこと" );
CERR( "- スタートとゴールを分断する壁のパスの非存在性" );
CERR( "などから翻訳できます。" );
}
} else if( num == num_temp++ ){
ASK_DETAILS(
"２点最短径路（迷路）問題。" ,
"多点最短経路（スタンプラリー）問題。" ,
"木の問題。" ,
"連結性問題。"
);
if( num == num_temp++ ){
CERR( "特定の経路を進むと思い込んで考察漏れをする可能性があります。" );
CERR( "なるべく全ての経路を許した探索アルゴリズムを適用した方が無難です。" );
CERR( "- 特定の２点のみを考える場合、BFSやDijkstra" );
CERR( " \\Utility\\Search\\BreadthFirst" );
CERR( " \\Utility\\Search\\Dijkstra" );
CERR( "- 全ての２点の組み合わせを考える場合、" );
CERR( " - 一般のモノイド演算を考えておりO(V^3)が通りそうならば、FloydWarshall" );
CERR( " \\Utility\\Search\\FloydWarshall" );
CERR( " - max演算を考えておりO(E(log_2 E + α(V)))が通りそうならば、UnionFind" );
CERR( " \\Utility\\VLTree\\UnionFindForest" );
CERR( "を検討しましょう。" );
} else if( num == num_temp++ ){
CERR( "HeldKarpや、移動方法を分類するパラメータの全探策などを検討しましょう。" );
} else if( num == num_temp++ ){
CERR( "深さ優先探索や動的木を検討しましょう。" );
CERR( "\\Utility\\Search\\DepthFirst" );
CERR( "\\Utility\\VLTree" );
} else if( num == num_temp++ ){
CERR( "- 0次の連結性はUnionFind" );
CERR( " \\Utility\\VLTree\\UnionFindForest" );
CERR( "- 高次の連結性は深さ優先探索" );
CERR( " \\Utility\\Search\\DepthFirst" );
CERR( "- マルチテストケースの場合は座標圧縮との併用" );
CERR( "を検討しましょう。" );
}
} else if( num == num_temp++ ){
ASK_DETAILS(
"最大化問題。" ,
"数え上げ問題。"
);
if( num == num_temp++ ){
CERR( "項数N、重みの総和の上限W（多変数も可）、価値（和を取る値）の上限Vとします。" );
CERR( "- B=∞ならば、通常のナップサック問題と同様の動的計画法" );
CERR( "- B<∞でO(2^N)が通りそうならば愚直に全探策" );
CERR( "- B<∞でO(2^{N/2} N)が通りそうならば半分全列挙" );
CERR( "- B<∞で重みと価値が等しくO(NV)が通りそうならば[B-V,B+V]での実現可能性を" );
CERR( " 遷移する動的計画法" );
CERR( " https://stackoverflow.com/a/18949218" );
CERR( "- Wが10^5オーダーで重みと価値が等しくO((N+W)log_2 W)が通りそうならば" );
CERR( " 適当な法での畳み込み（確率的解法）" );
CERR( " \\Mathematics\\Polynomial" );
} else if( num == num_temp++ ){
CERR( "項数N、重みの総和の上限Wとします。" );
CERR( "- 重みと価値が異なりO(2^N)が通りそうならば愚直な２変数多項式乗算" );
CERR( "- 重みと価値が等しくO(2^N)が通りそうならば愚直な１変数多項式乗算" );
CERR( "- 重みと価値が等しくO(2^{N/2}N)が通りそうならば半分ずつ多項式乗算を" );
CERR( " して最後にそれらの積の１係数のみの計算" );
CERR( "- 重みと価値が等しくWが10^5オーダーで法が与えられていてO((N+W)log_2 W)が" );
CERR( " 通りそうならば畳み込み" );
CERR( " \\Mathematics\\Polynomial" );
}
CERR( "を検討しましょう。" );
} else if( num == num_temp++ ){
CERR( "確率計算は" );
CERR( "- 余事象や包除原理（高速ゼータ変換）" );
CERR( "- 同様に確からしい事象の特定" );
CERR( "- ベイズの定理" );
CERR( "を検討しましょう。" );
CERR( "" );
CERR( "期待値計算は" );
CERR( "- 上記方法での確率計算" );
CERR( "- 対象を独立な和で表して線形性" );
CERR( "を検討しましょう。" );
} else if( num == num_temp++ ){
CERR( "ゲームの和に分解できる場合は最小単位で考察をし、グランディ数を実装しましょう。" );
CERR( "これ以上分解できないゲームには整礎構造を探し、順序数の小さい順に実験をしましょう。" );
} else if( num == num_temp++ ){
CERR( "複数の命題に対する" );
CERR( "- 区間処理は各種代数的データ構造" );
CERR( " \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace" );
CERR( "- 充足可能性判定は頂点数を２倍してUnionFind" );
CERR( " \\Utility\\VLTree\\UnionFindForest" );
CERR( "を検討しましょう。" );
} else if( num == num_temp++ ){
ASK_DETAILS(
"関数の計算問題。" ,
"降下／上昇列に関する問題。"
);
if( num == num_temp++ ){
CERR( "ゼータ変換を検討しましょう。" );
} else if( num == num_temp++ ){
CERR( "順序関係の定める有向グラフの問題に帰着させましょう。" );
return LibrarySearch( num = num_graph );
}
} else if( num == num_temp++ ){
ASK_DETAILS(
"反復合成の計算問題。" ,
"反復合成による到達可能性問題。"
);
if( num == num_temp++ ){
CERR( "定義域の要素数N、テストケース数T、反復回数の上限Kとします。" );
CERR( "- O((N + T)log_2 K)が通りそうならばダブリング" );
CERR( " \\Mathematics\\Function\\Iteration\\Doubling" );
CERR( "- O(TN)が通りそうならばループ検出" );
CERR( " \\Mathematics\\Function\\Iteration\\LoopDetection" );
CERR( "- O(N)すら通らなさそうならば関数の規則性を見付けるための実験" );
CERR( "を検討しましょう。" );
} else if( num == num_temp++ ){
CERR( "関数による遷移が定める有向グラフの問題に帰着させましょう。" );
return LibrarySearch( num = num_graph );
}
} else if( num == num_temp++ ){
ASK_DETAILS(
"数や配列や文字列の構築。" ,
"経路の構築。" ,
"戦略の構築。"
);
if( num == num_temp++ ){
CERR( "p進法を検討しましょう。" );
} else if( num == num_temp++ ){
CERR( "可能な経路の定める無向グラフの問題に帰着させましょう。" );
return LibrarySearch( num = num_graph );
} else if( num == num_temp++ ){
CERR( "ゲームの問題に帰着させましょう。" );
return LibrarySearch( num = num_game );
}
}
CERR( "" );
CERR( "ライブラリー探索は以上です。終了します。" );
CERR( "" );
return -1;
}
#define DECRALATION_OF_INTERVAL_MAX_BIT( MAX ) \
template <typename T , int N> \
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 <typename T , int N> inline Interval ## MAX ## BIT<T,N>::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 <typename T , int N> inline Interval ## MAX ## BIT<T,N>::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 <typename T , int N> inline Interval ## MAX ## BIT<T,N>::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 <typename T , int N> inline const T& Interval ## MAX ## BIT<T,N>::operator[]( const int& i ) const { return m_a[i]; } \
template <typename T , int N> inline const T& Interval ## MAX ## BIT<T,N>::Get( const int& i ) const { return m_a[i]; } \
\
template <typename T , int N> \
T Interval ## MAX ## BIT<T,N>::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 <typename T , int N> \
void Interval ## MAX ## BIT<T,N>::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 <typename T , int N> \
void Interval ## MAX ## BIT<T,N>::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<is_convertible_v<U,T> >* DEFAULT
template <typename T , typename U , U f(const T&) , int size_max , int digit>
class DoublingBody
{
private:
int m_length;
map<T,int> m_memory;
vector<T> m_memory_inv;
protected:
int m_size;
int m_doubling[digit][size_max];
public:
template <SFINAE_FOR_DOUBLING_BODY( = nullptr )> inline DoublingBody( const int& size );
// n < 2のdigit乗 の場合のみサポート
template <typename INT> T IteratedComposition( T t , INT n );
private:
virtual T e( const int& i );
virtual int e_inv( const T& t );
};
template <int f(const int&) , int size_max , int digit = 64>
class Doubling :
public DoublingBody<int,int,f,size_max,digit>
{
public:
inline Doubling( const int& size );
private:
inline int e( const int& i );
inline int e_inv( const int& t );
};
template <typename T, typename U , U f(const T&) , int size_max , int digit = 64>
class MemorisationDoubling :
public DoublingBody<T,U,f,size_max,digit>
{
public:
inline MemorisationDoubling( const int& size );
};
template <typename T, typename U , U f(const T&) , int size_max , T enum_T(const int&) , int enum_T_inv(const T&) , int digit = 64>
class EnumerationDoubling :
public DoublingBody<T,U,f,size_max,digit>
{
public:
inline EnumerationDoubling( const int& size );
private:
inline T e( const int& i );
inline int e_inv( const T& t );
};
template <typename T, typename U , U f(const T&) , int size_max , int digit> template <SFINAE_FOR_DOUBLING_BODY()> inline DoublingBody<T,U,f,size_max,digit>::DoublingBody( const int& size ) : m_length() , m_memory() , m_memory_inv() , m_size( size ) , m_doubling() {}
template <int f(const int&) , int size_max , int digit> inline Doubling<f,size_max,digit>::Doubling( const int& size ) :
DoublingBody<int,int,f,size_max,digit>( size )
{
using base = DoublingBody<int,int,f,size_max,digit>;
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 <typename T, typename U , U f(const T&) , int size_max , int digit> inline MemorisationDoubling<T,U,f,size_max,digit>::MemorisationDoubling( const int& size ) :
DoublingBody<T,U,f,size_max,digit>( size )
{
using base = DoublingBody<T,U,f,size_max,digit>;
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 <typename T, typename U , U f(const T&) , int size_max , T enum_T(const int&) , int enum_T_inv(const T&) , int digit> inline EnumerationDoubling<T,U,f,size_max,enum_T,enum_T_inv,digit>::EnumerationDoubling( const int& size ) :
DoublingBody<T,U,f,size_max,digit>( size )
{
using base = DoublingBody<T,U,f,size_max,digit>;
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 <typename T, typename U , U f(const T&) , int size_max , int digit> template <typename INT>
T DoublingBody<T,U,f,size_max,digit>::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 <typename T, typename U , U f(const T&) , int size_max , int digit>
T DoublingBody<T,U,f,size_max,digit>::e( const int& i )
{
assert( i < m_length );
return m_memory_inv[i];
}
template <int f(const int&) , int size_max , int digit> inline int Doubling<f,size_max,digit>::e( const int& i ) { return i; }
template <typename T, typename U , U f(const T&) , int size_max , T enum_T(const int&) , int enum_T_inv(const T&) , int digit> inline T EnumerationDoubling<T,U,f,size_max,enum_T,enum_T_inv,digit>::e( const int& i ) { return enum_T( i ); }
template <typename T, typename U , U f(const T&) , int size_max , int digit>
int DoublingBody<T,U,f,size_max,digit>::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 <int f(const int&) , int size_max , int digit> inline int Doubling<f,size_max,digit>::e_inv( const int& t ) { return t; }
template <typename T, typename U , U f(const T&) , int size_max , T enum_T(const int&) , int enum_T_inv(const T&) , int digit> inline int EnumerationDoubling<T,U,f,size_max,enum_T,enum_T_inv,digit>::e_inv( const T& t ) { return enum_T_inv( t ); }
inline CEXPR( int , bound_N , 200000 ); // 0が5個
IntervalMaxBIT<int,bound_N + 1>* 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<int,int,int> 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<int,bound_N + 1> T{ 0 , move( T_prep ) };
pT = &T;
FOREQ( i , 1 , N ){
get<0>( HTI[get<2>( HTI[i] )] ) = i;
}
static Doubling<f,bound_N + 1,20> 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;
}