結果
| 問題 |
No.3153 probability max K
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2025-05-20 23:43:02 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 812 bytes |
| コンパイル時間 | 1,432 ms |
| コンパイル使用メモリ | 163,608 KB |
| 実行使用メモリ | 7,848 KB |
| 最終ジャッジ日時 | 2025-05-20 23:43:08 |
| 合計ジャッジ時間 | 5,000 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | WA * 20 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<ll,ll> P;
#define REP(i,n) for(ll i=0;i<ll(n);i++)
const ll MOD=998244353;
ll f(ll x,ll n){
ll res=1;
while(n!=0LL){
if(n%2!=0LL) res=(res*x)%MOD;
x=(x*x)%MOD;
n/=2;
}
return res;
}
ll g(ll a,ll b){
if(b==0) return a;
else return g(b,a%b);
}
int main(void){
cin.tie(nullptr); ios_base::sync_with_stdio(false);
ll i,j;
ll N,K;
cin >> N >> K;
vector<ll> a(N+1);
for(i=1;i<=N;i++){
cin >> a[i];
}
ll x=1,y=1;
for(i=1;i<=N;i++){
if(a[i]<K) continue;
ll e=g(a[i],K);
ll q=a[i]/e,p=K/e;
ll qq=f(q,MOD-2);
x*=qq*p;
x%=MOD;
e=g(a[i],K-1);
q=a[i]/e,p=(K-1)/e;
qq=f(q,MOD-2);
y*=qq*p;
y%=MOD;
}
cout << (x-y+MOD)%MOD << endl;
return 0;
}