結果
問題 | No.1036 Make One With GCD 2 |
ユーザー |
![]() |
提出日時 | 2024-02-29 19:33:27 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 391 ms / 2,000 ms |
コード長 | 1,739 bytes |
コンパイル時間 | 2,756 ms |
コンパイル使用メモリ | 203,000 KB |
最終ジャッジ日時 | 2025-02-19 22:01:56 |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 41 |
ソースコード
#include <bits/stdc++.h>using namespace std;using ll = long long;template<class S, S (*op)(S, S), S (*e)()>struct SegTree {private:vector<S> seg; int N = 1;S _prod (int l, int r, int idx, int bitl, int bitr) const{if (r < bitl || l > bitr) return e();if (l <= bitl && bitr <= r) return seg[idx];ll bitm = (bitl+bitr)/2;return op(_prod(l, r, idx*2+1, bitl, bitm), _prod(l, r, idx*2+2, bitm+1, bitr));}public:SegTree (ll n) : SegTree(vector<S>(n, e())) {}SegTree (const vector<S> &v){ll n = v.size(); while (N < n) N *= 2; seg.resize(N*2-1, e());for (ll i=0; i<n; i++) seg[i+N-1] = v[i];for (ll i=N-2; i>=0; i--) seg[i] = op(seg[i*2+1], seg[i*2+2]);}void set(int loc, S val){loc += N-1; seg[loc] = val;while (loc != 0) loc = (loc-1)/2, seg[loc] = op(seg[loc*2+1], seg[loc*2+2]);}S prod (int l, int r) const {return _prod(l, r, 0, 0, N-1);}S all_prod() const {return seg[0];}S get (int i) const{return seg[i+N-1];}void show() const{for (int i=N-1; i<N*2-1; i++) cout << seg[i] << " ";cout << endl;}};//Range GCD Queryusing S = ll;S op(S a, S b){return gcd(a, b);}S e(){return 0;}int main(){cin.tie(nullptr);ios_base::sync_with_stdio(false);ll N, ans=0, r=0;cin >> N;vector<S> v(N);for (int i=0; i<N; i++) cin >> v[i];SegTree<S, op, e> tree(v);for (int l=0; l<N; l++){while(tree.prod(l, r) != 1 && r < N) r++;ans += N-r;if (l == r) r++;}cout << ans << endl;return 0;}