結果
| 問題 | No.2952 Invision of Multiples |
| ユーザー |
 沙耶花
|
| 提出日時 | 2024-10-25 22:46:32 |
| 言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,374 bytes |
| コンパイル時間 | 5,015 ms |
| コンパイル使用メモリ | 268,236 KB |
| 実行使用メモリ | 39,456 KB |
| 最終ジャッジ日時 | 2024-10-25 22:46:51 |
| 合計ジャッジ時間 | 19,173 ms |
|
ジャッジサーバーID （参考情報） |
judge1 / judge3 |
（要ログイン）
テストケース
テストケース表示| 入力 | 結果 | 実行時間 実行使用メモリ |
|---|---|---|
| testcase_00 | AC | 2 ms
39,456 KB |
| testcase_01 | AC | 2 ms
6,816 KB |
| testcase_02 | AC | 3,319 ms
33,152 KB |
| testcase_03 | WA | - |
| testcase_04 | WA | - |
| testcase_05 | WA | - |
| testcase_06 | WA | - |
| testcase_07 | WA | - |
| testcase_08 | WA | - |
| testcase_09 | WA | - |
| testcase_10 | WA | - |
| testcase_11 | TLE | - |
| testcase_12 | -- | - |
| testcase_13 | -- | - |
| testcase_14 | -- | - |
| testcase_15 | -- | - |
| testcase_16 | -- | - |
| testcase_17 | -- | - |
| testcase_18 | -- | - |
| testcase_19 | -- | - |
| testcase_20 | -- | - |
| testcase_21 | -- | - |
| testcase_22 | -- | - |
| testcase_23 | -- | - |
| testcase_24 | -- | - |
| testcase_25 | -- | - |
| testcase_26 | -- | - |
| testcase_27 | -- | - |
| testcase_28 | -- | - |
| testcase_29 | -- | - |
| testcase_30 | -- | - |
| testcase_31 | -- | - |
| testcase_32 | -- | - |
| testcase_33 | -- | - |
| testcase_34 | -- | - |
| testcase_35 | -- | - |
| testcase_36 | -- | - |
| testcase_37 | -- | - |
| testcase_38 | -- | - |
| testcase_39 | -- | - |
| testcase_40 | -- | - |
| testcase_41 | -- | - |
| testcase_42 | -- | - |
| testcase_43 | -- | - |
ソースコード
#include <stdio.h>
#include <atcoder/all>
#include <bits/stdc++.h>
using namespace std;
using namespace atcoder;
using mint = modint998244353;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf32 1000000001
#define Inf64 1000000000000000000LL
int main(){
int N,M;
cin>>N>>M;
int B = min(M,150);
vector cnt(B+1,vector<mint>(M+1));
for(int i=1;i<=B;i++){
for(int j=1;j<=M;j++){
int a0 = i-1;
int m = j;
int n = (M-a0+i-1)/i;
cnt[i][j] = floor_sum(n,m,i,a0);
cnt[i][j] /= mint(M/i) * mint(M/j);
}
}
vector<mint> num(B+1);
vector<int> a(N);
rep(i,N)scanf("%d",&a[i]);
fenwick_tree<mint> fw(M+1);
mint ans = 0;
rep(i,N){
//small, *
for(int j=1;j<=B;j++){
ans += num[j] * cnt[j][a[i]];
}
//big big
for(int j=a[i];j<=M;j+=a[i]){
ans += fw.sum(j+1,M+1)/mint(M/a[i]);
}
if(a[i]<=B)num[a[i]]++;
else{
for(int j=a[i];j<=M;j+=a[i]){
fw.add(j,mint(M/a[i]).inv());
}
}
}
num.assign(B+1,0);
reverse(a.begin(),a.end());
for(int i=1;i<=B;i++){
for(int j=1;j<=M;j++){
int a0 = j-1;
int m = i;
int n = (M-a0+j-1)/j;
cnt[i][j] = floor_sum(n,m,j,a0);
cnt[i][j] /= mint(M/j) * mint(M/i);
}
}
rep(i,N){
//big, small
if(a[i]>B){
for(int j=1;j<=B;j++){
ans += num[j] * cnt[j][a[i]];
}
}
else num[a[i]]++;
}
rep(i,N)ans *= mint(M/a[i]);
cout<<ans.val()<<endl;
return 0;
}