結果
問題 |
No.868 ハイパー部分和問題
|
ユーザー |
![]() |
提出日時 | 2019-08-17 01:31:08 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 1,781 ms / 7,000 ms |
コード長 | 1,741 bytes |
コンパイル時間 | 846 ms |
コンパイル使用メモリ | 102,932 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-09-24 09:11:01 |
合計ジャッジ時間 | 23,052 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 38 |
ソースコード
#include <cstdio> #include <cstring> #include <iostream> #include <string> #include <cmath> #include <bitset> #include <vector> #include <map> #include <set> #include <queue> #include <deque> #include <algorithm> #include <complex> #include <unordered_map> #include <unordered_set> #include <random> #include <cassert> #include <fstream> #include <utility> #include <functional> #define popcount __builtin_popcount using namespace std; typedef long long int ll; typedef pair<int, int> P; int MOD[3]; bool isprime(int x){ for(int i=2; i*i<=x; i++){ if(x%i==0) return false; } return true; } random_device rnd; mt19937 mt(rnd()); uniform_int_distribution<> rndp(700000000, 1000000000); int main() { int n, k; cin>>n>>k; for(int i=0; i<3; i++){ int x=rndp(mt); for(int j=x; ; j++){ if(isprime(j)){ MOD[i]=j; //cout<<j<<endl; break; } } } int a[15010], a0[15010]; for(int i=0; i<n; i++) cin>>a[i], a0[i]=a[i]; int q; cin>>q; int ans[15010]={}; int x[15010], v[15010]; for(int i=0; i<q; i++){ cin>>x[i]>>v[i]; x[i]--; } for(int t=0; t<3; t++){ int dp[15010]={}; dp[0]=1; for(int i=0; i<n; i++){ a[i]=a0[i]; if(a[i]==0) continue; for(int j=k-a[i]; j>=0; j--){ dp[j+a[i]]+=dp[j]; if(dp[j+a[i]]>=MOD[t]) dp[j+a[i]]-=MOD[t]; } } for(int i=0; i<q; i++){ if(a[x[i]]!=0){ for(int j=0; j<=k-a[x[i]]; j++){ dp[j+a[x[i]]]+=(MOD[t]-dp[j]); if(dp[j+a[x[i]]]>=MOD[t]) dp[j+a[x[i]]]-=MOD[t]; } } a[x[i]]=v[i]; if(a[x[i]]!=0){ for(int j=k-a[x[i]]; j>=0; j--){ dp[j+a[x[i]]]+=dp[j]; if(dp[j+a[x[i]]]>=MOD[t]) dp[j+a[x[i]]]-=MOD[t]; } } if(dp[k]!=0) ans[i]=1; } } for(int i=0; i<q; i++) printf("%d\n", ans[i]); return 0; }