#884581 (C++17(gcc12)) No.2359 A in S ?

提出ソース
結果

問題	No.2359 A in S ?
コンテスト
ユーザー	👑 p-adic
提出日時	2023-06-24 16:29:08
言語	C++17(gcc12) (gcc 12.3.0 + boost 1.87.0)
結果	AC
実行時間	564 ms / 2,000 ms
コード長	47,491 bytes
コンパイル時間	13,270 ms
コンパイル使用メモリ	296,124 KB
最終ジャッジ日時	2025-02-15 02:04:59
ジャッジサーバーID （参考情報）	judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン	結果
sample	AC * 2
other	AC * 18
権限があれば一括ダウンロードができます
ソースコード

raw source code
#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 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 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 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 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 FOR_ITR( ARRAY , ITR , END ) for( auto ITR = ARRAY .begin() , END = ARRAY .end() ; ITR != END ; ITR ++ )
#define REPEAT( HOW_MANY_TIMES ) FOR( VARIABLE_FOR_REPEAT , 0 , HOW_MANY_TIMES )
#define QUIT return 0
#define COUT( ANSWER ) cout << ( ANSWER ) << "\n"
#define RETURN( ANSWER ) COUT( ANSWER ); QUIT
#define SET_PRECISION( PRECISION ) cout << fixed << setprecision( PRECISION )
#define DOUBLE( PRECISION , ANSWER ) SET_PRECISION << ( ANSWER ) << "\n"; 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&

int QuitLibrarySearch( const int& problems_size ){
  cerr << "返答は" << problems_size - 1 << "以下の非負整数にしてください。";
  CERR( "終了します。" );
  CERR( "" );
  return -1;
}

