#710814 (C++17) No.137 貯金箱の焦り

• 問題: No.137 貯金箱の焦り
• ユーザー: 沙耶花
• 提出日時: 2021-10-28 06:41:00
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 3,116 ms / 5,000 ms
• コード長: 1,547 bytes
• コンパイル時間: 4,151 ms
• コンパイル使用メモリ: 267,232 KB
• 実行使用メモリ: 7,168 KB
• 最終ジャッジ日時: 2024-05-17 05:17:16
• 合計ジャッジ時間: 58,468 ms
• ジャッジサーバーID （参考情報）: judge3 / judge2

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	211 ms; 7,160 KB
testcase_01	AC	576 ms; 7,036 KB
testcase_02	AC	3,094 ms; 7,156 KB
testcase_03	AC	54 ms; 7,156 KB
testcase_04	AC	1,025 ms; 7,164 KB
testcase_05	AC	354 ms; 7,164 KB
testcase_06	AC	950 ms; 7,164 KB
testcase_07	AC	950 ms; 7,160 KB
testcase_08	AC	1,008 ms; 7,160 KB
testcase_09	AC	1,008 ms; 7,168 KB
testcase_10	AC	961 ms; 7,156 KB
testcase_11	AC	3,042 ms; 7,160 KB
testcase_12	AC	3,032 ms; 7,164 KB
testcase_13	AC	410 ms; 7,164 KB
testcase_14	AC	2,988 ms; 7,160 KB
testcase_15	AC	3,116 ms; 7,164 KB
testcase_16	AC	3,091 ms; 7,032 KB
testcase_17	AC	3,021 ms; 7,160 KB
testcase_18	AC	3,035 ms; 7,164 KB
testcase_19	AC	3,054 ms; 7,036 KB
testcase_20	AC	2,847 ms; 7,160 KB
testcase_21	AC	3,042 ms; 7,164 KB
testcase_22	AC	3,024 ms; 7,028 KB
testcase_23	AC	2,882 ms; 7,032 KB
testcase_24	AC	2,895 ms; 7,164 KB
testcase_25	AC	2,998 ms; 7,164 KB

ソースコード

#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){
	vector<long long> a1(a.size()),a2(a.size()),b1(b.size()),b2(b.size());
	rep(i,a.size()){
		a1[i] = a[i].val()>>15;
		a2[i] = a[i].val()&((1<<15)-1);
	}
	rep(i,b.size()){
		b1[i] = b[i].val()>>15;
		b2[i] = b[i].val()&((1<<15)-1);
	}
	
	vector<long long> A = convolution_ll(a1,b1), C = convolution_ll(a2,b2);
	
	rep(i,a1.size())a1[i] += a2[i];
	rep(i,b1.size())b1[i] += b2[i];
	
	vector<long long> B = convolution_ll(a1,b1);
	
	rep(i,B.size()){
		B[i] -= A[i] + C[i];
	}
	
	
	vector<T> ret(A.size(),0);
	rep(i,A.size()){
		ret[i] += (long long)T(A[i]).val()<<30;
		ret[i] += (long long)T(B[i]).val()<<15;
		ret[i] += T(C[i]).val();
	}
	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;
}