結果
| 問題 |
No.2075 GCD Subsequence
|
| コンテスト | |
| ユーザー |
kwm_t
|
| 提出日時 | 2022-09-16 22:32:57 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 2,188 bytes |
| コンパイル時間 | 4,429 ms |
| コンパイル使用メモリ | 252,560 KB |
| 最終ジャッジ日時 | 2025-02-07 10:06:31 |
|
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 4 WA * 24 |
ソースコード
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
//using mint = modint1000000007;
//const int mod = 1000000007;
using mint = modint998244353;
const int mod = 998244353;
//const int INF = 1e9;
//const long long LINF = 1e18;
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define rep2(i,l,r)for(int i=(l);i<(r);++i)
#define rrep(i, n) for (int i = (n-1); i >= 0; --i)
#define rrep2(i,l,r)for(int i=(r-1);i>=(l);--i)
#define all(x) (x).begin(),(x).end()
#define allR(x) (x).rbegin(),(x).rend()
#define endl "\n"
#define P pair<int,int>
template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; }
template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; }
// 素因数分解
vector<int> PrimeFactorization(long long n) {
vector<int> ret;
for (long long i = 2; i*i <= n; i++) {
int c = 0;
while (n%i == 0) { // 素数で割り切れなくなるまで割っていく
c++;//割った個数を配列に足す
n /= i;
}
if (0 != c) ret.push_back(i);
if (1 == n) break;
}
if (1 != n) ret.push_back(n);
return ret;
}
int popcount(long long n) {
n = (n & 0x5555555555555555) + (n >> 1 & 0x5555555555555555);
n = (n & 0x3333333333333333) + (n >> 2 & 0x3333333333333333);
n = (n & 0x0f0f0f0f0f0f0f0f) + (n >> 4 & 0x0f0f0f0f0f0f0f0f);
n = (n & 0x00ff00ff00ff00ff) + (n >> 8 & 0x00ff00ff00ff00ff);
n = (n & 0x0000ffff0000ffff) + (n >> 16 & 0x0000ffff0000ffff);
n = (n & 0x00000000ffffffff) + (n >> 32 & 0x00000000ffffffff);
return n;
}
int main() {
ios::sync_with_stdio(false);
cin.tie(nullptr);
mint ans = 0;
vector<mint>memo(1000006);
memo[1] = 1;
int n; cin >> n;
rep(i, n) {
int a; cin >> a;
auto v = PrimeFactorization(a);
mint val = 1;
vector<int>list;
rep2(j, 1,1 << v.size()) {
int cnt = 0;
int x = 1;
rep(k, v.size()) {
if (1 & (j >> k)) {
x *= v[k];
cnt++;
}
}
list.push_back(x);
if (1 == cnt) val += memo[x];
else val -= memo[x];
}
rep(j, list.size())memo[list[j]] += val;
ans += val;
}
cout << ans.val() << endl;
return 0;
}
kwm_t