結果

問題 No.1554 array_and_me
ユーザー momoyuu
提出日時 2024-08-07 02:05:38
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,379 bytes
コンパイル時間 1,525 ms
コンパイル使用メモリ 112,676 KB
実行使用メモリ 6,400 KB
最終ジャッジ日時 2024-08-07 02:05:43
合計ジャッジ時間 5,179 ms
ジャッジサーバーID
(参考情報)
judge3 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 1
other AC * 10 WA * 31
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;
using ll = long long;
#include<atcoder/modint>
using mint = atcoder::modint998244353;
void solve(){
ll n,k;
cin>>n>>k;
vector<ll> a(n);
for(int i = 0;i<n;i++) cin>>a[i];
ll sum = 0;
for(int i = 0;i<n;i++) sum += a[i];
vector<ll> use(n),res(n);
for(int i = 0;i<n;i++) {
use[i] = (a[i]*k) / sum;
res[i] = (a[i]*k) % sum;
}
vector<int> idx(n);
for(int i = 0;i<n;i++) idx[i] = i;
sort(idx.begin(),idx.end(),[&](int i,int j){
return a[i] > a[j];
});
ll now = 0;
for(int i = 0;i<n;i++) now += use[i];
now = k - now;
for(int i = 0;i<now;i++){
int ni = idx[i];
use[ni]++;
}
mint ans = 1;
ll all = k;
vector<mint> fac(k+1,0);
fac[0] = 1;
for(int i = 1;i<=k;i++) fac[i] = fac[i-1] * i;
for(int i = 0;i<n;i++){
mint tmp = fac[all] * (fac[all-use[i]].inv()) * (fac[use[i]].inv());
mint p = mint(a[i]) * (mint(sum).inv());
p = p.pow(use[i]);
tmp *= p;
ans *= tmp;
all -= use[i];
}
cout<<ans.val()<<endl;
return;
for(int i = 0;i<n;i++) cout<<use[i]<<" ";
cout<<endl;
}
int main(){
cin.tie(nullptr);
ios::sync_with_stdio(false);
int t;
cin>>t;
while(t--){
solve();
}
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0