結果
問題 | No.1036 Make One With GCD 2 |
ユーザー |
![]() |
提出日時 | 2020-04-25 00:28:26 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 620 ms / 2,000 ms |
コード長 | 3,121 bytes |
コンパイル時間 | 2,874 ms |
コンパイル使用メモリ | 204,300 KB |
最終ジャッジ日時 | 2025-01-10 00:45:56 |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 41 |
ソースコード
#include <bits/stdc++.h>using namespace std;#define rep2(i, m, n) for(int i=int(m); i<int(n); ++i)#define drep2(i, m, n) for(int i=int(m-1); i>=int(n); --i)#define rep(i, n) rep2(i, 0, n)#define drep(i, n) drep2(i, n, 0)#define all(a) a.begin(), a.end()using ll = long long;using ld = long double;using V = vector<int>;using Vll = vector<ll>;using Vld = vector<ld>;using Vbo = vector<bool>;using VV = vector<V>;using VVll = vector<Vll>;using VVld = vector<Vld>;using VVbo = vector<Vbo>;using P = pair<int, int>;using Pll = pair<ll, ll>;using Pld = pair<ld, ld>;struct fast_ios { fast_ios(){ cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(20); }; } fast_ios_;template<typename T> inline bool chmin(T &a, T b) { if (a > b) { a = b; return true; } return false; }template<typename T> inline bool chmax(T &a, T b) { if (a < b) { a = b; return true; } return false; }template<typename T1, typename T2> istream &operator>>(istream &is, pair<T1, T2> &p) { is >> p.first >> p.second; return is; }template<typename T1, typename T2> ostream &operator<<(ostream &os, const pair<T1, T2> &p) { os << "(" << p.first << ", " << p.second << ")"; returnos; }template<typename T> istream &operator>>(istream &is, vector<T> &v) { for (auto &e : v) is >> e; return is; }template<typename T> ostream &operator<<(ostream &os, const vector<T> &v) { for (auto &e : v) os << e << " "; return os; }template<typename T> inline int count_between(vector<T> &a, T l, T r) { return lower_bound(all(a), r) - lower_bound(all(a), l); } // [l, r)inline int Log2(ll x) { int k; for (k = 0; x > 0; ++k) x >>= 1; return k; } // number of binary digitsconst int INF = 1<<30;const ll INFll = 1ll<<62;const ld EPS = 1e-10;const int MOD = int(1e9)+7;template<typename T>struct SegmentTree {int m, n; // m: size of original vector, n: size of segment treevector<T> d;T idt = 0;T op(T &a, T &b) {// return a + b;return gcd(a, b);}SegmentTree(vector<T> v) : m(v.size()) {n = 1; while(n < m) n <<= 1;d.assign(2*n-1, 0);rep(i, m) d[i+n-1] = v[i];drep(i, n-1) d[i] = op(d[i*2+1], d[i*2+2]);}void add(int k, T val) {k += n-1, d[k] += val;while (k) k = (k-1)/2, d[k] = op(d[2*k+1], d[2*k+2]);}void set(int k, T val) {k += n-1, d[k] = val;while (k) k = (k-1)/2, d[k] = op(d[2*k+1], d[2*k+2]);}T op_range(int a, int b, int k=0, int l=0, int r=-1) {if (r < 0) r = n;if (b <= l || r <= a) return idt;if (a <= l && r <= b) return d[k];T vl = op_range(a, b, 2*k+1, l, (l+r)/2);T vr = op_range(a, b, 2*k+2, (l+r)/2, r);return op(vl, vr);}T get(int a) { return d[n-1+a]; }vector<T> get_all() {vector<T> res(m);copy(d.begin()+n-1, d.begin()+n+m-1, res.begin());return res;}};int main() {ll n; cin >> n;Vll a(n); cin >> a;ll ans = 0;SegmentTree<ll> st(a);int j = 0;rep(i, n) {ll y = st.op_range(i, j+1);while (j < n && st.op_range(i, j+1) > 1) ++j;ans += n-j;if (i == j) ++j;}cout << ans << endl;return 0;}