結果
問題 |
No.1036 Make One With GCD 2
|
ユーザー |
![]() |
提出日時 | 2020-04-25 21:55:14 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 318 ms / 2,000 ms |
コード長 | 1,684 bytes |
コンパイル時間 | 2,307 ms |
コンパイル使用メモリ | 202,656 KB |
最終ジャッジ日時 | 2025-01-10 01:41:55 |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 41 |
ソースコード
#include <bits/stdc++.h> using namespace std; using lint = long long int; struct fast_ios { fast_ios(){ cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(20); }; } fast_ios_; // Static sequence sparse table // Complexity: O(NlogN) for precalculation, O(1) per query template<typename T, typename F> struct SparseTable { int N, lgN; T defaultT; F func; std::vector<std::vector<T>> data; std::vector<int> lgx_table; SparseTable(F func) : func(func) {} SparseTable(const std::vector<T> &sequence, T defaultT, F func) : N(sequence.size()), defaultT(defaultT), func(func) { lgx_table.resize(N + 1); for (int i = 2; i < N + 1; i++) lgx_table[i] = lgx_table[i >> 1] + 1; lgN = lgx_table[N] + 1; data.assign(lgN, std::vector<T>(N, defaultT)); data[0] = sequence; for (int d = 1; d < lgN; d++) { for (int i = 0; i + (1 << d) <= N; i++) { data[d][i] = func(data[d - 1][i], data[d - 1][i + (1 << (d - 1))]); } } } T get(int l, int r) { // [l, r), 0-indexed assert(l >= 0 and r <= N); if (l >= r) return defaultT; int d = lgx_table[r - l]; return func(data[d][l], data[d][r - (1 << d)]); } }; int main() { int N; cin >> N; vector<lint> A(N); for (auto &a : A) cin >> a; auto f = [](lint l, lint r) -> lint { return gcd(l, r); }; SparseTable<lint, decltype(f)> st(A, 0, f); lint ret = 0; int r = 1; for (int l = 0; l < N; l++) { if (l == r) r++; while (r <= N and st.get(l, r) > 1) r++; ret += N - r + 1; } cout << ret << '\n'; }