結果

問題 No.2075 GCD Subsequence
ユーザー SSRSSSRS
提出日時 2022-09-16 21:56:42
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,321 ms / 4,000 ms
コード長 1,061 bytes
コンパイル時間 1,536 ms
コンパイル使用メモリ 173,644 KB
実行使用メモリ 95,392 KB
最終ジャッジ日時 2023-08-23 15:24:46
合計ジャッジ時間 43,343 ms
ジャッジサーバーID
(参考情報)
judge12 / judge14
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1,177 ms
94,432 KB
testcase_01 AC 1,172 ms
94,472 KB
testcase_02 AC 1,175 ms
94,476 KB
testcase_03 AC 1,195 ms
94,432 KB
testcase_04 AC 1,201 ms
94,596 KB
testcase_05 AC 1,184 ms
94,468 KB
testcase_06 AC 1,193 ms
94,484 KB
testcase_07 AC 1,199 ms
94,476 KB
testcase_08 AC 1,276 ms
95,016 KB
testcase_09 AC 1,321 ms
95,228 KB
testcase_10 AC 1,286 ms
95,068 KB
testcase_11 AC 1,299 ms
95,392 KB
testcase_12 AC 1,277 ms
95,012 KB
testcase_13 AC 1,261 ms
95,132 KB
testcase_14 AC 1,301 ms
95,008 KB
testcase_15 AC 1,275 ms
94,952 KB
testcase_16 AC 1,271 ms
94,948 KB
testcase_17 AC 1,317 ms
95,268 KB
testcase_18 AC 1,315 ms
95,204 KB
testcase_19 AC 1,319 ms
95,152 KB
testcase_20 AC 1,288 ms
95,160 KB
testcase_21 AC 1,315 ms
95,272 KB
testcase_22 AC 1,301 ms
95,272 KB
testcase_23 AC 1,283 ms
95,384 KB
testcase_24 AC 1,296 ms
95,312 KB
testcase_25 AC 1,281 ms
95,224 KB
testcase_26 AC 1,288 ms
95,312 KB
testcase_27 AC 1,290 ms
95,184 KB
testcase_28 AC 1,287 ms
95,276 KB
testcase_29 AC 1,226 ms
94,436 KB
testcase_30 AC 1,196 ms
94,476 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
const int MAX = 1000000;
const long long MOD = 998244353;
int main(){
  int N;
  cin >> N;
  vector<int> A(N);
  for (int i = 0; i < N; i++){
    cin >> A[i];
  }
  vector<bool> prime(MAX + 1, true);
  vector<int> mobius(MAX + 1, 1);
  for (int i = 2; i <= MAX; i++){
    if (prime[i]){
      for (int j = 1; j * i <= MAX; j++){
        if (j > 1){
          prime[i * j] = false;
        }
        mobius[i * j] *= -1;
        if (j % i == 0){
          mobius[i * j] = 0;
        }
      }
    }
  }
  vector<vector<int>> F(MAX + 1);
  for (int i = 2; i <= MAX; i++){
    if (mobius[i] != 0){
      for (int j = i; j <= MAX; j += i){
        F[j].push_back(i);
      }
    }
  }
  vector<long long> dp(MAX + 1, 0);
  long long ans = 0;
  for (int i = 0; i < N; i++){
    long long sum = 1;
    for (int j : F[A[i]]){
      sum += dp[j] * (-mobius[j]) + MOD;
    }
    sum %= MOD;
    ans += sum;
    ans %= MOD;
    for (int j : F[A[i]]){
      dp[j] += sum;
      dp[j] %= MOD;
    }
  }
  cout << ans << endl;
}
0