結果

問題 No.1079 まお
ユーザー snukesnuke
提出日時 2020-06-12 22:22:35
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,537 ms / 2,000 ms
コード長 4,452 bytes
コンパイル時間 2,967 ms
コンパイル使用メモリ 225,896 KB
実行使用メモリ 38,248 KB
最終ジャッジ日時 2024-06-24 05:15:25
合計ジャッジ時間 32,180 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 5 ms
5,376 KB
testcase_05 AC 14 ms
5,376 KB
testcase_06 AC 14 ms
5,376 KB
testcase_07 AC 42 ms
5,376 KB
testcase_08 AC 41 ms
5,376 KB
testcase_09 AC 43 ms
5,376 KB
testcase_10 AC 43 ms
5,376 KB
testcase_11 AC 41 ms
5,376 KB
testcase_12 AC 723 ms
14,928 KB
testcase_13 AC 843 ms
17,272 KB
testcase_14 AC 1,024 ms
19,684 KB
testcase_15 AC 1,164 ms
22,068 KB
testcase_16 AC 1,302 ms
24,388 KB
testcase_17 AC 1,466 ms
24,296 KB
testcase_18 AC 1,411 ms
24,424 KB
testcase_19 AC 1,503 ms
24,420 KB
testcase_20 AC 1,484 ms
24,420 KB
testcase_21 AC 1,440 ms
24,424 KB
testcase_22 AC 1,486 ms
27,492 KB
testcase_23 AC 1,481 ms
27,496 KB
testcase_24 AC 1,491 ms
27,496 KB
testcase_25 AC 1,523 ms
27,496 KB
testcase_26 AC 1,537 ms
27,236 KB
testcase_27 AC 1,347 ms
17,768 KB
testcase_28 AC 1,370 ms
25,572 KB
testcase_29 AC 1,366 ms
38,248 KB
testcase_30 AC 1,360 ms
27,368 KB
testcase_31 AC 1,340 ms
18,276 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define fi first
#define se second
#define rep(i,n) for(int i = 0; i < (n); ++i)
#define rrep(i,n) for(int i = 1; i <= (n); ++i)
#define drep(i,n) for(int i = (n)-1; i >= 0; --i)
#define srep(i,s,t) for (int i = s; i < t; ++i)
#define rng(a) a.begin(),a.end()
#define rrng(a) a.rbegin(),a.rend()
#define maxs(x,y) (x = max(x,y))
#define mins(x,y) (x = min(x,y))
#define isin(x,l,r) ((l) <= (x) && (x) < (r))
#define pb push_back
#define eb emplace_back
#define sz(x) (int)(x).size()
#define pcnt __builtin_popcountll
#define uni(x) x.erase(unique(rng(x)),x.end())
#define snuke srand((unsigned)clock()+(unsigned)time(NULL));
#define show(x) cout<<#x<<" = "<<x<<endl;
#define PQ(T) priority_queue<T,v(T),greater<T> >
#define bn(x) ((1<<x)-1)
#define dup(x,y) (((x)+(y)-1)/(y))
#define newline puts("")
#define v(T) vector<T>
#define vv(T) v(v(T))
using namespace std;
typedef long long int ll;
typedef unsigned uint;
typedef unsigned long long ull;
typedef pair<int,int> P;
typedef tuple<int,int,int> T;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<ll> vl;
typedef vector<P> vp;
typedef vector<T> vt;
inline int getInt() { int x; scanf("%d",&x); return x;}
template<typename T>inline istream& operator>>(istream&i,v(T)&v)
{rep(j,sz(v))i>>v[j];return i;}
template<typename T>string join(const v(T)&v)
{stringstream s;rep(i,sz(v))s<<' '<<v[i];return s.str().substr(1);}
template<typename T>inline ostream& operator<<(ostream&o,const v(T)&v)
{if(sz(v))o<<join(v);return o;}
template<typename T1,typename T2>inline istream& operator>>(istream&i,pair<T1,T2>&v)
{return i>>v.fi>>v.se;}
template<typename T1,typename T2>inline ostream& operator<<(ostream&o,const pair<T1,T2>&v)
{return o<<v.fi<<","<<v.se;}
template<typename T>inline ll suma(const v(T)& a) { ll res(0); for (auto&& x : a) res += x; return res;}
const double eps = 1e-10;
const ll LINF = 1001002003004005006ll;
const int INF = 1001001001;
#define dame { puts("-1"); return 0;}
#define yn {puts("Yes");}else{puts("No");}
const int MX = 200005;
typedef v(vl) vvl;

