結果

問題 No.332 数列をプレゼントに
ユーザー rickythetarickytheta
提出日時 2015-12-25 01:45:27
言語 C++11
(gcc 11.4.0)
結果
WA  
実行時間 -
コード長 4,175 bytes
コンパイル時間 2,430 ms
コンパイル使用メモリ 172,772 KB
実行使用メモリ 19,200 KB
最終ジャッジ日時 2024-09-18 23:38:54
合計ジャッジ時間 6,974 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 5 ms
5,248 KB
testcase_01 AC 4 ms
5,376 KB
testcase_02 AC 3 ms
5,376 KB
testcase_03 AC 4 ms
5,376 KB
testcase_04 AC 3 ms
5,376 KB
testcase_05 AC 30 ms
12,800 KB
testcase_06 WA -
testcase_07 AC 78 ms
18,560 KB
testcase_08 AC 53 ms
14,976 KB
testcase_09 AC 8 ms
5,760 KB
testcase_10 AC 9 ms
7,424 KB
testcase_11 AC 37 ms
12,416 KB
testcase_12 AC 77 ms
18,816 KB
testcase_13 WA -
testcase_14 AC 24 ms
10,240 KB
testcase_15 AC 73 ms
16,640 KB
testcase_16 AC 66 ms
17,152 KB
testcase_17 AC 79 ms
18,816 KB
testcase_18 WA -
testcase_19 AC 78 ms
17,408 KB
testcase_20 AC 73 ms
16,384 KB
testcase_21 AC 80 ms
17,408 KB
testcase_22 AC 70 ms
15,872 KB
testcase_23 AC 77 ms
19,072 KB
testcase_24 AC 79 ms
19,200 KB
testcase_25 AC 74 ms
19,200 KB
testcase_26 AC 79 ms
19,200 KB
testcase_27 AC 73 ms
19,072 KB
testcase_28 AC 73 ms
19,200 KB
testcase_29 AC 76 ms
19,200 KB
testcase_30 AC 76 ms
19,200 KB
testcase_31 AC 78 ms
19,200 KB
testcase_32 AC 81 ms
19,072 KB
testcase_33 AC 78 ms
19,200 KB
testcase_34 AC 89 ms
19,072 KB
testcase_35 AC 138 ms
19,200 KB
testcase_36 AC 2 ms
5,376 KB
testcase_37 AC 84 ms
19,072 KB
testcase_38 AC 81 ms
19,200 KB
testcase_39 AC 81 ms
19,200 KB
testcase_40 AC 89 ms
19,200 KB
testcase_41 AC 13 ms
19,200 KB
testcase_42 AC 91 ms
19,200 KB
testcase_43 AC 88 ms
19,200 KB
testcase_44 AC 88 ms
19,072 KB
testcase_45 AC 1,166 ms
19,200 KB
testcase_46 AC 41 ms
10,496 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef complex<double> P;
typedef pair<int,int> pii;
#define REP(i,n) for(ll i=0;i<n;++i)
#define REPR(i,n) for(ll i=1;i<n;++i)
#define FOR(i,a,b) for(ll i=a;i<b;++i)

#define DEBUG(x) cout<<#x<<": "<<x<<endl
#define DEBUG_VEC(v) cout<<#v<<":";REP(i,v.size())cout<<" "<<v[i];cout<<endl
#define ALL(a) (a).begin(),(a).end()

#define MOD (ll)(1e9+7)
#define ADD(a,b) a=((a)+(b))%MOD
#define FIX(a) ((a)%MOD+MOD)%MOD

