結果

問題 No.1857 Gacha Addiction
コンテスト
ユーザー 蜜蜂
提出日時 2021-12-31 19:16:20
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 990 ms / 6,000 ms
コード長 1,657 bytes
コンパイル時間 4,557 ms
コンパイル使用メモリ 250,468 KB
実行使用メモリ 31,208 KB
最終ジャッジ日時 2024-07-01 21:45:02
合計ジャッジ時間 34,741 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 43
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'int main()':
main.cpp:56:10: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
   56 |     auto [sz1,A,B]=que.top();
      |          ^
main.cpp:58:10: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
   58 |     auto [sz2,a,b]=que.top();
      |          ^
main.cpp:78:8: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
   78 |   auto [sz,q,r]=que.top();
      |        ^

ソースコード

diff #
raw source code 

//g++ 2.cpp -std=c++14 -O2 -I .
#include <bits/stdc++.h>
using namespace std;

#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 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())

// size pi/sx+1 (pi/s)^2x
using P = tuple<int,vll,vll>;

constexpr ll mod = 998244353;

int main(){
  ios::sync_with_stdio(false);
  cin.tie(nullptr);

  int n;
  ll s;
  cin>>n>>s;

  ll inv_s=pow_mod(s,mod-2,mod);

  pq_small(P) que;

  rep(i,0,n){
    ll p;
    cin>>p;
    p*=inv_s;
    p%=mod;
    que.push({1,{1,p},{0,(p*p)%mod}});
  }

  while(que.size()>1){
    auto [sz1,A,B]=que.top();
    que.pop();
    auto [sz2,a,b]=que.top();
    que.pop();

    vll x=convolution(A,a);
    vll y=convolution(A,b);
    vll z=convolution(B,a);

    rep(i,0,z.size()){
      if(i<y.size()){
        y[i]+=z[i];
        y[i]%=mod;
      }
      else{
        y.emplace_back(z[i]);
      }
    }

    que.push({sz1+sz2,x,y});
  }

  auto [sz,q,r]=que.top();

  ll ans=0;
  ll fractal=2;
  rep(i,1,r.size()){
    ans+=r[i]*fractal;
    ans%=mod;
    fractal*=i+2;
    fractal%=mod;
  }

  cout<<ans<<endl;
}
0