結果
問題 | No.2876 Infection |
ユーザー | 蜜蜂 |
提出日時 | 2024-09-06 23:12:53 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 255 ms / 2,000 ms |
コード長 | 2,308 bytes |
コンパイル時間 | 6,861 ms |
コンパイル使用メモリ | 284,908 KB |
実行使用メモリ | 15,576 KB |
最終ジャッジ日時 | 2024-09-06 23:13:04 |
合計ジャッジ時間 | 10,629 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 16 ms
15,404 KB |
testcase_01 | AC | 17 ms
15,448 KB |
testcase_02 | AC | 17 ms
15,396 KB |
testcase_03 | AC | 16 ms
15,472 KB |
testcase_04 | AC | 15 ms
15,464 KB |
testcase_05 | AC | 254 ms
15,576 KB |
testcase_06 | AC | 253 ms
15,440 KB |
testcase_07 | AC | 254 ms
15,484 KB |
testcase_08 | AC | 255 ms
15,520 KB |
testcase_09 | AC | 255 ms
15,432 KB |
testcase_10 | AC | 53 ms
15,496 KB |
testcase_11 | AC | 50 ms
15,428 KB |
testcase_12 | AC | 41 ms
15,424 KB |
testcase_13 | AC | 73 ms
15,488 KB |
testcase_14 | AC | 218 ms
15,428 KB |
testcase_15 | AC | 155 ms
15,492 KB |
testcase_16 | AC | 23 ms
15,440 KB |
testcase_17 | AC | 17 ms
15,448 KB |
testcase_18 | AC | 62 ms
15,516 KB |
testcase_19 | AC | 17 ms
15,380 KB |
testcase_20 | AC | 123 ms
15,472 KB |
testcase_21 | AC | 19 ms
15,524 KB |
testcase_22 | AC | 136 ms
15,472 KB |
testcase_23 | AC | 74 ms
15,560 KB |
testcase_24 | AC | 56 ms
15,436 KB |
testcase_25 | AC | 131 ms
15,424 KB |
testcase_26 | AC | 18 ms
15,476 KB |
testcase_27 | AC | 61 ms
15,452 KB |
testcase_28 | AC | 20 ms
15,440 KB |
testcase_29 | AC | 199 ms
15,568 KB |
ソースコード
// g++-13 1.cpp -std=c++17 -O2 -I . #include <bits/stdc++.h> using namespace std; #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include<ext/pb_ds/assoc_container.hpp> #include<ext/pb_ds/tree_policy.hpp> #include<ext/pb_ds/tag_and_trait.hpp> using namespace __gnu_pbds; #include <atcoder/all> using namespace atcoder; using ll = long long; using ld = long double; using vi = vector<int>; using vvi = vector<vi>; using vll = vector<ll>; using vvll = vector<vll>; using vld = vector<ld>; using vvld = vector<vld>; using vst = vector<string>; using vvst = vector<vst>; #define fi first #define se second #define pb push_back #define eb emplace_back #define pq_big(T) priority_queue<T,vector<T>,less<T>> #define pq_small(T) priority_queue<T,vector<T>,greater<T>> #define all(a) a.begin(),a.end() #define rep(i,start,end) for(ll i=start;i<(ll)(end);i++) #define per(i,start,end) for(ll i=start;i>=(ll)(end);i--) #define uniq(a) sort(all(a));a.erase(unique(all(a)),a.end()) using mint = modint998244353; random_device seed; mt19937_64 randint(seed()); ll grr(ll mi, ll ma) { // [mi, ma) return mi + randint() % (ma - mi); } const int MAX=510000; const long long MOD=998244353; long long fac[MAX],finv[MAX],inv[MAX]; void COMinit(){ fac[0]=fac[1]=1; finv[0]=finv[1]=1; inv[1]=1; for(int i=2;i<MAX;i++){ fac[i]=fac[i-1]*i%MOD; inv[i]=MOD-inv[MOD%i]*(MOD/i)%MOD; finv[i]=finv[i-1]*inv[i]%MOD; } } int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); COMinit(); int n,x;cin>>n>>x; vector<mint> dp(n,0),coef(n,0); mint p=100-x,q=x; p/=100,q/=100; vector<mint> powp(n+10,0); dp[n-1]=1; coef[0]=1; powp[0]=1; rep(i,1,n){ coef[i]=coef[i-1]*q; } rep(i,0,n){ coef[i]/=fac[i]; powp[i+1]=powp[i]*p; } mint ans=0; rep(i,1,n+1){ // rep(j,0,n){ // cout<<dp[j].val()<<" "; // } // cout<<endl; rep(j,0,n+1-i){ ans+=dp[j]; // dp[j]*=fac[j]; } rep(j,n+1-i,n){ dp[j]=0; } rep(j,0,n){ dp[j]*=fac[j]; } reverse(all(dp)); vector<mint> ndp=convolution(dp,coef); while(ndp.size()>n){ ndp.pop_back(); } reverse(all(ndp)); rep(j,0,n){ ndp[j]*=powp[j]; ndp[j]/=fac[j]; } dp=ndp; } cout<<ans.val()<<endl; }