結果

問題 No.15 カタログショッピング
コンテスト
ユーザー plasma_e
提出日時 2017-06-06 00:04:47
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 12 ms / 5,000 ms
コード長 4,965 bytes
コンパイル時間 1,179 ms
コンパイル使用メモリ 104,912 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-22 07:06:46
合計ジャッジ時間 1,718 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 10
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<iostream>
#include<vector>
#include<map>
#include<set>
#include<algorithm>
#include<numeric>
#include<limits>
#include<bitset>
#include<functional>
#include<type_traits>
#include<array>
#include<queue>
#include<memory>

namespace lib
{
	template<class T, class Comp>using p_queue = std::priority_queue<T, std::vector<T>, Comp>;

	template<std::uint64_t Mod>struct modnum;
	template<class T>constexpr T pow(T base, std::size_t p)
	{
		if (p == 0)
		{
			return T(1);
		}
		else if (p == 1)
		{
			return base;
		}
		else if (p == 2)
		{
			return base*base;
		}
		else if (p % 2 == 0)
		{
			return pow(pow(base, p / 2), 2);
		}
		else
		{
			return pow(pow(base, p / 2), 2)*base;
		}
	}
	template<std::uint64_t Mod>constexpr auto inverse(modnum<Mod> const&);

	template<std::uint64_t Mod>struct modnum
	{
		static constexpr auto mod = Mod;
		std::uint64_t val;
		modnum() = default;
		constexpr modnum(std::uint64_t v):val(v%Mod)
		{

		}

		constexpr modnum& operator+=(modnum const& v)
		{
			val += v.val;
			val %= mod;
			return *this;
		}
		constexpr modnum& operator-=(modnum const& v)
		{
			val += mod - v.val;
			val %= mod;
			return *this;
		}
		constexpr modnum& operator*=(modnum const& v)
		{
			val *= v.val;
			val %= mod;
			return *this;
		}
		constexpr modnum& operator/=(modnum const& v)
		{
			return operator*=(inverse(v));
		}
	};
	template<std::uint64_t Mod>constexpr auto operator+(modnum<Mod> lhs, modnum<Mod>const& rhs)
	{
		return lhs += rhs;
	}
	template<std::uint64_t Mod>constexpr auto operator+(modnum<Mod> lhs, int rhs)
	{
		return lhs += rhs;
	}
	template<std::uint64_t Mod>constexpr auto operator+(int lhs, modnum<Mod> rhs)
	{
		return rhs += lhs;
	}
	template<std::uint64_t Mod>constexpr auto operator-(modnum<Mod> lhs, modnum<Mod>const& rhs)
	{
		return lhs -= rhs;
	}
	template<std::uint64_t Mod>constexpr auto operator*(modnum<Mod> lhs, modnum<Mod>const& rhs)
	{
		return lhs *= rhs;
	}
	template<std::uint64_t Mod>constexpr auto operator*(modnum<Mod> lhs, int rhs)
	{
		return lhs *= rhs;
	}
	template<std::uint64_t Mod>constexpr auto operator*(int lhs, modnum<Mod> rhs)
	{
		return rhs *= lhs;
	}

	template<std::uint64_t Mod>constexpr auto operator/(modnum<Mod> lhs, modnum<Mod>const& rhs)
	{
		return lhs /= rhs;
	}

	template<std::uint64_t Mod>constexpr auto inverse(modnum<Mod>const& base)
	{
		return pow(base, Mod - 2);
	}

	template<class T>constexpr auto clamp(T v)
	{
		return std::max(v, T());
	}
	template<class T>void sort(std::vector<T>& vec)
	{
		std::sort(vec.begin(), vec.end());
	}
	template<class T, class Comp>void sort(std::vector<T>& vec, Comp&& comp)
	{
		std::sort(vec.begin(), vec.end(), std::forward<Comp>(comp));
	}

	template<class T>void vector_input(std::vector<T>& vec)
	{
		for (auto& v : vec)
		{
			std::cin >> v;
		}
	}

	template<class T>auto lower_bound(std::vector<T>const& vec, T v)
	{
		return std::distance(vec.begin(), std::lower_bound(vec.begin(), vec.end(), v));
	}
	template<class T>auto upper_bound(std::vector<T>const& vec, T v)
	{
		return std::distance(vec.begin(), std::upper_bound(vec.begin(), vec.end(), v));
	}

	struct scope_exit
	{
		std::function<void(void)> func;
		scope_exit(std::function<void(void)> f):func(f)
		{

		}
		~scope_exit()
		{
			func();
		}
	};

	template<int val>using int_tag = std::integral_constant<int, val>;

	int gcd(int x, int y)
	{
		if (x > y)
		{
			return gcd(y, x);
		}
		if (x == 0)
		{
			return y;
		}
		return gcd(y%x, x);
	}
}

void Main();
int main()
{
	std::cin.tie(nullptr);
	std::cin.sync_with_stdio(false);
	Main();
}

void Main()
{
	int N, S;
	std::cin >> N >> S;
	std::vector<int> data(N);
	for (auto& v : data)
	{
		std::cin >> v;
	}
	std::vector<std::pair<int, int>> prefix(1 << (N / 2));
	std::vector<std::pair<int, int>> suffix(1 << (N - N / 2));
	std::vector<int> result;
	prefix[0] = std::make_pair(0, 0);
	suffix[0] = std::make_pair(0, 0);
	for (int i{}; i < N / 2; ++i)
	{
		for (int j{}; j < (1 << i); ++j)
		{
			prefix[(1 << i) + j] = std::make_pair(prefix[j].first + data[i], prefix[j].second | (1 << (N - i - 1)));
		}
	}
	for (int i = N / 2; i < N; ++i)
	{
		for (int j{}; j < (1 << (i - N / 2)); ++j)
		{
			suffix[(1 << (i - N / 2)) + j] = std::make_pair(suffix[j].first + data[i], suffix[j].second | (1 << (N - i - 1)));
		}
	}
	lib::sort(prefix);
	lib::sort(suffix);
	int start = (1 << (N - N / 2)) - 1;
	for (int i{}; i < (1 << N / 2); ++i)
	{
		for (; start >= 0; --start)
		{
			if (prefix[i].first + suffix[start].first <= S)
			{
				break;
			}
		}
		if (prefix[i].first + suffix[start].first != S)
		{
			continue;
		}
		for (int j = start; j >= 0; --j)
		{
			if (prefix[i].first + suffix[j].first != S)
			{
				break;
			}
			result.emplace_back(prefix[i].second + suffix[j].second);
		}
	}
	lib::sort(result, std::greater<>());
	for (auto v : result)
	{
		for (int i = N - 1; i >= 0; --i)
		{
			if (v&(1 << i))
			{
				std::cout << N - i << " ";
			}
		}
		std::cout << std::endl;
	}
}
0