// coordinate compression
struct X {
  typedef int T;
  vector<T> d;
  X() {}
  X(vector<T>& a): d(a) {
    init();
    for (T& na : a) na = (*this)(na);
  }
  void add(const T& x) { d.pb(x);}
  void init() {
    sort(rng(d));
    d.erase(unique(rng(d)), d.end());
  }
  int size() const { return sz(d);}
  T operator[](int i) const { return d[i];}
  int operator()(const T& x) const { return upper_bound(rng(d),x)-d.begin()-1;}
  int get(int x) {
    int i = upper_bound(rng(d),x)-d.begin()-1;
    if (i == sz(d)) return -1;
    if (d[i] != x) return -1;
    return i;
  }
} xs;
//

// sparse table
int bsr(uint x) { return 31^__builtin_clz(x);}
struct Sparse {
  typedef P Int;
  int n;
  vv(Int) x;
  Sparse(int n=0):n(n),x(1,v(Int)(n)) {}
  Sparse(v(Int)& a):n(sz(a)),x(1,a) { init();}
  void init() {
    int i = 0;
    for (int w = 2; w <= n; w <<= 1) {
      x.pb(v(Int)(n-w+1));
      rep(j,sz(x.back())) {
        x.back()[j] = min(x[i][j], x[i][j+(w>>1)]);
      }
      ++i;
    }
  }
  Int get(int l, int r) {
    if (l >= r) return P(INF,l); // INF!!
    int i = bsr(r-l);
    return min(x[i][l], x[i][r-(1<<i)]);
  }
} s;

vvi is;
vvl iss;
typedef pair<int,ll> LP;
LP get(int x, int r) {
  int i = lower_bound(rng(is[x]), r) - is[x].begin();
  return LP(i,iss[x][i]);
}
LP get(int x, int l, int r) {
  LP a = get(x,l);
  LP b = get(x,r);
  return LP(b.fi-a.fi, b.se-a.se);
}

int n;
int k;
vi a;
ll ans;
void dfs(int l, int r) {
  if (l == r) return;
  int c = s.get(l,r).se;
  int rc = s.get(c+1,r).se;
  if (c+1 == r || rc == r || a[c] != a[rc]) rc = r;
  if (c-l < rc-c) {
    srep(i,l,c+1) {
      int ni = xs.get(k-xs[a[i]]);
      if (ni == -1) continue;
      LP p = get(ni,c,rc);
      ans += p.se;
      ans -= ll(i-1)*p.fi;
    }
  } else {
    srep(i,c,rc) {
      int ni = xs.get(k-xs[a[i]]);
      if (ni == -1) continue;
      LP p = get(ni,l,c+1);
      ans += ll(i+1)*p.fi;
      ans -= p.se;
    }
  }
  cerr<<l<<" "<<r<<" "<<c<<" "<<rc<<" "<<ans<<endl;
  dfs(l,c);
  dfs(c+1,r);
}

int main() {
  scanf("%d%d",&n,&k);
  a = vi(n);
  cin>>a;
  xs = X(a);
  s = Sparse(n);
  rep(i,n) s.x[0][i] = P(a[i],i);
  is = vvi(sz(xs));
  iss = vvl(sz(xs),vl(1));
  rep(i,n) is[a[i]].pb(i);
  rep(i,n) iss[a[i]].pb(i);
  rep(i,sz(xs)) {
    rep(j,sz(iss[i])-1) iss[i][j+1] += iss[i][j];
  }
  s.init();
  dfs(0,n);
  cout<<ans<<endl;
  return 0;
}
0