結果

問題 No.137 貯金箱の焦り
ユーザー 沙耶花沙耶花
提出日時 2021-10-29 00:20:42
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,048 ms / 5,000 ms
コード長 1,747 bytes
コンパイル時間 5,033 ms
コンパイル使用メモリ 266,148 KB
実行使用メモリ 6,316 KB
最終ジャッジ日時 2024-04-16 13:18:31
合計ジャッジ時間 24,192 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 74 ms
6,188 KB
testcase_01 AC 193 ms
6,188 KB
testcase_02 AC 1,037 ms
6,060 KB
testcase_03 AC 21 ms
6,168 KB
testcase_04 AC 348 ms
6,060 KB
testcase_05 AC 126 ms
6,116 KB
testcase_06 AC 332 ms
6,188 KB
testcase_07 AC 333 ms
6,056 KB
testcase_08 AC 350 ms
6,192 KB
testcase_09 AC 349 ms
6,060 KB
testcase_10 AC 334 ms
6,192 KB
testcase_11 AC 1,041 ms
6,188 KB
testcase_12 AC 1,042 ms
6,064 KB
testcase_13 AC 143 ms
6,060 KB
testcase_14 AC 1,020 ms
6,184 KB
testcase_15 AC 1,035 ms
6,188 KB
testcase_16 AC 1,040 ms
6,188 KB
testcase_17 AC 1,016 ms
6,060 KB
testcase_18 AC 1,023 ms
6,188 KB
testcase_19 AC 1,045 ms
6,316 KB
testcase_20 AC 967 ms
6,188 KB
testcase_21 AC 1,034 ms
6,060 KB
testcase_22 AC 1,041 ms
6,188 KB
testcase_23 AC 1,001 ms
6,060 KB
testcase_24 AC 983 ms
6,192 KB
testcase_25 AC 1,048 ms
6,064 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
using mint = modint;
using namespace std;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf 1000000000000000000

template <class T>
vector<T> convolution_mint(vector<T> a,vector<T> b){
	
	static constexpr unsigned long long M0 = 998244353;
	static constexpr unsigned long long M1 = 754974721;
	static constexpr unsigned long long M2 = 469762049;
	
	vector<long long> aa(a.size()),bb(b.size());
	rep(i,a.size())aa[i] = a[i].val();
	rep(i,b.size())bb[i] = b[i].val();
	
	auto c0 = convolution<M0>(aa,bb);
	auto c1 = convolution<M1>(aa,bb);
	auto c2 = convolution<M2>(aa,bb);
	
	vector<mint> ret(c0.size());
	vector<long long> m = {M0,M1,M2};
	rep(i,c0.size()){		
		
		ret[i] += c0[i];
		
		{
			long long cur = c0[i];
			cur %= M1;
			cur = c1[i]-cur;
			if(cur<0)cur += M1;
			cur *= 416537774;
			cur %= M1;
			
			mint m = M0;
			ret[i] += m*cur;
			
			cur *= M0;
			cur += c0[i];
			cur %= M2;
			cur = c2[i] - cur;
			if(cur<0)cur += M2;
			cur *= 429847628;
			cur %= M2;
			m *= M1;
			ret[i] += m * cur;
		}
	}
	
	return ret;
}

int main(){
	
	mint::set_mod(1234567891);
	
	int n;
	cin>>n;
	
	long long m;
	cin>>m;
	
	vector<int> a(n);
	rep(i,n)cin>>a[i];
	
	vector<mint> dp(25005,0);
	dp[0] = 1;
	
	rep(i,n){
		vector<mint> ndp(25005,0);
		rep(j,dp.size()){
			if(dp[j]==0)continue;
			ndp[j] += dp[j];
			
			ndp[j+a[i]] += dp[j];
		}
		swap(dp,ndp);
	}
	
	vector<mint> ddp(25005,0);
	ddp[0] = 1;
	
	while(m!=0){
		vector<mint> ndp = convolution_mint(dp,ddp);
		vector<mint> temp(25005,0);
		rep(i,ndp.size()){
			if(i%2!=m%2)continue;
			temp[i>>1] = ndp[i];
		}
		swap(temp,ddp);
		m>>=1;
	}
	
	cout<<ddp[0].val()<<endl;
	
	return 0;
}
0