結果
問題 |
No.1036 Make One With GCD 2
|
ユーザー |
|
提出日時 | 2020-04-25 17:57:52 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 1,438 bytes |
コンパイル時間 | 2,081 ms |
コンパイル使用メモリ | 172,300 KB |
実行使用メモリ | 86,912 KB |
最終ジャッジ日時 | 2024-11-07 11:21:57 |
合計ジャッジ時間 | 20,406 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 WA * 2 |
other | AC * 28 WA * 1 RE * 12 |
ソースコード
#include <bits/stdc++.h> using namespace std; using ll = long long; #define SPEED cin.tie(0);ios::sync_with_stdio(false); template<class Operator> class SparseTable{ public: Operator Op; using typeNode = decltype(Op.unitNode); vector<typeNode> node; vector<int> idx; size_t depth; size_t length; SparseTable(const vector<typeNode>& vec) { for(depth = 0;(1<<depth)<=vec.size();++depth); length = (1<<depth); node.resize(depth*length); for(int i = 0; i < vec.size(); ++i) node[i] = vec[i]; for(int i = 1; i < depth; ++i) for(int j = 0; j + (1<<i) < (1<<depth); ++j) node[i*length+j] = Op.funcNode(node[(i-1)*length+j],node[(i-1)*length+j + (1 << (i-1))]); idx.resize(vec.size()+1); for(int i = 2; i < vec.size()+1; ++i) idx[i] = idx[i>>1] + 1; } //[l,r) typeNode get(int l,int r) { int bit = idx[r-l]; return Op.funcNode(node[bit*length+l],node[bit*length+r - (1<<bit)]); } }; template<class typeNode> struct nodeGCD { typeNode unitNode = 0; typeNode funcNode(typeNode l,typeNode r){return ((r == 0) ? l : funcNode(r, l % r));} }; // solution by binary search in arbitary range on disjn sparse table int main() { SPEED ll N; cin >> N; vector<ll> A(N+1,1); for(int i = 0; i < N; ++i) cin >> A[i]; SparseTable<nodeGCD<ll>> st(A); ll ans=0; int r=1; for(int i=0; i<N; i++){ while(r<=N){ if(st.get(i, r)==1) break; r++; } ans+=N+1-r; } cout << ans << endl; }