結果
問題 | No.1036 Make One With GCD 2 |
ユーザー |
![]() |
提出日時 | 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 |
ソースコード
#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; }