#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 ) ) #define AUTO_CHECK bool auto_checked = true; AutoCheck( auto_checked ); if( auto_checked ){ return 0; }; #define START_WATCH( PROCESS_NAME ) StartWatch( PROCESS_NAME ) #define STOP_WATCH( HOW_MANY_TIMES ) StopWatch( HOW_MANY_TIMES ) #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 ) ) #define AUTO_CHECK #define START_WATCH( PROCESS_NAME ) #define STOP_WATCH( HOW_MANY_TIMES ) #endif // #define RANDOM_TEST #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 ) TYPE_OF( MAX ) A; SET_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 SET_PRECISION( DECIMAL_DIGITS ) cout << fixed << setprecision( DECIMAL_DIGITS ) #define QUIT goto END_MAIN #define TEST_CASE_NUM( BOUND ) DEXPR( int , bound_T , BOUND , min( BOUND , 100 ) ); int T = 1; if constexpr( bound_T > 1 ){ SET_ASSERT( T , 1 , bound_T ); } #define START_MAIN REPEAT( T ){ if constexpr( bound_T > 1 ){ CERR( "testcase " << VARIABLE_FOR_REPEAT_T << ":" ); } #define FINISH_MAIN QUIT; } END_MAIN: CERR( "" ); #ifdef DEBUG inline void AlertAbort( int n ) { CERR( "abort関数が呼ばれました。assertマクロのメッセージが出力されていない場合はオーバーフローの有無を確認をしてください。" ); } void StartWatch( const string& process_name = "nothing" ); void StopWatch( const int& how_many_times = 1 ); #endif #if defined( DEBUG ) && defined( RANDOM_TEST ) ll GetRand( const ll& Rand_min , const ll& Rand_max ); #define SET_ASSERT( A , MIN , MAX ) CERR( #A << " = " << ( A = GetRand( MIN , MAX ) ) ) #define RETURN( ANSWER ) if( ( ANSWER ) == guchoku ){ CERR( ( ANSWER ) << " == " << guchoku ); goto END_MAIN; } else { CERR( ( ANSWER ) << " != " << guchoku ); QUIT; } #else #define SET_ASSERT( A , MIN , MAX ) cin >> A; ASSERT( A , MIN , MAX ) #define RETURN( ANSWER ) COUT( ( ANSWER ) ); QUIT #endif // グリッド問題用関数 int H , W , H_minus , W_minus , HW; inline pair EnumHW( const int& v ) { return { v / W , v % W }; } inline int EnumHW_inv( const int& h , const int& w ) { return h * W + w; } const string direction[4] = {"U","R","D","L"}; // (i,j)->(k,h)の方向番号を取得 inline int DirectionNumberOnGrid( const int& i , const int& j , const int& k , const int& h ){return ik?0:jh?3:(assert(false),-1);} // v->wの方向番号を取得 inline int DirectionNumberOnGrid( const int& v , const int& w ){auto [i,j]=EnumHW(v);auto [k,h]=EnumHW(w);return DirectionNumberOnGrid(i,j,k,h);} // 方向番号の反転U<->D、R<->L inline int ReverseDirectionNumberOnGrid( const int& n ){assert(0<=n&&n<4);return(n+2)%4;} // 圧縮用 #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_NUMBER( ... ) \ CERR( "" ); \ CERR( "何番を選択しますか?" ); \ problems = { __VA_ARGS__ }; \ problems_size = problems.size(); \ FOR( i , 0 , problems_size ){ \ CERR( i << ": " << problems[i] ); \ } \ cin >> num; \ CERR( "" ); \ if( num < 0 || num >= problems_size ){ \ CERR( "返答は" << problems_size - 1 << "以下の非負整数にしてください。" ); \ CERR( "終了します。" ); \ CERR( "" ); \ return; \ } \ num_temp = 0; \ #define ASK_YES_NO( QUESTION ) \ CERR( "" ); \ CERR( QUESTION << "[y/n]" ); \ cin >> reply; \ if( reply != "y" && reply != "n" ){ \ CERR( "y/nのいずれかで答えてください。" ); \ CERR( "終了します。" ); \ CERR( "" ); \ return; \ } \ CERR( "" ); \ #define AC( BRANCH ) void BRANCH ## AutoCheck( int& num , int& num_temp , string& reply , vector& problems , int& problems_size ) #define CALL_AC( BRANCH ) BRANCH ## AutoCheck( num , num_temp , reply , problems , problems_size ) AC( DebugHint ); AC( LibrarySearch ); void AutoCheck( bool& auto_checked ) { int num = 0; int num_temp = 0; string reply{}; vector problems{}; int problems_size = 0; ASK_NUMBER( "提出用ファイルを実行する。" , "提出済みファイルをデバッグする。" , "ライブラリーを探索する。" ); if( num == num_temp++ ){ auto_checked = false; CERR( "提出用ファイルを実行します。" ); } else if( num == num_temp++ ){ CALL_AC( DebugHint ); } else if( num == num_temp++ ){ CALL_AC( LibrarySearch ); } CERR( "" ); } AC( DebugHint ) { ASK_NUMBER( "WAのデバッグ" , "REのデバッグ" , "TLEのデバッグ" , "MLEのデバッグ" ); if( num == num_temp++ ){ CERR( "- operator<<()をint型に適用していませんか?" ); CERR( "- イベント管理にsetを使っていませんか?" ); CERR( " - keyに狭義全順序を用いる場合、" ); CERR( " - 同一keyの処理を同時に行うならばmapを検討しましょう。" ); CERR( " - 同一keyの処理を同時に行わないならばmultisetを検討しましょう。" ); CERR( " - keyに狭義全順序の積順序を用いる場合、mapの入れ子を検討しましょう。" ); CERR( "- 動的計画法においてindexのswapやmodを用いてメモリ削減を行う場合、" ); CERR( " 各ループの最初にdpテーブルの初期化をし忘れていませんか?" ); CERR( "- マルチテストケースで配列にstaticをつけて値が持ち越されていませんか?" ); CERR( "- cLay使用畤にsetmod以前に10^9+7より大きい値を代入していませんか?" ); CERR( "" ); CERR( "原因に心当たりがない場合はランダムテストを検討しましょう。" ); } else if( num == num_temp++ ){ CERR( "- 制約を間違えてassertしていませんか?" ); CERR( "- 配列外参照していませんか?" ); CERR( " - 座標圧縮畤に圧縮前の添字操作と圧縮後の添字操作を混同していませんか?" ); CERR( " - グラフや座標圧縮の操作で存在しない座標-1の場合分けを忘れていませんか?" ); CERR( "- 二分探索マクロで探索値が存在しない場合を回避し忘れていませんか?" ); CERR( "- イテレータを範囲外まで--/++していませんか?" ); CERR( "- 0で割っていませんか?" ); } else if( num == num_temp++ ){ CERR( "- マルチテストケースでデータ構造の初期化が反復されていませんか?" ); CERR( "- whileループで添字等の更新忘れによる無限ループはありませんか?" ); CERR( "- 重過ぎる定数倍を考慮し忘れてませんか?" ); CERR( " - bool値の処理はbit演算による並列化を検討しましょう。" ); CERR( " - 変数の和や積に上限があるグリッド問題は動的配列を検討しましょう。" ); CERR( "- リアクティブ問題でflushと改行をし忘れていませんか?" ); } else if( num == num_temp++ ){ CERR( "- グラフで辺を持ち過ぎていませんか?" ); CERR( "- (多重)配列でサイズが10^6を大幅に上回ってませんか?" ); CERR( " - データ構造は座標圧縮を検討しましょう。" ); CERR( " - 動的計画法は初期化忘れに注意して添字のswapやmodを検討しましょう。" ); CERR( "- Python使用畤に変数名の重複による型エラーが起きてませんか?" ); CERR( " - rangeのRと[L,R)のRの重複は後者に[l,r)を用いて解消しましょう。" ); CERR( " - 素数のPと順列のPの重複は後者にAやQを用いて解消しましょう。" ); CERR( "- Python使用畤に桁数が大きくなり過ぎていませんか?" ); } CERR( "" ); CERR( "提出済みファイルのデバッグは以上です。終了します。" ); } AC( ExplicitExpression ); AC( ExplicitExpressionUnary ); AC( ExplicitExpressionUnaryLinearRecursion ); AC( ExplicitExpressionUnaryIteration ); AC( ExplicitExpressionUnaryOther ); AC( ExplicitExpressionArraySum ); AC( ExplicitExpressionFunctionOnPermutation ); AC( ExplicitExpressionFunctionOnTree ); AC( FunctionOnTree ); AC( ExplicitExpressionFunctionOnNonTreeGraph ); AC( ExplicitExpressionOrder ); AC( ExplicitExpressionProbability ); AC( ExplicitExpressionOther ); AC( Maximisation ); AC( MaximisationFunctionOnAffineSpace ); AC( MaximisationFunctionOnArray ); AC( MaximisationSubArraySum ); AC( MaximisationArrayFunction ); AC( MaximisationArrayLength ); AC( MaximisationFunctionOnTree ); AC( MinimisationMovingCost ); AC( MaximisationStringMatching ); AC( MaximisationBipartiteMatching ); AC( MaximisationProbability ); AC( Counting ); AC( CountingExplicitExpression ); AC( CountingArray ); AC( CountingSubArray ); AC( CountingSumFixedSubArray ); AC( CountingRestrctedSubArray ); AC( CountingRestrctedContinuousSubArray ); AC( CountingRestrctedDiscontinuousSubArray ); AC( CountingRestrctedSubPermutation ); AC( CountingArbitraryArray ); AC( CountingSubArrayImage ); AC( CountingPartitionOfTree ); AC( CountingString ); AC( CountingMatchingSubString ); AC( Solving ); AC( Query ); AC( QueryArray ); AC( QueryGraph ); AC( Decision ); AC( DecisionConnectedness ); AC( DecisionHigherConnectedness ); AC( DecisionGame ); AC( DecisionAccessibility ); AC( DecisionSatisfiability ); AC( Construction ); AC( LibrarySearch ) { ASK_NUMBER( "明示式の計算問題" , "最大/最小化問題" , "数え上げ問題" , "求解問題" , "クエリ処理問題" , "真偽判定問題" , "構築問題" ); if( num == num_temp++ ){ CALL_AC( ExplicitExpression ); } else if( num == num_temp++ ){ CALL_AC( Maximisation ); } else if( num == num_temp++ ){ CALL_AC( Counting ); } else if( num == num_temp++ ){ CALL_AC( Solving ); } else if( num == num_temp++ ){ CALL_AC( Query ); } else if( num == num_temp++ ){ CALL_AC( Decision ); } else if( num == num_temp++ ){ CALL_AC( Construction ); } ASK_YES_NO( "マルチテストケースですか?" ); if( reply == "y" ){ CERR( "テストケースを跨ぐ前計算が可能か否かを優先的に考察しましょう。" ); CERR( "" ); CERR( "テストケース全体でのNの総和に直接上限が与えられている問題では、" ); CERR( "ライブラリーの使用時は配列の初期化が各テストケースに必要となる場合に" ); CERR( "TLEとなる可能性が高いです。" ); CERR( "- 動的配列への置き換え" ); CERR( "- 座標圧縮" ); CERR( " \\Mathematics\\SetTheory\\DirectProduct\\CoordinateCompress" ); CERR( "を検討しましょう。" ); CERR( "" ); } CERR( "ライブラリー探索は以上です。終了します。" ); } AC( ExplicitExpression ) { ASK_NUMBER( "1変数関数/数列の計算問題" , "配列上の関数の総和の計算問題" , "順列上の関数の計算問題" , "木上の関数の総和の計算問題" , "木以外のグラフ上の関数の計算問題" , "序数の計算問題" , "確率/期待値の計算問題" , "その他の明示式の計算問題" ); if( num == num_temp++ ){ CALL_AC( ExplicitExpressionUnary ); } else if( num == num_temp++ ){ CALL_AC( ExplicitExpressionArraySum ); } else if( num == num_temp++ ){ CALL_AC( ExplicitExpressionFunctionOnPermutation ); } else if( num == num_temp++ ){ CALL_AC( ExplicitExpressionFunctionOnTree ); } else if( num == num_temp++ ){ CALL_AC( ExplicitExpressionFunctionOnNonTreeGraph ); } else if( num == num_temp++ ){ CALL_AC( ExplicitExpressionOrder ); } else if( num == num_temp++ ){ CALL_AC( ExplicitExpressionProbability ); } else if( num == num_temp++ ){ CALL_AC( ExplicitExpressionOther ); } } AC( ExplicitExpressionUnary ) { ASK_NUMBER( "線形漸化式の問題" , "1変数関数の反復合成の問題" , "その他の1変数関数の計算問題" ); if( num == num_temp++ ){ CALL_AC( ExplicitExpressionUnaryLinearRecursion ); } else if( num == num_temp++ ){ CALL_AC( ExplicitExpressionUnaryIteration ); } else if( num == num_temp++ ){ CALL_AC( ExplicitExpressionUnaryOther ); } CERR( "を検討しましょう。" ); } AC( ExplicitExpressionUnaryLinearRecursion ) { CERR( "L項間線形漸化式a_n=sum(i,0,L){c_i a_{n-L+i}}が与えられているとします。" ); CERR( "" ); CERR( "適宜Z/B Zの型" ); CERR( "\\Mathematics\\Arithmetic\\Mod\\ConstexprModulo" ); CERR( "\\Mathematics\\Arithmetic\\Mod\\QuotientRing" ); CERR( "を用いて" ); CERR( "- O(LN)が間に合いそうならば直近L+1項をメモ化する動的計画法" ); CERR( "- O(L^2 log N)が間に合いそうならば行列累乗" ); CERR( " \\Mathematics\\LinearAlgebra" ); CERR( "- 係数列(c_i)_iが定数列でO(N log N)が間に合いそうならば区間加算BIT" ); CERR( " \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\BIT\\IntervalAddo" ); CERR( "- 係数列(c_i)_iが単調でD階差分が定数列でO((N+L)D)が間に合いそうならば" ); CERR( " 「係数列のd階Δ^d(c_i)_i差分と(a_i)_{i=n-L+d+1}^{n}の内積dp[n][d]」" ); CERR( " を管理する動的計画法" ); } AC( ExplicitExpressionUnaryIteration ) { 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)すら通らなさそうならば関数の規則性を見付けるための実験" ); } AC( ExplicitExpressionUnaryOther ) { CERR( "まずは小さい入力の場合を愚直に計算し、OEISで検索しましょう。" ); CERR( "https://oeis.org/?language=japanese" ); CERR( "" ); CERR( "次に出力の定義と等価な式を考察しましょう。" ); CERR( "- 単調ならば、冪乗や階乗" ); CERR( "- 定義にp進法が使われていれば、各種探索アルゴリズム" ); CERR( "- 入力が素数に近い場合に規則性があれば、p進付値、p進法、" ); CERR( " オイラー関数、約数の個数など" ); CERR( "を検討しましょう。" ); CERR( "" ); CERR( "前計算の候補としては" ); CERR( "- 素数列挙" ); CERR( "- 1つまたは複数の整数の約数列挙" ); CERR( "- オイラー関数の値の列挙" ); CERR( "- サブゴールとなる関係式を満たす解の列挙" ); } AC( ExplicitExpressionArraySum ) { ASK_NUMBER( "成分を受け取る関数の総和の計算問題" , "部分列を受け取る関数の総和の計算問題" ); if( num == num_temp++ ){ CERR( "成分を受け取る関数fが与えられているとします。" ); CERR( "fが一次式の場合、実質内積と定数の和となります。" ); CERR( "内積は片方の添え字を反転させることで畳み込みに帰着させることができます。" ); CERR( "配列への操作がシフトである場合は繰り返し内積を求めることになるので、" ); CERR( "適当な法での高速フーリエ変換" ); CERR( "\\Mathematics\\Arithmetic\\Mod" ); CERR( "\\Mathematics\\Polynoial" ); CERR( "を検討しましょう。" ); } else if( num == num_temp++ ){ ASK_NUMBER( "連続部分列への分割に関する関数の総和の計算問題" , "連続とは限らない部分列への分割に関する関数の総和の計算問題" ); if( num == num_temp++ ){ CERR( "配列の添字集合は全順序集合なので、木の分割の問題に一般化されます。" ); CALL_AC( ExplicitExpressionFunctionOnTree ); CERR( "" ); CERR( "更にfが部分列の長さに関する再帰的な構造を持つ場合、全ての連続部分列に" ); CERR( "対しfの値を前計算することを検討しましょう。" ); } else if( num == num_temp++ ){ CERR( "配列の並び換えによって答えが変わらないので、適切にソートしてから" ); CERR( "計算することを検討しましょう。" ); } } CERR( "" ); CERR( "入力が大きい場合と小さい場合で解法を変える考察を忘れないようにしましょう。" ); } AC( ExplicitExpressionFunctionOnPermutation ) { CERR( "- 符号そのものの計算問題は" ); CERR( " - O(N log_2 N)やO(N^2)が間に合いそうなら転倒数の計算" ); CERR( " - O(N log_2 N)が間に合わなさそうなら互換表示(O(N))" ); CERR( "- 符号と何かの積の和は行列式に帰着させ、" ); CERR( " - 行列式そのものなら行基本変形(O(N^3))" ); CERR( " - 余因子展開の途中の値が必要ならメモ化再帰(O(N 2^N))" ); CERR( "を検討しましょう。" ); CERR( "" ); CERR( "1つの順列の転倒数は、" ); 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=3かつ(n,L)=(2,2)で" ); CERR( "- Q(i)=「i=3」" ); CERR( "- R_0(b_0,b_1)=「b_0b_1」" ); CERR( "と表されます。" ); CERR( "" ); if( num == num_temp++ ){ CALL_AC( CountingRestrctedContinuousSubArray ); } else if( num == num_temp++ ){ CALL_AC( CountingRestrctedDiscontinuousSubArray ); } else if( num == num_temp++ ){ CALL_AC( CountingRestrctedSubPermutation ); } CERR( "を検討しましょう。" ); CERR( "" ); CERR( "特にR_l(B)たちがgcdやmaxなどの羃等演算に関する等式で与えられる場合は、" ); CERR( "不等式の方が扱いやすいのでゼータ変換/メビウス変換" ); CERR( "\\Mathematics\\Combinatorial\\ZetaTransform" ); CERR( "を検討しましょう。" ); } AC( CountingRestrctedContinuousSubArray ) { CERR( "P(B)を満たすAの連続部分列Bの数え上げは、" ); CERR( "- R_lたちがlに依存しないならば尺取り法O(N)" ); CERR( "- R_lたちがlに依存する場合、" ); CERR( " - O(N^2)が通りそうなら左端を固定した愚直探索" ); CERR( " - O(N^2)が通らなさそうならR_lたちの読み替え" ); CERR( "を検討しましょう。" ); } AC( CountingRestrctedDiscontinuousSubArray ) { CERR( "P(B)を満たすAの連続とは限らない部分列Bの数え上げは、" ); CERR( "- n-1<=i<=max{j<=N|Q(j)}を満たす各i" ); CERR( "- (0,1,...,N-1)の長さn-1の各部分列s" ); CERR( "に対する" ); CERR( "「長さiで、任意の0<=l<=i-nに対しR_l(B)を満たし、" ); CERR( " 末尾n-1項がsに対応するAの部分列Bの個数dp[i][s]」" ); CERR( "を管理するi,sに関する動的計画法" ); CERR( "\\Mathematics\\Combinatorial\\Counting\\IncreasingSubsequence" ); CERR( "\\Mathematics\\Combinatorial\\Counting\\IncreasingSubsequence\\Subwalk" ); CERR( "を検討しましょう。" ); } AC( CountingRestrctedSubPermutation ) { CERR( "P(B)を満たすAの部分順列Bの数え上げは、" ); CERR( "- n-1<=|S|<=max{j<=N|Q(j)}を満たす(0,1,...,N-1)の部分集合S" ); CERR( "- Sの長さn-1の各部分順列s" ); CERR( "に対する" ); CERR( "「任意の0<=l<=|S|-nに対しR_l(B)を満たし、末尾n-1項がsに対応し、」" ); CERR( " 全体がSに対応するAの部分順列Bの個数dp[S][s]」" ); CERR( "を管理するS,sに関する動的計画法を検討しましょう。" ); CERR( "" ); CERR( "sの網羅は[0,N)^{n-1}の全探策でもS内の順列探索と定数倍しか変わらないので" ); CERR( "実装の速さを優先しましょう。" ); CERR( "" ); CERR( "Nが小さい場合は概算値" ); CERR( "7! ≒ 5×10^3" ); CERR( "8! ≒ 4×10^4" ); CERR( "9! ≒ 4×10^5" ); CERR( "10! ≒ 4×10^6" ); CERR( "11! ≒ 4×10^7" ); CERR( "12! ≒ 5×10^8" ); CERR( "を参考に順列の全列挙" ); CERR( "\\Mathematics\\Combinatorial\\Permutation" ); CERR( "を検討しましょう。" ); } AC( CountingArbitraryArray ) { ASK_NUMBER( "配列を受け取る関数の値の数え上げ問題" , "隣接成分間関係式を満たす配列の数え上げ問題" ); if( num == num_temp++ ){ CERR( "- 配列の種類が少ない場合は、全ての配列に対する関数の値の前計算" ); CERR( "- 取り得る値が少なく関数が長さに関して再帰的構造を持つ場合は、" ); CERR( " 「長さiの時に値vである配列の総数dp[i][v]」" ); CERR( " を管理するi,vに関する動的計画法" ); } else if( num == num_temp++ ){ CERR( "- いくつかの条件の重ね合わせの時は包除原理" ); CERR( "- 全順序の場合は数の分割方法などへの翻訳" ); CERR( "- 疎な半順序の場合はグラフの前計算" ); } CERR( "を検討しましょう。" ); } AC( CountingSubArrayImage ) { CERR( "入力で与えられる配列をAと置きます。" ); CERR( "配列として等しいAの部分列のうち辞書式順序最小のものを数え上げる" ); CERR( "部分列DPを検討しましょう。具体的には" ); CERR( "「末尾が第i成分由来で辞書順最小なAの部分列の個数dp[i]」" ); CERR( "「s=A[j]を満たすj ( &e )[bound_HW] , const char& walkable = '.' ){FOR(j,0,W){if(Si[j]==walkable){int v = EnumHW_inv(i,j);if(i>0){e[EnumHW_inv(i-1,j)].push_back(v);}if(i+10){e[EnumHW_inv(i,j-1)].push_back(v);}if(j+1; // list e[bound_M]; // bound_Mのデフォルト値は10^5 // // list e[bound_HW]; // bound_HWのデフォルト値は10^6 // list E( const int& i ) // { // list answer = e[i]; // // 入力によらない処理 // return answer; // } int main() { UNTIE; AUTO_CHECK; TEST_CASE_NUM( 1 ); START_MAIN; CIN( ll , N ); // CIN( ll , M ); // CIN( ll , K ); // // CIN_ASSERT( N , 1 , bound_N ); // 基本不要、上限のデフォルト値は10^5 // // CIN_ASSERT( M , 1 , bound_M ); // 基本不要、上限のデフォルト値は10^5 // CIN( string , S ); // CIN( string , T ); // cin >> H >> W; // // SET_ASSERT( H , 1 , bound_H ); // 基本不要、上限のデフォルト値は10^3 // // SET_ASSERT( W , 1 , bound_W ); // 基本不要、上限のデフォルト値は10^3 // H_minus = H - 1; // W_minus = W - 1; // HW = H * W; // // assert( HW <= bound_HW ); // 基本不要、上限のデフォルト値は10^6 // string S[H]; // // bool non_wall[bound_H+1][bound_W+1]={}; // FOR( i , 0 , H ){ // cin >> S[i]; // // SetEdgeOnGrid( S[i] , i , e ); // eの宣言のコメントアウト必要 // // SetWallOnGrid( S[i] , i , non_wall ); // } // // (i,j)->(k,h)の方向番号を取得: DirectionNumberOnGrid( i , j , k , h ); // // v->wの方向番号を取得: DirectionNumberOnGrid( v , w ); // // 方向番号の反転U<->D、R<->L: ReverseDirectionNumberOnGrid( n ); ll A[N+1]; // ll B[N]; // ll A[bound_N]; // 基本不要、長さのデフォルト値は10^5 // ll B[bound_N]; // 基本不要、長さのデフォルト値は10^5 FOR( i , 0 , N ){ cin >> A[i]; // cin >> B[i]; } if( N == 0 ){ RETURN( 0 ); } A[N] = A[N-1]; // FOR( i , 0 , M ){ // CIN_ASSERT( ui , 1 , N ); // CIN_ASSERT( vi , 1 , N ); // ui--; // vi--; // e[ui].push_back( vi ); // e[vi].push_back( ui ); // } // CIN( int , Q ); // tuple query[Q]; // DEXPR( int , bound_Q , 100000 , 100 ); // 基本不要 // CIN_ASSERT( Q , 1 , bound_Q ); // 基本不要 // FOR( q , 0 , Q ){ // CIN( int , type ); // if( type == 1 ){ // CIN( int , x ); // CIN( int , y ); // // query[q] = { type , x , y }; // } else if( type == 2 ){ // CIN( int , x ); // CIN( int , y ); // // query[q] = { type , x , y }; // } else { // CIN( int , x ); // CIN( int , y ); // // query[q] = { type , x , y }; // } // } // // ll guchoku = Guchoku(); ll answer = 0; set S{}; DEXPR( int , bound_Ai , 1000000 , 100 ); int A_inv[bound_Ai+1]; int L = 0; FOREQ( i , 0 , N ){ if( S.count( A[i] ) == 1 ){ int temp = i - L; if( answer < temp ){ answer = temp; } FOR( j , L , A_inv[A[i]] ){ S.erase( A[j] ); } L = A_inv[A[i]] + 1; } else { S.insert( A[i] ); } A_inv[A[i]] = i; } // COUT( ( answer ) ); RETURN( answer ); FINISH_MAIN; }