結果
問題 | 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;}