結果

問題 No.1036 Make One With GCD 2
コンテスト
ユーザー kappybar
提出日時 2020-04-25 17:03:23
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,074 ms / 2,000 ms
コード長 1,686 bytes
コンパイル時間 1,815 ms
コンパイル使用メモリ 173,404 KB
実行使用メモリ 19,328 KB
最終ジャッジ日時 2024-09-16 13:47:04
合計ジャッジ時間 22,058 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#define rep(i,n) for(int i=0;i<n;i++)
using namespace std;
using ll =  long long ;
using P = pair<int,int> ;
const int INF = 1e9;
const int MOD = 1000000007;

ll gcd(ll i,ll j){
        if(j == 0) return i;
        return gcd(j,i%j);
}

struct SegmentTree{
        private:
        ll n;
        vector<ll> node;
        
        public:
        SegmentTree(vector<ll> v){
                ll sz = v.size();
                n = 1;
                while(n < sz) n *= 2;
                node.resize(2*n-1,0);
                for(ll i=0;i<sz;i++) node[i+n-1] = v[i];
                for(ll i=n-2;i>=0;i--) node[i] = gcd(node[2*i+1],node[2*i+2]);
        }
        
        void update(ll x,ll val){
                x += n-1;
                node[x] = val;
                while(x > 0){
                        x = (x-1)/2;
                        node[x] = gcd(node[2*x+1],node[2*x+2]);
                }
        }
        
        ll getgcd(ll a,ll b,ll k=0,ll l=0,ll r=-1){
                if(r < 0) r = n;
                if(r <= a || b <= l) return 0;
                if(a <= l && r <= b) return node[k];
                ll vl = getgcd(a,b,2*k+1,l,(l+r)/2);
                ll vr = getgcd(a,b,2*k+2,(l+r)/2,r);
                return gcd(vl,vr);
        }
};

int main(){
    int n;
    cin >> n;
    vector<ll> a(n);
    rep(i,n) cin >> a[i];
    SegmentTree seg(a);
    ll ans = 0;
    int j = 0;
    ll res = 0;
    rep(i,n){
        while(j < n && gcd(res,a[j]) != 1){
            res = gcd(res,a[j]);
            ++j;
        }
        ans += n - j;
        if(i == j) ++j;
        res = seg.getgcd(i+1,j+1);
    }
    cout << ans << endl;
    return 0;
}
0