結果

問題 No.409 ダイエット
ユーザー umezo
提出日時 2025-05-14 02:59:00
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 39 ms / 2,000 ms
コード長 1,327 bytes
コンパイル時間 3,623 ms
コンパイル使用メモリ 280,772 KB
実行使用メモリ 7,844 KB
最終ジャッジ日時 2025-05-14 02:59:09
合計ジャッジ時間 8,685 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 92
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define ALL(v) v.begin(), v.end()
typedef long long ll;

#include <bits/stdc++.h>
using namespace std;
template <class T> using V=vector<T>;
template <class T> using VV=V<V<T>>;

template<typename T>
class ConvexHullTrick{
  struct F{
    T a,b;
    F(T a,T b):a(a),b(b){}
  };

  deque<F> deq;

  bool check(F &f1,F &f2,F &f3){
    return (f2.a-f1.a)*(f3.b-f2.b)>=(f2.b-f1.b)*(f3.a-f2.a);
  }

  ll f(F &f1,T x){return f1.a*x+f1.b;}

public:
  // a_{prev} >= a
  void add_line(T a,T b){
    F f1=F(a,b);
    while(deq.size()>=2 && check(deq[deq.size()-2],deq[deq.size()-1],f1)){
      deq.pop_back();
    }
    deq.push_back(f1);
  }

  // x_{prev} <= x
  ll query(T x){
    while(deq.size()>=2 && f(deq[0],x)>=f(deq[1],x)){
      deq.pop_front();
    }
    return f(deq[0],x);
  }
  
  deque<F> get_deq(){
    return deq;
  }
};

int main(){
  ios::sync_with_stdio(false);
  std::cin.tie(nullptr);
  
  ll n,a,b,w;
  cin>>n>>a>>b>>w;
  V<ll> D(n+1);
  rep(i,n) cin>>D[i+1];
  
  ConvexHullTrick<ll> cht;
  cht.add_line(0,0);
  V<ll> dp(n+1);
  for(ll i=1;i<=n;i++){
    dp[i]=cht.query(i)+a*(1-i)+i*(i-1)/2*b+D[i];
    cht.add_line(-b*i,dp[i]+a*i+i*(i+1)/2*b);
  }
  ll mi=1e18;
  rep(i,n+1) mi=min(mi,dp[i]-(n-i)*a+(n-i)*(n-i+1)/2*b+w);
  cout<<mi<<endl;
  
  return 0;
}
0