結果

問題 No.1036 Make One With GCD 2
ユーザー SPD_9X2
提出日時 2020-04-24 22:39:23
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 708 ms / 2,000 ms
コード長 1,345 bytes
コンパイル時間 2,630 ms
コンパイル使用メモリ 201,376 KB
最終ジャッジ日時 2025-01-10 00:16:47
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
#include <iostream>
#include <limits>
#include <numeric>
#include <type_traits>

using namespace std;

#define rep(i,n,m) for(ll (i)=(n);(i)<(m);(i)++)
#define rrep(i,n,m) for(ll (i)=(n);(i)>(m);(i)--)
using ll = long long;
const ll mod = 998244353;


template< typename T >
struct SparseTable {
  vector< vector< T > > st;
  vector< ll > lookup;

  SparseTable(const vector< T > &v) {
    ll b = 0;
    while((1 << b) <= v.size()) ++b;
    st.assign(b, vector< T >(1 << b));
    for(ll i = 0; i < v.size(); i++) {
      st[0][i] = v[i];
    }
    for(ll i = 1; i < b; i++) {
      for(ll j = 0; j + (1 << i) <= (1 << b); j++) {
        st[i][j] = gcd(st[i - 1][j], st[i - 1][j + (1 << (i - 1))]);
      }
    }
    lookup.resize(v.size() + 1);
    for(ll i = 2; i < (ll)lookup.size(); i++) {
      lookup[i] = lookup[i >> 1] + 1;
    }
  }

  inline T rmq(ll l, ll r) {
    ll b = lookup[r - l];
    return gcd(st[b][l], st[b][r - (1 << b)]);
  }
};

int main(){

    ll N;
    cin >> N;
    vector<ll> A(N);
    rep(i,0,N) cin >> A[i];

    SparseTable<ll> spt(A);

    ll l = 0 , ans = N*(N-1)/2 + N;
    rep(r,0,N){
        if (A[r] == 1) continue;

        while (spt.rmq(l,r+1) == 1) l++;
        ans -= (r-l+1);

        //cout << l << " " << r << " " << spt.rmq(l,r+1) << endl;
    }

    cout << ans << endl;
}
0