結果

問題 No.15 カタログショッピング
コンテスト
ユーザー Chanyuh
提出日時 2019-10-03 18:09:03
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 26 ms / 5,000 ms
コード長 2,491 bytes
コンパイル時間 1,232 ms
コンパイル使用メモリ 128,408 KB
実行使用メモリ 7,268 KB
最終ジャッジ日時 2024-10-03 06:33:34
合計ジャッジ時間 1,799 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 10
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<iostream>
#include<string>
#include<cstdio>
#include<vector>
#include<cmath>
#include<algorithm>
#include<functional>
#include<iomanip>
#include<queue>
#include<ciso646>
#include<random>
#include<map>
#include<set>
#include<complex>
#include<bitset>
#include<stack>
#include<unordered_map>
#include<utility>
#include<tuple>
using namespace std;
typedef long long ll;
typedef unsigned int ui;
const ll mod = 1000000007;
const ll INF = (ll)1000000007 * 1000000007;
typedef pair<ll, vector<int>> P;
#define stop char nyaa;cin>>nyaa;
#define rep(i,n) for(int i=0;i<n;i++)
#define per(i,n) for(int i=n-1;i>=0;i--)
#define Rep(i,sta,n) for(int i=sta;i<n;i++)
#define Per(i,sta,n) for(int i=n-1;i>=sta;i--)
#define rep1(i,n) for(int i=1;i<=n;i++)
#define per1(i,n) for(int i=n;i>=1;i--)
#define Rep1(i,sta,n) for(int i=sta;i<=n;i++)
typedef long double ld;
typedef complex<ld> Point;
const ld eps = 1e-8;
const ld pi = acos(-1.0);
typedef pair<ll, ll> LP;

int n;
ll s,p[32];

map<ll,vector<int>> knapsack(vector<ll> x){
    int m=x.size();
    int M=(1 << m);
    int k=0;
    vector<ll> dp(M,0);
    map<ll,vector<int>> res;
    Rep(i,1,M){
        if (i==(1 << (k+1))) k++;
        dp[i]=dp[i^(1<<k)]+x[k];
    }
    rep(i,M){
        if(res.find(dp[i])==res.end()){
            res[dp[i]]={i};
        }
        else res[dp[i]].push_back(i);
    }
    return res;
}

void solve(){
    cin >> n >> s;
    rep(i,n){
        cin >> p[i];
    }
    vector<ll> s1={},s2={};
    map<ll,vector<int>> t1,t2;
    rep(i,n/2){
        s1.push_back(p[i]);
    }
    t1=knapsack(s1);
    vector<vector<int>> ans;
    int n2=n-(n/2),m2=(1 << n2);
    vector<ll> dp(m2,0);
    int k=0;
    rep(i,m2){
        if (i!=0){
        if (i==(1 << (k+1))) k++;
        dp[i]=dp[i^(1<<k)]+p[n/2+k];
        }
        if (s-dp[i]<0) continue;
        if (t1.find(s-dp[i])==t1.end()) continue;
        for(auto a:t1[s-dp[i]]){
            vector<int> vec={};
                rep(l,n/2){
                    if(a & (1<<l)) vec.push_back(l+1);
                }
                Rep(l,n/2,n){
                    if(i & (1<<(l-(n/2)))) vec.push_back(l+1);
                }
                ans.push_back(vec);
                }
        
    }

    sort(ans.begin(),ans.end());

    for(auto a:ans){
        for(auto c:a){
            cout << c << " ";
        }
        cout << "" << endl;
    }
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout << fixed << setprecision(50);
    solve();
}
0