結果

問題 No.913 木の燃やし方
ユーザー beetbeet
提出日時 2019-10-18 22:29:26
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,825 bytes
コンパイル時間 2,145 ms
コンパイル使用メモリ 208,988 KB
実行使用メモリ 18,704 KB
最終ジャッジ日時 2024-06-25 17:34:33
合計ジャッジ時間 7,993 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 WA -
testcase_04 WA -
testcase_05 AC 84 ms
5,376 KB
testcase_06 WA -
testcase_07 AC 23 ms
5,376 KB
testcase_08 TLE -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
using Int = long long;
template<typename T1,typename T2> inline void chmin(T1 &a,T2 b){if(a>b) a=b;}
template<typename T1,typename T2> inline void chmax(T1 &a,T2 b){if(a<b) a=b;}


namespace MonotoneMinima{
  template<typename T,typename F>
  void induce(Int l,Int r,Int a,Int b,vector<Int> &dp,F dist){
    if(l==r) return;
    Int m=(l+r)>>1;
    Int &idx=(dp[m]=a);
    T res=dist(idx,m);
    for(Int i=a;i<b;i++){
      T tmp=dist(i,m);
      if(tmp<res) res=tmp,idx=i;
    }
    induce<T>(l,m,a,idx+1,dp,dist);
    induce<T>(m+1,r,idx,b,dp,dist);
  }

  template<typename T,typename F>
  vector<Int> args(Int n,F dist){
    vector<Int> dp(n,-1);
    induce<T>(0,n,0,n,dp,dist);
    return dp;
  }
}


template <typename E>
struct SegmentTree{
  using H = function<E(E,E)>;
  Int n,height;
  H h;
  E ei;
  vector<E> laz;
  SegmentTree(H h,E ei):h(h),ei(ei){}
  void init(Int n_){
    n=1;height=0;
    while(n<n_) n<<=1,height++;
    laz.assign(2*n,ei);
  }
  inline void eval(Int k){
    if(laz[k]==ei) return;
    laz[(k<<1)|0]=h(laz[(k<<1)|0],laz[k]);
    laz[(k<<1)|1]=h(laz[(k<<1)|1],laz[k]);
    laz[k]=ei;
  }
  inline void thrust(Int k){
    for(Int i=height;i;i--) eval(k>>i);
  }
  void update(Int a,Int b,E x){
    thrust(a+=n);
    thrust(b+=n-1);
    for(Int l=a,r=b+1;l<r;l>>=1,r>>=1){
      if(l&1) laz[l]=h(laz[l],x),l++;
      if(r&1) --r,laz[r]=h(laz[r],x);
    }
  }
  E get_val(Int a){
    thrust(a+=n);
    return laz[a];
  }
  void set_val(Int a,E x){
    thrust(a+=n);
    laz[a]=x;
  }
};


struct FastIO{
  FastIO(){
    cin.tie(0);
    ios::sync_with_stdio(0);
  }
}fastio_beet;

//INSERT ABOVE HERE
signed main(){
  Int n;
  cin>>n;
  vector<Int> as(n);
  for(Int i=0;i<n;i++) cin>>as[i];

  const Int INF = 1e18;
  auto g=[&](Int a,Int b){return min(a,b);};
  SegmentTree<Int> seg(g,INF);
  seg.init(n);

  for(Int uku=0;uku<2;uku++){
    vector<Int> sm(n+1,0);
    for(Int i=0;i<n;i++) sm[i+1]=sm[i]+as[i];

    auto dist=
      [&](Int i,Int j)->Int{
        if(i<j) return INF;
        return (i-j+1)*(i-j+1)+(sm[i+1]-sm[j]);
      };
    auto bs=MonotoneMinima::args<Int>(n,dist);

    auto latte=[&](Int k){return uku?n-k:k;};
    for(Int i=0;i<n;i++){
      Int l=min(latte(i),latte(bs[i]+1));
      Int r=max(latte(i),latte(bs[i]+1));
      seg.update(l,r,dist(bs[i],i));
    }
    reverse(as.begin(),as.end());
  }

  {
    vector<Int> sm(n+1,0);
    for(Int i=0;i<n;i++) sm[i+1]=sm[i]+as[i];

    auto dist=
      [&](Int i,Int j)->Int{
        if(i<j) return INF;
        return (i-j+1)*(i-j+1)+(sm[i+1]-sm[j]);
      };

    const int BS = 300;
    for(int i=0;i<n;i++)
      for(int l=1;l<BS&&i+l<=n;l++)
        seg.update(i,i+l,dist(i+l-1,i));
  }

  for(Int i=0;i<n;i++)
    cout<<seg.get_val(i)<<"\n";
  cout<<flush;
  return 0;
}
0