結果

問題 No.1017 Reiwa Sequence
ユーザー HIR180
提出日時 2020-04-03 21:48:51
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 44 ms / 2,000 ms
コード長 2,832 bytes
コンパイル時間 3,500 ms
コンパイル使用メモリ 199,684 KB
実行使用メモリ 22,684 KB
最終ジャッジ日時 2024-07-03 02:08:56
合計ジャッジ時間 36,179 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 50
権限があれば一括ダウンロードができます

ソースコード

diff # 

//Let's join Kaede Takagaki Fan Club !!
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
typedef pair<int,P> P1;
typedef pair<P,P> P2;
#define pu push
#define pb push_back
#define mp make_pair
#define eps 1e-7
#define INF 1000000000
#define fi first
#define sc second
#define rep(i,x) for(int i=0;i<x;i++)
#define repn(i,x) for(int i=1;i<=x;i++)
#define SORT(x) sort(x.begin(),x.end())
#define ERASE(x) x.erase(unique(x.begin(),x.end()),x.end())
#define POSL(x,v) (lower_bound(x.begin(),x.end(),v)-x.begin())
#define POSU(x,v) (upper_bound(x.begin(),x.end(),v)-x.begin())
#define all(x) x.begin(),x.end()
template<class T>
void dmp(T a){
	rep(i,a.size()) cout << a[i] << " ";
	cout << endl;
}
template<class T>
bool chmax(T&a, T b){
	if(a < b){
		a = b;
		return 1;
	}
	return 0;
}
template<class T>
bool chmin(T&a, T b){
	if(a > b){
		a = b;
		return 1;
	}
	return 0;
}
template<class T>
void g(T &a){
	cin >> a;
}
template<class T>
void o(const T &a,bool space=false){
	cout << a << (space?' ':'\n');
}
//ios::sync_with_stdio(false);
const ll mod = 1000000007;//998244353
template<class T>
void add(T&a,T b){
	a+=b;
	if(a >= mod) a-=mod;
}
#define SZ 1000000
bitset<SZ*2+5>cur, nxt, nxt2;
int n,a[150005],aa[150005],cr[150005];
P fr[SZ*2+5];
int main(){
	ios::sync_with_stdio(false);
	g(n);
	repn(i,n) g(aa[i]);
	vector<P>vec;
	repn(i,n) vec.pb(mp(aa[i], i));
	mt19937 mt(614);
	shuffle(vec.begin(), vec.end(), mt);
	repn(i,n){
		a[i] = vec[i-1].fi;
		cr[i] = vec[i-1].sc;
	}
	
	cur[SZ] = 1;
	fr[SZ] = mp(0, 0);
	auto output=[&](vector<P>vi){
		vector<int>ans(n+1);
		repn(i,n) ans[i] = 0;
		rep(i,vi.size()) ans[cr[vi[i].fi]] = a[vi[i].fi] * vi[i].sc;
		repn(i,n) o(ans[i], 1);
		cout<<'\n';
		return;
	};
	repn(i,n){
	     //cout << i << endl;
		if(cur[SZ-a[i]]){ //cout << i << endl;
			vector<P>ans;
			ans.pb(mp(i, 1));
			int cur = SZ-a[i];
			//cout << fr[cur].fi << " " << fr[cur].sc << endl;
			while(cur != SZ){
				P k = fr[cur];
				ans.pb(mp(k.fi, k.sc));
				cur -= k.sc * a[k.fi];
			}
			o("Yes");
			output(ans);
			return 0;
		}
		if(cur[SZ+a[i]]){
			vector<P>ans;
			ans.pb(mp(a[i], -1));
			int cur = SZ+a[i];
			while(cur != SZ){
				P k = fr[cur];
				ans.pb(mp(k.fi, k.sc));
				cur -= k.sc * a[k.fi];
			}
			o("Yes");
			output(ans);
			return 0;
		}
		
		nxt = cur | (cur << a[i]);
		nxt ^= cur;
		for(int ii=nxt._Find_first();ii< nxt.size();ii = nxt._Find_next(ii)){
			fr[ii]  = mp(i, 1); //cout << ii << " " << i << " " << 1 << endl;
		}
		nxt ^= cur;
		nxt2 = nxt | (cur >> a[i]);
		nxt2 ^= nxt;
		for(int ii=nxt2._Find_first();ii<nxt2.size();ii = nxt2._Find_next(ii)){
			fr[ii]  = mp(i, -1); //cout << ii << " " << i << " " << -1 << endl;
		}
		cur = nxt2^nxt;
	}
	o("No");
}
0