結果
問題 |
No.1036 Make One With GCD 2
|
ユーザー |
|
提出日時 | 2020-04-24 23:32:45 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,260 ms / 2,000 ms |
コード長 | 1,625 bytes |
コンパイル時間 | 3,277 ms |
コンパイル使用メモリ | 201,488 KB |
最終ジャッジ日時 | 2025-01-10 00:39:14 |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 41 |
ソースコード
#include <bits/stdc++.h> using namespace std; template<typename T> struct Segtree { int n, n_org; T e; vector<T> dat; typedef function<T(T a, T b)> Func; Func f; Segtree(){} Segtree(int n_input, Func f_input, T e_input){ initialize(n_input, f_input, e_input); } void initialize(int n_input, Func f_input, T e_input){ n_org = n_input; f = f_input; e = e_input; n = 1; while(n < n_input) n <<= 1; dat.assign(2*n-1, e); } void build(vector<T>& A){ for(int k=0; k<int(A.size()); k++) dat[k+n-1] = A[k]; for(int k=n-2; k>=0; k--) dat[k] = f(dat[2*k+1], dat[2*k+2]); } void update(int k, T a){ k += n - 1; dat[k] = a; while(k > 0){ k = (k - 1)/2; dat[k] = f(dat[2*k+1], dat[2*k+2]); } } T get(int k){ return dat[k+n-1]; } T between(int a, int b){ return query(a, b+1, 0, 0, n); } T query(int a, int b, int k, int l, int r){ if(r<=a || b<=l) return e; if(a<=l && r<=b) return dat[k]; T vl = query(a, b, 2*k+1, l, (l+r)/2); T vr = query(a, b, 2*k+2, (l+r)/2, r); return f(vl, vr); } }; int main(){ int N; cin >> N; vector<int64_t> A(N); for(int i=0; i<N; i++) cin >> A[i]; Segtree<int64_t> st(N, [](auto a, auto b){ return gcd(a, b); }, 0); st.build(A); int r = 0; int64_t ans = 0; for(int l=0; l<N; l++){ while(r < N && st.between(l, r) != 1) r++; ans += N-r; } cout << ans << endl; return 0; }