結果
問題 | No.332 数列をプレゼントに |
ユーザー |
![]() |
提出日時 | 2015-12-25 13:07:44 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 4,272 bytes |
コンパイル時間 | 1,608 ms |
コンパイル使用メモリ | 171,816 KB |
実行使用メモリ | 19,456 KB |
最終ジャッジ日時 | 2024-09-18 23:46:13 |
合計ジャッジ時間 | 11,097 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 35 WA * 6 TLE * 1 |
ソースコード
#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)%MODint 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++なら復元適当でも間に合うでしょw// よく考えたら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> 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;// 大きい値の全探索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(0<=x-s && x-s<=mxval){X_S[x-s]=true;}else if(0<s && s<=mxval){// 小さい値の方に突っ込むdp[0][s] = true;pred[0][s] = 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);// 大きいやつREPR(mask,1<<more.size()){ll s = 0;REP(i,more.size()){if(mask>>i&1)s+=more[i];}if(x-s==ansid || s==cst){REP(i,more.size()){if(mask>>i&1)result.insert(more[i]);}}}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;}