int main(){
  ll n,x;
  cin>>n>>x;
  vl a(n);
  REP(i,n)cin>>a[i];
  // ウルトラスーパー考察ターイム

  // 総積条件が怪しい
  // ウェイト置くナップサックはX<=1e12で間に合わず
  // 当然全部殴ろうにも2^100パターン出て無理

  // Xは1e12
  // A1~ANをソートしたとして、総積が1e100(オーダー)だとしておく
  // 最悪の場合ANを1e12オーダーとして、次は1e8ぐらい
  // N=100なのでO(NT)としてTが1e6ぐらいなら間に合う
  // 最大の重みが1e6のナップサックに落とせれば行けそう
  // 1e6を超えるAは、(1e6)^18=1e108で、たかだか18個
  // -> 最大2^18パターンを保持しておく
  // -> X-S<=1e6となるようなSだけを取っておく(|S|<=1e6)
  // -> 残ったAで1e6ナップサック
  // -> 和を求める
  // -> あとはこの値を元に復元する
  // C++なら復元適当でも間に合うでしょｗ

  // よく考えたら1e5より大きい数はたかだか20個
  // 閾値を10^Aとする
  // 閾値より大きい数は(10^A)^N=10^100 => N=100/A
  // 組み合わせ全探索ループは2^N
  // ナップサックは10^A*100
  // ぐーぐるせんせー！に 2^(100/x)-10^(x+2) を計算してもらったところx=4.5ぐらいが最適っぽい
  // ので閾値は10^4.5=30000 を閾値にしよう

  // 待って待って
  // 1e6を100個積んだら1e8になるな
  // つまり
  // 閾値を10^Aにする
  // 閾値以上の数は N=100/A
  // 組み合わせ全探索は 2^N
  // ナップサックは 10^A *100 *100
  // 2^(100/x) = 10^(x+4) を解く、だいたい3.8
  // 6400らしい
  // 本当かなぁ・・・

  ll limit = 6400;
  ll mxval = 640000;

  vl more,less;
  REP(i,n){
    if(a[i]>x)continue;
    if(a[i]>=limit)more.push_back(a[i]);
    else less.push_back(a[i]);
  }

  // 大きい値の全探索
  vector<bool> X_S(mxval+10,false);
  if(x>mxval) REP(mask,1<<more.size()){
    ll s = 0;
    REP(i,more.size()){
      if(mask>>i&1)s+=more[i];
    }
    if(x-s<0)continue;
    if(x-s>mxval)continue;
    X_S[x-s] = true;
  }

  // 小さい値のナップサック
  // vector<bool> dp(mxval+10,false);
  vector< vector<bool> > dp(n+1,vector<bool>(mxval+10,false));
  vector< vector<bool> > pred(n+1,vector<bool>(mxval+10,false));
  dp[0][0]=true;

  REP(i,less.size()){
    ll cost = less[i];
    REP(p,mxval+10){
      if(!dp[i][p])continue;
      dp[i+1][p] = true;
      pred[i+1][p] = false;
      if(p+cost<mxval+10){
        dp[i+1][p+cost] = true;
        pred[i+1][p+cost] = true;
      }
    }
  }

  // 合体
  ll ansid = -1;
  if(x<=mxval){
    // 合体しなくてよい
    if(dp[less.size()][x])ansid = x;
  }
  if(ansid==-1)REP(i,mxval+1){
    if(dp[less.size()][i]&&X_S[i]){
      ansid = i;
      break;
    }
  }
  if(ansid==-1){
    cout << "No" << endl;
    return 0;
  }

  // 復元
  multiset<ll> result; // ここに使う値をぶっこんでいく

  // dp,predから
  ll cst = ansid;
  REP(_i,less.size()){
    ll i = less.size()-_i;
    if(pred[i][cst]){
      result.insert(less[i-1]);
      cst -= less[i-1];
    }
  }
  assert(cst==0);

  // 大きいやつ
  if(x>mxval) REP(mask,1<<more.size()){
    ll s = 0;
    REP(i,more.size()){
      if(mask>>i&1)s+=more[i];
    }
    if(x-s==ansid){
      REP(i,more.size()){
        if(mask>>i&1)result.insert(more[i]);
      }
      break;
    }
  }
  multiset<ll>::iterator iter;
  REP(i,n){
    iter = result.find(a[i]);
    if(iter!=result.end()){
      cout<<"o";
      result.erase(iter);
    }else{
      cout<<"x";
    }
  }
  cout<<endl;

  return 0;
}
0