int LibrarySearch( int num = -1 )
{
  vector<string> problems =
    {
      "数に関する問題。" ,
      "配列に関する問題。" ,
      "文字列に関する問題。" ,
      "順列に関する問題。" ,
      "矩形領域に関する問題。" ,
      "グラフに関する問題。" ,
      "部分和問題。" ,
      "確率／期待値に関する問題。" ,
      "ゲームに関する問題。" ,
      "論理に関する問題。" ,
      "順序に関する問題。" ,
      "関数適用に関する問題。" ,
      "構築問題。"
    };
  CEXPR( int , num_graph , 5 );
  CEXPR( int , num_subsequence_sum , 6 );
  CEXPR( int , num_game , 8 );
  int problems_size = problems.size();
  string reply{};
  if( num == -1 ){
    CERR( "ライブラリーを探索しますか？[y/n]" );
    CIN( string , reply );
    if( reply == "n" ){
      CERR( "ライブラリーを探索せずに続行します。" );
      CERR( "" );
      return 0;
    } else if( reply != "y" ){
      CERR( "y/nのいずれかで答えてください。" );
      CERR( "終了します。" );
      CERR( "" );
      return -1;
    }
    CERR( "" );
    CERR( "ライブラリーを探索します。" );
    CERR( "問題の種類を番号で指定してください；" );
    FOR( i , 0 , problems_size ){
      CERR( i << ": " << problems[i] );
    }
    cin >> num;
  }
  CERR( "" );
  int num_temp = 0;
  if( num < 0 || num >= problems_size ){
      return QuitLibrarySearch( problems_size );
  } else 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" ){
      CERR( "このケースのライブラリー探索は未実装です。" );
    } else {
      CERR( "y/nのいずれかで答えてください。" );
      CERR( "終了します。" );
      CERR( "" );
      return -1;
    }
    CERR( "" );
    CERR( "マルチテストケースの場合は以下の前計算を検討しましょう；" );
    CERR( "素数列挙、約数列挙、サブゴールとなる関係式を満たす解列挙。" );
  } else if( num == num_temp++ ){
    CERR( "より詳細な問題の種類を番号で指定してください；" );
    problems =
      {
	"区間処理問題。" ,
	"最大化問題。" ,
	"最長部分列問題。" ,
	"数え上げ問題。" ,
	"部分和問題。"
      };
    problems_size = problems.size();
    FOR( i , 0 , problems_size ){
      CERR( i << ": " << problems[i] );
    }
    CIN( int , num );
    CERR( "" );
    num_temp = 0;
    if( num < 0 || num >= problems_size ){
      return QuitLibrarySearch( problems_size );
    } else if( num == num_temp++ ){
      CERR( "クエリ処理を整理し、代数構造を抽出しましょう。" );
      CERR( "- 可換群構造+しか使わない場合は可換群BIT" );
      CERR( "  \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\BIT\\Template")
      CERR( "- 可換羃等モノイド構造∨しか使わない場合は可換羃等モノイドBIT" );
      CERR( "  \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\BIT\\IntervalMax\\Template" );
      CERR( "- モノイド構造*しか使わない場合はモノイドBITかモノイド平方分割" );
      CERR( "  \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\BIT\\Template\\Monoid" );
      CERR( "  \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\SqrtDecomposition\\Template\\Monoid" );
      CERR( "- 集合へのマグマ作用(*,\\cdot)しか使わない場合は双対平方分割" );
      CERR( "  \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\SqrtDecomposition\\Template\\Dual" );
      CERR( "- モノイドへのマグマ作用(+,\\cdot)しか使わない場合は遅延評価平方分割" );
      CERR( "  \\Mathematics\\SetTheory\\DirectProduct\\AffineSpace\\SqrtDecomposition\\Template\\LazyEvaluation" );
      CERR( "- 非結合演算が使われている時は写像のコード化" );
      CERR( "  \\Mathematics\\Function\\Encoder" );
      CERR( "- (+,max)が使われている時は場合分けで表し、配列の各成分の場合分けが切り替わる位置（イベント）と" );
      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)で処理できます。" );
      CERR( "を検討しましょう。" );
    } 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)が通りそうで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 {
      return LibrarySearch( num = num_subsequence_sum );
    }
  } else if( num == num_temp++ ){
    CERR( "より詳細な問題の種類を番号で指定しましょう；" );
    problems =
      {
	"部分列マッチング問題。" ,
	"回文探索問題。" ,
      };
    problems_size = problems.size();
    FOR( i , 0 , problems_size ){
      CERR( i << ": " << problems[i] );
    }
    CIN( int , num );
    CERR( "" );
    num_temp = 0;
    if( num < 0 || num >= problems_size ){
      return QuitLibrarySearch( problems_size );
    } else if( num == num_temp++ ){
      CERR( "基本的には丁寧にループを回して解きましょう。" );
      CERR( "- 比較対象が少ない場合、前または後ろから順に探索（貪欲法／動的計画法）" );
      CERR( "- ワイルドカードを含む場合、" );
      CERR( "  - 前または後ろから順に場合分けをしてO(N)で処理できるか" );
      CERR( "  - 可能な代入方法を絞り込んでO(N)種類に落せるか" );
      CERR( "- 比較回数が多い場合、ローリングハッシュ" );
      CERR( "- マッチングする文字列の最長化をする場合、Zアルゴリズム" );
      CERR( "  https://qiita.com/Pro_ktmr/items/16904c9570aa0953bf05" );
      CERR( "を検討しましょう。" );
    } else {
      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++ ){
    CERR( "より詳細な問題の種類を番号で指定してください；" );
    problems =
      {
	"探索問題。" ,
	"絶対値の最小／最大化問題。" ,
	"その他の最小／最大化問題。"
      };
    problems_size = problems.size();
    FOR( i , 0 , problems_size ){
      CERR( i << ": " << problems[i] );
    }
    CIN( int , num );
    CERR( "" );
    num_temp = 0;
    if( num < 0 || num >= problems_size ){
      return QuitLibrarySearch( problems_size );
    } else 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++ ){
    CERR( "より詳細な問題の種類を番号で指定してください；" );
    problems =
      {
	"２点最短径路（迷路）問題。" ,
	"多点最短経路（スタンプラリー）問題。" ,
	"木の問題。" ,
	"連結性問題。"
      };
    problems_size = problems.size();
    FOR( i , 0 , problems_size ){
      CERR( i << ": " << problems[i] );
    }
    CIN( int , num );
    CERR( "" );
    num_temp = 0;
    if( num < 0 || num >= problems_size ){
      return QuitLibrarySearch( problems_size );
    } else 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 {
      CERR( "- 0次の連結性はUnionFind" );
      CERR( "  \\Utility\\VLTree\\UnionFindForest" );
      CERR( "- 高次の連結性は深さ優先探索" );
      CERR( "  \\Utility\\Search\\DepthFirst" );
      CERR( "- マルチテストケースの場合は座標圧縮との併用" );
      CERR( "を検討しましょう。" );
    }
  } else if( num == num_temp++ ){
    CERR( "より詳細な問題の種類を番号で指定してください；" );
    problems =
      {
	"部分和の最大化問題。" ,
	"部分和の数え上げ問題。"
      };
    problems_size = problems.size();
    FOR( i , 0 , problems_size ){
      CERR( i << ": " << problems[i] );
    }
    CIN( int , num );
    CERR( "" );
    num_temp = 0;
    if( num < 0 || num >= problems_size ){
      return QuitLibrarySearch( problems_size );
    } else 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 {
      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++ ){
    CERR( "より詳細な問題の種類を番号で指定してください；" );
    problems =
      {
	"半順序集合上の関数の計算問題。" ,
	"半順序集合上の降下／上昇列に関する問題。"
      };
    problems_size = problems.size();
    FOR( i , 0 , problems_size ){
      CERR( i << ": " << problems[i] );
    }
    CIN( int , num );
    CERR( "" );
    num_temp = 0;
    if( num < 0 || num >= problems_size ){
      return QuitLibrarySearch( problems_size );
    } else if( num == num_temp++ ){
      CERR( "ゼータ変換を検討しましょう。" );
    } else if( num == num_temp++ ){
      CERR( "順序関係の定める有向グラフの問題に帰着させましょう。" );
      return LibrarySearch( num = num_graph );
    }
  } else if( num == num_temp++ ){
    CERR( "より詳細な問題の種類を番号で指定してください；" );
    problems =
      {
	"反復合成の計算問題。" ,
	"反復合成による到達可能性問題。"
      };
    problems_size = problems.size();
    FOR( i , 0 , problems_size ){
      CERR( i << ": " << problems[i] );
    }
    CIN( int , num );
    CERR( "" );
    num_temp = 0;
    if( num < 0 || num >= problems_size ){
      return QuitLibrarySearch( problems_size );
    } else 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 {
      CERR( "関数による遷移が定める有向グラフの問題に帰着させましょう。" );
      return LibrarySearch( num = num_graph );
    }
  } else if( num == num_temp++ ){
    CERR( "より詳細な問題の種類を番号で指定してください；" );
    problems =
      {
	"数や配列や文字列の構築。" ,
	"経路の構築。" ,
	"戦略の構築。"
      };
    problems_size = problems.size();
    FOR( i , 0 , problems_size ){
      CERR( i << ": " << problems[i] );
    }
    CIN( int , num );
    CERR( "" );
    num_temp = 0;
    if( num < 0 || num >= problems_size ){
      return QuitLibrarySearch( problems_size );
    } else 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;
}

template <typename V , typename T , typename U , int size_max>
class DifferenceSequenceBody
{

protected:
  int m_size;
  U m_a[size_max];
  U m_lazy_addition[size_max];
  bool m_updated;

public:
  inline DifferenceSequenceBody( const int& size );

  // 作用の遅延評価を解消してから値を代入する。
  inline void Set( const V& v , const U& u );
  // 作用の遅延評価を解消してから値を参照する。
  inline const U& Get( const V& v );
  // 作用の遅延評価を解消せずに元々の値を参照する。
  inline U& Ref( const V& v );

  // vでの値にnを遅延評価せずに加算する。
  inline void Add( const V& v , const U& u );

  // 要件
  // (1) v_startからv_finalの各要素への有向パスが存在する。
  // (2) v_finalの要素間の有向パスが存在しない。
  // を満たす場合にのみサポート。
  // 条件
  // (1) v_startからvへの有向パスが存在する。
  // (2) v_finalのいずれの要素からもvへの有向パスが存在しない。
  // を満たすVの各要素vに対し、vでの値にnを遅延評価で加算する。
  void SubTreeAdd( const V& v_start , const list<V>& v_final , const U& u );

  inline DifferenceSequenceBody<V,T,U,size_max>& operator+=( const DifferenceSequenceBody<V,T,U,size_max>& a );

  // 作用の遅延評価を解消する。
  void Update();

private:
  virtual list<V> Edge( const V& v ) const = 0;
  virtual V Edge_inv( const V& v ) const = 0;
  virtual V e( const int& i ) const = 0;
  virtual int e_inv( const V& v ) const = 0;
  virtual const T& Unit() const = 0;
  virtual T& Addition( const T& t0 , T& t1 ) const = 0;
  virtual U Inverse( const U& u ) const = 0;
  virtual T& Action( const U& u , T& t ) const = 0;

};

// 入力の範囲内で要件
// (6) (T,operator+:T^2->T,T(),operator-:T->T)は可換群である。
// が成り立つ場合にのみサポート。

// 0による初期化O(size_max)
// 配列による初期化O(size_max)

// 一点代入O(size)（作用の遅延評価を解消する。元々作用の遅延評価がない場合はO(1)）
// 一点取得O(size)（作用の遅延評価を解消する。元々作用の遅延評価がない場合はO(1)）

// 一点加算O(1)（作用を遅延評価しない）
// 始切片加算O(1)（作用を遅延評価する）
// 終切片加算O(1)（作用を遅延評価する）
// 区間加算O(1)（作用を遅延評価する）
// 加法O(size)（作用の遅延評価を解消する）
template <typename T , int size_max>
class DifferenceSequence :
  public DifferenceSequenceBody<int,T,T,size_max>
{

public:
  inline DifferenceSequence( const int& size );
  inline DifferenceSequence( const int& size , const T ( &a )[size_max] );
  inline DifferenceSequence( const int& size , T ( &&a )[size_max] );

  inline void InitialSegmentAdd( const int& v_start , const T& u );
  inline void FinalSegmentAdd( const int& v_final , const T& u );
  inline void IntervalAdd( const int& v_start , const int& v_final , const T& u );
  inline void IntervalSubtract( const int& v_start , const int& v_final , const T& u );

private:
  inline list<int> Edge( const int& v ) const;
  inline int Edge_inv( const int& v ) const;
  inline int e( const int& i ) const;
  inline int e_inv( const int& v ) const;
  inline const T& Unit() const;
  inline T& Addition( const T& t0 , T& t1 ) const;
  inline T Inverse( const T& u ) const;
  inline T& Action( const T& u , T& t ) const;

};

// 入力の範囲内で要件
// (1) グラフ(V,E)は有向グラフとして各連結成分が根つき木である。
// (2) (V,E)の各辺v0->v1に対しv0 = E_inv(v1)が成り立つ。
// (3) (V,E)の各連結成分の根vに対しE_inv(v) = vが成り立つ。
// (4) (enum_V:N->V,enum_V_inv:V->N)は互いに逆写像である。
// (5) (V,E)の各辺v0->v1に対しe_inv(v0) < e_inv(v1)が成り立つ。
// (6) (T,m_T:T^2->T,e_T:1->T)は可換モノイドである。
// (7) (U,inv_U:U->U)は羃等演算つき集合である。
// (8) o_U:U×T->Tは(U,inv_U)から(T,m_T,e_T)の可逆元全体の群への準同型とm_Tの合成である。
// が成り立つ場合にのみサポート。

// e_T()による初期化O(size_max)

// 一点代入O(size)（作用の遅延評価を解消する。元々作用の遅延評価がない場合はO(1)）
// 一点取得O(size)（作用の遅延評価を解消する。元々作用の遅延評価がない場合はO(1)）

// 一点加算O(1)（作用を遅延評価しない）
// 部分木加算O(1+t_finalのサイズ)（作用を遅延評価する）
// 加法O(size)（作用の遅延評価を解消する）
template <typename V , list<V> E(const V&) , V E_inv(const V&) , V enum_V(const int&) , int enum_V_inv(const V&) , typename T , T& m_T(const T&,T&) , const T& e_T() , typename U , U inv_U(const U&) , T& o_U(const U&,T&) , int size_max>
class EnumerationDifferenceSequence :
  public DifferenceSequenceBody<V,T,U,size_max>
{

public:
  inline EnumerationDifferenceSequence( const int& size );
  
private:
  inline list<V> Edge( const V& v ) const;
  inline V Edge_inv( const V& v ) const;
  inline V e( const int& i ) const;
  inline int e_inv( const V& v ) const;
  inline const T& Unit() const;
  inline T& Addition( const T& t0 , T& t1 ) const;
  inline U Inverse( const U& u ) const;
  inline T& Action( const U& u , T& t ) const;

};

template <typename T , int size_max>
class MultiDimensionalDifferenceSequence :
  public DifferenceSequence<T,size_max>
{

public:
  inline MultiDimensionalDifferenceSequence( const int& size );

  template <typename... Args> inline void IntervalAdd( const int& v_start , const int& v_final , const Args&... args );
  template <typename... Args> inline void IntervalSubtract( const int& v_start , const int& v_final , const Args&... args );

private:
  inline T Inverse( const T& u ) const = delete;

};

template <typename V , typename T , typename U , int size_max> inline DifferenceSequenceBody<V,T,U,size_max>::DifferenceSequenceBody( const int& size ) :
  m_size( size ) , m_a() , m_lazy_addition() , m_updated( false ) { assert( m_size <= size_max ); }
template <typename T , int size_max> inline DifferenceSequence<T,size_max>::DifferenceSequence( const int& size ) :
  DifferenceSequenceBody<int,T,T,size_max>( size )
{

  using base = DifferenceSequenceBody<int,T,T,size_max>;
  const T& unit = Unit();
  
  if( base::m_lazy_addition[0] != unit ){

    for( int i = 0 ; i < size ; i++ ){

      base::m_a[i] = base::m_lazy_addition[i] = unit;

    }

  }

}

template <typename T , int size_max> inline DifferenceSequence<T,size_max>::DifferenceSequence( const int& size , const T ( &a )[size_max] ) :
  DifferenceSequenceBody<int,T,T,size_max>( size )
{

  using base = DifferenceSequenceBody<int,T,T,size_max>;

  for( int i = 0 ; i < size ; i++ ){

    base::m_a[i] = a[i];

  }

  const T& unit = Unit();
  
  if( base::m_lazy_addition[0] != unit ){

    for( int i = 0 ; i < size ; i++ ){

      base::m_lazy_addition[i] = unit;

    }

  }
  
}

template <typename T , int size_max> inline DifferenceSequence<T,size_max>::DifferenceSequence( const int& size , T ( &&a )[size_max] ) :
  DifferenceSequenceBody<int,T,T,size_max>( size )
{

  using base = DifferenceSequenceBody<int,T,T,size_max>;
  swap( base::m_a , a );
  const T& unit = Unit();
  
  if( base::m_lazy_addition[0] != unit ){

    for( int i = 0 ; i < size ; i++ ){

      base::m_lazy_addition[i] = unit;

    }

  }

}

template <typename V , list<V> E(const V&) , V E_inv(const V&) , V enum_V(const int&) , int enum_V_inv(const V&) , typename T , T& m_T(const T&,T&) , const T& e_T() , typename U , U inv_U(const U&) , T& o_U(const U&,T&) , int size_max> inline EnumerationDifferenceSequence<V,E,E_inv,enum_V,enum_V_inv,T,m_T,e_T,U,inv_U,o_U,size_max>::EnumerationDifferenceSequence( const int& size ) :
  DifferenceSequenceBody<V,T,U,size_max>( size )
{

  using base = DifferenceSequenceBody<V,T,U,size_max>;
  const T& unit = Unit();
  
  if( base::m_lazy_addition[0] != unit ){

    for( int i = 0 ; i < size ; i++ ){

      base::m_a[i] = base::m_lazy_addition[i] = unit;

    }

  }

}

template <typename T , int size_max> inline MultiDimensionalDifferenceSequence<T,size_max>::MultiDimensionalDifferenceSequence( const int& size ) :
  DifferenceSequence<T,size_max>( size ) {}

template <typename V , typename T , typename U , int size_max> inline void DifferenceSequenceBody<V,T,U,size_max>::Set( const V& v , const U& u ) { Update(); m_a[e_inv(v)] = u; }
template <typename V , typename T , typename U , int size_max> inline const U& DifferenceSequenceBody<V,T,U,size_max>::Get( const V& v ) { Update(); return m_a[e_inv(v)]; }
template <typename V , typename T , typename U , int size_max> inline U& DifferenceSequenceBody<V,T,U,size_max>::Ref( const V& v ) { return m_a[e_inv(v)]; }

template <typename V , typename T , typename U , int size_max> inline void DifferenceSequenceBody<V,T,U,size_max>::Add( const V& v , const U& u ) { Addition( m_a[e_inv(v)] , u ); }

template <typename V , typename T , typename U , int size_max>
void DifferenceSequenceBody<V,T,U,size_max>::SubTreeAdd( const V& v_start , const list<V>& v_final , const U& u )
{

  if( u == Unit() ){

    return;

  }

  m_updated = true;
  Addition( u , m_lazy_addition[e_inv(v_start)] );
  const U u_inv = Inverse( u );
  
  for( auto itr_v_final = v_final.begin() , end_v_final = v_final.end() ; itr_v_final != end_v_final ; itr_v_final++ ){

    auto edge = Edge( *itr_v_final );

    for( auto itr_edge = edge.begin() , end_edge = edge.end() ; itr_edge != end_edge ; itr_edge++ ){

      Addition( u_inv , m_lazy_addition[ e_inv( *itr_edge ) ] );

    }

  }

  return;
  
}

template <typename T , int size_max> inline void DifferenceSequence<T,size_max>::InitialSegmentAdd( const int& v_start , const T& u ) { SubTreeAdd( v_start , {} , u ); }
template <typename T , int size_max> inline void DifferenceSequence<T,size_max>::FinalSegmentAdd( const int& v_final , const T& u ) { IntervalAdd( 0 , v_final , u ); }
template <typename T , int size_max> inline void DifferenceSequence<T,size_max>::IntervalAdd( const int& v_start , const int& v_final , const T& u ) { SubTreeAdd( v_start , { v_final } , u ); }
template <typename T , int size_max> template <typename... Args> inline void MultiDimensionalDifferenceSequence<T,size_max>::IntervalAdd( const int& v_start , const int& v_final , const Args&... args ) { using base = DifferenceSequence<T,size_max>; base::m_a[v_start].IntervalAdd( args... ); int w = v_final + 1; if( w < base::m_size ){ base::m_a[w].IntervalSubtract( args... ); } }

template <typename T , int size_max> inline void DifferenceSequence<T,size_max>::IntervalSubtract( const int& v_start , const int& v_final , const T& u ) { SubTreeAdd( v_start , { v_final } , -u ); }
template <typename T , int size_max> template <typename... Args> inline void MultiDimensionalDifferenceSequence<T,size_max>::IntervalSubtract( const int& v_start , const int& v_final , const Args&... args ) { using base = DifferenceSequence<T,size_max>; base::m_a[v_start].IntervalSubtract( args... ); int w = v_final + 1; if( w < base::m_size ){ base::m_a[w].IntervalAdd( args... ); } }


template <typename V , typename T , typename U , int size_max>
DifferenceSequenceBody<V,T,U,size_max>& DifferenceSequenceBody<V,T,U,size_max>::operator+=( const DifferenceSequenceBody<V,T,U,size_max>& a )
{

  for( int i = 0 ; i < m_size ; i++ ){

    Addition( a.m_a[i] , m_a[i] );
    Addition( a.m_lazy_addition[i] , m_lazy_addition[i] );
    
  }

  Update();
  return *this;
  
}

template <typename V , typename T , typename U , int size_max>
void DifferenceSequenceBody<V,T,U,size_max>::Update()
{

  if( ! m_updated ){

    return;

  }

  for( int i = 0 ; i < m_size ; i++ ){

    const int j = e_inv( Edge_inv( e( i ) ) );
    
    if( j != i ){

      assert( j < i );
      Addition( m_lazy_addition[j] , m_lazy_addition[i] );

    }

  }

  const T& unit = Unit();

  for( int i = 0 ; i < m_size ; i++ ){

    U& m_lazy_addition_i = m_lazy_addition[i];
    Action( m_lazy_addition_i , m_a[i] );
    m_lazy_addition_i = unit;

  }

  m_updated = false;
  return;
  
}


template <typename T , int size_max> inline list<int> DifferenceSequence<T,size_max>::Edge( const int& v ) const { using base = DifferenceSequenceBody<int,T,T,size_max>; const int w = v + 1; return w < base::m_size ? list<int>{ w } : list<int>{}; }
template <typename V , list<V> E(const V&) , V E_inv(const V&) , V enum_V(const int&) , int enum_V_inv(const V&) , typename T , T& m_T(const T&,T&) , const T& e_T() , typename U , U inv_U(const U&) , T& o_U(const U&,T&) , int size_max> inline list<V> EnumerationDifferenceSequence<V,E,E_inv,enum_V,enum_V_inv,T,m_T,e_T,U,inv_U,o_U,size_max>::Edge( const V& v ) const { return E( v ); }

template <typename T , int size_max> inline int DifferenceSequence<T,size_max>::Edge_inv( const int& v ) const { const int w = v - 1; return w < 0 ? v : w; }
template <typename V , list<V> E(const V&) , V E_inv(const V&) , V enum_V(const int&) , int enum_V_inv(const V&) , typename T , T& m_T(const T&,T&) , const T& e_T() , typename U , U inv_U(const U&) , T& o_U(const U&,T&) , int size_max> inline V EnumerationDifferenceSequence<V,E,E_inv,enum_V,enum_V_inv,T,m_T,e_T,U,inv_U,o_U,size_max>::Edge_inv( const V& v ) const { return E_inv( v ); }

template <typename T , int size_max> inline int DifferenceSequence<T,size_max>::e( const int& i ) const { return i; }
template <typename V , list<V> E(const V&) , V E_inv(const V&) , V enum_V(const int&) , int enum_V_inv(const V&) , typename T , T& m_T(const T&,T&) , const T& e_T() , typename U , U inv_U(const U&) , T& o_U(const U&,T&) , int size_max> inline V EnumerationDifferenceSequence<V,E,E_inv,enum_V,enum_V_inv,T,m_T,e_T,U,inv_U,o_U,size_max>::e( const int& i ) const { return enum_V( i ); }

template <typename T , int size_max> inline int DifferenceSequence<T,size_max>::e_inv( const int& v ) const { return v; }
template <typename V , list<V> E(const V&) , V E_inv(const V&) , V enum_V(const int&) , int enum_V_inv(const V&) , typename T , T& m_T(const T&,T&) , const T& e_T() , typename U , U inv_U(const U&) , T& o_U(const U&,T&) , int size_max> inline int EnumerationDifferenceSequence<V,E,E_inv,enum_V,enum_V_inv,T,m_T,e_T,U,inv_U,o_U,size_max>::e_inv( const V& v ) const { return enum_V_inv( v ); }

template <typename T , int size_max> inline const T& DifferenceSequence<T,size_max>::Unit() const { static const T unit{}; return unit; }
template <typename V , list<V> E(const V&) , V E_inv(const V&) , V enum_V(const int&) , int enum_V_inv(const V&) , typename T , T& m_T(const T&,T&) , const T& e_T() , typename U , U inv_U(const U&) , T& o_U(const U&,T&) , int size_max> inline const T& EnumerationDifferenceSequence<V,E,E_inv,enum_V,enum_V_inv,T,m_T,e_T,U,inv_U,o_U,size_max>::Unit() const { return e_T(); }

template <typename T , int size_max> inline T& DifferenceSequence<T,size_max>::Addition( const T& t0 , T& t1 ) const { return t1 += t0; }
template <typename V , list<V> E(const V&) , V E_inv(const V&) , V enum_V(const int&) , int enum_V_inv(const V&) , typename T , T& m_T(const T&,T&) , const T& e_T() , typename U , U inv_U(const U&) , T& o_U(const U&,T&) , int size_max> inline T& EnumerationDifferenceSequence<V,E,E_inv,enum_V,enum_V_inv,T,m_T,e_T,U,inv_U,o_U,size_max>::Addition( const T& t0 , T& t1 ) const { return m_T( t0 , t1 ); }

template <typename T , int size_max> inline T DifferenceSequence<T,size_max>::Inverse( const T& u ) const { return -u; }
template <typename V , list<V> E(const V&) , V E_inv(const V&) , V enum_V(const int&) , int enum_V_inv(const V&) , typename T , T& m_T(const T&,T&) , const T& e_T() , typename U , U inv_U(const U&) , T& o_U(const U&,T&) , int size_max> inline U EnumerationDifferenceSequence<V,E,E_inv,enum_V,enum_V_inv,T,m_T,e_T,U,inv_U,o_U,size_max>::Inverse( const U& u ) const { return inv_U( u ); }

template <typename T , int size_max> inline T& DifferenceSequence<T,size_max>::Action( const T& u , T& t ) const { return t += u; }
template <typename V , list<V> E(const V&) , V E_inv(const V&) , V enum_V(const int&) , int enum_V_inv(const V&) , typename T , T& m_T(const T&,T&) , const T& e_T() , typename U , U inv_U(const U&) , T& o_U(const U&,T&) , int size_max> inline T& EnumerationDifferenceSequence<V,E,E_inv,enum_V,enum_V_inv,T,m_T,e_T,U,inv_U,o_U,size_max>::Action( const U& u , T& t ) const { return o_U( u , t ); }


inline DEXPR( int , bound , 100000 , 100 ); // 0が5個
using V = pair<int,int>;
inline list<V> E( const V& v )
{
  const int& x = get<0>( v );
  const int A = get<1>( v ) + x;
  return A < bound ? list<V>{ { x , A } } : list<V>{};
}
inline V E_inv( const V& v )
{
  const int& x = get<0>( v );
  const int& A = get<1>( v );
  return { x , A < x ? A : A - x };
}
inline V enum_V( const int& i ) { return { i / bound , i % bound }; }
inline int enum_V_inv( const V& v ) { return get<0>( v ) * bound + get<1>( v ); }
inline int& add( const int& t0 , int& t1 ) { return t1 += t0; }
inline const int& zero() { static const int& z = 0; return z; }
inline int neg( const int& t ) { return -t; }

int main()
{
  UNTIE;
  LIBRARY_SEARCH;
  CIN_ASSERT( N , 1 , bound );
  CIN_ASSERT( M , 1 , bound );
  DEXPR( int , sqrt_bound , 316 , 10 );
  int C[bound] = {};
  static EnumerationDifferenceSequence<V,E,E_inv,enum_V,enum_V_inv,int,add,zero,int,neg,add,sqrt_bound * bound> DS{ sqrt_bound * bound };
  int L,R,X,Y,Z;
  REPEAT( N ){
    SET_ASSERT( L , 1 , bound );
    SET_ASSERT( R , L , bound );
    SET_ASSERT( X , 1 , bound );
    SET_ASSERT( Y , 0 , X - 1 );
    if( X < sqrt_bound ){
      ( ( Y < ( Z = L % X ) ? L += X : L ) -= Z ) += Y;
      ( ( Y <= ( Z = R % X ) ? R : R -= X ) -= Z ) += Y;
      DS.SubTreeAdd( { X , --L } , { { X , --R } } , 1 );
    } else {
      ( ( Y < ( Z = L % X ) ? L : L -= X ) -= Z ) += Y;
      --L;
      --R;
      while( ( L += X ) <= R ){
	++C[L];
      }
    }
  }
  int A;
  REPEAT( M ){
    SET_ASSERT( A , 1 , bound );
    int answer = C[--A];
    FOR( X , 1 , sqrt_bound ){
      answer += DS.Get( { X , A } );
    }
    COUT( answer );
  }
  QUIT;
}
yukicoder

結果

ソースコード
