結果

問題 No.705 ゴミ拾い Hard
ユーザー beetbeet
提出日時 2018-06-21 14:31:06
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 392 ms / 1,500 ms
コード長 2,420 bytes
コンパイル時間 2,384 ms
コンパイル使用メモリ 187,888 KB
実行使用メモリ 12,672 KB
最終ジャッジ日時 2024-06-30 17:36:21
合計ジャッジ時間 9,951 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 1 ms
6,940 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,944 KB
testcase_04 AC 2 ms
6,944 KB
testcase_05 AC 2 ms
6,944 KB
testcase_06 AC 2 ms
6,944 KB
testcase_07 AC 2 ms
6,944 KB
testcase_08 AC 2 ms
6,944 KB
testcase_09 AC 2 ms
6,940 KB
testcase_10 AC 1 ms
6,940 KB
testcase_11 AC 2 ms
6,940 KB
testcase_12 AC 1 ms
6,944 KB
testcase_13 AC 1 ms
6,940 KB
testcase_14 AC 2 ms
6,940 KB
testcase_15 AC 3 ms
6,944 KB
testcase_16 AC 2 ms
6,944 KB
testcase_17 AC 3 ms
6,944 KB
testcase_18 AC 2 ms
6,944 KB
testcase_19 AC 3 ms
6,944 KB
testcase_20 AC 3 ms
6,940 KB
testcase_21 AC 3 ms
6,940 KB
testcase_22 AC 3 ms
6,940 KB
testcase_23 AC 3 ms
6,940 KB
testcase_24 AC 383 ms
12,672 KB
testcase_25 AC 388 ms
12,644 KB
testcase_26 AC 392 ms
12,544 KB
testcase_27 AC 386 ms
12,544 KB
testcase_28 AC 380 ms
12,552 KB
testcase_29 AC 383 ms
12,580 KB
testcase_30 AC 363 ms
12,544 KB
testcase_31 AC 389 ms
12,544 KB
testcase_32 AC 368 ms
12,544 KB
testcase_33 AC 253 ms
12,544 KB
testcase_34 AC 265 ms
12,544 KB
testcase_35 AC 297 ms
12,544 KB
testcase_36 AC 340 ms
12,596 KB
testcase_37 AC 296 ms
12,544 KB
testcase_38 AC 349 ms
12,600 KB
testcase_39 AC 364 ms
12,604 KB
testcase_40 AC 2 ms
6,940 KB
testcase_41 AC 2 ms
6,940 KB
testcase_42 AC 2 ms
6,940 KB
testcase_43 AC 385 ms
12,544 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
using Int = long long;
//BEGIN CUT HERE
template<typename T>
struct MonotoneMinima{
  using F = function<T(int, int)>;
  int n;
  F dist;
  T INF;
  vector<T> dp;

  // dist(i,l) + dist(j,k) >= dist(i,k) + dist(j,l)
  // if( i<=j && k<=l )
  MonotoneMinima(int n,F dist,T INF):
    n(n),dist(dist),INF(INF),dp(n+1,INF){}

  void dfs(int l,int r,int a,int b){
    if(l==r) return;
    int m=(l+r)>>1;
    T ma=INF;
    int mi=-1;
    for(int i=a;i<b;i++){
      T tmp=dp[i+1]+dist(i,m);
      if(tmp<ma) ma=tmp,mi=i;
    }
    dp[m]=min(dp[m],ma);
    dfs(l,m,a,mi+1);
    dfs(m+1,r,mi,b);
  }
  
  void solve(int l,int r){
    if(l+1==r){
      dp[l]=min(dp[l],dp[r]+dist(l,l));
      return;
    }
    int m=(l+r)>>1;
    solve(m,r);
    dfs(l,m,m,r);
    solve(l,m);
  }

  T solve(){
    dp[n]=0;
    solve(0,n);
    return dp[0];
  }
};
//END CUT HERE
//INSERT ABOVE HERE

signed YUKI_703(){
  int n;
  cin>>n;
  vector<Int> a(n),x(n),y(n);
  for(int i=0;i<n;i++) cin>>a[i];
  for(int i=0;i<n;i++) cin>>x[i];
  for(int i=0;i<n;i++) cin>>y[i];

  auto dist=
    [&](int i,int j)->Int{
      Int s=abs(a[i]-x[j]);
      Int t=abs(y[j]);
      return s*s+t*t;
    };

  const Int INF = 1e18;
  MonotoneMinima<Int> monge(n,dist,INF);
  
  cout<<monge.solve()<<endl;
  return 0;
}
/*
  verified on 2018/06/21
  https://yukicoder.me/problems/no/703
*/

signed YUKI_704(){
  int n;
  cin>>n;
  vector<Int> a(n),x(n),y(n);
  for(int i=0;i<n;i++) cin>>a[i];
  for(int i=0;i<n;i++) cin>>x[i];
  for(int i=0;i<n;i++) cin>>y[i];

  auto dist=
    [&](int i,int j)->Int{
      Int s=abs(a[i]-x[j]);
      Int t=abs(y[j]);
      return s+t;
    };

  const Int INF = 1e18;
  MonotoneMinima<Int> monge(n,dist,INF);
  
  cout<<monge.solve()<<endl;
  return 0;
}
/*
  verified on 2018/06/21
  https://yukicoder.me/problems/no/704
*/

signed YUKI_705(){
  int n;
  cin>>n;
  vector<Int> a(n),x(n),y(n);
  for(int i=0;i<n;i++) cin>>a[i];
  for(int i=0;i<n;i++) cin>>x[i];
  for(int i=0;i<n;i++) cin>>y[i];

  auto dist=
    [&](int i,int j)->Int{
      Int s=abs(a[i]-x[j]);
      Int t=abs(y[j]);
      return s*s*s+t*t*t;
    };

  const Int INF = 1e18;
  MonotoneMinima<Int> monge(n,dist,INF);
  
  cout<<monge.solve()<<endl;
  return 0;
}
/*
  verified on 2018/06/21
  https://yukicoder.me/problems/no/705
*/


signed main(){
  //YUKI_703();
  //YUKI_704();
  YUKI_705();
  return 0;
}
0